Mark P. Miller wrote: >Robert: Just a thought on this topic: > >Would it be possible for the Scipy folks to add a new module based >solely off your old mtrand code (pre-broadcast)? I have to say that the >mtrand code from numpy 0.9.8 has some excellent advantages over the core >python random number generators. > > I'm not convinced that the broadcasting is causing the slow-downs. Currently, the code has two path-ways. One gets called when the inputs are scalars which is equivalent to the old code and the second gets called when broadcasting is desired. Their should be only a very small speed hit because of the check. So, it should be possible to do both without great differences in speed.
Perhaps the speed hit is due to something else (like creating 0-d arrays from scalars, or something similar). So, in other words, I think the problem is fixable, but we need someone who can track down where the slowness is coming from. -Travis _______________________________________________ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion
