Re: [Numpy-discussion] speeding up y[:, i] for y.shape = (big number, small number)

2006-10-04 Thread Bruce Southey
Hi, I think what you are after is the multivariate normal distribution. Assuming I have it correctly, it is clearer to see (and probably more accurate to compute) in the log form as: -(N/2)*log(2*PI) - 0.5*log(determinant of V) - 0.5*(transpose of (x-mu))*inverse(V)*(x-mu) where N is the number o

Re: [Numpy-discussion] Vectorizing code, for loops, and all that

2006-10-04 Thread Bruce Southey
Hi, On 10/2/06, Albert Strasheim <[EMAIL PROTECTED]> wrote: > I recently started looking at David Cournapeau's PyEM package, > specifically his implementation of the K-Means algorithm. He > implemented part of this algorithm with in pure Python version and > also provided a Pyrex alternative that

Re: [Numpy-discussion] RandomArray module

2006-07-13 Thread Bruce Southey
Hi, As Francesc meintioned, Robert Kern did a great job of replacing ranlib: "numpy.random uses the Mersenne Twister. RANLIB is dead! Long live MT19937!" So throw away ran1! Regards Bruce On 7/12/06, Sebastian Żurek <[EMAIL PROTECTED]> wrote: > Hi All, > > Has anyone worked with the RandomArray

Re: [Numpy-discussion] Call for a vote on .M .A .T .H attributes

2006-07-07 Thread Bruce Southey
On 7/6/06, Travis Oliphant <[EMAIL PROTECTED]> wrote: > > This is a call for a vote on each of the math attributes. Please post > your vote as > > +1 : support > +0 : don't care so go ahead > -0 : don't care so why do it > -1 : against > > Vote on the following issues separately: > > > > 1) .T

Re: [Numpy-discussion] Time for beta1 of NumPy 1.0

2006-07-02 Thread Bruce Southey
Hi, Linux Weekly News (http://lwn.net) had an very interesting article on KDE's switch on June 19, 2006 by Alexander Neundorf: http://lwn.net/Articles/187923/ The full article is at: http://lwn.net/Articles/188693/ This should be freely available to all. Also, the current US Linux Magazine (Jun

Re: [Numpy-discussion] Time for beta1 of NumPy 1.0

2006-06-30 Thread Bruce Southey
Hi, You should be encouraged by the trend from Numeric to numarray because the tar users clearly are prepared to upgrade. In terms of the education program, the 1.0 release is the best starting point as there is a general phobia for pre-1.0 releases (and dot zero releases). Also, Python 2.5 is co

Re: [Numpy-discussion] Distance Matrix speed

2006-06-16 Thread Bruce Southey
Hi, Please run the exact same code in Matlab that you are running in NumPy. Many of Matlab functions are very highly optimized so these are provided as binary functions. I think that you are running into this so you are not doing the correct comparison So the ways around it are to write an extensi

Re: [Numpy-discussion] Any Numeric or numarray users on this list?

2006-05-31 Thread Bruce Southey
Hi, On 5/30/06, Travis Oliphant <[EMAIL PROTECTED]> wrote: Please help the developers by responding to a few questions. 1) Have you transitioned or started to transition to NumPy (i.e. import numpy)? Yes and No 2) Will you transition within the next 6 months? (if you answered No to #1)