[Numpy-discussion] Enhancements for NumPy's FFTs

Charles R Harris charlesr.harris@gmail....
Sat Mar 14 20:36:07 CDT 2009


On Sat, Mar 14, 2009 at 7:02 PM, Sturla Molden <sturla@molden.no> wrote:

> > On Sat, Mar 14, 2009 at 3:58 PM, Sturla Molden <sturla@molden.no> wrote:
>
> > There is also a ticket (#579) to add an implementation of the Bluestein
> > algorithm for doing prime order fft's.  This could also be used for zoom
> > type fft's. There is lots of fft stuff to be done. I wonder if some of it
> > shouldn't go in Scipy? I think David added some dcts to Scipy.
>
> I am not changing or adding algorithms for now. This is just to prevent
> NumPy from locking up the interpreter while doing FFTs.
>

Well, I was hoping to get you sucked into doing some work here ;)


>
> The loops that are worth multithreading are done in C in
> fftpack_litemodule.c, not in Python in fftpack.py. I have added OpenMP
> pragmas around them. When NumPy gets a build process that supports OpenMP,
> they will execute in parallel. On GCC 4.4 is means compiling with -fopenmp
> and linking -lgomp -lpthread (that goes for mingw/cygwin as well).
>
> The init function seems to be thread safe. cffti and rffti work on arrays
> created in the callers (fftpack_cffti and fftpack_rffti), no global
> objects are touched.
>
> I'm attaching a version of fftpack_litemodule.c that fixes most of what I
> mentioned.
>

Can you put it somewhere for review? I don't think this should go into 1.3
at this late date but 1.4 is a good chance. Hopefully we will get the next
release out a bit faster than this one.

Chuck
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/numpy-discussion/attachments/20090314/244e33e6/attachment.html 


More information about the Numpy-discussion mailing list