| 查看: 2400 | 回复: 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 ] |
回复此楼» 猜你喜欢
时间戳变了，能看出什么问题？
已经有18人回复
-
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
言之凿凿，8月21号(明天)放榜
已经有14人回复
-
只有每年这种时候来逛逛小木虫
已经有21人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有125人回复
-
今天基金会出结果吗？20260819
已经有15人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
Fermi
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 3979.4
- 散金: 105
- 红花: 1
- 帖子: 352
- 在线: 90.6小时
- 虫号: 166962
- 注册: 2006-01-14
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2010-11-08 20:26:34
