Other Parts Discussed in Thread: FFTLIB
Hi,
I'm evaluating the FFTLIB FFT 1D complrex-to-complex basic routine (no OMP, single core), but re result of the computation are wrong.
I execute the test side-by-side with DSPLIB FFT (DSPF_sp_fftSPxSP): with same input (including twiddle) the DSPLIB results are correct while the FFLIB are not correct (I also check the results against Python numpy).
For instance, for N=16, I invoke FFLIB as:
_fft_rad_2p2_next_sm1_notwid_br(in, in, 1, svg, 16, wd, out)
Where svg in a static int[64], wd is float[32] (twiddle), in and out are float[16][2] (cast to double).
The 1st and 9th points are correct (out[0] and out[7]) while all the other are wrong.
As test vector I use a simple ramp over real part:
for(int i=0; i<16; ++i) { in[i][0]=i; in[i][1]=1;}
I test from 16 to 4096 points, with the same problems (first point is always good).
Any advice?
Thanks