test_bessel_j.cpp 4.5 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119
  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 https://www.boost.org/LICENSE_1_0.txt
  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 "table_type.hpp"
  12. #include <boost/math/special_functions/math_fwd.hpp>
  13. #include "libs/math/test/test_bessel_j.hpp"
  14. void expected_results()
  15. {
  16. //
  17. // Define the max and mean errors expected for
  18. // various compilers and platforms.
  19. //
  20. add_expected_result(
  21. ".*", // compiler
  22. ".*", // stdlib
  23. ".*", // platform
  24. ".*mpfr_float_backend<18>.*", // test type(s)
  25. ".*J0.*Tricky.*", // test data group
  26. ".*", 10000000000L, 5000000000L); // test function
  27. add_expected_result(
  28. ".*", // compiler
  29. ".*", // stdlib
  30. ".*", // platform
  31. ".*", // test type(s)
  32. ".*J0.*Tricky.*", // test data group
  33. ".*", 400000000, 400000000); // test function
  34. add_expected_result(
  35. ".*", // compiler
  36. ".*", // stdlib
  37. ".*", // platform
  38. ".*mpfr_float_backend<18>.*", // test type(s)
  39. ".*J1.*Tricky.*", // test data group
  40. ".*", 1000000000L, 500000000L); // test function
  41. add_expected_result(
  42. ".*", // compiler
  43. ".*", // stdlib
  44. ".*", // platform
  45. ".*", // test type(s)
  46. ".*J1.*Tricky.*", // test data group
  47. ".*", 10000000, 5000000); // test function
  48. add_expected_result(
  49. ".*", // compiler
  50. ".*", // stdlib
  51. ".*", // platform
  52. ".*cpp_bin_float.*", // test type(s)
  53. ".*JN.*Integer.*", // test data group
  54. ".*", 500000, 150000); // test function
  55. add_expected_result(
  56. ".*", // compiler
  57. ".*", // stdlib
  58. ".*", // platform
  59. ".*", // test type(s)
  60. ".*JN.*Integer.*", // test data group
  61. ".*", 50000, 15000); // test function
  62. add_expected_result(
  63. ".*", // compiler
  64. ".*", // stdlib
  65. ".*", // platform
  66. ".*mpfr_float_backend<18>.*", // test type(s)
  67. ".*(JN|j).*|.*Tricky.*", // test data group
  68. ".*", 70000, 7000); // test function
  69. add_expected_result(
  70. ".*", // compiler
  71. ".*", // stdlib
  72. ".*", // platform
  73. ".*cpp_bin_float.*", // test type(s)
  74. ".*(JN|j).*|.*Tricky.*", // test data group
  75. ".*", 500000, 200000); // test function
  76. add_expected_result(
  77. ".*", // compiler
  78. ".*", // stdlib
  79. ".*", // platform
  80. ".*", // test type(s)
  81. ".*(JN|j).*|.*Tricky.*", // test data group
  82. ".*", 7000, 3000); // test function
  83. add_expected_result(
  84. ".*", // compiler
  85. ".*", // stdlib
  86. ".*", // platform
  87. ".*", // test type(s)
  88. ".*", // test data group
  89. ".*", 40, 20); // test function
  90. //
  91. // Finish off by printing out the compiler/stdlib/platform names,
  92. // we do this to make it easier to mark up expected error rates.
  93. //
  94. std::cout << "Tests run with " << BOOST_COMPILER << ", "
  95. << BOOST_STDLIB << ", " << BOOST_PLATFORM << std::endl;
  96. }
  97. template <class T>
  98. void test(T t, const char* p)
  99. {
  100. test_bessel(t, p);
  101. }
  102. BOOST_AUTO_TEST_CASE(test_main)
  103. {
  104. using namespace boost::multiprecision;
  105. expected_results();
  106. //
  107. // Test at:
  108. // 18 decimal digits: tests 80-bit long double approximations
  109. // 30 decimal digits: tests 128-bit long double approximations
  110. // 35 decimal digits: tests arbitrary precision code
  111. //
  112. ALL_TESTS
  113. }