Hello all,
I am trying to do 3-D FFT complex to complex. First i converted Q to complex QC then reshape 3D to 1D. Then I did FFT. Could anybody tell, if i did any mistake?
! convert Q to complex QC QC=cmplx(Q,0.d0,16) ! Reshape 3D to 1D QC_1d=reshape(QC,(/nfftdim1*nfftdim2*nfftdim3/)) ! Initializing DFFT stuff FFT_length(1)=nfftdim1 FFT_length(2)=nfftdim2 FFT_length(3)=nfftdim3 ! Initializing Forward DFFT status=DftiCreateDescriptor(forward_desc,DFTI_DOUBLE,DFTI_COMPLEX, 3, FFT_length) status=DftiCommitDescriptor(forward_desc) ! Doing Forward FFT status=DftiComputeForward(forward_desc, QC_1d)