Use pointers instead of reshape
Use subroutine
Get rid of spread
Avoid copying p=tmp
Do the loops explicitly
Use slow DFT for N <= 4
Merge pull request #20 from certik/fft2

Speedup the FFT implementation