[SciPy-Dev] Accuracy of single-precision FFT

Sturla Molden sturla at molden.no
Fri Jun 25 08:20:32 EDT 2010


Den 25.06.2010 03:47, skrev David:
> Yep, the bluestein algorithm. I thought fftpack in scipy did this, but
> it actually does not (I tested size 4096 vs 4093, one is 100 times
> slower than the other).
>    

Bluestein's algorithm is just a few lines of Python code using the FFT 
we already have.

(It requires that we factor N into primes though, but we can memoize 
with a dict.)

Sturla



More information about the SciPy-Dev mailing list