[SciPy-User] Float32 FFT

[David Baddeley]

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


----- Original Message ----
From: Sebastian Haase <seb.haase at gmail.com>
To: David Baddeley <david_baddeley at yahoo.com.au>; SciPy Users List <scipy-user at 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 at 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