x86与x86_64平台下浮点数转整数的截断误差差异问题

为什么x86-64和x86-32结果不同？ #

核心原因是两种平台默认使用的浮点运算单元（FPU）行为存在差异：

• x86-32默认采用x87 FPU，内部以80位扩展精度执行计算。1.234的双精度值（1.2339999999999999857...）在80位精度下会被更精确存储，乘以1000.0后得到的结果略小于1234，强制转换为整数时直接截断小数部分，最终得到1233。
• x86-64默认使用SSE/SSE2指令集，全程保持64位双精度运算。1.2339999999999999857... * 1000.0的计算结果，因双精度的舍入规则，最终会被表示为恰好的1234.0，所以强制转换后得到1234。

即使禁用优化，两种架构的默认浮点运算路径不同，直接导致了结果差异。

转换前加0.5是不是正确解决方案？ #

这个方案仅适用于正数且仅需四舍五入到最近整数的场景，存在明显局限性：

1. 负数处理错误：若输入为负数（如-1.234），-1.234*1000 +0.5 = -1233.5，强制转换会得到-1233，但符合预期的四舍五入结果应为-1234。
2. 边缘值舍入歧义：对于恰好为x.5的数值，加0.5后强制转换的行为，可能与预期的舍入逻辑（如银行家舍入）不一致。

更可靠的替代方案 #

如果需求是将秒数转换为毫秒并按标准规则四舍五入到最近整数，推荐使用标准库的舍入函数：

#include <cmath>

int64_t deltatime(double numSeconds) {
    return static_cast<int64_t>(std::round(numSeconds * 1000.0));
}

std::round遵循标准的"远离零舍入到最近整数"规则，能正确处理正负数值，且跨平台行为一致。

如果需求是向零截断（即直接丢弃小数部分），则可以通过编译选项统一运算精度：在x86-32环境下添加-mfpmath=sse，强制使用SSE指令集执行浮点运算，与x86-64平台保持一致。

内容的提问来源于stack exchange，提问作者umläute