如何借助Boost多精度库精准将双精度数按指定精度取整？

需求 #

输入一个double类型数值和一个整数精度值，精准将该double数按指定精度取整，最终输出精度更可靠的double类型结果。

现有实现代码 #

用户基于Boost多精度库编写了以下代码（已补充缺失的分号与必要头文件）：

#include <boost/multiprecision/cpp_dec_float.hpp>
#include <cmath>
#include <limits>
#include <iostream>
#include <iomanip>

int main() {
    double input = 7.7778;
    int precision = 3;

    using my_dec_100 = boost::multiprecision::cpp_dec_float_50;

    double intpart, fracpart;
    fracpart = modf(input, &intpart);
    int decimalLength = 0;
    if (fracpart > 0){
        decimalLength = std::ceil(-std::log10(std::pow(10, -std::numeric_limits<double>::digits10) * fracpart));
    }

    double rounded = std::floor(input * std::pow(10, precision) + 0.5) / std::pow(10, precision);
    my_dec_100 f11(my_dec_100(rounded * std::pow(10, decimalLength)) / std::pow(10, decimalLength));
    double output = f11.convert_to<double>();
    std::cout << std::setprecision(precision) << std::fixed << output << std::endl;  
    return 0;
}

用户疑问 #

1. 不确定上述代码是否能适配所有数值场景，比如包含超长小数部分的数值（如7.1346745428482672454818495）。
2. 当前仅用Boost库创建了f11对象做处理，想了解是否有其他基于Boost多精度库的实现方法。

解答 #

现有代码的局限性 #

现有代码依赖std::pow和std::floor组合实现取整，但double的二进制存储特性会导致部分十进制小数无法精确表示（比如0.1），用pow(10, precision)计算时易引入额外误差；面对超长小数输入时，这种误差会被放大，导致取整结果不准确。另外，计算decimalLength的逻辑复杂冗余，没必要绕此弯路。

基于Boost多精度库的更简洁实现方法 #

直接利用Boost高精度浮点类型处理全流程，避免原生double的精度损失，步骤如下：

1. 将输入double转换为Boost高精度浮点类型，保留完整数值信息；
2. 用Boost提供的取整函数或精度设置完成指定精度的取整；
3. 将处理后的高精度数值转换回double输出。

示例代码：

#include <boost/multiprecision/cpp_dec_float.hpp>
#include <iostream>
#include <iomanip>

int main() {
    double input = 7.1346745428482672454818495;
    int precision = 3;

    // 定义50位精度的浮点类型，足够覆盖绝大多数场景
    using DecimalFloat = boost::multiprecision::cpp_dec_float_50;

    // 转换为高精度浮点，规避原生double的精度丢失
    DecimalFloat high_precision_input(input);

    // 计算10^precision的缩放因子
    DecimalFloat scale = boost::multiprecision::pow(DecimalFloat(10), precision);

    // 四舍五入：放大后取整，再缩小回原精度
    DecimalFloat rounded_value = boost::multiprecision::round(high_precision_input * scale) / scale;

    // 转换回double类型
    double output = rounded_value.convert_to<double>();

    // 按指定精度输出
    std::cout << std::setprecision(precision) << std::fixed << output << std::endl;
    return 0;
}

适配超长小数的说明 #

像7.1346745428482672454818495这类超长小数，原生double本身最多只能存储约15-17位有效数字，输入时就会丢失部分精度。但用Boost多精度处理时，能基于double提供的最大有效数字精准取整，避免后续计算额外引入误差，比原生方法更可靠。

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