| 查看: 2299 | 回复: 1 | |||
Fermi木虫 (正式写手)
|
[交流]
【求助】Fortran调用FFTW的问题【已完结】
|
|
最近要用到fftw，按照fftw的说明文档写了下面这个简单的测试程序： --------------------------------------------------------------------------- program fftwtest integer*8,parameter::N=2,M=2 double precision in dimension in(0:M,0:N) double complex out dimension out(0:M/2, 0:N) integer*8 plan,i,j open(10,file='in.txt') do i=0,M do j=0,N read(10,*) i,j,in(i,j) end do end do close(10) open(15,file='in_1.txt') do i=0,M do j=0,N write(15,*) i,j,in(i,j) end do end do close(15) call dfftw_plan_dft_r2c_2d(plan,M,N,in,out,FFTW_FORWARD,FFTW_ESTIMATE) call dfftw_execute_dft_r2c(plan, in, out) call dfftw_destroy_plan(plan) open(20,file='out.txt') do i=0,M/2 do j=0,N write(20,*) i,j,out(i,j) end do end do close(20) print*, 'Program finished' end program —————————————————————————————————————————— 我用pgi编译如下： pgf90 -fast -o fftw fftwtest1.f90 -L/opt/fftw/lib -lfftw3_f77 -lm 编译的时候只是说我在读in（i,j）的时候有两个警告，其他没有问题。编译后的程序能够顺利运行，但是不管输入什么数据，结果都是0，不知道是哪里出了问题？请各位兄弟帮小弟看看。 以下是我的输入文件： 0 0 1 0 1 2 0 2 4 1 0 4 1 1 5 1 2 6 2 0 7 2 1 9 2 2 3 下面是运算结果： 0 0 (0.000000000000000,0.000000000000000) 0 1 (0.000000000000000,0.000000000000000) 0 2 (0.000000000000000,0.000000000000000) 1 0 (0.000000000000000,0.000000000000000) 1 1 (0.000000000000000,0.000000000000000) 1 2 (0.000000000000000,0.000000000000000) [ Last edited by nono2009 on 2010-11-11 at 07:17 ] |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
Fermi
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 3969.4
- 散金: 105
- 红花: 1
- 帖子: 350
- 在线: 90.6小时
- 虫号: 166962
- 注册: 2006-01-14
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2010-11-08 20:26:34