I am applying a complex-to-complex FFT on some data, first ‘FORWARD’, then ‘INVERSE’ via code
CUFFT_CALL(cufftExecZ2Z(dev_mem->nufft_fftplan, dev_mem->d_nufft_I_tau_comp.data(), dev_mem->d_nufft_fftdata_ext.data(), CUFFT_FORWARD));
CUFFT_CALL(cufftExecZ2Z(dev_mem->nufft_fftplan, dev_mem->d_nufft_fftdata_ext.data(), dev_mem->d_nufft_I_tau_comp.data(), CUFFT_INVERSE));
I would have expected, as I do not modify the data inbetween, that after inverse FFT I get the original data again - like in Matlab. But, it is completely different. Do I miss a scaling factor or something like that, any clue from someone ?