Fortran Coder

查看: 325|回复: 3

[数学库] 使用MKL_DFTI进行快速傅里叶变换时,输出复数实虚部皆为0

[复制链接]

5

帖子

1

主题

0

精华

新人

F 币
16 元
贡献
7 点
发表于 2022-12-6 15:56:02 | 显示全部楼层 |阅读模式
5F 币
运行代码如下:
[Fortran] 纯文本查看 复制代码
program test2
  use mkl_dfti
  complex :: x_2D(2,4),y_2D(4,8)
  complex :: x(8),y(32)
  type(DFTI_DESCRIPTOR), POINTER :: my_desc1
  integer :: status,l(2)
  equivalence(x_2D,x)
  equivalence(y_2D,y)
  l(1)=2
  l(2)=4
  x_2D=reshape([1,2,3,4,5,6,7,8],[2,4])
  print *, x
  status = DftiCreateDescriptor(my_desc1,DFTI_DOUBLE,DFTI_COMPLEX,2,l)
  status = DftiCommitDescriptor(my_desc1)
  status = DftiComputeForward(my_desc1,x,y)
  status = DftiFreeDescriptor(my_desc1)
  print *, x
  print *, y
end program test2

运行结果如图:
结果.png
运行网页:如何配置MKL的MKL_DFTI来计算FFT(f90接口)-编程工具交流-专业Fortran论坛 - (fcode.cn)http://bbs.fcode.cn/forum.php?mo ... ghlight=MKL%5C_DFTI)的代码时输出复数实虚部也为0。。。
结果如下:
结果2.png

最佳答案

查看完整内容

1. 你要判断一下返回值 status 是否为0,如果不是0,说明这个步骤出了错。 2. DFTI_DOUBLE 代表双精度,你需要把 complex 改为 complex(8) 才能匹配 DFTI_DOUBLE 双精度。 当然你也可以改为 DFTI_SINGLE 这样表示单精度。总之要和传入的数据匹配。 3. DftiCreateDescriptor 函数,你传入 2,l 是要做二维DFT吗? 4. 如果你不写 status = DftiSetValue( my_desc1, DFTI_PLACEMENT , DFTI_NOT_INPLACE ) 这个语句,代表输出y存入x,并 ...

1820

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
760 元
贡献
305 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2022-12-6 15:56:03 | 显示全部楼层
1. 你要判断一下返回值 status 是否为0,如果不是0,说明这个步骤出了错。
2. DFTI_DOUBLE 代表双精度,你需要把 complex 改为 complex(8) 才能匹配 DFTI_DOUBLE 双精度。
当然你也可以改为 DFTI_SINGLE 这样表示单精度。总之要和传入的数据匹配。
3. DftiCreateDescriptor 函数,你传入 2,l 是要做二维DFT吗?
4. 如果你不写
status = DftiSetValue( my_desc1, DFTI_PLACEMENT , DFTI_NOT_INPLACE )
这个语句,代表输出y存入x,并且你不能在 status = DftiComputeForward(my_desc1,x,y) 里给定 y
如果你写了上述语句,则可以指定 y。
以下代码,假定你想要计算双精度,复数域,一维DFT,结果存入y数组(而不是x数组)
如果这不是你的目的,则要根据你的目的决定如何修改。

[Fortran] 纯文本查看 复制代码
program test2
  use mkl_dfti
  complex(8) :: x(8),y(8)
  type(DFTI_DESCRIPTOR), POINTER :: my_desc1
  integer :: status
  x = [1,2,3,4,5,6,7,8]
  print *, x
  status = DftiCreateDescriptor(my_desc1,DFTI_DOUBLE,DFTI_COMPLEX,1,size(x))
  status = DftiSetValue( my_desc1, DFTI_PLACEMENT , DFTI_NOT_INPLACE )
  status = DftiCommitDescriptor(my_desc1)
  status = DftiComputeForward(my_desc1,x,y)
  status = DftiFreeDescriptor(my_desc1)
  print *, y
end program test2

回复

使用道具 举报

5

帖子

1

主题

0

精华

新人

F 币
16 元
贡献
7 点
 楼主| 发表于 2022-12-7 11:15:00 | 显示全部楼层
fcode 发表于 2022-12-6 20:03
1. 你要判断一下返回值 status 是否为0,如果不是0,说明这个步骤出了错。
2. DFTI_DOUBLE 代表双精度,你 ...

非常感谢您的回答。
原代码是做的二维DFT,但是在做一维DFT时出现了同样的问题,输出的参数实虚部皆为0.
运行了一下您所附代码,输出结果y时实虚部同样为0。根据建议对status进行了判断,的确不为0。请问这个大概是什么原因引起的呢?
回复

使用道具 举报

5

帖子

1

主题

0

精华

新人

F 币
16 元
贡献
7 点
 楼主| 发表于 2022-12-7 11:23:53 | 显示全部楼层
th0rn4 发表于 2022-12-7 11:15
非常感谢您的回答。
原代码是做的二维DFT,但是在做一维DFT时出现了同样的问题,输出的参数实虚部皆为0.
...

问题解决了,太感谢您了,调用mkl库的时候除了一点小问题。。
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2023-2-4 17:07

Powered by Tencent X3.4

© 2013-2023 Tencent

快速回复 返回顶部 返回列表