Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

scaling factor for CUFFT

Tags:

c++

cuda

fft

fftw

I wrote a simple test program, where I was doing Complex to Complex FT's and I just generated some data 1..50 and stuck it in for the real and imaginary part for each index of the array.

When I do an operation like this IFFT(FFT(A)) = A

to test them out, I am getting different results for each library.

FFTW, I have to divide the output by len(A) to get back the original A

However, doing this forward then reverse FFT with CUFFT, it looks like I have to divide by (sqrt(2)*50) to get back to the original data.

Where is this extra square root factor coming from?

according to the CUFFT documentation: CUFFT performs un normalized FFTs; that is, performing a forward FFT on an input data set followed by an inverse FFT on the resulting set yields data that is equal to the input scaled by the number of elements. Scaling either transform by the reciprocal of the size of the data set is left for the user to perform as seen fit.

Thanks in Advance

like image 939
Derek Avatar asked Mar 09 '11 21:03

Derek


2 Answers

CUFFT has the same behavior as FFTW, it computes unnormalized FFTs. IFFT(FFT(A))=n A where n is the length of the vector. The length n is in number of samples (not floats or bytes). There are some padding differences between FFTW and CUFFT with C2R and R2C that can screw up a simple comparison, but not for C2C. I would double-check your data setup and length calculations, and verify your plan in both FFTW and CUFFT.

like image 65
Nathan Whitehead Avatar answered Oct 06 '22 18:10

Nathan Whitehead


This ended up being a problem with the way the absolute values of the complex number was being calculated. in the std::complex library, it was computing the distance of the vector.

like image 27
Derek Avatar answered Oct 06 '22 17:10

Derek