浮点收缩仅允许在单个表达式内？IEEE754/FP_CONTRACT行为确认

以下是触发断言失败的不稳定代码：

int x = ~~~;
const double s = 1.0 / x;
const double r = 1.0 - x * s;

assert(r >= 0); // fail

这段代码的问题在于编译器可能将1.0 - x * s优化为FMA（融合乘加）指令——FMA会直接计算1.0 - (x*s)而不单独对x*s进行舍入，导致结果可能为负数。

而拆分后的代码则不会触发断言失败：

int x = ~~~;
const double s = 1.0 / x;
const double tmp = x * s;
const double r = 1.0 - tmp;

assert(r >= 0); // OK

拆分表达式是否被标准保证能阻止浮点收缩？ #

是的，这种拆分行为符合C/C++标准，能够可靠阻止浮点收缩（如FMA优化）：

1. C/C++标准中的FP_CONTRACT规则：
   C11、C++11及后续标准通过FP_CONTRACT宏控制编译器是否允许将多个浮点运算合并为单个收缩指令。默认情况下编译器可选择允许收缩，但收缩仅允许在单个表达式内部进行。当用中间变量（如tmp）存储中间结果时，标准要求该结果必须按对应浮点类型（此处为double）完成舍入后再存储，编译器无法跨语句合并运算——x * s的结果必须先舍入为double，再参与后续的1.0 - tmp运算，无法直接合并为FMA。

2. IEEE754的角色：
   IEEE754标准定义了浮点运算的精度和行为，但并未直接规定编程语言的编译优化规则。编程语言对浮点运算的约束由C/C标准规范，只要编译器遵循C/C标准，拆分后的代码就不会被合并为FMA。

是否需要采用更具数值稳定性的方案？ #

虽然拆分表达式的方式是可靠的，但从长期维护和数值健壮性角度，建议优先采用数值稳定性更好的实现方案：

• 避免依赖浮点运算的舍入特性保证结果范围，尽量使用数学上等价且精度更可控的计算形式。
• 如果必须保证r >= 0，可以采用更直观的范围约束（如std::max(0.0, 1.0 - x*s)），虽然看似不够优雅，但能更明确地表达代码意图，避免后续维护者对拆分逻辑产生误解。

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