I am using Boost.Multiprecision for wrappers around the mpfr backend, and I am having some trouble creating pi (and e or any other math constant) to my desired precision. I feel like what I want to do should be possible, because of the use of Boost.Math for constants on a tutorial page for Boost.Multiprecision. In the tutorial, they use fixed-precision numbers of types such as cpp_dec_float_50
-- I want to do it with variable_precision mpfr_float
. Check out the following code:
#include <boost/multiprecision/mpfr.hpp>
#include <boost/math/constants/constants.hpp>
#include <iostream>
...
int main() {
boost::multiprecision::mpfr_float::default_precision(1000);
boost::multiprecision::mpfr_float pi = boost::math::constants::pi<boost::multiprecision::mpfr_float>();
std::cout << std::fixed;
std::cout.precision(1000);
std::cout << pi.precision() << " " << pi << std::endl;
}
The result is a number, pi
, which has precision 1000, but only has ~165 digits of pi, as evidenced by the result of the output statements, which confirm that my pi
is of precision 1000, and prints about 165 correct digits, and ~835 zeros. this is clearly wrong.
is it possible to make a boost::multiprecision::mpfr_float at high precision populated with constants from boost::math::constants?
Note that I need to use a variable precision type, and that other high-but-fixed-precision types are not an option. I must be able to change precision on the fly at runtime.