[Numpy-discussion] strange problem compiling against ndarrayobject.h

2007-10-31 Thread Vincent Schut
Hmm, it seems my original message did not come through? Not in gmane, at least... Well, here's again: Hi numpy and/or gdal guru's, I'm suddenly getting into trouble compiling gdal's python extension, when it includes ndarrayobject.h from numpy. First it complains about my python not being unicode

[Numpy-discussion] how to get/set column

2007-10-31 Thread dev new
is there a method for numpy arrays and matrices to get/set a particular column i know that a row can be fetched by mymat[1,:] etc can this be done for a column dn ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mail

Re: [Numpy-discussion] numpy FFT memory accumulation

2007-10-31 Thread Anne Archibald
On 31/10/2007, Ray S <[EMAIL PROTECTED]> wrote: > I am using > fftRes = abs(fft.rfft(data_array[end-2**15:end])) > to do running analysis on streaming data. The N never changes. > It sucks memory up at ~1MB/sec with 70kHz data rate and 290 ffts/sec. > (Interestingly, Numeric FFT accumulates much sl

Re: [Numpy-discussion] numpy FFT memory accumulation

2007-10-31 Thread Travis E. Oliphant
Ray S wrote: > I am using > fftRes = abs(fft.rfft(data_array[end-2**15:end])) > to do running analysis on streaming data. The N never changes. > It sucks memory up at ~1MB/sec with 70kHz data rate and 290 ffts/sec. > (Interestingly, Numeric FFT accumulates much slower..) > (Commenting out that one

[Numpy-discussion] numpy FFT memory accumulation

2007-10-31 Thread Ray S
I am using fftRes = abs(fft.rfft(data_array[end-2**15:end])) to do running analysis on streaming data. The N never changes. It sucks memory up at ~1MB/sec with 70kHz data rate and 290 ffts/sec. (Interestingly, Numeric FFT accumulates much slower..) (Commenting out that one line stops memory growth.

Re: [Numpy-discussion] vectorizing loops

2007-10-31 Thread Timothy Hochberg
On Oct 31, 2007 3:18 AM, Francesc Altet <[EMAIL PROTECTED]> wrote: [SNIP] > Incidentally, all the improvements of the PyTables flavor of numexpr > have been reported to the original authors, but, for the sake of > keeping numexpr simple, they decided to implement only some of them. > However, peo

Re: [Numpy-discussion] matrix problem: float to matrix power

2007-10-31 Thread Robert Kern
Charles R Harris wrote: > On 10/31/07, Robert Kern <[EMAIL PROTECTED]> wrote: >> Charles R Harris wrote: >>> On 10/31/07, Alan G Isaac <[EMAIL PROTECTED]> wrote: >>> 1.0**numpy.array([1,2,3]) array([ 1., 1., 1.]) >>> 1.0**numpy.mat([1,2,3]) Traceback (most recent call last): >>>

Re: [Numpy-discussion] matrix problem: float to matrix power

2007-10-31 Thread Charles R Harris
On 10/31/07, Robert Kern <[EMAIL PROTECTED]> wrote: > Charles R Harris wrote: > > On 10/31/07, Alan G Isaac <[EMAIL PROTECTED]> wrote: > > 1.0**numpy.array([1,2,3]) > >> array([ 1., 1., 1.]) > > 1.0**numpy.mat([1,2,3]) > >> Traceback (most recent call last): > >> File "", line 1, in >

Re: [Numpy-discussion] matrix problem: float to matrix power

2007-10-31 Thread Robert Kern
Charles R Harris wrote: > On 10/31/07, Alan G Isaac <[EMAIL PROTECTED]> wrote: > 1.0**numpy.array([1,2,3]) >> array([ 1., 1., 1.]) > 1.0**numpy.mat([1,2,3]) >> Traceback (most recent call last): >> File "", line 1, in >> TypeError: unsupported operand type(s) for ** or pow(): 'float' a

Re: [Numpy-discussion] matrix problem: float to matrix power

2007-10-31 Thread Timothy Hochberg
On Oct 30, 2007 11:40 PM, Alan G Isaac <[EMAIL PROTECTED]> wrote: > >>> 1.0**numpy.array([1,2,3]) > array([ 1., 1., 1.]) > >>> 1.0**numpy.mat([1,2,3]) > Traceback (most recent call last): > File "", line 1, in > TypeError: unsupported operand type(s) for ** or pow(): 'float' and 'matrix' > > W

[Numpy-discussion] strange problem compiling against ndarrayobject.h

2007-10-31 Thread Vincent Schut
Btw forgot some info: gentoo linux, amd64 numpy 1.0.3.1 gdal from svn, updated today python 2.5.1 Cheers, Vincent. ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion

Re: [Numpy-discussion] vectorizing loops

2007-10-31 Thread Francesc Altet
A Wednesday 31 October 2007, Mathew Yeates escrigué: [...] > I also took a look at NumExpr. While it wasn't something I needed for > vectorizing, it still looks very interesting. What kinds of > performance improvements would be expected using this? Well, for some speedup figures you can always ch