Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Robert Kern
On Tue, Sep 30, 2008 at 01:34, Charles R Harris <[EMAIL PROTECTED]> wrote: > > On Tue, Sep 30, 2008 at 12:24 AM, Robert Kern <[EMAIL PROTECTED]> wrote: >> OTOH, Python 2.6 and up will be following to the C99 standard as >> closely as possible. I would prefer to keep up with them. It's true >> that

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Tue, Sep 30, 2008 at 12:24 AM, Robert Kern <[EMAIL PROTECTED]> wrote: > On Tue, Sep 30, 2008 at 00:46, Nathan Bell <[EMAIL PROTECTED]> wrote: > > On Tue, Sep 30, 2008 at 1:20 AM, Robert Kern <[EMAIL PROTECTED]> > wrote: > >> > >> F.9.9.2 The fmax functions > >> 1 If just one argument is a NaN,

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 11:56 PM, David Cournapeau < [EMAIL PROTECTED]> wrote: > Nathan Bell wrote: > > > > +1 for NaN-always-propagates since we have explicit variants for the > > alternative semantics. > > > > Users are more likely to remember that "NaNs always propagate" than > > "as stated in

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Robert Kern
On Tue, Sep 30, 2008 at 00:46, Nathan Bell <[EMAIL PROTECTED]> wrote: > On Tue, Sep 30, 2008 at 1:20 AM, Robert Kern <[EMAIL PROTECTED]> wrote: >> >> F.9.9.2 The fmax functions >> 1 If just one argument is a NaN, the fmax functions return the other >> argument (if both arguments are NaNs, the funct

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread David Cournapeau
Nathan Bell wrote: > > +1 for NaN-always-propagates since we have explicit variants for the > alternative semantics. > > Users are more likely to remember that "NaNs always propagate" than > "as stated in the C99 standard...". I don't know. I would like to agree with you, but OTOH, starting to go

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread David Cournapeau
David Cournapeau wrote: > > when run tells me: > > signbit(NAN) = 0 > signbit(-NAN) = -2147483648 > signbit((-1) * NAN) = 0 > signbit(- NAN + NAN) = -2147483648 > signbit(NAN - NAN) = 0 > > Does not this indicate that signbit(NAN) is undefined ? I guess I am > afraid that the glibc NAN is just one

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Nathan Bell
On Tue, Sep 30, 2008 at 1:20 AM, Robert Kern <[EMAIL PROTECTED]> wrote: > > F.9.9.2 The fmax functions > 1 If just one argument is a NaN, the fmax functions return the other > argument (if both arguments are NaNs, the functions return a NaN). > 2 The body of the fmax function might be > {return (is

Re: [Numpy-discussion] nonuniform scatter operations

2008-09-29 Thread Nathan Bell
On Sun, Sep 28, 2008 at 4:15 PM, Geoffrey Irving <[EMAIL PROTECTED]> wrote: > > Thanks. That works great. A slightly cleaner version is > >X += coo_matrix((Y, (K, zeros_like(K.sum(axis=1) > > The next question is: is there a similar way that generalizes to the > case where X is n by 3 and

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 11:20 PM, Robert Kern <[EMAIL PROTECTED]> wrote: > On Mon, Sep 29, 2008 at 20:22, Charles R Harris > <[EMAIL PROTECTED]> wrote: > > > > On Mon, Sep 29, 2008 at 4:28 PM, Robert Kern <[EMAIL PROTECTED]> > wrote: > >> > >> On Mon, Sep 29, 2008 at 17:13, Charles R Harris > >> <

Re: [Numpy-discussion] Are there command similar as Matlab find command?

2008-09-29 Thread Anne Archibald
2008/9/30 frank wang <[EMAIL PROTECTED]>: > I really do not know the difference of debug mode and the pdb debugger. To > me, it seems that pdb is only way to debug the python code. How do the > expert of numpy/python debug their code? Are there any more efficient way to > debug the code in python

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Robert Kern
On Mon, Sep 29, 2008 at 20:22, Charles R Harris <[EMAIL PROTECTED]> wrote: > > On Mon, Sep 29, 2008 at 4:28 PM, Robert Kern <[EMAIL PROTECTED]> wrote: >> >> On Mon, Sep 29, 2008 at 17:13, Charles R Harris >> <[EMAIL PROTECTED]> wrote: >> > >> > On Mon, Sep 29, 2008 at 3:52 PM, Charles R Harris >> >

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 10:50 PM, Robert Kern <[EMAIL PROTECTED]> wrote: > On Mon, Sep 29, 2008 at 23:02, Charles R Harris > <[EMAIL PROTECTED]> wrote: > > > > > > On Mon, Sep 29, 2008 at 9:02 PM, David Cournapeau > > <[EMAIL PROTECTED]> wrote: > >> > >> Charles R Harris wrote: > >> > > >> > So th

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread David Cournapeau
Robert Kern wrote: > > numpy.signbit() should work like C99 signbit() (where possible), IMO. > It can only return (integer) 0 or 1, and it does differentiate between > NAN and -NAN. I don't think we should invent new semantics if we can > avoid it. Agreed, but for signbit case, I can't find what t

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread David Cournapeau
Charles R Harris wrote: > > umathmodule.c is cleaned/done. Or was until I put in some of the nan > handling. But some of the nan comparisons were just odd, in particular > the sign function, which was effectively undefined. And since the > current behavior is an accident of comparison choices and o

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 10:09 PM, David Cournapeau < [EMAIL PROTECTED]> wrote: > Charles R Harris wrote: > > > > Yes, that too. But I was thinking of the ufunc returning nan when needed > > > > I think this is better for consistency, yes. NaN is not a number, so it > has no sign :) More serious

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Robert Kern
On Mon, Sep 29, 2008 at 23:02, Charles R Harris <[EMAIL PROTECTED]> wrote: > > > On Mon, Sep 29, 2008 at 9:02 PM, David Cournapeau > <[EMAIL PROTECTED]> wrote: >> >> Charles R Harris wrote: >> > >> > So the proposition is, sign, max, min return nan when any of the >> > arguments is nan. >> >> Note

Re: [Numpy-discussion] Are there command similar as Matlab find command?

2008-09-29 Thread David Cournapeau
frank wang wrote: > Hi, Robert, > > Thank you very much for clarification of this. > > I really do not know the difference of debug mode and the pdb debugger. That's totally different in principle. Debug mode means your code is compiled/interpreted differently to keep as much information as po

Re: [Numpy-discussion] Are there command similar as Matlab find command?

2008-09-29 Thread frank wang
Hi, Robert, Thank you very much for clarification of this. I really do not know the difference of debug mode and the pdb debugger. To me, it seems that pdb is only way to debug the python code. How do the expert of numpy/python debug their code? Are there any more efficient way to debug the

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread David Cournapeau
Charles R Harris wrote: > > Yes, that too. But I was thinking of the ufunc returning nan when needed I think this is better for consistency, yes. NaN is not a number, so it has no sign :) More seriously, I think those features should be clearly listed and thought out before being implemented, part

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 9:53 PM, David Cournapeau < [EMAIL PROTECTED]> wrote: > Charles R Harris wrote: > > > > Do any of the MS compilers handle these things correctly? > > Don't know. To be 100 % honest, one of the problem for MS compilers is > the /Ox flag (for IEEE-754 rules). This should not

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 9:54 PM, David Cournapeau < [EMAIL PROTECTED]> wrote: > Charles R Harris wrote: > > > > > > On Mon, Sep 29, 2008 at 9:02 PM, David Cournapeau > > <[EMAIL PROTECTED] > > > wrote: > > > > Charles R Harris wrote: > > > > > > So the proposi

Re: [Numpy-discussion] Proposal: scipy.spatial

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 9:24 PM, Anne Archibald <[EMAIL PROTECTED]>wrote: > Hi, > > Once again there has been a thread on the numpy/scipy mailing lists > requesting (essentially) some form of spatial data structure. Pointers > have been posted to ANN (sadly LGPLed and in C++) as well as a handful

Re: [Numpy-discussion] Proposal: scipy.spatial

2008-09-29 Thread Christopher Barker
Anne Archibald wrote: > I suggest the creation of > a new submodule of scipy, scipy.spatial, +1 Here's one to consider: http://pypi.python.org/pypi/Rtree and perhaps other stuff from: http://trac.gispython.org/spatialindex/wiki which I think is LGPL -- can scipy use that? By the way, a gener

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread David Cournapeau
Charles R Harris wrote: > > > On Mon, Sep 29, 2008 at 9:02 PM, David Cournapeau > <[EMAIL PROTECTED] > > wrote: > > Charles R Harris wrote: > > > > So the proposition is, sign, max, min return nan when any of the > > arguments is nan. > > Note that inte

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread David Cournapeau
Charles R Harris wrote: > > Do any of the MS compilers handle these things correctly? Don't know. To be 100 % honest, one of the problem for MS compilers is the /Ox flag (for IEEE-754 rules). This should not be used for numpy, period (I am sure you could break numpy with gcc and -ffast-math and co

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 9:02 PM, David Cournapeau < [EMAIL PROTECTED]> wrote: > Charles R Harris wrote: > > > > So the proposition is, sign, max, min return nan when any of the > > arguments is nan. > > Note that internally, signbit (the C function) returns an integer. > That is the signature of

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
Hi David, On Mon, Sep 29, 2008 at 9:07 PM, David Cournapeau < [EMAIL PROTECTED]> wrote: > Charles R Harris wrote: > > > > Thanks, that was very helpful. I wonder how widespread the less, > > lessequal, etc. macros are? > > If it is C99, count on some platforms (MS in particular) to not do it. > A

[Numpy-discussion] Proposal: scipy.spatial

2008-09-29 Thread Anne Archibald
Hi, Once again there has been a thread on the numpy/scipy mailing lists requesting (essentially) some form of spatial data structure. Pointers have been posted to ANN (sadly LGPLed and in C++) as well as a handful of pure-python implementations of kd-trees. I suggest the creation of a new submodul

Re: [Numpy-discussion] Numpy Nu-bee: "forward fill" function

2008-09-29 Thread Pierre GM
On Monday 29 September 2008 22:31:47 Joshua Ford wrote: > I would like to conduct the equivalent of a "forward fill" function on an > array, that repeats the last non-null value until the next non-null value > is reached. I saw that the scipy TimeSeries has a "forward_fill" > function: > > http:/

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread David Cournapeau
Charles R Harris wrote: > > Thanks, that was very helpful. I wonder how widespread the less, > lessequal, etc. macros are? If it is C99, count on some platforms (MS in particular) to not do it. Also, when doing things in C, beware that some compilers break most reasonable expectations about floati

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread David Cournapeau
Charles R Harris wrote: > > So the proposition is, sign, max, min return nan when any of the > arguments is nan. Note that internally, signbit (the C function) returns an integer. > > Complex numbers are more complicated because we first compare the real > parts, then the imaginary. Arguably 1 >

[Numpy-discussion] Numpy Nu-bee: "forward fill" function

2008-09-29 Thread Joshua Ford
Hello all! I'm a python and numpy newbie and I'm having difficulty doing something that seems to be relatively straightforward. So I figured I would ask the experts! I would like to conduct the equivalent of a "forward fill" function on an array, that repeats the last non-null value until t

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 4:28 PM, Robert Kern <[EMAIL PROTECTED]> wrote: > On Mon, Sep 29, 2008 at 17:13, Charles R Harris > <[EMAIL PROTECTED]> wrote: > > > > On Mon, Sep 29, 2008 at 3:52 PM, Charles R Harris > > <[EMAIL PROTECTED]> wrote: > >> > >> Hi All, > >> > >> I've been cleaning up the ufun

Re: [Numpy-discussion] Are there command similar as Matlab find command?

2008-09-29 Thread Robert Kern
On Mon, Sep 29, 2008 at 20:16, frank wang <[EMAIL PROTECTED]> wrote: > Thanks for the help. > > It seems that the where command has problem when I tried to run it in the > debug mode. It does not return any thing such as: > > (Pdb) aa=array([1,2,3,4] > (Pdb) where(aa>2) > (1)() >> c:\dhg\docsis\l

Re: [Numpy-discussion] Are there command similar as Matlab find command?

2008-09-29 Thread frank wang
Thanks for the help. It seems that the where command has problem when I tried to run it in the debug mode. It does not return any thing such as: (Pdb) aa=array([1,2,3,4] (Pdb) where(aa>2) (1)()> c:\dhg\docsis\lab_test\parseadc.py(70)parsempeg()-> bb=array(fid).astype('int') (Pdb) It doe

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 5:16 PM, Robert Kern <[EMAIL PROTECTED]> wrote: > On Mon, Sep 29, 2008 at 18:10, Charles R Harris > <[EMAIL PROTECTED]> wrote: > > > > On Mon, Sep 29, 2008 at 4:54 PM, Charles R Harris > > <[EMAIL PROTECTED]> wrote: > >> > >> On Mon, Sep 29, 2008 at 4:40 PM, Charles R Harri

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Robert Kern
On Mon, Sep 29, 2008 at 18:10, Charles R Harris <[EMAIL PROTECTED]> wrote: > > On Mon, Sep 29, 2008 at 4:54 PM, Charles R Harris > <[EMAIL PROTECTED]> wrote: >> >> On Mon, Sep 29, 2008 at 4:40 PM, Charles R Harris >> <[EMAIL PROTECTED]> wrote: >>> >>> On Mon, Sep 29, 2008 at 4:28 PM, Robert Kern <[

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 4:54 PM, Charles R Harris <[EMAIL PROTECTED] > wrote: > > > On Mon, Sep 29, 2008 at 4:40 PM, Charles R Harris < > [EMAIL PROTECTED]> wrote: > >> >> >> On Mon, Sep 29, 2008 at 4:28 PM, Robert Kern <[EMAIL PROTECTED]>wrote: >> >>> On Mon, Sep 29, 2008 at 17:13, Charles R Harr

Re: [Numpy-discussion] problem calling SD.SDim.setname in pyhdf

2008-09-29 Thread Robert Kern
On Mon, Sep 29, 2008 at 17:39, Catherine Moroney <[EMAIL PROTECTED]> wrote: > Hi, > > I recently upgraded pyhdf from 0.7-3 to 0.8-2, and have noticed some > problems with the SD.SDim interface. > > Specifically, the short test code attached works fine with version 0.7-3 > and fails with version 0.8

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 4:40 PM, Charles R Harris <[EMAIL PROTECTED] > wrote: > > > On Mon, Sep 29, 2008 at 4:28 PM, Robert Kern <[EMAIL PROTECTED]>wrote: > >> On Mon, Sep 29, 2008 at 17:13, Charles R Harris >> <[EMAIL PROTECTED]> wrote: >> > >> > On Mon, Sep 29, 2008 at 3:52 PM, Charles R Harris

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 4:28 PM, Robert Kern <[EMAIL PROTECTED]> wrote: > On Mon, Sep 29, 2008 at 17:13, Charles R Harris > <[EMAIL PROTECTED]> wrote: > > > > On Mon, Sep 29, 2008 at 3:52 PM, Charles R Harris > > <[EMAIL PROTECTED]> wrote: > >> > >> Hi All, > >> > >> I've been cleaning up the ufun

[Numpy-discussion] problem calling SD.SDim.setname in pyhdf

2008-09-29 Thread Catherine Moroney
Hi, I recently upgraded pyhdf from 0.7-3 to 0.8-2, and have noticed some problems with the SD.SDim interface. Specifically, the short test code attached works fine with version 0.7-3 and fails with version 0.8-2 when it's trying to assign a name to the dimension. I get the following error (below

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Robert Kern
On Mon, Sep 29, 2008 at 17:13, Charles R Harris <[EMAIL PROTECTED]> wrote: > > On Mon, Sep 29, 2008 at 3:52 PM, Charles R Harris > <[EMAIL PROTECTED]> wrote: >> >> Hi All, >> >> I've been cleaning up the ufunc loops and the sign function currently >> doesn't have a defined behavior for nans. This m

Re: [Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 3:52 PM, Charles R Harris <[EMAIL PROTECTED] > wrote: > Hi All, > > I've been cleaning up the ufunc loops and the sign function currently > doesn't have a defined behavior for nans. This makes the results depend on > the order/type of comparisons in the code, which looks fr

[Numpy-discussion] What is the sign of nan?

2008-09-29 Thread Charles R Harris
Hi All, I've been cleaning up the ufunc loops and the sign function currently doesn't have a defined behavior for nans. This makes the results depend on the order/type of comparisons in the code, which looks fragile to me. So what should it return? I vote for nan but am open for suggestions. Chuc

Re: [Numpy-discussion] Are there command similar as Matlab find command?

2008-09-29 Thread Travis E. Oliphant
frank wang wrote: > Hi, > > I am trying to find a command in numpy or python that perform similar > function as Matlab find command. It will return the indexes of array > that satisfy a condition. So far I have not found anything. There are several ways to do this, but what are you trying to d

Re: [Numpy-discussion] Are there command similar as Matlab find command?

2008-09-29 Thread Nathan Bell
On Mon, Sep 29, 2008 at 5:32 PM, frank wang <[EMAIL PROTECTED]> wrote: > Hi, > > I am trying to find a command in numpy or python that perform similar > function as Matlab find command. It will return the indexes of array that > satisfy a condition. So far I have not found anything. > If you're fa

[Numpy-discussion] Are there command similar as Matlab find command?

2008-09-29 Thread frank wang
Hi, I am trying to find a command in numpy or python that perform similar function as Matlab find command. It will return the indexes of array that satisfy a condition. So far I have not found anything. Thanks Frank> Date: Thu, 25 Sep 2008 09:51:46 +0200> From: [EMAIL PROTECTED]> To: nump

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Anne Archibald
2008/9/29 Francesc Alted <[EMAIL PROTECTED]>: > But I agree with Robert Kern that tests *should* never be run under -OO > because they are *tests*, not production code, and using the -OO for > running the test units is simply a wrong practice. I have to disagree with this. What if some piece of c

Re: [Numpy-discussion] should tests use assert?

2008-09-29 Thread Alan G Isaac
> On Mon, Sep 29, 2008 at 12:43 PM, Francesc Alted wrote:` >> But I agree with Robert Kern that tests *should* never be >> run under -OO >> because they are *tests*, not production code, and using >> the -OO for >> running the test units is simply a wrong practice. Charles R Harris wrote: > Wha

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Francesc Alted
A Monday 29 September 2008, Charles R Harris escrigué: > On Mon, Sep 29, 2008 at 12:43 PM, Francesc Alted <[EMAIL PROTECTED]>wrote: > > A Monday 29 September 2008, Charles R Harris escrigué: > > > On Mon, Sep 29, 2008 at 11:50 AM, Robert Kern > > > <[EMAIL PROTECTED]> > > > > wrote: > > > > On Mon

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Francesc Alted
A Monday 29 September 2008, Robert Kern escrigué: > On Mon, Sep 29, 2008 at 13:43, Francesc Alted <[EMAIL PROTECTED]> wrote: > > Plain assert clauses used to be the most common way to check for > > test units in the original "unittest" module. Later on, the > > authors started to introduce things

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 12:43 PM, Francesc Alted <[EMAIL PROTECTED]>wrote: > A Monday 29 September 2008, Charles R Harris escrigué: > > On Mon, Sep 29, 2008 at 11:50 AM, Robert Kern <[EMAIL PROTECTED]> > wrote: > > > On Mon, Sep 29, 2008 at 12:40, Charles R Harris > > > > > > <[EMAIL PROTECTED]> w

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Robert Kern
On Mon, Sep 29, 2008 at 13:43, Francesc Alted <[EMAIL PROTECTED]> wrote: > Plain assert clauses used to be the most common way to check for test > units in the original "unittest" module. Later on, the authors started > to introduce things like self.assert_() and family and warned that the > use

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 12:32 PM, Robert Kern <[EMAIL PROTECTED]> wrote: > On Mon, Sep 29, 2008 at 13:22, Charles R Harris > <[EMAIL PROTECTED]> wrote: > > > > On Mon, Sep 29, 2008 at 12:17 PM, Robert Kern <[EMAIL PROTECTED]> > wrote: > >> > >> On Mon, Sep 29, 2008 at 13:07, Charles R Harris > >>

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Francesc Alted
A Monday 29 September 2008, Charles R Harris escrigué: > On Mon, Sep 29, 2008 at 11:50 AM, Robert Kern <[EMAIL PROTECTED]> wrote: > > On Mon, Sep 29, 2008 at 12:40, Charles R Harris > > > > <[EMAIL PROTECTED]> wrote: > > > On Mon, Sep 29, 2008 at 11:35 AM, Robert Kern > > > <[EMAIL PROTECTED]> > >

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Robert Kern
On Mon, Sep 29, 2008 at 13:22, Charles R Harris <[EMAIL PROTECTED]> wrote: > > On Mon, Sep 29, 2008 at 12:17 PM, Robert Kern <[EMAIL PROTECTED]> wrote: >> >> On Mon, Sep 29, 2008 at 13:07, Charles R Harris >> <[EMAIL PROTECTED]> wrote: >> > >> > On Mon, Sep 29, 2008 at 11:50 AM, Robert Kern <[EMAIL

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 12:17 PM, Robert Kern <[EMAIL PROTECTED]> wrote: > On Mon, Sep 29, 2008 at 13:07, Charles R Harris > <[EMAIL PROTECTED]> wrote: > > > > On Mon, Sep 29, 2008 at 11:50 AM, Robert Kern <[EMAIL PROTECTED]> > wrote: > >> > >> On Mon, Sep 29, 2008 at 12:40, Charles R Harris > >>

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Robert Kern
On Mon, Sep 29, 2008 at 13:07, Charles R Harris <[EMAIL PROTECTED]> wrote: > > On Mon, Sep 29, 2008 at 11:50 AM, Robert Kern <[EMAIL PROTECTED]> wrote: >> >> On Mon, Sep 29, 2008 at 12:40, Charles R Harris >> <[EMAIL PROTECTED]> wrote: >> > >> > On Mon, Sep 29, 2008 at 11:35 AM, Robert Kern <[EMAIL

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 11:50 AM, Robert Kern <[EMAIL PROTECTED]> wrote: > On Mon, Sep 29, 2008 at 12:40, Charles R Harris > <[EMAIL PROTECTED]> wrote: > > > > On Mon, Sep 29, 2008 at 11:35 AM, Robert Kern <[EMAIL PROTECTED]> > wrote: > >> > >> On Mon, Sep 29, 2008 at 12:30, Charles R Harris > >>

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Robert Kern
On Mon, Sep 29, 2008 at 12:40, Charles R Harris <[EMAIL PROTECTED]> wrote: > > On Mon, Sep 29, 2008 at 11:35 AM, Robert Kern <[EMAIL PROTECTED]> wrote: >> >> On Mon, Sep 29, 2008 at 12:30, Charles R Harris >> <[EMAIL PROTECTED]> wrote: >> > >> > On Mon, Sep 29, 2008 at 11:27 AM, Robert Kern <[EMAIL

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Pauli Virtanen
Mon, 29 Sep 2008 11:40:26 -0600, Charles R Harris wrote: [clip] > Let's put it this way, if the tests were consistent with assert then > they couldn't be run with the -OO option. If they *can* be run with the > -OO option then they had better run correctly. Perhaps the Nose author would have input

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 11:35 AM, Robert Kern <[EMAIL PROTECTED]> wrote: > On Mon, Sep 29, 2008 at 12:30, Charles R Harris > <[EMAIL PROTECTED]> wrote: > > > > On Mon, Sep 29, 2008 at 11:27 AM, Robert Kern <[EMAIL PROTECTED]> > wrote: > >> > >> On Mon, Sep 29, 2008 at 12:23, Charles R Harris > >>

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Robert Kern
On Mon, Sep 29, 2008 at 12:30, Charles R Harris <[EMAIL PROTECTED]> wrote: > > On Mon, Sep 29, 2008 at 11:27 AM, Robert Kern <[EMAIL PROTECTED]> wrote: >> >> On Mon, Sep 29, 2008 at 12:23, Charles R Harris >> <[EMAIL PROTECTED]> wrote: >> > >> > >> > On Mon, Sep 29, 2008 at 11:14 AM, Pauli Virtanen

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 11:27 AM, Robert Kern <[EMAIL PROTECTED]> wrote: > On Mon, Sep 29, 2008 at 12:23, Charles R Harris > <[EMAIL PROTECTED]> wrote: > > > > > > On Mon, Sep 29, 2008 at 11:14 AM, Pauli Virtanen <[EMAIL PROTECTED]> wrote: > >> > >> Mon, 29 Sep 2008 11:05:09 -0600, Charles R Harri

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Robert Kern
On Mon, Sep 29, 2008 at 12:23, Charles R Harris <[EMAIL PROTECTED]> wrote: > > > On Mon, Sep 29, 2008 at 11:14 AM, Pauli Virtanen <[EMAIL PROTECTED]> wrote: >> >> Mon, 29 Sep 2008 11:05:09 -0600, Charles R Harris wrote: >> [clip: assert in tests] >> > Using assert for this sort of thing is a bug, a

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 11:14 AM, Pauli Virtanen <[EMAIL PROTECTED]> wrote: > Mon, 29 Sep 2008 11:05:09 -0600, Charles R Harris wrote: > [clip: assert in tests] > > Using assert for this sort of thing is a bug, assert goes away with the > > -OO options. This needs to be fixed. > > "assert FOO" is

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 11:09 AM, Robert Kern <[EMAIL PROTECTED]> wrote: > On Mon, Sep 29, 2008 at 12:05, Charles R Harris > <[EMAIL PROTECTED]> wrote: > > > > > > On Mon, Sep 29, 2008 at 2:00 AM, Matthieu Brucher > > <[EMAIL PROTECTED]> wrote: > >> > >> I've tested the latest Numpy on a RHEL 4, a

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Pauli Virtanen
Mon, 29 Sep 2008 11:05:09 -0600, Charles R Harris wrote: [clip: assert in tests] > Using assert for this sort of thing is a bug, assert goes away with the > -OO options. This needs to be fixed. "assert FOO" is used in quite many of Numpy's test cases, and it appears Nose endorses using assert (an

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Robert Kern
On Mon, Sep 29, 2008 at 12:05, Charles R Harris <[EMAIL PROTECTED]> wrote: > > > On Mon, Sep 29, 2008 at 2:00 AM, Matthieu Brucher > <[EMAIL PROTECTED]> wrote: >> >> I've tested the latest Numpy on a RHEL 4, and I got this result : >> >> Running unit tests for numpy >> NumPy version 1.2.0 >> NumPy

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Charles R Harris
On Mon, Sep 29, 2008 at 2:00 AM, Matthieu Brucher < [EMAIL PROTECTED]> wrote: > I've tested the latest Numpy on a RHEL 4, and I got this result : > > Running unit tests for numpy > NumPy version 1.2.0 > NumPy is installed in /.../BRUCHER/local/lib/python2.5/site-packages/numpy > Python version 2.5

Re: [Numpy-discussion] numpy installation as non root using site.cfg

2008-09-29 Thread pgmoses
Thank you Robert It is now possible for me to install numpy. I get an error with the defaults libraries, which I have reported in a pending mail to the mailen list. Then I tried to use the intel mkl libraries in this way: [mkl] library_dirs = /opt/intel/mkl70/lib/32 mkl_libs = mkl, guide lapack_l

Re: [Numpy-discussion] profiling line by line

2008-09-29 Thread Arnar Flatberg
On Sat, Sep 27, 2008 at 2:36 AM, Robert Kern <[EMAIL PROTECTED]> wrote: > On Wed, Sep 17, 2008 at 07:12, Arnar Flatberg <[EMAIL PROTECTED]> > wrote: > > > > On Wed, Sep 17, 2008 at 3:56 AM, Robert Kern <[EMAIL PROTECTED]> > wrote: > > There is now! Excellent work! I've had no problem with using

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Pauli Virtanen
Mon, 29 Sep 2008 12:33:32 +0200, Matthieu Brucher wrote: > 2008/9/29 Pauli Virtanen <[EMAIL PROTECTED]>: >> Mon, 29 Sep 2008 12:07:53 +0200, Matthieu Brucher wrote: >> >>> np.log(-1 + 0j) 3.1415926535897931j >>> np.log(-1 - 1e-99j) -3.1415926535897931j >>> np.log(-1 + 1e-99j)

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Matthieu Brucher
2008/9/29 Pauli Virtanen <[EMAIL PROTECTED]>: > Mon, 29 Sep 2008 12:07:53 +0200, Matthieu Brucher wrote: > >> np.log(-1 + 0j) >>> 3.1415926535897931j >> np.log(-1 - 1e-99j) >>> -3.1415926535897931j >> np.log(-1 + 1e-99j) >>> 3.1415926535897931j >>> >>> I'd guess this is typically harmle

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Pauli Virtanen
Mon, 29 Sep 2008 12:07:53 +0200, Matthieu Brucher wrote: > np.log(-1 + 0j) >> 3.1415926535897931j > np.log(-1 - 1e-99j) >> -3.1415926535897931j > np.log(-1 + 1e-99j) >> 3.1415926535897931j >> >> I'd guess this is typically harmless, but may lead to problems if your >> code relies on th

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Matthieu Brucher
np.log(-1 + 0j) > 3.1415926535897931j np.log(-1 - 1e-99j) > -3.1415926535897931j np.log(-1 + 1e-99j) > 3.1415926535897931j > > I'd guess this is typically harmless, but may lead to problems if your > code relies on the choice of branch at the branch cut. > > Could you check if this i

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Pauli Virtanen
Mon, 29 Sep 2008 17:19:10 +0900, David Cournapeau wrote: > On Mon, 2008-09-29 at 10:00 +0200, Matthieu Brucher wrote: >> I've tested the latest Numpy on a RHEL 4, and I got this result : >> >> > Hi Matthieu, > > Are you on 32 or 64 bits ? As you can see, those errors are test > errors, no

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Matthieu Brucher
2008/9/29 David Cournapeau <[EMAIL PROTECTED]>: > On Mon, 2008-09-29 at 10:21 +0200, Matthieu Brucher wrote: >> >> Hi, I compiled Python with a 64bits compiler (icc 10.1.018) as well as >> Numpy (that reminds me that I should link against the MKL as well). I >> know it's only tolerance errors, but

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread David Cournapeau
On Mon, 2008-09-29 at 10:21 +0200, Matthieu Brucher wrote: > > Hi, I compiled Python with a 64bits compiler (icc 10.1.018) as well as > Numpy (that reminds me that I should link against the MKL as well). I > know it's only tolerance errors, but they still show up ;) Did you compile with the -fast

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Matthieu Brucher
2008/9/29 David Cournapeau <[EMAIL PROTECTED]>: > On Mon, 2008-09-29 at 10:00 +0200, Matthieu Brucher wrote: >> I've tested the latest Numpy on a RHEL 4, and I got this result : >> > > Hi Matthieu, > >Are you on 32 or 64 bits ? As you can see, those errors are test > errors, not errors in t

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread David Cournapeau
On Mon, 2008-09-29 at 10:00 +0200, Matthieu Brucher wrote: > I've tested the latest Numpy on a RHEL 4, and I got this result : > Hi Matthieu, Are you on 32 or 64 bits ? As you can see, those errors are test errors, not errors in the function themselves (tolerance too low; maybe influence

Re: [Numpy-discussion] ANN: NumPy 1.2.0

2008-09-29 Thread Matthieu Brucher
I've tested the latest Numpy on a RHEL 4, and I got this result : Running unit tests for numpy NumPy version 1.2.0 NumPy is installed in /.../BRUCHER/local/lib/python2.5/site-packages/numpy Python version 2.5.2 (r252:60911, Sep 29 2008, 09:47:20) [GCC Intel(R) C++ gcc 3.4 mode] nose version 0.10.3