Java的Math.exp()与fdlibm的C实现计算结果不一致的原因探究

你遇到的结果差异（两者仅相差1个ULP，即双精度浮点数最后一位的最小单位），主要源于以下几个实际情况，和文档描述的“委托关系”并不矛盾：

• JVM的StrictMath并非完全照搬原始fdlibm代码：虽然Oracle文档说明Math.exp()会委托StrictMath，而StrictMath基于fdlibm实现，但主流JVM（如OpenJDK）会对fdlibm的参考代码做适配修改——比如调整舍入步骤、优化计算逻辑，以满足Java的可移植性要求或性能需求，并非直接使用原始的e_exp.c代码。
• 浮点运算环境的差异：C语言实现的fdlibm依赖系统的浮点运算单元（FPU），不同平台或编译器的FPU中间精度（比如x86平台默认用80位扩展精度存储中间值）、舍入模式可能和Java虚拟机的处理逻辑不同，这些细微差异会累积导致最终结果差1个ULP。
• 平台特定的替代实现：部分JVM会根据运行平台选择系统级数学库（如Linux的glibc、Windows的msvcrt）来实现StrictMath的函数，而非严格使用fdlibm，这些系统库的exp函数采用了不同的近似算法，结果自然会有差异。

这种1个ULP的差异符合IEEE 754双精度浮点数的精度规范，属于正常的计算误差范围。

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