PPC64架构下long double类型机器epsilon计算异常问题咨询

问题根因分析 #

1. 环境long double类型确认 #

你所使用的PPC64 POWER8环境默认long double为double-double实现，符合你观测到的LDBL_MANT_DIG=106的特征：该类型由两个连续的IEEE 754 binary64（double）拼接而成，高位存储数值的主近似值，低位存储误差修正值，有效尾数为两个double的53位尾数减去重叠的高位符号位，总长度刚好为106位，不属于IEEE 754 float128标准实现。

2. 异常epsilon输出的核心原因 #

你得到的4.94066e-324不是long double的真实机器epsilon，是输出环节的精度丢失导致的错误结果：

• C++标准库的std::cout默认不会对long double做全精度输出，绝大多数实现会在输出前将long double隐式转换为普通double类型
• 4.94066e-324本身是binary64（double）类型能表示的最小非规格化值，是long double的真实epsilon转换为double时发生下溢得到的截断值，和实际值完全不符。
• 哪怕你直接打印LDBL_EPSILON或者std::numeric_limits<long double>::epsilon()，只要用默认配置的std::cout输出，都会触发同样的转换截断问题。

3. 正确验证方法 #

替换输出逻辑为以下两种方案之一即可得到正确结果：

方案1：使用printf直接输出long double

#include <cstdio>
int main()
{
    long double eps = 1.0l;
    while (1.0l + 0.5l * eps != 1.0l)
        eps = 0.5l * eps;
    printf("long double epsilon: %Lg\n", eps);
    return 0;
}

方案2：给std::cout配置long double全精度输出

#include <iostream>
#include <iomanip>
int main()
{
    long double eps = 1.0l;
    while (1.0l + 0.5l * eps != 1.0l)
        eps = 0.5l * eps;
    std::cout << std::setprecision(40) << eps << std::endl;
    return 0;
}

正确输出结果约为1.23e-32（对应2^-106），确实比IEEE 754 float128的epsilon（1.93e-34）精度更低，和公开资料对double-double类型的精度描述完全吻合。

补充：切换为IEEE 128位浮点的方法 #

如果需要在POWER8环境下使用符合IEEE 754标准的float128作为long double实现，可以在编译时添加参数-mfloat128 -mabi=ieeelongdouble，此时LDBL_MANT_DIG会变为112，epsilon也会匹配IEEE float128的标准值。

内容的提问来源于stack exchange，提问作者anymous.asker