[SciPy-dev] Cleaning and fixing fft in scipy ?

[David Cournapeau]

Hi,

    When fixing the fft problem with prime number in numpy, I came again 
across the fft implementation in scipy, and related problems (suboptimal 
solution for fftw, see eg scipy ticket #1). I would like to improve the 
situation:
    - First, I think the module needs serious cleaning, as for now, it 
is a bunch of C files with many #define all accross the code, making it 
difficult to track things. I propose to split the sources by 
implementation (fft_fftw3.c, fft_fft2, fft_mkl, etc...); this can be 
done without any consequence on the implementation.
    - Then, improving fft where it is needed.

    Does this sound ok with scipy developers ?

    I should come up with a patch for zfft within today; even if I try 
to keep changes minimal, I would need people willing to test (I can 
easily test on linux with fftw2 and 3, and if necessary with mkl on the 
same platform. I would prefer avoiding testing on windows myself, as I 
have only a Japanese windows, making things extremely painful for me on 
this already painful platform).

    cheers,

    David