libs/math/test/float128/powm1_sqrtp1m1_test.cpp

   1 ///////////////////////////////////////////////////////////////
   2 //  Copyright Christopher Kormanyos 2002 - 2011.
   3 //  Copyright 2011 John Maddock. Distributed under the Boost
   4 //  Software License, Version 1.0. (See accompanying file
   5 //  LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_
   6 //
   7 // This work is based on an earlier work:
   8 // "Algorithm 910: A Portable C++ Multiple-Precision System for Special-Function Calculations",
   9 // in ACM TOMS, {VOL 37, ISSUE 4, (February 2011)} (C) ACM, 2011. http://doi.acm.org/10.1145/1916461.1916469
  10
  11 #include "setup.hpp"
  12
  13 #include <boost/test/tools/floating_point_comparison.hpp>
  14 #include <boost/math/special_functions/powm1.hpp>
  15 #include <boost/math/special_functions/sqrt1pm1.hpp>
  16
  17 #include "table_type.hpp"
  18
  19 #include "libs/math/test/powm1_sqrtp1m1_test.hpp"
  20
  21 //
  22 // DESCRIPTION:
  23 // ~~~~~~~~~~~~
  24 //
  25 // This file tests the functions log1p and expm1.  The accuracy tests
  26 // use values generated with NTL::RR at 1000-bit precision
  27 // and our generic versions of these functions.
  28 //
  29 // Note that when this file is first run on a new platform many of
  30 // these tests will fail: the default accuracy is 1 epsilon which
  31 // is too tight for most platforms.  In this situation you will
  32 // need to cast a human eye over the error rates reported and make
  33 // a judgement as to whether they are acceptable.  Either way please
  34 // report the results to the Boost mailing list.  Acceptable rates of
  35 // error are marked up below as a series of regular expressions that
  36 // identify the compiler/stdlib/platform/data-type/test-data/test-function
  37 // along with the maximum expected peek and RMS mean errors for that
  38 // test.
  39 //
  40
  41 void expected_results()
  42 {
  43    //
  44    // Define the max and mean errors expected for
  45    // various compilers and platforms.
  46    //
  47
  48    //
  49    // Catch all cases come last:
  50    //
  51    add_expected_result(
  52       ".*",                          // compiler
  53       ".*",                          // stdlib
  54       ".*",                          // platform
  55       ".*",                          // test type(s)
  56       ".*",                          // test data group
  57       ".*",                          // test function
  58       15,                             // Max Peek error
  59       5);                            // Max mean error
  60
  61    //
  62    // Finish off by printing out the compiler/stdlib/platform names,
  63    // we do this to make it easier to mark up expected error rates.
  64    //
  65    std::cout << "Tests run with " << BOOST_COMPILER << ", "
  66       << BOOST_STDLIB << ", " << BOOST_PLATFORM << std::endl;
  67 }
  68
  69 template <class T>
  70 void test(T t, const char* p)
  71 {
  72    test_powm1_sqrtp1m1(t, p);
  73 }
  74
  75 BOOST_AUTO_TEST_CASE( test_main )
  76 {
  77    expected_results();
  78    ALL_TESTS
  79 }
  80
  81
  82