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