使用FFTW Guru接口计算实输出逆傅里叶变换时计划创建失败排查

核心问题：输入尺寸不匹配FFTW c2r变换的要求 #

FFTW的fftw_plan_guru_dft_c2r是用于从Hermitian对称的压缩复数组恢复实数组的逆变换（对应实信号DFT的逆过程），它对输入复数组的尺寸有严格要求：

• 若输出是N×M的实矩阵，输入复矩阵的最后一维长度必须是M//2 + 1（无论M是奇数还是偶数）。
• 对于你3×3的实输出需求，输入复矩阵的正确尺寸应该是3×2，而非你当前使用的3×3。

你当前传入3×3的复数组作为输入，FFTW无法识别符合c2r变换的Hermitian压缩格式，因此返回NULL计划。而fftw_plan_guru_dft（复转复）不要求输入对称，且输入输出尺寸一致，所以能正常工作。

其他可能的错误点 #

• Hermitian对称性验证：即使尺寸正确，输入复数组必须严格满足Hermitian对称性：对于输入复数组X[k][l]，需满足X[k][l] = conjugate(X[N−k mod N][M−l mod M])，且X[0][0]必须为实数。你的矩阵看起来符合对称，但建议用代码验证浮点精度下的对称性（避免微小误差导致计划创建失败）。
• Guru接口参数设置错误：调用fftw_plan_guru_dft_c2r时，需正确指定维度参数：输出实数组的维度是{3,3}，输入复数组的维度应设为{3, 2}（最后一维为3//2 +1=2），同时要确保dims和howmany_dims等参数正确对应输入输出的维度结构，避免维度顺序或数量错误。

修正步骤 #

1. 将你的3×3复矩阵压缩为3×2的Hermitian格式：
   保留第0列全部元素，第1列保留前2个元素（因为3//2+1=2），丢弃第2列（可通过对称关系推导）。示例压缩后的输入：

   [-1.2190 + 0.0000i   1.0162 - 0.6787i;
    -3.3477 + 1.3537i   -1.1505 + 3.4269i;
    -3.3477 - 1.3537i   -2.8059 + 0.5052i]
   

2. 调整fftw_plan_guru_dft_c2r的参数，指定输入为3×2复数组，输出为3×3实数组。
3. 确保输入数组的内存布局符合FFTW要求（建议用fftw_malloc分配内存）。

验证建议 #

如果不想手动压缩数组，可以先用fftw_plan_guru_dft（复转复）执行逆变换得到实数组，再用fftw_plan_guru_dft_r2c执行正变换，查看输出的压缩复数组格式，对比你当前的输入格式，确认差异。

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