program main
!
integer :: N
parameter (N=80)
real in, var
real out
dimension in(N), out(N)
dimension var(N)
integer :: plan,plan2
integer :: i
real :: ran
!
do i=1,N
call random_number(ran)
in(i-1)=ran
end do
!
call rfftw_f77_create_plan(plan,N,FFTW_FORWARD,FFTW_ESTIMATE)
call rfftw_f77_one(plan,in,var)
call rfftw_f77_destroy_plan(plan)
call rfftw_f77_create_plan(plan2,N,FFTW_BACKWARD,FFTW_ESTIMATE)
call rfftw_f77_one(plan2,var,out)
call rfftw_f77_destroy_plan(plan2)
!
open(18,file='testft.dat')
do i=1,N
write(18,*)i,real(in(i)),real(out(i))
end do
close(18)
end program main
楚香饭 发表于 2014-10-28 19:26
使用 fftw2 的话,你得首先编译 wrappers 接口库。
位于 MKL 安装目录下的 interfaces\fftw2xf 文件夹。
11.83 KB, 下载次数: 5
楚香饭 发表于 2014-10-29 16:44
icl 找不到。
icl 是 Intel C++ 编译器,找不到是自然的,你没安装。
11.84 KB, 下载次数: 4
楚香饭 发表于 2014-10-29 18:51
先给你做个single的lib,代码结果不对,稍候在调。
楚香饭 发表于 2014-10-29 18:51
先给你做个single的lib,代码结果不对,稍候在调。
欢迎光临 Fortran Coder (http://bbs.fcode.cn/) | Powered by Discuz! X3.2 |