log1p_expm1_test.cpp 3.1 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293
  1. // Copyright John Maddock 2005.
  2. // Copyright Paul A. Bristow 2010
  3. // Use, modification and distribution are subject to the
  4. // Boost Software License, Version 1.0. (See accompanying file
  5. // LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
  6. #include <pch_light.hpp>
  7. // Requires MS extensions permitted or fails to link.
  8. #define BOOST_MATH_OVERFLOW_ERROR_POLICY ignore_error
  9. #include <boost/math/concepts/real_concept.hpp>
  10. #define BOOST_TEST_MAIN
  11. #include <boost/test/unit_test.hpp>
  12. #include <boost/test/tools/floating_point_comparison.hpp>
  13. #include <boost/math/special_functions/math_fwd.hpp>
  14. #include "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. 4, // Max Peek error
  51. 3); // 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. BOOST_MATH_CONTROL_FP;
  63. test(float(0), "float");
  64. test(double(0), "double");
  65. //
  66. // The long double version of these tests fails on some platforms
  67. // due to poor std lib support (not enough digits returned from
  68. // std::log and std::exp):
  69. //
  70. #if !defined(BOOST_MATH_NO_LONG_DOUBLE_MATH_FUNCTIONS)
  71. test((long double)(0), "long double");
  72. #ifndef BOOST_MATH_NO_REAL_CONCEPT_TESTS
  73. test((boost::math::concepts::real_concept)(0), "real_concept");
  74. #endif
  75. #else
  76. std::cout << "<note>The long double tests have been disabled on this platform "
  77. "either because the long double overloads of the usual math functions are "
  78. "not available at all, or because they are too inaccurate for these tests "
  79. "to pass.</note>" << std::endl;
  80. #endif
  81. }