C++计算函数(1-exp(-x))/x时极小x输入（如10^-30）输出为0的问题排查

• 异常的核心是双精度浮点数的有效精度限制：双精度double类型在1附近的最小可分辨数值增量约为2.22e-16，简单来说只要两个数的差值小于这个量级，double就会判定两者完全相等。
• 当输入x=1e-30时，exp(-x)的理论值是1减去一个1e-30量级的极小值，这个差值远小于double的分辨阈值。此时代码中直接调用exp(-x)得到的结果会被直接截断为精确的1.0，分子部分1 - exp(-x)的计算结果直接变为0，最终fx输出自然为0，和理论上趋近于1的预期完全不符。
• 原代码附带的小隐患：循环终止条件依赖动态计算出的x值判断范围，而非根据迭代的指数步长判断，容易因浮点累计精度误差导致实际循环范围偏移；输出时未正确配置固定精度，也可能造成数值显示偏差。

• 放弃直接硬算1 - exp(-x)的写法，使用C数学库专门针对小值精度损失设计的expm1()函数计算。该函数从设计之初就瞄准了这类精度丢失场景，它不会先计算exp的值再做减法，而是直接高精度输出exp(y) - 1的结果，哪怕y小到1e-30量级，也能完整保留有效差值。
• 计算式等价替换：原计算式fx = (1 - exp(-x))/x可等价替换为fx = -expm1(-x) / x，这套写法可以完整覆盖1e-30到1e9的测试范围：
  • 极小x场景：expm1(-x)会高精度返回-x + x²/2 - ...的结果，与x作除后自然得到趋近于1的正确值
  • 极大x场景：exp(-x)会下溢为0，expm1(-x)返回-1，此时fx计算结果为1/x，也完全匹配理论值

#include <math.h>
#include <iostream>
#include <iomanip>

int main()
{
    std::cout << "fx\t\t\t\t\tx" << std::endl;
    // 按指数步长循环，范围从-30到9，对应x取值1e-30到1e9
    for(double exp_idx = -30; exp_idx <= 9; exp_idx += 0.01)
    {
        double x = pow(10.0, exp_idx);
        // 用expm1替换原计算逻辑，消除极小值下的精度损失
        double fx = -expm1(-x) / x;
        std::cout << std::setprecision(15) << fx << "\t\t\t\t" << x << std::endl;
    }
    return 0;
}

注意：部分老旧编译器需要开启对应C标准支持、正确关联数学库，才能正常使用expm1函数。

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