// Boost pow_test.cpp test file
// Tests the pow function

//  (C) Copyright Bruno Lalande 2008.
//  Distributed under the Boost Software License, Version 1.0.
//  (See accompanying file LICENSE_1_0.txt or copy at
//  http://www.boost.org/LICENSE_1_0.txt)


#include <cmath>
#include <iostream>
#include <boost/test/included/test_exec_monitor.hpp>
#include <boost/test/floating_point_comparison.hpp>
#include <boost/math/pow.hpp>


template <int N, class T>
void test_pow(T base)
{
    if (base == 0 && N < 0) return;

    BOOST_CHECK_CLOSE(::boost::math::pow<N>(base),
                      ::pow(base, static_cast<double>(N)),
                      0.001);
}

template <int N, class T>
void test_with_big_bases()
{
    for (T base = T(); base < T(1000); ++base)
        test_pow<N>(base);
}

template <int N>
void test_with_small_bases()
{
    double base = 0.9;
    for (int i = 0; i < 10; ++i)
    {
        base += base/50;
        test_pow<N>(base);
    }
}

template <class T, int Factor>
void test_with_small_exponents()
{
    test_with_big_bases<0, T>();
    test_with_big_bases<Factor*1, T>();
    test_with_big_bases<Factor*2, T>();
    test_with_big_bases<Factor*3, T>();
    test_with_big_bases<Factor*5, T>();
    test_with_big_bases<Factor*6, T>();
    test_with_big_bases<Factor*7, T>();
    test_with_big_bases<Factor*8, T>();
    test_with_big_bases<Factor*9, T>();
    test_with_big_bases<Factor*10, T>();
    test_with_big_bases<Factor*11, T>();
    test_with_big_bases<Factor*12, T>();
}

template <int Factor>
void test_with_big_exponents()
{
    test_with_small_bases<Factor*50>();
    test_with_small_bases<Factor*100>();
    test_with_small_bases<Factor*150>();
    test_with_small_bases<Factor*200>();
    test_with_small_bases<Factor*250>();
    test_with_small_bases<Factor*300>();
    test_with_small_bases<Factor*350>();
    test_with_small_bases<Factor*400>();
    test_with_small_bases<Factor*450>();
    test_with_small_bases<Factor*500>();
}


int test_main(int, char* [])
{
    using namespace std;

    cout << "Testing with integral bases and positive small exponents" << endl;
    test_with_small_exponents<int, 1>();
    cout << "Testing with integral bases and negative small exponents" << endl;
    test_with_small_exponents<int, -1>();

    cout << "Testing with double precision bases and positive small exponents" << endl;
    test_with_small_exponents<double, 1>();
    cout << "Testing with double precision bases and negative small exponents" << endl;
    test_with_small_exponents<double, -1>();

    cout << "Testing with positive big exponents" << endl;
    test_with_big_exponents<1>();
    cout << "Testing with negative big exponents" << endl;
    test_with_big_exponents<-1>();

    return 0;
}