re: fftw3 positive, real output in dB SPL
I tried the r2c approach and did 20 * log10(fabs(out[i])), and seems to be coming out more of what I would expect to see now if this is what is supposed to happen.
It could probably use a window, and moving on to see if I can figure that out now.
out = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * (N/2)+1);
p= fftw_plan_dft_r2c_1d(N,input,out, FFTW_ESTIMATE);
outf1 = fopen("output.dat", "w");
fprintf(outf1, "%d\t%5.6f\t%5.6f \n", i, 20 * log10(fabs(out[i])), 20 * log10(fabs(out[i])));
The New Busy is not the old busy. Search, chat and e-mail from your inbox.