[Numpy-discussion] Quick Question about Optimization

Charles R Harris charlesr.harris@gmail....
Mon May 19 17:27:32 CDT 2008

On Mon, May 19, 2008 at 2:53 PM, Christopher Barker <Chris.Barker@noaa.gov>

> Anne Archibald wrote:
> > 2008/5/19 James Snyder <jbsnyder@gmail.com>:
> >> I can provide the rest of the code if needed, but it's basically just
> >> filling some vectors with random and empty data and initializing a few
> >> things.
> >
> > It would kind of help, since it would make it clearer what's a scalar
> > and what's an array, and what the dimensions of the various arrays
> > are.
> It would also help if you provided a complete example (as little code as
> possible), so we could try out and time our ideas before suggesting them.
> >> np.random.standard_normal(size=(1,self.naff))
> Anyone know how fast this is compared to Matlab? That could be the
> difference right there.

The latest versions of Matlab use the ziggurat method to generate random
normals and it is faster than the method used in numpy. I have ziggurat code
at hand, but IIRC, Robert doesn't trust the method ;) I don't know if it
would actually speed things up, though.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://projects.scipy.org/pipermail/numpy-discussion/attachments/20080519/3a70c94e/attachment.html 

More information about the Numpy-discussion mailing list