g++ 5.4.0中std::exp处理小float值返回NaN是否为标准库Bug？

你碰到的这个问题确实是g++ 5.4.0配套的libstdc标准库的一个已知实现Bug，完全不是你对std::exp的理解有误——按照C标准的规定，std::exp对于负的有限输入，哪怕数值极小，也应该返回趋近于0的正浮点数；如果输入是负无穷，返回+0，绝对不应该返回NaN。

问题原因 #

g++ 5.x系列的libstdc++中，float版本的std::exp实现存在边界处理缺陷。float类型的指数范围是-126到127，当输入的x小于约-87时，exp(x)的结果已经接近float能表示的最小正正规数（约1.175e-38），此时标准库的实现错误地触发了NaN，而没有正确处理下溢情况（返回极小值或+0）。

完整复现代码 #

你提供的代码可以补全为以下形式，用来验证这个问题：

#include <cmath>
#include <fenv.h>
#include <iostream>

int main() {
    feenableexcept(FE_ALL_EXCEPT & ~FE_INEXACT);
    for (float x = 30; x > -1000; --x) {
        float result = std::exp(x);
        if (std::isnan(result)) {
            std::cout << "std::exp(" << x << ") = NaN" << std::endl;
            break;
        }
    }
    return 0;
}

用g++ 5.4.0编译运行这段代码，会输出类似std::exp(-88) = NaN的结果，这明显不符合标准行为。

解决办法 #

• 优先升级编译器：这个Bug在g++ 6及以上的版本中已经被修复，升级后std::exp会正确返回趋近于0的极小值或+0。
• 临时 workaround：如果暂时无法升级编译器，可以手动添加边界判断，提前处理会触发Bug的输入：

float safe_exp(float x) {
    // 阈值基于float的最小正正规数对应的exp输入
    const float EXP_UNDERFLOW_THRESHOLD = -87.0f;
    if (x < EXP_UNDERFLOW_THRESHOLD) {
        return 0.0f;
    }
    return std::exp(x);
}

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