x86-64下浮点数自加与乘2运算的精度对比

标准x86-64环境下，只要浮点运算遵循IEEE 754语义（默认编译配置，未开启-ffast-math这类破坏标准合规性的激进优化），两个取值完全相同的float/double相加，结果是完全精确的，和乘以2的运算没有任何区别，不存在精度损失。

IEEE 754下的所有浮点数（包括规格化数、非规格化数、特殊值）都可以统一表示为 $s \times m \times 2^e$ 的形式：s是符号位，m是尾数，e是指数。两个值完全相等的浮点数相加，本质就是计算原值的2倍，和乘2的逻辑完全一致：

• 如果计算2m之后，m仍然落在对应类型的尾数合法区间，只需要调整尾数位即可，不需要做任何舍入；
• 如果2m超出了当前尾数的范围（比如规格化数算完2m≥2，非规格化数算完2m≥1），只需要把尾数逻辑左移一位、指数加1，整个过程也不会产生任何舍入误差；
• 特殊值场景也不存在精度问题：同符号无穷大相加结果还是对应符号的无穷大，NaN加自身还是NaN，结果溢出到无穷大的行为也和乘2完全一致；
• 哪怕是用传统x87 FPU做运算（内部用80位扩展精度计算），float/double加载到80位寄存器的转换是完全无损的，两个相同值相加后再转回原类型的结果，和乘2的结果也完全一致，不会产生额外误差。

如果你开启了编译器的激进浮点优化（比如GCC/Clang的-ffast-math、MSVC的/fp:fast），浮点运算不再严格遵守IEEE 754规则，编译器可能随意调整运算顺序、降低中间计算精度，极端场景下可能出现a+a和2*a结果不一致的情况，但这不属于标准x86-64环境的默认行为。

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