David Cournapeau wrote: > Here we are (the difference may not seem that much between the two > multiple_ga, but in reality, _diag_gauss_den is an internal function > which is done in C, and is much faster... By writing this example, I've > just realized that the function _diag_gauss_den may be slow for exactly > the same reasons): > > I checked my assumption about the _diag_gauss_den function being slow because of the same problem, and this is indeed the case:
If I replace X = randn(nframes, d) by X = randn(d, nframes); X = X.T, the function is mode than twice as fast ! This seems way to much of a difference to me.... David ------------------------------------------------------------------------- Take Surveys. Earn Cash. Influence the Future of IT Join SourceForge.net's Techsay panel and you'll get the chance to share your opinions on IT & business topics through brief surveys -- and earn cash http://www.techsay.com/default.php?page=join.php&p=sourceforge&CID=DEVDEV _______________________________________________ Numpy-discussion mailing list Numpy-discussion@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpy-discussion