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

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

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

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

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

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

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 :: 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