Fix FFT normalization

This commit is contained in:
Gerald Fournier 2026-02-16 00:02:20 +02:00 committed by dwarning
parent 053363f85d
commit ce3f3c6221
1 changed files with 2 additions and 2 deletions

View File

@ -181,7 +181,7 @@ com_fft(wordlist *wl)
fftw_execute(plan_forward);
scale = (double) fpts - 1.0;
scale = (double) (length*1.0/2.0);
fdvec[i][0].cx_real = out[0][0]/scale/2.0;
fdvec[i][0].cx_imag = 0.0;
for (j = 1; j < fpts; j++) {
@ -215,7 +215,7 @@ com_fft(wordlist *wl)
rffts(in, M, 1);
fftFree();
scale = (double) fpts - 1.0;
scale = (double) (length*1.0/2.0);
/* Re(x[0]), Re(x[N/2]), Re(x[1]), Im(x[1]), Re(x[2]), Im(x[2]), ... Re(x[N/2-1]), Im(x[N/2-1]). */
fdvec[i][0].cx_real = in[0]/scale/2.0;
fdvec[i][0].cx_imag = 0.0;