//
//  boost_multiprecision_power_demo.cpp
//
//  Created by Jan Bouwer on 9/7/13.


#include <cassert>
#include <cmath>
#include <iostream>
#include <boost/multiprecision/cpp_dec_float.hpp>

template<class T>
void test_power_float(double e)
{
    const T xb(-1), xe(e);
    assert( xe > 0 );
    
    std::cout << "pow(-1, " << e << ") = " << pow(xb, xe) << std::endl;
}

void test_boost_mp()
{
    using namespace boost::multiprecision;
    
    test_power_float<cpp_dec_float_50>(0x8000000000000000u);
    test_power_float<cpp_dec_float_50>(0x8000000000000001u);
    test_power_float<cpp_dec_float_50>(0xFFFFFFFFFFFFFFFEu);
    test_power_float<cpp_dec_float_50>(0xFFFFFFFFFFFFFFFFu);
    test_power_float<cpp_dec_float_100>(0x8000000000000000u);
    test_power_float<cpp_dec_float_100>(0x8000000000000001u);
    test_power_float<cpp_dec_float_100>(0xFFFFFFFFFFFFFFFEu);
    test_power_float<cpp_dec_float_100>(0xFFFFFFFFFFFFFFFFu);
}