Hello Davide !

Our cfft/icfft are there ... so fast.
Just can't yet augment the 2D... shit !

Jean

Fourier al_fftc1d.sm (91,28 KiB) scaricato 1594 volta(e).

Hello Jean

The problem seems that inside the functions stack/augment works in symbolical mode, and the input in symbolical mode is the one below, that is not managed by stack/augment (not sure if correct for symbolical calculations, bug or missing feature).



You have to place an eval() inside fft2D to force a numerical result.

Fourier al_fftc1d_[ALGLIB][DEBUG].sm (101,29 KiB) scaricato 1346 volta(e).
Fourier al_[FFTPACK].sm (52,49 KiB) scaricato 1363 volta(e).

Thanks Davide,

Your coming FFTPACK looks great/effective.
In the mean time, interested users can use
cfft ... spectrum added.

Cheers, Jean

Fourier al_fftc1d_ALGLIB.sm (88,06 KiB) scaricato 1575 volta(e).