使用MKL_DFTI进行快速傅里叶变换时，输出复数实虚部皆为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] 纯文本查看 复制代码

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