I tried the following:
#include <cmath>
#include <iostream>
#include <boost/multiprecision/cpp_bin_float.hpp>
using float256 = boost::multiprecision::number<boost::multiprecision::backends::cpp_bin_float<192,
boost::multiprecision::backends::digit_base_2, void, long long int>, boost::multiprecision::et_off>;
void Test()
{
double a = exp(-0.5);
double b = boost::multiprecision::exp(float256(-0.5)).convert_to<double>();
cout << a << endl << b << endl;
}
Unfortunately this doesn't work. I am getting 0.606531 for a but 1 for b.
I'm running: boost 1.58, Linux, gcc 5.4.0.
Any idea what I am doing wrong?