如何实现Intel与AMD平台间Fortran数值计算结果的可复现性？

核心原因 #

• 硬件指令集差异：Intel Xeon E5-2670（Sandy Bridge架构）仅支持AVX指令集，AMD 7700X（Zen 4）则支持AVX-512、AVX2等更高级指令集。即便用-O0编译，GCC仍会根据目标硬件生成对应指令集的浮点运算代码，不同指令集的运算单元在舍入精度、底层实现上存在细微差异，直接导致结果末尾数值不同。
• 二进制兼容性问题：Intel机器编译的可执行文件针对Sandy Bridge架构优化，直接放到AMD Zen4平台运行时，部分指令无法兼容，触发异常执行逻辑，结果自然与本地编译运行的结果都不匹配。
• 浮点运算顺序差异：编译器会根据硬件特性调整运算顺序（如自动向量化、循环展开），即使-O0也无法完全禁止这类优化。浮点运算不满足结合律，运算顺序的细微变化会导致最终结果的末尾数值出现差异。
• 系统依赖库差异：两台机器的Linux Mint版本不同（21.2 vs 21.3），对应的glibc、libgfortran等系统库版本存在差异，Fortran程序调用的数学库（如libm）实现细节不同，也会影响计算结果。

解决方法 #

• 统一编译指令集：编译时指定两台机器都支持的最低指令集，强制生成兼容代码。例如使用通用x86-64指令集：

  FFLAGS = -O0 -g -march=x86-64 -mtune=generic
  

  或者指定与Intel机器匹配的指令集-march=corei7，AMD 7700X向下兼容该指令集。
• 固定浮点运算环境：
  • 添加编译选项-ffloat-store，禁止将浮点值缓存在寄存器中（避免x87 80位精度与SSE 64位精度的差异）。
  • 在Fortran程序开头添加代码，强制设置IEEE浮点舍入模式为“四舍五入到最近值”：

    use, intrinsic :: iso_fortran_env
    use, intrinsic :: ieee_exceptions
    use, intrinsic :: ieee_arithmetic
    call ieee_set_rounding_mode(ieee_nearest)
    

• 统一系统依赖库：
  • 在AMD机器上安装与Intel机器相同版本的glibc、libgfortran库。
  • 使用静态编译选项-static，将所有依赖库打包进可执行文件，消除系统库差异的影响。
• 禁止自动优化导致的运算顺序变化：添加编译选项-fno-tree-vectorize，禁止编译器自动向量化循环，确保运算顺序与代码逻辑完全一致。
• 避免跨机器直接运行本地编译的二进制文件：使用统一指令集编译的可执行文件在两台机器上测试，确保二进制兼容性。

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