登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
PROGRAM MAIN
include 'fftw3.f'
INTEGER I,J,M,N,M1,N1,M2,N2,PLAN
PARAMETER(M=128,N=M,M2=M*2,N2=N*2)
REAL*8 PI,HP,PH,F_ETA0,RPM,R1,R2,Rm_PAD
REAL*8 DX,DY,XL,XR,YL,YR,D_PAD,D_PAD1,E_PAD
REAL*8 E_RUNNER,EM_RUNNER,E_PTFE,E_STEEL,EM_PTFE,EM_STEEL
REAL*8 CE(0:M,0:N)
REAL*8, ALLOCATABLE::D0(:,:), ZZ(:,:)
DOUBLE PRECISION, ALLOCATABLE::Z(:,:)
DOUBLE COMPLEX ZFFT(M2,N2)
。。。。。。
CALL dfftw_plan_dft_r2c_2d(PLAN,M2,N2,Z,ZFFT,FFTW_FORWARD,FFTW_ESTIMATE)
。。。。。
DEALLOCATE (D0, ZZ, Z)
END PROGRAM MAIN
win7 64位下的IntelFortran程序。
输入数组z(实双精度型)没有问题,但变换结果zfft都是零。哪位大牛帮看看,问题出在哪了。谢谢!
|