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

typedef boost::multiprecision::number<boost::multiprecision::backends::cpp_bin_float<128, boost::multiprecision::backends::digit_base_2> > ext_float_t;

int main(int, char**) {
  ext_float_t a0("1e500000000");
  ext_float_t a = a0*a0;  // +inf
  ext_float_t b = a*0;   // nan

  std::cout
    << a << " => " << a.convert_to<double>() << "\n"
    << b << " => " << b.convert_to<double>() << "\n"
    ;

  return 0;
}