From 483eb2f56657e8e7f419ab1a4fab8dce9ade8609 Mon Sep 17 00:00:00 2001 From: Daniel Baumann Date: Sat, 27 Apr 2024 20:24:20 +0200 Subject: Adding upstream version 14.2.21. Signed-off-by: Daniel Baumann --- src/boost/libs/math/test/test_hermite.hpp | 106 ++++++++++++++++++++++++++++++ 1 file changed, 106 insertions(+) create mode 100644 src/boost/libs/math/test/test_hermite.hpp (limited to 'src/boost/libs/math/test/test_hermite.hpp') diff --git a/src/boost/libs/math/test/test_hermite.hpp b/src/boost/libs/math/test/test_hermite.hpp new file mode 100644 index 00000000..0b00677e --- /dev/null +++ b/src/boost/libs/math/test/test_hermite.hpp @@ -0,0 +1,106 @@ +// Copyright John Maddock 2006. +// Copyright Paul A. Bristow 2007, 2009 +// 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) + +#ifdef _MSC_VER +# pragma warning(disable : 4127) // conditional expression is constant +# pragma warning(disable : 4512) // assignment operator could not be generated +# pragma warning(disable : 4756) // overflow in constant arithmetic +// Constants are too big for float case, but this doesn't matter for test. +#endif + +#include +#define BOOST_TEST_MAIN +#include +#include +#include +#include +#include +#include "functor.hpp" + +#include "handle_test_result.hpp" +#include "table_type.hpp" + +#ifndef SC_ +#define SC_(x) static_cast::type>(BOOST_JOIN(x, L)) +#endif + +template +void do_test_hermite(const T& data, const char* type_name, const char* test_name) +{ +#if !(defined(ERROR_REPORTING_MODE) && !defined(HERMITE_FUNCTION_TO_TEST)) + typedef Real value_type; + + typedef value_type (*pg)(unsigned, value_type); +#ifdef HERMITE_FUNCTION_TO_TEST + pg funcp = HERMITE_FUNCTION_TO_TEST; +#elif defined(BOOST_MATH_NO_DEDUCED_FUNCTION_POINTERS) + pg funcp = boost::math::hermite; +#else + pg funcp = boost::math::hermite; +#endif + + boost::math::tools::test_result result; + + std::cout << "Testing " << test_name << " with type " << type_name + << "\n~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\n"; + + // + // test hermite against data: + // + result = boost::math::tools::test_hetero( + data, + bind_func_int1(funcp, 0, 1), + extract_result(2)); + handle_test_result(result, data[result.worst()], result.worst(), type_name, "hermite", test_name); + + std::cout << std::endl; +#endif +} + +template +void test_hermite(T, const char* name) +{ + // + // The actual test data is rather verbose, so it's in a separate file + // + // The contents are as follows, each row of data contains + // three items, input value a, input value b and erf(a, b): + // +# include "hermite.ipp" + + do_test_hermite(hermite, name, "Hermite Polynomials"); +} + +template +void test_spots(T, const char* t) +{ + std::cout << "Testing basic sanity checks for type " << t << std::endl; + // + // basic sanity checks, tolerance is 100 epsilon: + // These spots were generated by MathCAD, precision is + // 14-16 digits. + // + T tolerance = (std::max)(boost::math::tools::epsilon() * 100, static_cast(1e-14)); + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(0, static_cast(1)), static_cast(1.L), tolerance); + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(1, static_cast(1)), static_cast(2.L), tolerance); + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(1, static_cast(2)), static_cast(4.L), tolerance); + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(1, static_cast(10)), static_cast(20), tolerance); + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(1, static_cast(100)), static_cast(200), tolerance); + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(1, static_cast(1e6)), static_cast(2e6), tolerance); + if(std::numeric_limits::max_exponent >= std::numeric_limits::max_exponent) + { + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(1, static_cast(1e307)), static_cast(2e307), tolerance); + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(99, static_cast(100)), static_cast(4.967223743011310E+227L), tolerance); + } + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(10, static_cast(30)), static_cast(5.896624628001300E+17L), tolerance); + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(10, static_cast(1000)), static_cast(1.023976960161280E+33L), tolerance); + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(10, static_cast(10)), static_cast(8.093278209760000E+12L), tolerance); + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(10, static_cast(-10)), static_cast(8.093278209760000E+12L), tolerance); + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(3, static_cast(-10)), static_cast(-7.880000000000000E+3L), tolerance); + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(3, static_cast(-1000)), static_cast(-7.999988000000000E+9L), tolerance); + BOOST_CHECK_CLOSE_FRACTION(::boost::math::hermite(3, static_cast(-1000000)), static_cast(-7.999999999988000E+18L), tolerance); +} + -- cgit v1.2.3