FFTW3中RODFT00变换的归一化问题（Intel Fortran环境）

问题背景 #

在Visual Studio环境中使用Intel Fortran编译器调用FFTW3库，针对sinx-siny分量的谱空间与网格空间互转（RODFT00类型2D实变换），归一化操作后结果不符合预期：

• 已创建变换计划（因FFTW为C库行主序，Fortran为列主序，x、y维度已反转）：

  plan_spectral_to_grid_sin_sin = fftwf_plan_r2r_2d(Ny, Nx, dummy, dummy, FFTW_RODFT00, FFTW_RODFT00, FFTW_MEASURE)
  plan_grid_to_spectral_sin_sin = fftwf_plan_r2r_2d(Ny, Nx, dummy, dummy, FFTW_RODFT00, FFTW_RODFT00, FFTW_MEASURE)
  

• 当前归一化方案：
  计算归一化因子：

  sin_sin_norm = (2.0 * REAL(currentDims%IXMAX + 1) * 2.0 * REAL(currentDims%IYMAX + 1))
  

  操作流程：将谱分量A(k,l)sinkxsinly除以该因子后执行正变换，再执行逆变换，预期得到未归一化的原始谱分量，但结果不符。
• 补充信息：horzArr%G1是维度为currentDims%IXMAX、currentDims%IYMAX的二维数组

排查方向与解决方案 #

1. 修正归一化时机 #

FFTW的RODFT00（离散正弦变换类型1）正逆变换为同一操作，变换对满足：连续执行两次变换后，结果为原始数据乘以[2*(Nx+1)]*[2*(Ny+1)]（Nx、Ny为对应维度的变换点数）。
正确操作流程应为：

1. 对谱分量A(k,l)直接执行正变换（谱→网格）
2. 对网格数据执行逆变换（网格→谱）
3. 将逆变换结果除以归一化因子sin_sin_norm，即可还原原始A(k,l)

提前在正变换前除以归一化因子会导致最终结果幅度偏差，不符合预期。

2. 校验计划与数组维度的一致性 #

确保计划中的Nx、Ny与数组horzArr%G1的维度严格对应：

• 计划参数fftwf_plan_r2r_2d(Ny, Nx, ...)中，Ny对应Fortran数组的第二维度（currentDims%IYMAX），Nx对应第一维度（currentDims%IXMAX）
• 若维度不匹配，变换过程会出现数据错位，直接导致结果异常

3. 简化计划创建（可选） #

由于RODFT00的正逆变换逻辑完全相同，无需创建两个独立计划，复用同一个计划即可：

plan_sin_sin = fftwf_plan_r2r_2d(Ny, Nx, dummy, dummy, FFTW_RODFT00, FFTW_RODFT00, FFTW_MEASURE)
! 谱→网格正变换
call fftwf_execute_r2r(plan_sin_sin, horzArr%G1, horzArr%G1)
! 网格→谱逆变换
call fftwf_execute_r2r(plan_sin_sin, horzArr%G1, horzArr%G1)

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