diff options
Diffstat (limited to 'src/boost/libs/math/test/test_expint.cpp')
| -rw-r--r-- | src/boost/libs/math/test/test_expint.cpp | 147 |
1 files changed, 147 insertions, 0 deletions
diff --git a/src/boost/libs/math/test/test_expint.cpp b/src/boost/libs/math/test/test_expint.cpp new file mode 100644 index 000000000..3f44a8091 --- /dev/null +++ b/src/boost/libs/math/test/test_expint.cpp @@ -0,0 +1,147 @@ +// (C) Copyright John Maddock 2007. +// Use, modification and distribution are subject to 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 <pch_light.hpp> +#include "test_expint.hpp" + +// +// DESCRIPTION: +// ~~~~~~~~~~~~ +// +// This file tests the expint functions. There are two sets of tests: +// 1) Sanity checks: comparison to test values created with the +// online calculator at functions.wolfram.com +// 2) Accuracy tests use values generated with NTL::RR at +// 1000-bit precision and our generic versions of these functions. +// +// Note that when this file is first run on a new platform many of +// these tests will fail: the default accuracy is 1 epsilon which +// is too tight for most platforms. In this situation you will +// need to cast a human eye over the error rates reported and make +// a judgement as to whether they are acceptable. Either way please +// report the results to the Boost mailing list. Acceptable rates of +// error are marked up below as a series of regular expressions that +// identify the compiler/stdlib/platform/data-type/test-data/test-function +// along with the maximum expected peek and RMS mean errors for that +// test. +// + +void expected_results() +{ + // + // Define the max and mean errors expected for + // various compilers and platforms. + // + + // + // On MacOS X erfc has much higher error levels than + // expected: given that the implementation is basically + // just a rational function evaluation combined with + // exponentiation, we conclude that exp and pow are less + // accurate on this platform, especially when the result + // is outside the range of a double. + // + add_expected_result( + ".*", // compiler + ".*", // stdlib + "Mac OS", // platform + "float|double|long double", // test type(s) + ".*E1.*", // test data group + ".*", 30, 10); // test function + add_expected_result( + ".*", // compiler + ".*", // stdlib + "Mac OS", // platform + "float|double|long double|real_concept", // test type(s) + ".*Ei.*", // test data group + ".*", 300, 200); // test function + add_expected_result( + ".*", // compiler + ".*", // stdlib + "Mac OS", // platform + ".*", // test type(s) + ".*", // test data group + ".*", 40, 15); // test function + + add_expected_result( + ".*", // compiler + ".*", // stdlib + ".*", // platform + "float|double|long double", // test type(s) + ".*E1.*", // test data group + ".*", 2, 1); // test function + add_expected_result( + ".*", // compiler + ".*", // stdlib + ".*", // platform + "float|double|long double", // test type(s) + ".*Ei.*", // test data group + ".*", 6, 3); // test function + if(std::numeric_limits<long double>::digits > 100) + { + add_expected_result( + ".*", // compiler + ".*", // stdlib + ".*", // platform + "real_concept", // test type(s) + ".*Ei.*", // test data group + ".*", 150, 50); // test function + } + add_expected_result( + ".*", // compiler + ".*", // stdlib + ".*", // platform + "real_concept", // test type(s) + ".*Ei.*", // test data group + ".*", 150, 50); // test function + add_expected_result( + ".*", // compiler + ".*", // stdlib + ".*", // platform + ".*", // test type(s) + ".*", // test data group + ".*", 25, 5); // test function + + std::cout << "Tests run with " << BOOST_COMPILER << ", " + << BOOST_STDLIB << ", " << BOOST_PLATFORM << std::endl; +} + +BOOST_AUTO_TEST_CASE( test_main ) +{ + expected_results(); + BOOST_MATH_CONTROL_FP; + + boost::math::expint(114.7); + + test_spots(0.0f, "float"); + test_spots(0.0, "double"); +#ifndef BOOST_MATH_NO_LONG_DOUBLE_MATH_FUNCTIONS + test_spots(0.0L, "long double"); +#ifndef BOOST_MATH_NO_REAL_CONCEPT_TESTS + test_spots(boost::math::concepts::real_concept(0.1), "real_concept"); +#endif +#else + std::cout << "<note>The long double tests have been disabled on this platform " + "either because the long double overloads of the usual math functions are " + "not available at all, or because they are too inaccurate for these tests " + "to pass.</note>" << std::endl; +#endif + + test_expint(0.1F, "float"); + test_expint(0.1, "double"); +#ifndef BOOST_MATH_NO_LONG_DOUBLE_MATH_FUNCTIONS + test_expint(0.1L, "long double"); +#ifndef BOOST_MATH_NO_REAL_CONCEPT_TESTS + test_expint(boost::math::concepts::real_concept(0.1), "real_concept"); +#endif +#else + std::cout << "<note>The long double tests have been disabled on this platform " + "either because the long double overloads of the usual math functions are " + "not available at all, or because they are too inaccurate for these tests " + "to pass.</note>" << std::endl; +#endif +} + + |