Quantcast
Channel: Clusters and HPC Technology
Viewing all articles
Browse latest Browse all 927

3-D FFT via 1-D conversion

$
0
0

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)

 


Viewing all articles
Browse latest Browse all 927

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>