log1p_expm1_test.cpp 2.7 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879
  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. #include "setup.hpp"
  11. #include <boost/math/special_functions/log1p.hpp>
  12. #include <boost/math/special_functions/expm1.hpp>
  13. #include "table_type.hpp"
  14. #include "libs/math/test/log1p_expm1_test.hpp"
  15. //
  16. // DESCRIPTION:
  17. // ~~~~~~~~~~~~
  18. //
  19. // This file tests the functions log1p and expm1. The accuracy tests
  20. // use values generated with NTL::RR at 1000-bit precision
  21. // and our generic versions of these functions.
  22. //
  23. // Note that when this file is first run on a new platform many of
  24. // these tests will fail: the default accuracy is 1 epsilon which
  25. // is too tight for most platforms. In this situation you will
  26. // need to cast a human eye over the error rates reported and make
  27. // a judgement as to whether they are acceptable. Either way please
  28. // report the results to the Boost mailing list. Acceptable rates of
  29. // error are marked up below as a series of regular expressions that
  30. // identify the compiler/stdlib/platform/data-type/test-data/test-function
  31. // along with the maximum expected peek and RMS mean errors for that
  32. // test.
  33. //
  34. void expected_results()
  35. {
  36. //
  37. // Define the max and mean errors expected for
  38. // various compilers and platforms.
  39. //
  40. //
  41. // Catch all cases come last:
  42. //
  43. add_expected_result(
  44. ".*", // compiler
  45. ".*", // stdlib
  46. ".*", // platform
  47. ".*", // test type(s)
  48. ".*", // test data group
  49. ".*", // test function
  50. 8, // Max Peek error
  51. 5); // Max mean error
  52. //
  53. // Finish off by printing out the compiler/stdlib/platform names,
  54. // we do this to make it easier to mark up expected error rates.
  55. //
  56. std::cout << "Tests run with " << BOOST_COMPILER << ", "
  57. << BOOST_STDLIB << ", " << BOOST_PLATFORM << std::endl;
  58. }
  59. BOOST_AUTO_TEST_CASE( test_main )
  60. {
  61. expected_results();
  62. //
  63. // Test at:
  64. // 18 decimal digits: tests 80-bit long double approximations
  65. // 30 decimal digits: tests 128-bit long double approximations
  66. // 35 decimal digits: tests arbitrary precision code
  67. //
  68. ALL_TESTS
  69. }