Re: [Numpy-discussion] zoom FFT with numpy?

2007-03-15 Thread Stefan van der Walt
Here is another implementation of the Chirp-Z transform in Python, that handles complex numbers. Regards Stéfan On Thu, Mar 15, 2007 at 11:53:49AM +0200, Nadav Horesh wrote: > A long time ago I translated a free code of chirp z transform (zoom fft) into > python. > Attached here the source and tw

Re: [Numpy-discussion] rant against from numpy import * / from pylab import *

2007-03-15 Thread Eric Firing
Sebastian Haase wrote: > Hi! > I use the wxPython PyShell. > I like especially the feature that when typing a module and then the > dot "." I get a popup list of all available functions (names) inside > that module. > > Secondly, I think it really makes code clearer when one can see where > a fun

Re: [Numpy-discussion] nd_image.affine_transform edge effects

2007-03-15 Thread Stefan van der Walt
Hi James On Thu, Mar 15, 2007 at 06:01:55PM -0400, James Turner wrote: > >Code snippets to illustrate the problem would be welcome. > > OK. I have had a go at producing a code snippet. I apologize that > this is based on numarray rather than numpy, since I'm using STScI > Python, but I think it s

Re: [Numpy-discussion] numpy.distutils, windows dll vs lib

2007-03-15 Thread Albert Strasheim
Hello all - Original Message - From: "David Cournapeau" <[EMAIL PROTECTED]> To: "Discussion of Numerical Python" Sent: Friday, March 16, 2007 6:51 AM Subject: Re: [Numpy-discussion] numpy.distutils, windows dll vs lib > Robert Kern wrote: >> David Cournapeau wrote: >> Generally, you

Re: [Numpy-discussion] numpy.distutils, windows dll vs lib

2007-03-15 Thread Robert Kern
David Cournapeau wrote: > I don't use the library for linking: I only need to be able to load it > dynamically through ctypes. What I did is simply overriding the > calc_info method, in which I try to detect both library and header > files. For the library, I do the following: > > # Look for t

Re: [Numpy-discussion] zoom FFT with numpy? (Nadav Horesh)

2007-03-15 Thread Charles R Harris
On 3/15/07, Ray S <[EMAIL PROTECTED]> wrote: "Charles R Harris" wrote: >BTW, did the s/n you gave refer to the signal or to it's transform? That is the time domain s/n, approximately... Ah, good. The s/n for lines spectra in the transform will go up as the sqrt of the number of sample poi

[Numpy-discussion] rant against from numpy import * / from pylab import *

2007-03-15 Thread Sebastian Haase
Hi! I use the wxPython PyShell. I like especially the feature that when typing a module and then the dot "." I get a popup list of all available functions (names) inside that module. Secondly, I think it really makes code clearer when one can see where a function comes from. I have a default imp

Re: [Numpy-discussion] numpy.distutils, windows dll vs lib

2007-03-15 Thread David Cournapeau
Robert Kern wrote: > David Cournapeau wrote: >> Hi, >> >> recently, I got some problems detecting a dynamic library (dll) with >> numpy.distutils. Basically, I have a package which uses a class derived >> from system_info from numpy.distutils to detect a dll I use through ctypes. >> If on

Re: [Numpy-discussion] numpy.distutils, windows dll vs lib

2007-03-15 Thread Robert Kern
David Cournapeau wrote: > Hi, > > recently, I got some problems detecting a dynamic library (dll) with > numpy.distutils. Basically, I have a package which uses a class derived > from system_info from numpy.distutils to detect a dll I use through ctypes. > If only the dll is present, my

[Numpy-discussion] numpy.distutils, windows dll vs lib

2007-03-15 Thread David Cournapeau
Hi, recently, I got some problems detecting a dynamic library (dll) with numpy.distutils. Basically, I have a package which uses a class derived from system_info from numpy.distutils to detect a dll I use through ctypes. If only the dll is present, my numpy.distutils.system_info derived

Re: [Numpy-discussion] zoom FFT with numpy? (Nadav Horesh)

2007-03-15 Thread Ray S
"Charles R Harris" wrote: >BTW, did the s/n you gave refer to the signal or to it's transform? That is the time domain s/n, approximately... Thanks to all for the replies, some digestion is needed. Ray ___ Numpy-discussion mailing list Numpy-discus

Re: [Numpy-discussion] nd_image.affine_transform edge effects

2007-03-15 Thread Robert Kern
James Turner wrote: > Sorry I accidentally broke this thread in 2 (see thread of March 9). > > I tried manually adding the right reference in the mail header to > continue the same thread, but obviously got it wrong (I think because > I replied to the digest email instead of starting a new one). N

[Numpy-discussion] nd_image.affine_transform edge effects

2007-03-15 Thread James Turner
Sorry I accidentally broke this thread in 2 (see thread of March 9). I tried manually adding the right reference in the mail header to continue the same thread, but obviously got it wrong (I think because I replied to the digest email instead of starting a new one). Not sure whether there is a bet

Re: [Numpy-discussion] nd_image.affine_transform edge effects

2007-03-15 Thread James Turner
Hi Stefan, Thanks for the suggestions! > Is this related to > > http://projects.scipy.org/scipy/scipy/ticket/213 > > in any way? As far as I can see, the problems look different, but thanks for the example of how to document this. I did confirm that your example exhibits the same behaviour und

Re: [Numpy-discussion] correct way to specify type in array definition

2007-03-15 Thread David M. Cooke
On Thu, Mar 15, 2007 at 11:23:36AM +0100, Francesc Altet wrote: > El dj 15 de 03 del 2007 a les 06:01 -0400, en/na Brian Blais va > escriure: > > Hello, > > > > Can someone tell me what the preferred way to specify the type of an array? > > I want > > it to be a float array, no matter what is gi

Re: [Numpy-discussion] zoom FFT with numpy?

2007-03-15 Thread Charles R Harris
On 3/15/07, Warren Focke <[EMAIL PROTECTED]> wrote: On Wed, 14 Mar 2007, Charles R Harris wrote: > On 3/14/07, Ray Schumacher <[EMAIL PROTECTED]> wrote: > > > > What I had been doing is a 2048 N full real_FFT with a Hann window, and > > further analyzing the side lobe/bin energy (via linear i

Re: [Numpy-discussion] zoom FFT with numpy? (Nadav Horesh)

2007-03-15 Thread Charles R Harris
On 3/15/07, Ray S <[EMAIL PROTECTED]> wrote: Hi Nadev, >A long time ago I translated a free code of chirp z transform (zoom fft) >into python. Thanks, I'll try it out. I did, however read before on the differences: From Numerix http://www.numerix-dsp.com/zoomfft.html: "One common question is :

Re: [Numpy-discussion] zoom FFT with numpy?

2007-03-15 Thread Warren Focke
On Wed, 14 Mar 2007, Charles R Harris wrote: > On 3/14/07, Ray Schumacher <[EMAIL PROTECTED]> wrote: > > > > What I had been doing is a 2048 N full real_FFT with a Hann window, and > > further analyzing the side lobe/bin energy (via linear interp) to try to > > more precisely determine the f wit

Re: [Numpy-discussion] zoom FFT with numpy? (Nadav Horesh)

2007-03-15 Thread Ray S
Hi Nadev, >A long time ago I translated a free code of chirp z transform (zoom fft) >into python. Thanks, I'll try it out. I did, however read before on the differences: From Numerix http://www.numerix-dsp.com/zoomfft.html: "One common question is : Is the zoom FFT the same as the chirp z-tr

Re: [Numpy-discussion] Problem compiling numpy with Python 2.5 on Powerbook

2007-03-15 Thread Robert Kern
Roberto De Almeida wrote: > Hi, Steve. > > On 3/14/07, Steve Lianoglou <[EMAIL PROTECTED]> wrote: >> I'm not sure what the problem is exactly, but is it weird that >> there's something to do w/ 'i686' when you're running on a powerbook >> being that the pbook is PowerPC? > > I managed to compile

Re: [Numpy-discussion] Problem compiling numpy with Python 2.5 on Powerbook

2007-03-15 Thread Roberto De Almeida
Hi, Steve. On 3/14/07, Steve Lianoglou <[EMAIL PROTECTED]> wrote: > I'm not sure what the problem is exactly, but is it weird that > there's something to do w/ 'i686' when you're running on a powerbook > being that the pbook is PowerPC? I managed to compile numpy by first compiling Python 2.5 as

Re: [Numpy-discussion] arrayrange

2007-03-15 Thread Robert Kern
Miguel Oliveira, Jr. wrote: > Hi Eike! > > Yes, I did subscribe! Thanks once again for your reply... I'm sorry > but I don't know exactly what you mean... Do you think I should replace > t = arange(0, durSecs, 1.0/SRate) > by > t = linspace(0, durSecs, durSecs*SRate)? > > That won't work... *H

Re: [Numpy-discussion] arrayrange

2007-03-15 Thread Robert Kern
Miguel Oliveira, Jr. wrote: > Hi Robert, > > I'm expecting that the code generates a sine wave sound file in aiff > format... Apparently the code worked before when arrayrange() was > used... Then the problem lies somewhere else; a number of other things changed between Numeric and numpy. arr

[Numpy-discussion] Broadcasting for hstack, vstack etc

2007-03-15 Thread Bill Baxter
I just had a need to append a column of 1's to an array, and given how big numpy is on broadcasting I thought this might work: column_stack((m1,m2, 1)) But it doesn't. Is there any reason why that couldn't or shouldn't be made to work? --bb ___ Nump

Re: [Numpy-discussion] Matlab -> NumPy translation and indexing

2007-03-15 Thread David Koch
On 3/15/07, David Koch <[EMAIL PROTECTED]> wrote: ... NumPy equiv for Matlab B(A_idx, A_Idx) Ok, I did like this: A_Idx = array([1, 0]) B = random.randn(3,3) rowInd = kron(ones((1, len(A_Idx)), int), A_Idx[:, newaxis]) colInd = kron(ones((len(A_Idx), 1), int), A_Idx) B[rowInd, colInd] - re

Re: [Numpy-discussion] Matlab -> NumPy translation and indexing

2007-03-15 Thread David Koch
On 3/14/07, Gael Varoquaux <[EMAIL PROTECTED]> wrote: I definitely second this comment. Using arrays when you are trying to append a lot of data is using the wrong data format. And the code is so much more readable with lists. Thank you, I will consider it, Next thing, I have A_Idx = arr

Re: [Numpy-discussion] correct way to specify type in array definition

2007-03-15 Thread Francesc Altet
El dj 15 de 03 del 2007 a les 06:01 -0400, en/na Brian Blais va escriure: > Hello, > > Can someone tell me what the preferred way to specify the type of an array? > I want > it to be a float array, no matter what is given (say, integers). I can do: > > a=numpy.array([1,2,3],numpy.dtype('float'

[Numpy-discussion] Possible bug in PyArray_RemoveSmallest?

2007-03-15 Thread Albert Strasheim
Hello all I was poking around in the NumPy internals and I came across the following code in PyArray_RemoveSmallest in arrayobject.c: intp sumstrides[NPY_MAXDIMS]; ... for (i=0; ind; i++) { sumstrides[i] = 0; for (j=0; jnumiter; j++) { sumstrides[i] = multi->iters[j]->strides[i

[Numpy-discussion] correct way to specify type in array definition

2007-03-15 Thread Brian Blais
Hello, Can someone tell me what the preferred way to specify the type of an array? I want it to be a float array, no matter what is given (say, integers). I can do: a=numpy.array([1,2,3],numpy.dtype('float')) or a=numpy.array([1,2,3],type(1.0)) or perhaps many others. Is there a way that i

Re: [Numpy-discussion] zoom FFT with numpy?

2007-03-15 Thread Nadav Horesh
A long time ago I translated a free code of chirp z transform (zoom fft) into python. Attached here the source and two translations (probably one of the is right) Nadav. On Wed, 2007-03-14 at 14:02 -0800, Ray S wrote: > We'd like to do what most call a "zoom FFT"; we only are interested > in

Re: [Numpy-discussion] arrayrange

2007-03-15 Thread Miguel Oliveira, Jr.
Hi Robert, I'm expecting that the code generates a sine wave sound file in aiff format... Apparently the code worked before when arrayrange() was used... Thanks! M. On 14 Mar 2007, at 17:31, Robert Kern wrote: > Miguel Oliveira, Jr. wrote: >> Hello, >> >> I've got a few codes that use "arr

Re: [Numpy-discussion] arrayrange

2007-03-15 Thread Miguel Oliveira, Jr.
Hi Eike! Yes, I did subscribe! Thanks once again for your reply... I'm sorry but I don't know exactly what you mean... Do you think I should replace t = arange(0, durSecs, 1.0/SRate) by t = linspace(0, durSecs, durSecs*SRate)? That won't work... Maybe I am missing something?... The resulting f

Re: [Numpy-discussion] Nonblocking Plots with Matplotlib

2007-03-15 Thread Andrew Straw
Bill Baxter wrote: > On 3/15/07, Bill Baxter <[EMAIL PROTECTED]> wrote: > >> Thanks, Sebastian. I'll take a look at Pyro. Hadn't heard of it. >> I'm using just xmlrpclib with pickle right now. >> > > I took a look at Pyro -- it looks nice. > The only thing I couldn't find, though, is how

Re: [Numpy-discussion] [Matplotlib-users] Nonblocking Plots with Matplotlib

2007-03-15 Thread Andrew Straw
Bill, very cool. Also, thanks for showing me how Twisted can be used like Pyro, more-or-less, I think. (If I understand your code from my 1 minute perusal.) On Mac OS X, there's one issue I don't have time to follow any further: sys.executable points to /Library/Frameworks/Python.framework/Ve