| 查看: 2324 | 回复: 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 ] |
回复此楼» 猜你喜欢
化工调剂303分，过四级
已经有5人回复
-
071000生物学调剂
已经有3人回复
-
一志愿华东理工大学，080500学硕，317分，求调剂
已经有10人回复
-
327求调剂
已经有9人回复
-
0856初试324分求调剂
已经有6人回复
-
08开头看过来！！！
已经有6人回复
-
材料调剂
已经有7人回复
-
一志愿北京科技大学085601材料工程英一数二初试总分335求调剂
已经有8人回复
-
085900土木水利336分求调剂
已经有6人回复
-
319求调剂
已经有13人回复
Fermi
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 3974.4
- 散金: 105
- 红花: 1
- 帖子: 351
- 在线: 90.6小时
- 虫号: 166962
- 注册: 2006-01-14
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2010-11-08 20:26:34
