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


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

const double errorMargin = 0.00001;

template<class T>
void test_power_int(double b, int e)
{
    const T xb(b), xe(e);
    
    using namespace std;
    cout << "pow(" << b << ", " << e << ") = " << pow(xb, e)
         << (double(pow(xb, e)) - pow(b, e) < errorMargin ? "" : " - Unexpected!") << endl;
}

template<class T>
void test_power_float(double b, double e)
{
    const T xb(b), xe(e);
    
    pow(xb, xe);
    using namespace std;
    cout << "pow(" << b << ", " << e << ") = " << pow(xb, xe)
         << (double(pow(xb, xe)) - pow(b, e) < errorMargin ? "" : " - Unexpected!") << endl;
}

void test_boost_mp()
{
    using namespace boost::multiprecision;
    
using T = cpp_dec_float_50;
assert( pow(T(0), T(2)) == 0 ); // Failure
    test_power_int<checked_cpp_int>(0, 0);
    test_power_int<checked_cpp_int>(0, 1);
    test_power_int<checked_cpp_int>(0, 2);
    test_power_int<checked_cpp_int>(1, 0);
    test_power_int<checked_cpp_int>(1, 1);
    test_power_int<checked_cpp_int>(1, 2);
    test_power_int<checked_cpp_int>(2, 0);
    test_power_int<checked_cpp_int>(2, 1);
    test_power_int<checked_cpp_int>(2, 2);
    test_power_float<cpp_dec_float_50>(0, 0);
    test_power_float<cpp_dec_float_50>(0, 1);
    test_power_float<cpp_dec_float_50>(0, 2);
    test_power_float<cpp_dec_float_50>(1, 0);
    test_power_float<cpp_dec_float_50>(1, 1);
    test_power_float<cpp_dec_float_50>(1, 2);
    test_power_float<cpp_dec_float_50>(2, 0);
    test_power_float<cpp_dec_float_50>(2, 1);
    test_power_float<cpp_dec_float_50>(2, 2);
    test_power_float<cpp_dec_float_100>(0, 0);
    test_power_float<cpp_dec_float_100>(0, 1);
    test_power_float<cpp_dec_float_100>(0, 2);
    test_power_float<cpp_dec_float_100>(1, 0);
    test_power_float<cpp_dec_float_100>(1, 1);
    test_power_float<cpp_dec_float_100>(1, 2);
    test_power_float<cpp_dec_float_100>(2, 0);
    test_power_float<cpp_dec_float_100>(2, 1);
    test_power_float<cpp_dec_float_100>(2, 2);
    test_power_float<number<cpp_dec_float<0, long long>>>(0, 0);
    test_power_float<number<cpp_dec_float<0, long long>>>(0, 1);
    test_power_float<number<cpp_dec_float<0, long long>>>(0, 2);
    test_power_float<number<cpp_dec_float<0, long long>>>(1, 0);
    test_power_float<number<cpp_dec_float<0, long long>>>(1, 1);
    test_power_float<number<cpp_dec_float<0, long long>>>(1, 2);
    test_power_float<number<cpp_dec_float<0, long long>>>(2, 0);
    test_power_float<number<cpp_dec_float<0, long long>>>(2, 1);
    test_power_float<number<cpp_dec_float<0, long long>>>(2, 2);
}