c++ - FFTW : how to prevent breaking aliasing rules?

Question

我有一个使用 std::complex<double> 的代码类型。来自 FFTW Manual :

if you have a variable complex<double> *x, you can pass it directly to FFTW via reinterpret_cast<fftw_complex*>(x).

但是，当我在我的代码中这样做时:

tmp_spectrum = reinterpret_cast<std::complex<double>*>(fftw_alloc_complex(conf.spectrumSize()));
plan_bw_temp = fftw_plan_dft_c2r_1d(conf.FFTSize(), reinterpret_cast<fftw_complex*>(tmp_spectrum), tmp_out, FFTW_ESTIMATE);

我得到 dereferencing type-punned pointer might break strict-aliasing rules [-Wstrict-aliasing] .如何解决这个警告？谢谢!

Answer 1

这里有三个选项:

• 只需在需要时创建一个fftw_complex:fftw_plan_dft_c2r_1d(conf.FFTSize(), fftw_complex(tmp_spectrum.real(), tmp_spectrum.imag()), tmp_out, FFTW_ESTIMATE);
• 不要在代码中使用 C++ 语言的复杂类型，而只使用 fftw_complex 类型。
• 使用 -fno-strict-aliasing 在适当的翻译单元中禁用所有严格别名优化和强制执行。仅使警告静音不安全，因为它可能会导致代码损坏。