[SciPy-User] Float32 FFT

Jochen Schroeder cycomanic@gmail....
Fri Feb 5 23:11:31 CST 2010


On 02/03/10 02:43, David Baddeley wrote:
> Thanks David & Sebastian,
> 
> will look into the scipy trunk & priithon options. Alternatively might find the motivation to dust off some old fftw3 based code I've got floating round which should in theory be a little faster (I'm doing multiple FFTs of the same size, with the potential to do the memory allocation and plan generation just once at the start).
> 
> cheers,
> David
> 

Hi David,

I've written a python wrapper around fftw3 in ctypes. It has quite a
pythonish interface around the way of doing things in fftw3 (using plans).
I'm planning to rewrite the detection of the c-library, before I do a release
but you can find the code at http://launchpad.net/pyfftw

Cheers
Jochen


> 
> ----- Original Message ----
> From: Sebastian Haase <seb.haase@gmail.com>
> To: David Baddeley <david_baddeley@yahoo.com.au>; SciPy Users List <scipy-user@scipy.org>
> Sent: Wed, 3 February, 2010 9:42:54 PM
> Subject: Re: [SciPy-User] Float32 FFT
> 
> On Wed, Feb 3, 2010 at 3:32 AM, David Baddeley
> <david_baddeley@yahoo.com.au> wrote:
> > Does anyone know if there is a way of doing single precision ffts (more explicitly n-dimensional ffts using fftn) in scipy/numpy? The standard functions all seem to force a conversion to double precision and with large (~1000x1000x50) images the memory starts to add up.
> >
> > I'm running numpy 1.2.1, scipy 0.7.0, python 2.6 on linux x64.
> >
> > Thanks,
> > David
> >
> I have an fftw2 wrapper (SWIGged) in Priithon (see googlecode).
> It is based on Python 2.5.  A kind of "super package" binary version
> (also for 64 bit Linux) is downloadable.
> (FFTW is GPL). (The wrapper could be BSD).
> 
> Regards,
> Sebastian Haase
> 
> 
> 
>       
> _______________________________________________
> SciPy-User mailing list
> SciPy-User@scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user


More information about the SciPy-User mailing list