C++比较函数中FMA指令破坏严格弱序属性的相关问题咨询

为什么FMA特性默认开启 #

• 收益远高于普遍风险：FMA能同时提升浮点运算的性能和精度——单条FMA指令即可完成乘法+加法操作，吞吐是分开执行两条指令的2倍左右；同时省略了乘法结果的中间舍入步骤，计算结果精度比分开运算更高。绝大多数数值计算、普通业务场景下，这一特性能带来明确收益，且不会触发逻辑问题。
• 符合C标准规范：C语言标准默认允许浮点表达式的融合收缩，并未要求编译器必须保留乘加运算的中间舍入步骤，GCC默认开启FMA代码生成属于符合标准的行为。
• 触发问题的场景极罕见：本次出现的严格弱序破坏问题仅发生在「比较函数中使用可被FMA优化的浮点表达式」这一非常边缘的场景，大部分开发者很难碰到这类问题，编译器没有理由为了极低概率的边缘case牺牲全局的性能、精度收益。

除禁用全局FMA外的规避方案 #

• 改写比较函数逻辑：不要使用a*b - c*d < 0这类可被融合的表达式写法，直接改为比较两个乘积的大小return s1.a * s1.b < s2.a * s2.b;，自然不会触发FMA优化。如果需要保证中间计算不被融合，可以显式将乘积存入局部变量后再比较：

bool compare(const S& s1, const S& s2) {
    double p1 = s1.a * s1.b;
    double p2 = s2.a * s2.b;
    return p1 < p2;
}

如果还是担心被优化，可以给局部变量加volatile修饰强制编译器读写内存，阻断跨语句的融合优化。

• 局部禁用FMA优化：不需要全局关FMA，只针对出问题的比较函数做特殊处理即可：
  • GCC/Clang下可以给函数加属性：__attribute__((optimize("ffp-contract=off")))，仅关闭该函数的浮点融合优化
  • 用标准pragma局部控制：在比较函数前加#pragma STDC FP_CONTRACT OFF，函数结束后加#pragma STDC FP_CONTRACT ON，仅影响中间的代码段
• 显式控制浮点运算行为：如果需要严格保证浮点运算的舍入逻辑，可以使用C++标准库提供的显式运算函数，或者在比较时增加可接受的精度容差，实现符合严格弱序要求的浮点比较逻辑。

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