Re: [Numpy-discussion] A few comments (MDP conversion)

2006-07-13 Thread Robert Kern
Pietro Berkes wrote: > We wondered whether it would be possible to obtain SVN write access in > orderto be able to fix this kind of issues by ourselves in the future. > We could also contribute docstrings for some of the functions. The best way to get SVN privileges in any open source project i

Re: [Numpy-discussion] numpy.average(fooArray, axis=0): dropping nans from calculation

2006-07-13 Thread Eric Firing
Webb Sprague wrote: > Could someone recommend a way to average an array along the columns > without propagating the nans and without turning them into some weird > number which bias the result? I guess I can just keep using an > indexing array for fooArray, but if there is somehting more graceful,

[Numpy-discussion] numpy.vectorize performance

2006-07-13 Thread Nick Fotopoulos
Dear all, I often make use of numpy.vectorize to make programs read more like the physics equations I write on paper. numpy.vectorize is basically a wrapper for numpy.frompyfunc. Reading Travis's Scipy Book (mine is dated Jan 6 2005) kind of suggests to me that it returns a full- fledged

[Numpy-discussion] numpy.average(fooArray, axis=0): dropping nans from calculation

2006-07-13 Thread Webb Sprague
Could someone recommend a way to average an array along the columns without propagating the nans and without turning them into some weird number which bias the result? I guess I can just keep using an indexing array for fooArray, but if there is somehting more graceful, I would love to know. Boy

[Numpy-discussion] ANN: Python Enthought Edition Version 1.0.0.beta4 Released

2006-07-13 Thread Bryce Hendrix
Enthought is pleased to announce the release of Python Enthought Edition Version 1.0.0.beta4 (http://code.enthought.com/enthon/) -- a python distribution for Windows. 1.0.0.beta4 Release Notes: There are two known issues:   * No documentation is included due to problems

Re: [Numpy-discussion] [Fwd: Re: A few comments (MDP conversion)]

2006-07-13 Thread David M. Cooke
On Jul 13, 2006, at 12:39 , Pietro Berkes wrote: > On Thu, 13 Jul 2006, Nils Wagner wrote: > >> It seems to be line 281 instead of 269. I am using latest svn. >> BTW, in linalg.py in >> >> def pinv: >> >> there is another Complex with capital C. > > Well, the problem is not really the capital 'C'

Re: [Numpy-discussion] A few comments (MDP conversion)

2006-07-13 Thread David M. Cooke
On Jul 13, 2006, at 09:11 , Pietro Berkes wrote: > Dear numpys, > > a couple of weeks ago Tiziano and I completed the conversion of our > data > processing library MDP to numpy. We collected a few ideas and > questions: > > - we found the convertcode.py module quite useful to perform a first, >

Re: [Numpy-discussion] Finding indices for all rows which contain nan (or anything else)

2006-07-13 Thread Webb Sprague
On 7/13/06, Robert Kern <[EMAIL PROTECTED]> wrote: > Webb Sprague wrote: > > Does anyone have some vectorized code that pulls out all the row > > def is_row_nan(a): > return numpy.isnan(a).any(axis=-1) I knew there was a way, but I didn't know to check any() and all(). Thanks to all (I love f

Re: [Numpy-discussion] Finding indices for all rows which contain nan (or anything else)

2006-07-13 Thread Robert Kern
Webb Sprague wrote: > Does anyone have some vectorized code that pulls out all the row > indices for any row which has an nan (or a number less than 1 or > whatever). I want to subsequently be able to perform an operation > with all the good rows. See the imaginary code below. > > a = numpy.arra

[Numpy-discussion] Finding indices for all rows which contain nan (or anything else)

2006-07-13 Thread Webb Sprague
Does anyone have some vectorized code that pulls out all the row indices for any row which has an nan (or a number less than 1 or whatever). I want to subsequently be able to perform an operation with all the good rows. See the imaginary code below. a = numpy.array([[1,2],[nan,1], [2,3]]) is_row

Re: [Numpy-discussion] how to use argsort result?

2006-07-13 Thread Eric Firing
> > Would it be reasonable if argsort returned the complete tuple of > indices, so that > A[A.argsort(ax)] would work ? +1 This is the behavior one would naturally expect. Eric - Using Tomcat but need to do more? Need to

Re: [Numpy-discussion] RandomArray module

2006-07-13 Thread Robert Kern
Sebastian Żurek wrote: > Hi All, > > Has anyone worked with the RandomArray module? I wonder, > if it's OK to use its pseudo-random numbers generators, or > maybe I shall find more trusted methods (ie. ran1 from Numerical Recipes)? At this point in time, I don't recommend using RandomArray if you

Re: [Numpy-discussion] [Fwd: Re: A few comments (MDP conversion)]

2006-07-13 Thread Pietro Berkes
On Thu, 13 Jul 2006, Nils Wagner wrote: > It seems to be line 281 instead of 269. I am using latest svn. > BTW, in linalg.py in > > def pinv: > > there is another Complex with capital C. Well, the problem is not really the capital 'C', but rather the lack of quotation marks... ---

[Numpy-discussion] [Fwd: Re: A few comments (MDP conversion)]

2006-07-13 Thread Nils Wagner
--- Begin Message --- Pietro Berkes wrote: Dear numpys, a couple of weeks ago Tiziano and I completed the conversion of our data processing library MDP to numpy. We collected a few ideas and questions: - we found the convertcode.py module quite useful to perform a first, low-level conversio

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

[Numpy-discussion] A few comments (MDP conversion)

2006-07-13 Thread Pietro Berkes
Dear numpys, a couple of weeks ago Tiziano and I completed the conversion of our data processing library MDP to numpy. We collected a few ideas and questions: - we found the convertcode.py module quite useful to perform a first, low-level conversion. We had some problem when 'typecode' was us

[Numpy-discussion] generalized_inverse

2006-07-13 Thread Jon Peirce
There used to be a function generalized_inverse in the numpy.linalg module (certainly in 0.9.2). In numpy0.9.8 it seems to have been moved to the numpy.linalg.old subpackage. Does that mean it's being dropped? Did it have to move? Now i have to add code to my package to try both locations becau

[Numpy-discussion] Voce Recebeu uma Charge do Humortadela

2006-07-13 Thread Humortadela
Title: Olá   Olá! Alguém que não tinha nada para fazer, numa de suas visitas ao Humor Tadela não sei por que cargas d'água, lhe recomendou a seguinte p

Re: [Numpy-discussion] how to use argsort result?

2006-07-13 Thread Pau Gargallo
On 7/13/06, Travis Oliphant <[EMAIL PROTECTED]> wrote: > Pau Gargallo wrote: > > On 7/12/06, Victoria G. Laidler <[EMAIL PROTECTED]> wrote: > > > >> Hi, > >> > >> Pardon me if I'm reprising an earlier discussion, as I'm new to the list. > >> > >> But is there a reason that this obscure syntax > >>

Re: [Numpy-discussion] LAPACK Problem with "import numpy" on AIX

2006-07-13 Thread Adam Carter
> Lapack_lite is in the numpy sources. I think it's the fallback if you don't have a LAPACK. Thanks Bill, that's what I thought, but I don't understand why it can't find it. > You can check the config that got built with: > numpy.show_config() OK, so it appears that NumPy is imported then, a

Re: [Numpy-discussion] LAPACK Problem with "import numpy" on AIX

2006-07-13 Thread Bill Baxter
On 7/13/06, Adam Carter <[EMAIL PROTECTED]> wrote: Hi all,Do I need lapack_lite aswell? Where can I get it? If I need this otherLAPACK, how can I ensure that my code uses the optimised version of LAPACKalready on this system I'm using? Lapack_lite is in the numpy sources.  I think it's the fallback

[Numpy-discussion] LAPACK Problem with "import numpy" on AIX

2006-07-13 Thread Adam Carter
Hi all, I'm new to this list so apologies if this is a solved problem, but I haven't been able to find anything in the archives. I've just installed Python 2.4.3 and Numpy-0.9.8 on AIX, and the configure/make/install of Python and the install of Numpy _appeared_ to go smoothly. However, now when

Re: [Numpy-discussion] int 'up'casted to uint

2006-07-13 Thread Pau Gargallo
On 7/11/06, Pau Gargallo <[EMAIL PROTECTED]> wrote: > On 7/11/06, Travis Oliphant <[EMAIL PROTECTED]> wrote: > > Pau Gargallo wrote: > > > hi, > > > > > > looking at the upcasting table at > > > http://www.scipy.org/Tentative_NumPy_Tutorial#head-4c1d53fe504adc97baf27b65513b4b97586a4fc5 > > > I saw

Re: [Numpy-discussion] how to use argsort result?

2006-07-13 Thread Travis Oliphant
Pau Gargallo wrote: > On 7/12/06, Victoria G. Laidler <[EMAIL PROTECTED]> wrote: > >> Hi, >> >> Pardon me if I'm reprising an earlier discussion, as I'm new to the list. >> >> But is there a reason that this obscure syntax >> >> A[arange(2)[:,newaxis],indexes] >> >> A[arange(A.shape[0])[:,newaxi

Re: [Numpy-discussion] how to use argsort result?

2006-07-13 Thread Pau Gargallo
On 7/12/06, Victoria G. Laidler <[EMAIL PROTECTED]> wrote: > Hi, > > Pardon me if I'm reprising an earlier discussion, as I'm new to the list. > > But is there a reason that this obscure syntax > > A[arange(2)[:,newaxis],indexes] > > A[arange(A.shape[0])[:,newaxis],indexes] > > is preferable to the

Re: [Numpy-discussion] RandomArray module

2006-07-13 Thread Francesc Altet
A Dijous 13 Juliol 2006 01:07, Sebastian Żurek va escriure: > Hi All, > > Has anyone worked with the RandomArray module? I wonder, > if it's OK to use its pseudo-random numbers generators, or > maybe I shall find more trusted methods (ie. ran1 from Numerical Recipes)? I'm not an expert, but my und

[Numpy-discussion] RandomArray module

2006-07-13 Thread Sebastian Żurek
Hi All, Has anyone worked with the RandomArray module? I wonder, if it's OK to use its pseudo-random numbers generators, or maybe I shall find more trusted methods (ie. ran1 from Numerical Recipes)? Please, give some comments. Thanks. Sebastian