Re: [Numpy-discussion] Bug in frompyfunc starting at 10000 elements?

2010-05-21 Thread David Warde-Farley
Confirmed in NumPy 1.4.1, Py 2.6.5. David On Fri, 21 May 2010, James Bergstra wrote: > Hi all, I'm wondering if this is a bug... > > Something strange happens with my ufunc as soon as I use 1 elements. As > the test shows, the ufunc computes the correct result for either the first > or last

[Numpy-discussion] Bug in frompyfunc starting at 10000 elements?

2010-05-21 Thread James Bergstra
Hi all, I'm wondering if this is a bug... Something strange happens with my ufunc as soon as I use 1 elements. As the test shows, the ufunc computes the correct result for either the first or last elements, but both at the same time is no good. Turns out I'm only running numpy 1.3.0 with

Re: [Numpy-discussion] Summation of large float32/float64 arrays

2010-05-21 Thread Christoph Gohlke
On 5/21/2010 2:30 PM, Alan G Isaac wrote: > On 5/21/2010 4:13 PM, Matthew Turk wrote: >> a1 = numpy.random.random((512,512,512)).astype("float32") > > > This consistently gives me a "MemoryError". > I believe I have plenty of physical memory. > (That should require about 1.3G during creation, rig

Re: [Numpy-discussion] Summation of large float32/float64 arrays

2010-05-21 Thread Alan G Isaac
On 5/21/2010 4:13 PM, Matthew Turk wrote: > a1 = numpy.random.random((512,512,512)).astype("float32") This consistently gives me a "MemoryError". I believe I have plenty of physical memory. (That should require about 1.3G during creation, right? I have 8G.) It seems I'm hitting some kind of 1G me

Re: [Numpy-discussion] Summation of large float32/float64 arrays

2010-05-21 Thread Matthew Turk
Hi Robert, > It's not quite an overflow. > > In [1]: from numpy import * > > In [2]: x = float32(16777216.0) > > In [3]: x + float32(0.9) > Out[3]: 16777216.0 > > You are accumulating your result in a float32. With the a.sum() > approach, you eventually hit a level where the next number to add is

Re: [Numpy-discussion] Summation of large float32/float64 arrays

2010-05-21 Thread Robert Kern
On Fri, May 21, 2010 at 15:13, Matthew Turk wrote: > Hi all, > > I have a possibly naive question.  I don't really understand this > particular set of output: > > In [1]: import numpy > > In [2]: a1 = numpy.random.random((512,512,512)).astype("float32") > > In [3]: a1.sum(axis=0).sum(axis=0).sum(a

[Numpy-discussion] Summation of large float32/float64 arrays

2010-05-21 Thread Matthew Turk
Hi all, I have a possibly naive question. I don't really understand this particular set of output: In [1]: import numpy In [2]: a1 = numpy.random.random((512,512,512)).astype("float32") In [3]: a1.sum(axis=0).sum(axis=0).sum(axis=0) Out[3]: 67110312.0 In [4]: a1.sum() Out[4]: 16777216.0 I re

[Numpy-discussion] calling C function from Python via f2py

2010-05-21 Thread Matt Fearon
Hello, I am trying to use f2py to generate a wrapped C function that I can call from Python (passing arguments to and from). I have this almost working, but I receive trouble with "exp and pow" related to C and some "pos (2) error" with one of my passed variables. My f2py syntax is: f2py -c -lm F

Re: [Numpy-discussion] Build failure at rev8246

2010-05-21 Thread T J
On Fri, May 21, 2010 at 8:51 AM, Pauli Virtanen wrote: > Fri, 21 May 2010 08:09:55 -0700, T J wrote: >> I tried upgrading today and had trouble building numpy (after rm -rf >> build).  My full build log is here: >> >>     http://www.filedump.net/dumped/build1274454454.txt >> > > Your SVN checkout

Re: [Numpy-discussion] Build failure at rev8246

2010-05-21 Thread Pauli Virtanen
Fri, 21 May 2010 08:09:55 -0700, T J wrote: > I tried upgrading today and had trouble building numpy (after rm -rf > build). My full build log is here: > > http://www.filedump.net/dumped/build1274454454.txt > Your SVN checkout might be corrupted, containing a mix of old and new files. Try

[Numpy-discussion] Build failure at rev8246

2010-05-21 Thread T J
Hi, I tried upgrading today and had trouble building numpy (after rm -rf build). My full build log is here: http://www.filedump.net/dumped/build1274454454.txt If someone can point me in the right direction, I'd appreciate it very much. To excerpts from the log file: Running from numpy so