123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226 |
- // 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 <string>
- #include <iostream>
- #include <boost/math/concepts/real_concept.hpp>
- #include <boost/math/tools/test.hpp>
- #define BOOST_TEST_MAIN
- #include <boost/test/unit_test.hpp>
- #include <boost/test/tools/floating_point_comparison.hpp>
- #include <boost/typeof/typeof.hpp>
- #include <boost/type_traits/is_same.hpp>
- #include <boost/static_assert.hpp>
- #include <boost/math/special_functions/pow.hpp>
- #include BOOST_TYPEOF_INCREMENT_REGISTRATION_GROUP()
- BOOST_TYPEOF_REGISTER_TYPE(boost::math::concepts::real_concept)
- using namespace boost;
- using namespace boost::math;
- template <int N, class T>
- void test_pow(T base)
- {
- typedef typename tools::promote_args<T>::type result_type;
- BOOST_MATH_STD_USING
- if ((base == 0) && N < 0)
- {
- BOOST_MATH_CHECK_THROW(math::pow<N>(base), std::overflow_error);
- }
- else
- {
- BOOST_CHECK_CLOSE(math::pow<N>(base),
- pow(static_cast<result_type>(base), static_cast<result_type>(N)),
- boost::math::tools::epsilon<result_type>() * 100 * 400); // 400 eps as a %
- }
- }
- 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, class T>
- void test_with_small_bases()
- {
- T base = 0.9f;
- 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 <class T, int Factor>
- void test_with_big_exponents()
- {
- test_with_small_bases<Factor*50, T>();
- test_with_small_bases<Factor*100, T>();
- test_with_small_bases<Factor*150, T>();
- test_with_small_bases<Factor*200, T>();
- test_with_small_bases<Factor*250, T>();
- test_with_small_bases<Factor*300, T>();
- test_with_small_bases<Factor*350, T>();
- test_with_small_bases<Factor*400, T>();
- test_with_small_bases<Factor*450, T>();
- test_with_small_bases<Factor*500, T>();
- }
- void test_return_types()
- {
- BOOST_STATIC_ASSERT((is_same<BOOST_TYPEOF(pow<2>('\1')), double>::value));
- BOOST_STATIC_ASSERT((is_same<BOOST_TYPEOF(pow<2>(L'\2')), double>::value));
- BOOST_STATIC_ASSERT((is_same<BOOST_TYPEOF(pow<2>(3)), double>::value));
- BOOST_STATIC_ASSERT((is_same<BOOST_TYPEOF(pow<2>(4u)), double>::value));
- BOOST_STATIC_ASSERT((is_same<BOOST_TYPEOF(pow<2>(5ul)), double>::value));
- BOOST_STATIC_ASSERT((is_same<BOOST_TYPEOF(pow<2>(6.0f)), float>::value));
- BOOST_STATIC_ASSERT((is_same<BOOST_TYPEOF(pow<2>(7.0)), double>::value));
- #ifndef BOOST_MATH_NO_LONG_DOUBLE_MATH_FUNCTIONS
- BOOST_STATIC_ASSERT((is_same<BOOST_TYPEOF(pow<2>(7.0l)), long double>::value));
- #endif
- }
- namespace boost { namespace math { namespace policies {
- template <class T>
- T user_overflow_error(const char*, const char*, const T&)
- { return T(123.45); }
- }}}
- namespace boost { namespace math { namespace policies {
- template <class T>
- T user_indeterminate_result_error(const char*, const char*, const T&)
- { return T(456.78); }
- }}}
- void test_error_policy()
- {
- using namespace policies;
- BOOST_CHECK(pow<-2>(
- 0.0,
- policy< ::boost::math::policies::overflow_error<user_error> >()
- )
- == 123.45);
- BOOST_CHECK(pow<0>(
- 0.0,
- policy< ::boost::math::policies::indeterminate_result_error<user_error> >()
- )
- == 456.78);
- }
- BOOST_AUTO_TEST_CASE( test_main )
- {
- 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 float precision bases and positive small exponents" << endl;
- test_with_small_exponents<float, 1>();
- cout << "Testing with float precision bases and negative small exponents" << endl;
- test_with_small_exponents<float, -1>();
- cout << "Testing with float precision bases and positive big exponents" << endl;
- test_with_big_exponents<float, 1>();
- cout << "Testing with float precision bases and negative big exponents" << endl;
- test_with_big_exponents<float, -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 double precision bases and positive big exponents" << endl;
- test_with_big_exponents<double, 1>();
- cout << "Testing with double precision bases and negative big exponents" << endl;
- test_with_big_exponents<double, -1>();
- #ifndef BOOST_MATH_NO_LONG_DOUBLE_MATH_FUNCTIONS
- cout << "Testing with long double precision bases and positive small exponents" << endl;
- test_with_small_exponents<long double, 1>();
- cout << "Testing with long double precision bases and negative small exponents" << endl;
- test_with_small_exponents<long double, -1>();
- cout << "Testing with long double precision bases and positive big exponents" << endl;
- test_with_big_exponents<long double, 1>();
- cout << "Testing with long double precision bases and negative big exponents" << endl;
- test_with_big_exponents<long double, -1>();
- cout << "Testing with concepts::real_concept precision bases and positive small exponents" << endl;
- test_with_small_exponents<boost::math::concepts::real_concept, 1>();
- cout << "Testing with concepts::real_concept precision bases and negative small exponents" << endl;
- test_with_small_exponents<boost::math::concepts::real_concept, -1>();
- cout << "Testing with concepts::real_concept precision bases and positive big exponents" << endl;
- test_with_big_exponents<boost::math::concepts::real_concept, 1>();
- cout << "Testing with concepts::real_concept precision bases and negative big exponents" << endl;
- test_with_big_exponents<boost::math::concepts::real_concept, -1>();
- #endif
- test_return_types();
- test_error_policy();
- }
- /*
- Running 1 test case...
- Testing with integral bases and positive small exponents
- Testing with integral bases and negative small exponents
- Testing with float precision bases and positive small exponents
- Testing with float precision bases and negative small exponents
- Testing with float precision bases and positive big exponents
- Testing with float precision bases and negative big exponents
- Testing with double precision bases and positive small exponents
- Testing with double precision bases and negative small exponents
- Testing with double precision bases and positive big exponents
- Testing with double precision bases and negative big exponents
- Testing with long double precision bases and positive small exponents
- Testing with long double precision bases and negative small exponents
- Testing with long double precision bases and positive big exponents
- Testing with long double precision bases and negative big exponents
- Testing with concepts::real_concept precision bases and positive small exponents
- Testing with concepts::real_concept precision bases and negative small exponents
- Testing with concepts::real_concept precision bases and positive big exponents
- Testing with concepts::real_concept precision bases and negative big exponents
-
- *** No errors detected
- */
|