//
//   g++ -o jm.cpp -I/JM/SOFT/boost_1_77_0  jm.cpp -lquadmath
//

#include <boost/multiprecision/float128.hpp>

#include <iostream>

using boost::multiprecision::float128 ;

int main( )
{
    std::cout.setf( std::ios_base::scientific ) ;
    std::cout.setf( std::ios_base::showpoint  ) ;  // include any trailing zeros

    double      double_pi =  3.141592653589793  ;

    float128  float128_pi =  3.141592653589793Q ;

    std::cout << "  double_pi = " << std::endl << std::setprecision(   std::numeric_limits<double>::digits10 ) <<    double_pi  << std::endl << std::endl ;
    std::cout << "float128_pi = " << std::endl << std::setprecision( std::numeric_limits<float128>::digits10 ) <<  float128_pi  << std::endl << std::endl ;

    float128_pi = double_pi ;

    std::cout << "float128_pi = " << std::endl << std::setprecision( std::numeric_limits<float128>::digits10 ) <<  float128_pi  << std::endl << std::endl ;

    std::cout << "sizeof   double_pi = " << sizeof(   double_pi ) << std::endl              ;
    std::cout << "sizeof float128_pi = " << sizeof( float128_pi ) << std::endl << std::endl ;
}