101 lines
3.5 KiB
C++
101 lines
3.5 KiB
C++
///////////////////////////////////////////////////////////////
|
|
// Copyright Christopher Kormanyos 2002 - 2011.
|
|
// Copyright 2011 John Maddock. Distributed under the Boost
|
|
// Software License, Version 1.0. (See accompanying file
|
|
// LICENSE_1_0.txt or copy at https://www.boost.org/LICENSE_1_0.txt
|
|
//
|
|
// This work is based on an earlier work:
|
|
// "Algorithm 910: A Portable C++ Multiple-Precision System for Special-Function Calculations",
|
|
// in ACM TOMS, {VOL 37, ISSUE 4, (February 2011)} (C) ACM, 2011. http://doi.acm.org/10.1145/1916461.1916469
|
|
|
|
#include "setup.hpp"
|
|
//#include <boost/math/special_functions/log1p.hpp>
|
|
//#include <boost/math/special_functions/expm1.hpp>
|
|
#define BOOST_TEST_MAIN
|
|
#include <boost/test/unit_test.hpp>
|
|
#include <boost/test/floating_point_comparison.hpp>
|
|
#include <boost/math/special_functions/math_fwd.hpp>
|
|
|
|
#include "table_type.hpp"
|
|
|
|
#include "libs/math/test/log1p_expm1_test.hpp"
|
|
|
|
//
|
|
// DESCRIPTION:
|
|
// ~~~~~~~~~~~~
|
|
//
|
|
// This file tests the functions log1p and expm1. The 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.
|
|
//
|
|
|
|
//
|
|
// Catch all cases come last:
|
|
//
|
|
add_expected_result(
|
|
".*", // compiler
|
|
".*", // stdlib
|
|
".*", // platform
|
|
".*gmp_float<18>.*", // test type(s)
|
|
".*", // test data group
|
|
".*", // test function
|
|
500, // Max Peek error
|
|
100); // Max mean error
|
|
add_expected_result(
|
|
".*", // compiler
|
|
".*", // stdlib
|
|
".*", // platform
|
|
".*mpfr_float_backend<18>.*", // test type(s)
|
|
".*", // test data group
|
|
".*", // test function
|
|
500, // Max Peek error
|
|
100); // Max mean error
|
|
add_expected_result(
|
|
".*", // compiler
|
|
".*", // stdlib
|
|
".*", // platform
|
|
".*", // test type(s)
|
|
".*", // test data group
|
|
".*", // test function
|
|
8, // Max Peek error
|
|
5); // Max mean error
|
|
|
|
//
|
|
// Finish off by printing out the compiler/stdlib/platform names,
|
|
// we do this to make it easier to mark up expected error rates.
|
|
//
|
|
std::cout << "Tests run with " << BOOST_COMPILER << ", "
|
|
<< BOOST_STDLIB << ", " << BOOST_PLATFORM << std::endl;
|
|
}
|
|
|
|
BOOST_AUTO_TEST_CASE(test_main)
|
|
{
|
|
using namespace boost::multiprecision;
|
|
expected_results();
|
|
//
|
|
// Test at:
|
|
// 18 decimal digits: tests 80-bit long double approximations
|
|
// 30 decimal digits: tests 128-bit long double approximations
|
|
// 35 decimal digits: tests arbitrary precision code
|
|
//
|
|
ALL_TESTS
|
|
}
|