FFTW3单精度（FP32）R2C FFT输出全零问题排查求助

问题排查与解决思路 #

针对你遇到的FFTW3单精度R2C FFT输出全零的问题，这些是最可能的原因及修复方式：

• 输出指针转换错误：
  代码中reinterpret_cast<fftwf_complex *>(&out)是核心错误。out本身已是complexf*类型，取地址后变成complexf**，强制转换会让FFTW把结果写入指针变量out的内存地址（而非输出数组首地址），导致真正的输出数组完全未被写入，呈现全零。
  正确写法应直接将out转换为fftwf_complex*：

  plan = fftwf_plan_dft_r2c_2d(N0, N1, input, reinterpret_cast<fftwf_complex*>(out), FFTW_ESTIMATE);
  

  （注：std::complex<float>与fftwf_complex内存布局兼容，实部、虚部存储顺序一致，该转换安全）

• 输出数组大小计算错误：
  R2C二维FFT的输出维度并非N0*N1，而是N0*(N1/2 + 1)——因实输入的FFT结果共轭对称，仅需存储正频率分量。若size_of_output未按此规则计算，要么内存不足导致写入越界，要么分配内存远超实际需要，最终读取到未初始化的零值。
  正确的输出大小计算：

  size_t size_of_output = N0 * (N1 / 2 + 1);
  

• 内存对齐或类型兼容性问题：
  部分编译器对std::complex<float>的布局可能有特殊处理，导致与fftwf_complex对齐不匹配。可直接用fftwf_complex*定义输出指针，避免类型转换问题：

  fftwf_complex* out = (fftwf_complex*)fftwf_malloc(sizeof(fftwf_complex) * size_of_output);
  plan = fftwf_plan_dft_r2c_2d(N0, N1, input, out, FFTW_ESTIMATE);
  

  读取结果时，可转换为complexf*，或直接通过out[i][0]（实部）、out[i][1]（虚部）访问。

• 库链接与头文件匹配问题：
  确认链接的是你编译的单精度库：

  • 检查-L/path/to/fftw_f/.libs/路径是否正确，目录下是否存在libfftw3f.a或libfftw3f.so；
  • 确保编译时include的是对应单精度版本的fftw3.h，未引入系统默认的双精度头文件；
  • 用nm工具检查可执行文件，确认链接的是fftwf_前缀的函数，而非双精度的fftw_函数。

• 输入数据隐性问题：
  即便已转为FP32，仍需验证：

  • 输入数组大小确实为N0*N1，无长度不足；
  • 输入数据未被其他代码意外覆盖或清零；
  • 用简单测试输入（如全1数组）验证FFT输出是否符合预期，排除输入本身问题。

内容的提问来源于stack exchange，提问作者mr_jnrdve