c++ - 将 STL complex<double> 转换为 fftw_complex 的问题

Question

FFTW 手册 says它的fftw_complex类型与 std::complex<double> 位兼容STL 中的类。但这对我不起作用:

#include <complex>
#include <fftw3.h>
int main()
{
   std::complex<double> x(1,0);
   fftw_complex fx;
   fx = reinterpret_cast<fftw_complex>(x);
}

这给了我一个错误:

error: invalid cast from type ‘std::complex<double>’ to type ‘double [2]’

我做错了什么？

Answer 1

fftw_complex 与 C99 和 C++ 复杂类型的位兼容性背后的想法不是它们可以很容易地相互创建，而是 FFTW 中所有采用指向 fftw_complex 的指针的函数也可以采用指向 c++ std::complex 的指针.因此，最好的方法可能是在整个程序中使用 std::complex<> 并且仅在调用 FFTW 函数时将指针转换为这些值:

std::vector<std::complex<double> > a1, a2;
....
....
fftw_plan_dft(N, reinterpret_cast<fftw_complex*>(&a1[0]),
                 reinterpret_cast<fftw_complex*>(&a2[0]),
                 FFTW_FORWARD, FFTW_ESTIMATE);
....