Re: [Numpy-discussion] newbie question about boolean testing of array equality result

2011-03-23 Thread Jonathan Hartley
Thanks for the responses, folks. Robert, thanks for helping me come to terms with the situation, for pointing out the perils of comparing NaNs, and particularly for the "not isinstance(x, bool)" idea - some variant of this or attribute checking as you suggest sounds like what we'll use. Daniel

Re: [Numpy-discussion] ANN: Numpy 1.6.0 beta 1

2011-03-23 Thread Derek Homeier
On 24 Mar 2011, at 00:34, Derek Homeier wrote: > tests with the fink-installed pythons on MacOS X mostly succeeded, > with one failure in python2.4 and a couple of issues seemingly > related to PPC floating point accuracy, as below: > Probably last update for tonight: with the 'full' test suite,

Re: [Numpy-discussion] ANN: Numpy 1.6.0 beta 1

2011-03-23 Thread Derek Homeier
Hi again, On 23 Mar 2011, at 17:07, Ralf Gommers wrote: > I am pleased to announce the availability of the first beta of NumPy > 1.6.0. Due to the extensive changes in the Numpy core for this > release, the beta testing phase will last at least one month. Please > test this beta and report any pr

Re: [Numpy-discussion] [SciPy-Dev] ANN: Numpy 1.6.0 beta 1

2011-03-23 Thread Derek Homeier
Hi Bruce, > Sorry as I should I have paid more attention to you first email as I > raised this on the numpy list. (At least with gmail, the reply is to > the > scipy-dev list rather than the numpy list.) > > Based on Pauli's answer I got around it by commenting out a line in > the > setup.py f

Re: [Numpy-discussion] Make array uncopyable

2011-03-23 Thread Dag Sverre Seljebotn
On 03/23/2011 09:24 PM, Daniel Lepage wrote: > On Wed, Mar 23, 2011 at 3:54 PM, Dag Sverre Seljebotn > wrote: >> On 03/23/2011 08:05 PM, Daniel Lepage wrote: >>> Hi all, >>> Is there a way to mark an array as uncopyable? If you have very >>> large data arrays some numpy functions will cause

Re: [Numpy-discussion] [SciPy-Dev] ANN: Numpy 1.6.0 beta 1

2011-03-23 Thread Ralf Gommers
On Wed, Mar 23, 2011 at 8:32 PM, Pauli Virtanen wrote: > On Wed, 23 Mar 2011 14:17:20 -0500, Bruce Southey wrote: >> I can not figure out what is different between git and this version :-( >> >> All the paths appear to be the same. >> Further it continues onwards when I do: $python3.1 -m pdb setup

Re: [Numpy-discussion] Make array uncopyable

2011-03-23 Thread Daniel Lepage
On Wed, Mar 23, 2011 at 3:54 PM, Dag Sverre Seljebotn wrote: > On 03/23/2011 08:05 PM, Daniel Lepage wrote: >> Hi all, >>     Is there a way to mark an array as uncopyable? If you have very >> large data arrays some numpy functions will cause MemoryErrors because >> behind the scenes they e.g. tra

Re: [Numpy-discussion] np.histogramdd of empty data

2011-03-23 Thread Ralf Gommers
On Tue, Mar 22, 2011 at 11:28 AM, Nils Becker wrote: > Hi, > > I was wondering why histogram2d and histogramdd raise a ValueError when > fed with empty data of the correct dimensions. I came across this as a > corner case when calling histogram2d from my own specialized histogram > function. Many

Re: [Numpy-discussion] Make array uncopyable

2011-03-23 Thread Dag Sverre Seljebotn
On 03/23/2011 08:05 PM, Daniel Lepage wrote: > Hi all, > Is there a way to mark an array as uncopyable? If you have very > large data arrays some numpy functions will cause MemoryErrors because > behind the scenes they e.g. transpose the arrays into Fortran order to > feed them into Fortran lib

Re: [Numpy-discussion] [SciPy-Dev] ANN: Numpy 1.6.0 beta 1

2011-03-23 Thread Bruce Southey
On 03/23/2011 02:32 PM, Pauli Virtanen wrote: > On Wed, 23 Mar 2011 14:17:20 -0500, Bruce Southey wrote: >> I can not figure out what is different between git and this version :-( >> >> All the paths appear to be the same. >> Further it continues onwards when I do: $python3.1 -m pdb setup.py build

Re: [Numpy-discussion] numpy log2 has bug

2011-03-23 Thread Robert Kern
On Wed, Mar 23, 2011 at 13:51, wrote: > 2011/3/23 Dmitrey : > from numpy import log2, __version__ > > log2(2**63) >> Traceback (most recent call >> last): >>   File "", line 1, in >> >> AttributeError: log2 > __version__ >> '2.0.0.dev-1fe8136' >> (doesn't work with 1.3.0 as well)

Re: [Numpy-discussion] [SciPy-Dev] ANN: Numpy 1.6.0 beta 1

2011-03-23 Thread Pauli Virtanen
On Wed, 23 Mar 2011 14:17:20 -0500, Bruce Southey wrote: > I can not figure out what is different between git and this version :-( > > All the paths appear to be the same. > Further it continues onwards when I do: $python3.1 -m pdb setup.py build > > I added these lines to the start of "numpy/com

Re: [Numpy-discussion] [SciPy-Dev] ANN: Numpy 1.6.0 beta 1

2011-03-23 Thread Bruce Southey
On 03/23/2011 12:25 PM, Ralf Gommers wrote: > On Wed, Mar 23, 2011 at 6:00 PM, Bruce Southey wrote: >> On 03/23/2011 11:38 AM, Ralf Gommers wrote: >>> On Wed, Mar 23, 2011 at 5:29 PM, Bruce Southeywrote: On 03/23/2011 11:07 AM, Ralf Gommers wrote: > Hi, > > I am pleased to ann

[Numpy-discussion] Make array uncopyable

2011-03-23 Thread Daniel Lepage
Hi all, Is there a way to mark an array as uncopyable? If you have very large data arrays some numpy functions will cause MemoryErrors because behind the scenes they e.g. transpose the arrays into Fortran order to feed them into Fortran libraries (numpy.linalg.svd is an example of this). It woul

Re: [Numpy-discussion] numpy log2 has bug

2011-03-23 Thread josef . pktd
2011/3/23 Dmitrey : from numpy import log2, __version__ log2(2**63) > Traceback (most recent call > last): >   File "", line 1, in > > AttributeError: log2 __version__ > '2.0.0.dev-1fe8136' > (doesn't work with 1.3.0 as well) >>> np.array([2**63]) array([9223372036854775808],

Re: [Numpy-discussion] bug with numpy 2 ** N

2011-03-23 Thread Matthieu Brucher
Hi, I don't thnk this is a bug. You are playign with C integers, not Python integers, and the former are limited. It's a common "feature" in all processors (even DSPs). Matthieu 2011/3/23 Dmitrey > >>> 2**64 > 18446744073709551616L > >>> 2**array(64) > -9223372036854775808 > >>> 2**100 > 1267

[Numpy-discussion] bug with numpy 2 ** N

2011-03-23 Thread Dmitrey
>>> 2**64 18446744073709551616L >>> 2**array(64) -9223372036854775808 >>> 2**100 1267650600228229401496703205376L >>> 2**array(100) -9223372036854775808 ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.o

[Numpy-discussion] numpy log2 has bug

2011-03-23 Thread Dmitrey
>>> from numpy import log2, __version__ >>> log2(2**63) Traceback (most recent call last):

[Numpy-discussion] Scientific Software developer wanted in Seattle.

2011-03-23 Thread Christopher Barker
Scientific Software Developer NOAA Emergency Response Division Help us develop our next-generation oil spill transport model. Background: The Emergency Response Division (ERD) of NOAA's Office of Response and Restoration (OR&R) provides scientific expertise to support the response to oil and ch

Re: [Numpy-discussion] [SciPy-Dev] ANN: Numpy 1.6.0 beta 1

2011-03-23 Thread Ralf Gommers
On Wed, Mar 23, 2011 at 6:00 PM, Bruce Southey wrote: > On 03/23/2011 11:38 AM, Ralf Gommers wrote: >> >> On Wed, Mar 23, 2011 at 5:29 PM, Bruce Southey  wrote: >>> >>> On 03/23/2011 11:07 AM, Ralf Gommers wrote: Hi, I am pleased to announce the availability of the first beta o

Re: [Numpy-discussion] [SciPy-Dev] ANN: Numpy 1.6.0 beta 1

2011-03-23 Thread Ralf Gommers
On Wed, Mar 23, 2011 at 5:29 PM, Bruce Southey wrote: > On 03/23/2011 11:07 AM, Ralf Gommers wrote: >> Hi, >> >> I am pleased to announce the availability of the first beta of NumPy >> 1.6.0. Due to the extensive changes in the Numpy core for this >> release, the beta testing phase will last at le

Re: [Numpy-discussion] [SciPy-Dev] ANN: Numpy 1.6.0 beta 1

2011-03-23 Thread Bruce Southey
On 03/23/2011 11:07 AM, Ralf Gommers wrote: > Hi, > > I am pleased to announce the availability of the first beta of NumPy > 1.6.0. Due to the extensive changes in the Numpy core for this > release, the beta testing phase will last at least one month. Please > test this beta and report any problems

Re: [Numpy-discussion] newbie question about boolean testing of array equality result

2011-03-23 Thread Daniel Lepage
On Wed, Mar 23, 2011 at 10:29 AM, Jonathan Hartley wrote: > Hey people, > > I'm writing an application in which we evaluate user-supplied Python > expressions. > > Sometimes, we want to do an equality test on the result of the evaluation, > eg: > >     result = eval(...) >     if result == float('

[Numpy-discussion] ANN: Numpy 1.6.0 beta 1

2011-03-23 Thread Ralf Gommers
Hi, I am pleased to announce the availability of the first beta of NumPy 1.6.0. Due to the extensive changes in the Numpy core for this release, the beta testing phase will last at least one month. Please test this beta and report any problems on the Numpy mailing list. Sources and binaries can b

Re: [Numpy-discussion] newbie question about boolean testing of array equality result

2011-03-23 Thread Robert Kern
On Wed, Mar 23, 2011 at 09:29, Jonathan Hartley wrote: > Hey people, > > I'm writing an application in which we evaluate user-supplied Python > expressions. > > Sometimes, we want to do an equality test on the result of the evaluation, > eg: > >     result = eval(...) >     if result == float('nan

[Numpy-discussion] newbie question about boolean testing of array equality result

2011-03-23 Thread Jonathan Hartley
Hey people, I'm writing an application in which we evaluate user-supplied Python expressions. Sometimes, we want to do an equality test on the result of the evaluation, eg: result = eval(...) if result == float('nan'): ... If the result is a numpy array, then this raises a

Re: [Numpy-discussion] Is numpy/scipy linux apt or PYPI installation linked with ACML?

2011-03-23 Thread Sturla Molden
ACML has a C interface for BLAS and LAPACK, but it is different from cblas. For portability, stick with the Fortran 77 interface, particularly for LAPACK where no C interface is specified. One can probably get cblas from Netlib and link with libacml. cblas is just a wrapper for the Fortran r

Re: [Numpy-discussion] Numpy: controlling ATLAS multithreading

2011-03-23 Thread Mark Dixon
On Wed, 23 Mar 2011, Ralf Gommers wrote: ... >> Does Numpy have the capability to choose between the single-threaded and >> multi-threaded ATLAS libraries at runtime? Or do I need to make this >> decision at build time? > > You need to decide at build time. ... Thanks Ralf Mark -- --

Re: [Numpy-discussion] Numpy: controlling ATLAS multithreading

2011-03-23 Thread Ralf Gommers
On Wed, Mar 23, 2011 at 10:11 AM, Mark Dixon wrote: > Hi, > > I've built Numpy 1.5.1 against a copy of ATLAS 3.8.3 and LAPACK 3.3.0, > it's passing all tests and is looking good. Great. > > I note that Numpy preferentially links against the multi-threaded ATLAS > libraries; however, we sometimes w

[Numpy-discussion] SciPy 2011 Call for Papers

2011-03-23 Thread Amenity Applewhite
Hello, SciPy 2011 , the 10th Python in Science conference, will be held July 11 - 16, 2011, in Austin, TX. At this conference, novel applications and breakthroughs made in the pursuit of science using Python are presented. Attended by leading figur

Re: [Numpy-discussion] Appending a numpy array to binary file

2011-03-23 Thread Scott Sinclair
On 22 March 2011 17:22, Robert Kern wrote: > On Tue, Mar 22, 2011 at 05:37, Alessandro > wrote: >> Hi everybody, >> >> I'm trying to append, inside a loop, into a binary file some arrays using >> >> numpy.save(file, arr) but wvhen I open it there is only the last array I >> saved. >> >> If I use

Re: [Numpy-discussion] Problems building NumPy with GotoBLAS

2011-03-23 Thread Sturla Molden
Den 22.03.2011 23:18, skrev Paul Anton Letnes: I'm no expert, but I just pulled off the scipy+numpy+GotoBLAS2 installation. From what I gather, the Makefile for libgoto2 downloads and compiles the generic lapack from netlib. It also wraps lapack into libgoto2.so/.a . I be

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-23 Thread Sturla Molden
I've been busy lately, so I haven't had time to answer. The idea is of course to use Fortran with Python, not to build a framework in Fortran. I also find C useful to interface with the operating system, and Cython to write C extensions for Python. Also, if the computationally demaning parts a

Re: [Numpy-discussion] Fortran was dead ... [was Re: rewriting NumPy code in C or C++ or similar]

2011-03-23 Thread Sturla Molden
Den 16.03.2011 14:46, skrev Neal Becker: > Also: > * can it adopt external memory? Yes. Using CRAY pointers with libc malloc/free is e.g. a common way to get dynamic memory in Fortran 77. > * can it interwork with numpy? (kinda required for this audience) > Yes, that is why NumPy has f2py :-)

[Numpy-discussion] Numpy: controlling ATLAS multithreading

2011-03-23 Thread Mark Dixon
Hi, I've built Numpy 1.5.1 against a copy of ATLAS 3.8.3 and LAPACK 3.3.0, it's passing all tests and is looking good. Great. I note that Numpy preferentially links against the multi-threaded ATLAS libraries; however, we sometimes want to constrain a process to a single thread, and sometimes n

Re: [Numpy-discussion] git tag question/warning

2011-03-23 Thread Ralf Gommers
On Wed, Mar 23, 2011 at 5:34 AM, Charles R Harris wrote: > > > On Tue, Mar 22, 2011 at 11:29 AM, Ralf Gommers > wrote: >> >> Hi, >> >> After tagging 1.6.0b1 and pushing it I found some issues that only showed >> up with released=True in setup.py. According to the docs at >> http://www.kernel.org/