Boost可变精度数输出不遵循流精度标志的技术问询

问题描述 #

将基于MPFI自研OOP封装的代码迁移至Boost Multiprecision时，发现输出可变精度类型mpfi_float时会忽略输出流的精度参数——无论数值实际精度设置为多少，输出表现都与精度设为20时一致。若这是预期行为而非Bug，希望了解其设计依据。

补充背景 #

使用mpfi_float_100这类固定精度类型时，输出区间会显示流指定精度的所有数字；但将可变精度类型mpfi_float设为相同的100位精度时，输出区间仅显示少量数字，似乎忽略了流精度。值得注意的是，直接输出上下界（通过lower()/upper()）时会显示所有预期数字，与固定精度情况一致（精度低于20时除外，这似乎是可变精度mpfi_float的最低精度）。

示例代码 #

#include <boost/multiprecision/mpfi.hpp>
#include <iostream>

int main(int argc, char** argv)
{
    constexpr unsigned TESTPREC = 100;

    using namespace std;
    using namespace boost::multiprecision;
    cout.precision(5000);

    // 固定精度数输出时始终显示最大位数
    number<mpfi_float_backend<TESTPREC>> a = 2;
    number<mpfi_float_backend<TESTPREC>> root_a = sqrt(a);
    cout << root_a << endl; // 分别输出332和334位的上下界
    cout << lower(root_a) << endl; // 输出332位下界
    cout << upper(root_a) << endl; // 输出334位上界
    cout << endl;

    // 可变精度数输出可能不显示全部位数
    mpfi_float b = 2;
    mpfi_float::default_precision(TESTPREC);
    mpfi_float root_b = sqrt(b);
    cout << root_b << endl; // 仅输出68和67位的上下界
    cout << lower(root_b) << endl; // 输出332位下界
    cout << upper(root_b) << endl; // 输出334位上界
    cout << endl;

    return 0;
}

测试结果汇总 #

构建环境 #

• Windows 11 Pro（Intel Core i9-11900H）
• Visual Studio Community 22（17.10.1 64位）
  • 编译标志：/std:c++20、/permissive-、/fp:strict
• Boost 1.85.0（x64-windows）
• GNU MPFR 4.2.1（x64-windows）

解答：这是预期行为，设计依据如下 #

1. 可变精度类型的输出精度策略 #

mpfi_float作为可变精度区间类型，其输出流运算符采用了保守的默认显示策略：默认以20位十进制精度作为输出基准（对应测试结果中的19位有效位数，因区间上下界的差异需要额外位数区分），而非直接使用流的精度参数。这一设计是为了避免在未显式指定时，输出过多冗余数字——毕竟可变精度类型的实际精度可能动态调整，默认输出全精度可能导致不必要的性能开销和冗长输出。

而固定精度类型（如mpfi_float_100）的精度是编译期固定的，输出逻辑会直接绑定该精度值，因此会遵循流的精度设置，展示对应位数的数字。

2. 区间输出与上下界输出的逻辑差异 #

直接输出mpfi_float对象时，调用的是针对区间类型的格式化逻辑，该逻辑的核心是仅显示足够区分上下界的最小有效位数，而非全部存储精度。而调用lower()/upper()返回的是mpfr_float（或其可变精度版本），这类单精度浮点类型的输出逻辑会严格遵守流的精度参数，因此能显示全精度数字。

3. 实现全精度输出的方法 #

若需要让mpfi_float输出全精度，有两种可行方式：

• 使用精度操纵器：通过boost::multiprecision::set_precision显式指定输出精度，覆盖默认策略：

  cout << boost::multiprecision::set_precision(5000) << root_b << endl;
  

• 手动格式化上下界：直接输出lower(root_b)和upper(root_b)，如示例代码所示，这类输出会遵循流的精度设置。

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