基于IEEE双精度标准的精准向上取整实现方案问询

取整无疑是被频繁问及与解决的问题之一。若预设小数点后精度（如10^N），可轻松实现如下代码：

double round(double input) {
    static constexpr double DIVISOR = std::pow(10, N); 
    return std::round(value * DIVISOR) / DIVISOR;
}

但该方法可能得到类似15.1499999999999的结果，而预期值为15.15；这是因为该结果是最接近预期值的双精度浮点数。但我的应用要求结果必须严格大于15.15。

最简单的实现方式是将上述代码中的std::round()替换为std::ceil()，但需选定N值。基于IEEE双精度浮点数标准，是否有更优雅的方式实现向上取整，确保结果始终为类似15.1500000000000001的数值并达到最高精度？

我了解到gtest的ASSERT_DOUBLE_EQ规则：

验证两个双精度浮点数val1与val2近似相等，误差不超过4个ULP。

但不确定在4个ULP范围内能否确保得到上述预期结果。

解决方案 #

核心思路：基于ULP的精确调整 #

要实现严格大于目标十进制值（如15.15）的最高精度双精度浮点数，核心是找到刚好大于目标值的最小双精度浮点数，而非依赖固定的10^N缩放逻辑。

具体实现方案 #

方案1：利用标准库函数直接实现

C++11及以上标准提供了std::nextafter函数（定义在<cmath>中），可以直接获取某个浮点数的下一个可表示值。结合精度缩放逻辑，代码如下：

#include <cmath>

double ceil_to_strictly_greater(double input, int decimal_places) {
    double divisor = std::pow(10, decimal_places);
    // 先按指定精度向上取整得到目标值
    double target = std::ceil(input * divisor) / divisor;
    // 获取刚好大于target的最小双精度浮点数
    return std::nextafter(target, INFINITY);
}

方案2：手动操作浮点位模式（理解原理用）

如果需要理解底层逻辑，可以直接操作双精度浮点数的二进制位模式：

#include <cstdint>
#include <cmath>

double next_after(double x) {
    // 将double的位模式转换为无符号整数
    uint64_t bits = *reinterpret_cast<uint64_t*>(&x);
    // 正数场景下，位模式加1即可得到下一个更大的浮点数
    if (x >= 0) {
        bits += 1;
    } else {
        bits -= 1;
    }
    // 转换回double类型
    return *reinterpret_cast<double*>(&bits);
}

double ceil_to_strictly_greater(double input, int decimal_places) {
    double divisor = std::pow(10, decimal_places);
    double target = std::ceil(input * divisor) / divisor;
    return next_after(target);
}

关于ULP验证的说明 #

• 上述方法得到的结果与目标值的差值恰好是1个ULP，远小于gtestASSERT_DOUBLE_EQ允许的4个ULP误差范围，因此用该断言验证时会判定为近似相等。
• 同时该结果是严格大于理论十进制目标值的最小双精度数，完全满足你的应用要求。

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