long double与std::uintmax_t转换时整数精度丢失问题排查

我实现了一个SplitLD类，用于将long double拆分为整数值和二进制指数以进行精确计算。该类通常运行正常，但在WSL（Ubuntu）环境下，long double与std::uintmax_t互转时会丢失若干最低有效位，且仅在代码分文件编译时出现此问题，单文件编译则正常。

相关代码 #

SplitLD.hh #

#include <cstdint>

// Splits `long double` into an integer and an exponent.
class SplitLD
{
public: // Everything is public to make the example easier to test.
    std::uintmax_t integer;
    int exponent;

    SplitLD(const long double number);
    operator long double() const;
};

SplitLD.cc #

#include <cfloat>
#include <cmath>
#include <limits>
#include <climits>
#include "SplitLD.hh"

SplitLD::SplitLD(long double number) // For the sake of simplicity, we ignore negative numbers and various corner cases.
{
    static_assert(FLT_RADIX == 2);
    static_assert(sizeof(std::uintmax_t) * CHAR_BIT >= std::numeric_limits<long double>::digits);
    // The following two operations change the exponent to make the represented value a whole number.
    number = std::frexp(number, &exponent);
    number = std::ldexp(number, std::numeric_limits<long double>::digits);
    exponent -= std::numeric_limits<long double>::digits;
    integer = number; // cast from `long double` to `std::uintmax_t`
}

SplitLD::operator long double() const
{
    long double number = integer; // cast from `std::uintmax_t` to `long double`
    number = std::ldexp(number, exponent);
    return number;
}

main.cc #

#include "SplitLD.hh"

int main()
{
    const long double x = 12345.67890123456789l; // arbitrarily chosen number for the test
    const SplitLD y = x;
    const long double z = y;
    return z == x ? 0 : 1;
}

复现条件 #

• 仅使用long double类型时出现，float和double无此问题；
• GCC和Clang均能复现；
• 仅代码分文件编译时出错，单文件编译正常；
• 仅在WSL环境下出现，同一机器的实体Linux或虚拟机Linux运行正常。

通过打印二进制表示排查发现，精度丢失发生在long double转std::uintmax_t以及std::uintmax_t转long double的赋值步骤，std::frexp和std::ldexp未改变尾数。

问题成因与分析 #

核心原因：WSL1的ABI兼容性缺陷 #

这是WSL1（Windows Subsystem for Linux 1）特有的long double ABI实现问题。在标准x86_64 Linux环境中，long double是80位扩展精度浮点类型，通过x87浮点寄存器传递；但WSL1底层依赖Windows内核，而Windows的x86_64 ABI仅支持将long double作为64位double处理（Windows不支持x87寄存器跨调用边界的正确传递）。

当代码分文件编译时，编译器会遵循WSL1的ABI规则：跨编译单元传递long double参数（比如main.cc调用SplitLD构造函数）时，会自动将80位long double降级为64位double，直接导致精度丢失。而单文件编译时，编译器可进行全程序优化，直接在x87寄存器中处理long double，无需跨调用边界传递，因此不会触发降级，精度得以保留。

程序是否存在未定义行为？ #

你的代码本身没有未定义行为：

• 静态断言确保std::uintmax_t能容纳long double的所有有效位；
• std::frexp和std::ldexp的使用符合标准规范；
• 类型转换操作均为标准允许的显式/隐式转换。

问题完全来自WSL1的ABI与标准Linux ABI不一致，导致跨编译单元的long double处理异常。

解决方案 #

1. 切换到WSL2：WSL2基于完整Linux内核，完全遵循x86_64 Linux ABI，不会出现long double精度丢失问题；
2. 强制使用64位long double：编译时添加-mlong-double-64参数，让long double等同于double，虽失去80位精度，但可避免跨单元ABI问题；
3. 保持单文件编译：通过单文件编译让编译器优化掉跨单元调用的ABI兼容处理，保留80位精度。

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