[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