Fortran Coder

标题: 使用MKL_DFTI进行快速傅里叶变换时，输出复数实虚部皆为0 [打印本页]

作者: th0rn4 时间: 2022-12-6 15:56
标题: 使用MKL_DFTI进行快速傅里叶变换时，输出复数实虚部皆为0
运行代码如下：

[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

运行结果如图：

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

作者: fcode 时间: 2022-12-6 15:56
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

作者: th0rn4 时间: 2022-12-7 11:15

fcode 发表于 2022-12-6 20:03
1. 你要判断一下返回值 status 是否为0，如果不是0，说明这个步骤出了错。
2. DFTI_DOUBLE 代表双精度，你 ...

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

作者: th0rn4 时间: 2022-12-7 11:23

th0rn4 发表于 2022-12-7 11:15
非常感谢您的回答。
原代码是做的二维DFT，但是在做一维DFT时出现了同样的问题，输出的参数实虚部皆为0.
...

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

作者: woduiren 时间: 2023-6-8 00:22

th0rn4 发表于 2022-12-7 11:23
问题解决了，太感谢您了，调用mkl库的时候除了一点小问题。。

楼主，请问问题怎么解决的？

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)