| 查看: 2323 | 回复: 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 ] |
回复此楼» 猜你喜欢
材料化工340求调剂
已经有4人回复
-
一志愿厦门大学材料工程专硕354找调剂！！！
已经有6人回复
-
271求调剂
已经有20人回复
-
一志愿085600中科院宁波所276分求调剂
已经有24人回复
-
一志愿北交大材料工程总分358
已经有4人回复
-
0856材料与化工调剂，339
已经有11人回复
-
一志愿郑大材料工程290求调剂
已经有19人回复
-
266分，一志愿电气工程，本科材料，求材料专业调剂
已经有9人回复
-
372分材料与化工（085600）一志愿湖南大学求调剂
已经有3人回复
-
一志愿南师大0703化学 275求调剂
已经有6人回复
Fermi
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 3974.4
- 散金: 105
- 红花: 1
- 帖子: 351
- 在线: 90.6小时
- 虫号: 166962
- 注册: 2006-01-14
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2010-11-08 20:26:34
