使用MKL_DFTI进行快速傅里叶变换时，输出复数实虚部皆为0

th0rn4

运行代码如下：

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	program test2

02	use mkl_dfti

03	complex :: x_2D(2,4),y_2D(4,8)

04	complex :: x(8),y(32)

05	type(DFTI_DESCRIPTOR), POINTER :: my_desc1

06	integer :: status,l(2)

07	equivalence(x_2D,x)

08	equivalence(y_2D,y)

09

l(1)=2

10

l(2)=4

11	x_2D=reshape([1,2,3,4,5,6,7,8],[2,4])

12	print *, x

13	status = DftiCreateDescriptor(my_desc1,DFTI_DOUBLE,DFTI_COMPLEX,2,l)

14	status = DftiCommitDescriptor(my_desc1)

15	status = DftiComputeForward(my_desc1,x,y)

16	status = DftiFreeDescriptor(my_desc1)

17	print *, x

18	print *, y

19	end program test2

运行结果如图：

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

woduiren

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

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

th0rn4

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

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

th0rn4

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

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

fcode

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] 纯文本查看 复制代码

查看代码

打印?

01	program test2

02	use mkl_dfti

03	complex(8) :: x(8),y(8)

04	type(DFTI_DESCRIPTOR), POINTER :: my_desc1

05	integer :: status

06	x = [1,2,3,4,5,6,7,8]

07	print *, x

08	status = DftiCreateDescriptor(my_desc1,DFTI_DOUBLE,DFTI_COMPLEX,1,size(x))

09	status = DftiSetValue( my_desc1, DFTI_PLACEMENT , DFTI_NOT_INPLACE )

10	status = DftiCommitDescriptor(my_desc1)

11	status = DftiComputeForward(my_desc1,x,y)

12	status = DftiFreeDescriptor(my_desc1)

13	print *, y

14	end program test2