Fortran代码SIGFPE错误回溯指向WHERE循环，求定位方法

这种幽灵式的浮点错误在大型Fortran代码里真的很棘手——明明回溯指向的代码看起来没问题，错误还跟着无关代码改来改去，完全摸不着头脑。别着急，咱们一步步来定位根源：

• 用内存检查工具揪出隐形污染
  Valgrind是对付这类“错误不在报错行”问题的神器，尤其是内存越界、未初始化变量这类会污染后续数据的问题。运行命令：

  valgrind --leak-check=full --track-origins=yes ./your_program
  

  它会帮你追踪到所有非法内存访问，哪怕是之前某个数组越界写入了WHERE循环用到的数组内存，Valgrind都能精准揪出来。

• 加强编译时的初始化检查
  给编译命令加上-finit-real=snan -finit-integer=-999999，这样所有未初始化的实数会被设为NaN，整数会被设为明显的异常值。一旦代码用到这些未初始化的值，会立刻触发FPE或者出现明显的错误输出，帮你找到真正的问题源头，而不是等到WHERE循环才爆错。

• 获取更详细的浮点错误信息
  在编译命令里加上-ffpe-summary=all，这样触发SIGFPE时会输出具体是哪种浮点错误（比如是除以零？还是无效操作？），结合这个信息去检查WHERE循环用到的变量在进入循环前的状态——比如是不是某个数组元素被之前的代码改成了NaN，或者出现了零分母？

• 用调试器追踪变量的异常变化
  用gdb启动你的程序：

  gdb ./your_program
  

  在触发错误的WHERE循环前设置断点（比如break f_mod.f90:130），运行到断点后，用print命令检查所有相关数组、变量的值，看看有没有异常。还可以用watch命令监视WHERE循环里用到的数组，比如watch my_array(1:10)，这样只要有代码修改这个数组，程序就会停下来，你就能找到是谁在偷偷改数据。

• 排查编译器版本或优化的影响
  你用的gfortran 6.4.0是比较老的版本了，有可能是编译器的bug导致的。试试用较新的gfortran版本（比如9.x或10.x）重新编译测试，如果错误消失，那大概率是旧编译器的问题；如果问题还在，就回到代码层面继续查。另外，记得加上-O0关闭所有优化，有时候优化会让错误的表现更诡异。

• 重点检查数组边界和内存布局
  既然错误会随无关代码修改出现/消失，这大概率是内存布局相关的问题——比如某个数组越界写入，刚好覆盖到了WHERE循环用到的变量内存，而修改无关代码会改变内存布局，导致覆盖的位置变化。可以单独加上-fbounds-check（虽然-fcheck=all包含它，但有时候单独启用会有更详细的输出），严格检查所有数组的访问是否越界。

内容的提问来源于stack exchange，提问作者R. N