gfortran指定编译选项下仍对cos(atan(x))进行不安全数学优化的咨询

问题概述 #

在使用gfortran编译代码时，即使指定-fno-fast-math和-fno-unsafe-math-optimizations选项，不同优化级别下仍会出现结果差异。测试覆盖RedHat的gfortran-11.2.1、gfortran-11.4.1以及gfortran-14.2版本，且使用OpenLibM替代gcc数学库保证结果可复现。

差异源于y = cos(atan(x))的计算：

• 在-O1及更高优化级别下，该表达式会被替换为y = 1/sqrt(1+x²)；
• 在-O0级别下，会直接调用atan和cos函数。

测试代码如下：

double precision function foo(x)
implicit none
double precision, intent(in) :: x
foo = cos(atan(x))
end function foo

编译选项的行为差异 #

根据文档描述，-fno-fast-math或-fno-unsafe-math-optimizations应该禁用此类数学变换，但实际对gfortran无效。唯一能禁用该优化、让各级别结果与-O0一致的选项是-fno-tree-forwprop。

对比gcc C编译器的行为：相同逻辑的C代码在默认-O1下不会执行此替换，仅当指定-ffast-math时才会触发，符合预期。

解决方案 #

若要保证所有优化级别下结果一致，同时尽可能保留性能，推荐使用-fno-tree-forwprop而非-O0：

• -O0会关闭所有优化，导致性能大幅下降；
• -fno-tree-forwprop仅禁用树向前传播阶段的代数恒等变换，不会影响其他安全优化，能在保证结果一致性的前提下维持较好性能。

目前没有其他更精细的选项可以单独禁用这一特定变换，因为该优化被归类为标准的代数简化，不属于-ffast-math管控的"不安全"数学优化范畴。

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