ULP为1时的舍入规则及C++与Java exp()结果一致性问题

问题背景 #

我发现C++的std::exp(x)（x为double类型）与Java的Math.exp(x)计算结果存在差异，二者采用的exp()算法不同，但均保证舍入误差不超过1 ULP。

核心问题及解答 #

1. 1 ULP约束下的舍入逻辑 #

• ULP（Unit in the Last Place）指浮点数表示中最后一位的单位值，1 ULP误差意味着计算结果与真实数学值的差距不超过该浮点数所能表示的最小精度单位。
• 舍入方向遵循IEEE 754标准默认的“就近舍入（round to nearest, ties to even）”规则：当真实值恰好落在两个可表示浮点数的中间时，选择最后一位为偶数的那个值；否则选择最接近真实值的浮点数。
• 舍入的执行者取决于具体实现：可能是算法本身在计算过程中处理舍入，也可能是底层硬件指令（如x87或SSE指令集）完成舍入，但最终都需符合IEEE 754规则。

2. 复制fdlibm算法能否获得与Java一致的结果？ #

• 可行，但需注意细节：Java的Math.exp确实基于fdlibm实现，但Java调用fdlibm代码时可能存在额外的标准化处理（如输入范围预处理、舍入步骤严格对齐）。直接引入fdlibm的e_exp.c代码，需确保编译环境的浮点数模式（如舍入模式、精度设置）与Java完全一致，否则结果仍可能有差异。
• 1 ULP的误差约束不代表舍入选择完全任意：不同算法的计算路径会得到不同中间值，最终舍入后的结果可能落在真实值的两侧（一个向上舍入，一个向下舍入），但都符合1 ULP要求。因此要完全对齐Java结果，必须复用其完整实现逻辑（包括fdlibm代码+Java的封装处理）。

附加问题解答 #

使用MSVC 14.16.27023编译时，C++在x86与x64平台下std::exp(x)结果不同、有时与Java结果一致的原因：

• x86平台默认使用x87浮点单元，默认精度为80位扩展精度，计算过程中会用更高精度的中间值，最后截断为64位double；而x64平台默认使用SSE/SSE2指令集，直接用64位double进行计算。
• 这种计算精度差异会导致最终舍入结果不同：x87的高精度中间值可能让舍入后的结果更接近真实值，或触发不同的舍入方向，有时刚好与Java结果一致（比如Java的实现也用了类似的高精度中间步骤，或舍入方向巧合相同）。

更新补充说明 #

引入fdlibm的exp()实现后，结果与C++ std::exp、Java Math.exp均不同，可能的原因：

• 编译环境浮点数设置不同：MSVC的浮点模式（如/fp:precise或/fp:fast）会影响舍入和中间计算精度，与fdlibm原本基于IEEE 754默认模式的假设不一致。
• 缺少Java对fdlibm的封装处理：Java调用fdlibm的exp实现时，可能对输入做了额外的范围检查、归一化处理，或强制设置了特定舍入模式，直接使用原始fdlibm代码没有这些步骤，导致结果偏差。

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