咨询GNU C++标准库libstdc++中std::exp的指数计算算法

Great question! Let's break down exactly how std::exp works under the hood in GNU's libstdc++ standard library.

GNU libstdc的std::exp复用了经典的**FDlibm（Freely Distributable Math Library）**实现，这是多数C/C标准库的底层数学计算依赖。具体计算流程分为三步：

1. 参数预处理与范围缩小 #

• 先处理特殊输入：比如exp(+infinity)直接返回+infinity，exp(-infinity)返回0，exp(nan)返回nan。
• 对于普通实数x，利用指数函数的性质 exp(x) = 2^k * exp(x - k*ln2)，将x分解为整数部分k和小数部分r（保证|r| ≤ ln2/2 ≈ 0.3466）。这一步把计算范围压缩到一个更易精确近似的小区间，大幅降低后续近似的难度。

2. 有理多项式逼近（核心计算） #

在缩小后的区间内，exp(r)通过6次有理多项式来近似，形式大致为：

exp(r) ≈ 1 + r * P(r²) / Q(r²)

其中P和Q是经过高精度拟合得到的多项式系数，在|r| ≤ 0.3466的区间内，这个近似的误差可以控制在IEEE双精度浮点数的精度上限以内，完全满足日常计算需求。

3. 结果缩放 #

最后将得到的exp(r)乘以2^k——这里的幂次计算不需要实际做乘法，而是直接操作IEEE754浮点数的指数域，实现高效的数值缩放。

补充：为什么不用算术几何均值（AGM）算法？ #

AGM算法确实可以计算指数函数，但它的优势体现在超高精度场景（比如数百位小数计算）。对于普通双精度浮点数，多项式/有理函数逼近的速度更快、实现更简洁，精度也完全达标，因此libstdc++选择了后者作为std::exp的实现方案。

内容的提问来源于stack exchange，提问作者Herpes Free Engineer