test_bessel_y.cpp 5.3 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139
  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_y.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. ".*(Y[nv]|y).*Random.*", // test data group
  26. ".*", 200000, 10000); // test function
  27. add_expected_result(
  28. ".*", // compiler
  29. ".*", // stdlib
  30. ".*", // platform
  31. ".*cpp_bin_float.*", // test type(s)
  32. ".*(Y[nv]|y).*Random.*", // test data group
  33. ".*", 2000000, 1000000); // test function
  34. add_expected_result(
  35. ".*", // compiler
  36. ".*", // stdlib
  37. ".*", // platform
  38. ".*mpfr_float_backend<0>.*", // test type(s)
  39. ".*(Y[nv]|y).*Random.*", // test data group
  40. ".*", 150000, 15000); // test function
  41. add_expected_result(
  42. ".*", // compiler
  43. ".*", // stdlib
  44. ".*", // platform
  45. ".*", // test type(s)
  46. ".*(Y[nv]|y).*Random.*", // test data group
  47. ".*", 70000, 4000); // test function
  48. add_expected_result(
  49. ".*", // compiler
  50. ".*", // stdlib
  51. ".*", // platform
  52. ".*mpfr_float_backend<18>.*", // test type(s)
  53. ".*Y0.*", // test data group
  54. ".*", 3000, 2000); // test function
  55. add_expected_result(
  56. ".*", // compiler
  57. ".*", // stdlib
  58. ".*", // platform
  59. ".*cpp_bin_float.*", // test type(s)
  60. ".*Y0.*", // test data group
  61. ".*", 40000, 20000); // test function
  62. add_expected_result(
  63. ".*", // compiler
  64. ".*", // stdlib
  65. ".*", // platform
  66. ".*", // test type(s)
  67. ".*Y0.*", // test data group
  68. ".*", 800, 400); // test function
  69. add_expected_result(
  70. ".*", // compiler
  71. ".*", // stdlib
  72. ".*", // platform
  73. ".*mpfr_float_backend<18>.*", // test type(s)
  74. ".*Yn.*", // test data group
  75. ".*", 400000, 70000); // test function
  76. add_expected_result(
  77. ".*", // compiler
  78. ".*", // stdlib
  79. ".*", // platform
  80. ".*cpp_bin_float.*", // test type(s)
  81. ".*Yn.*", // test data group
  82. ".*", 400000, 200000); // test function
  83. add_expected_result(
  84. ".*", // compiler
  85. ".*", // stdlib
  86. ".*", // platform
  87. ".*", // test type(s)
  88. ".*Yn.*", // test data group
  89. ".*", 1700, 600); // test function
  90. add_expected_result(
  91. ".*", // compiler
  92. ".*", // stdlib
  93. ".*", // platform
  94. ".*mpfr_float_backend<18>.*", // test type(s)
  95. ".*", // test data group
  96. ".*", 15000, 4000); // test function
  97. add_expected_result(
  98. ".*", // compiler
  99. ".*", // stdlib
  100. ".*", // platform
  101. ".*cpp_bin_float.*", // test type(s)
  102. ".*", // test data group
  103. ".*", 50000, 20000); // test function
  104. add_expected_result(
  105. ".*", // compiler
  106. ".*", // stdlib
  107. ".*", // platform
  108. ".*", // test type(s)
  109. ".*", // test data group
  110. ".*", 150, 60); // test function
  111. //
  112. // Finish off by printing out the compiler/stdlib/platform names,
  113. // we do this to make it easier to mark up expected error rates.
  114. //
  115. std::cout << "Tests run with " << BOOST_COMPILER << ", "
  116. << BOOST_STDLIB << ", " << BOOST_PLATFORM << std::endl;
  117. }
  118. template <class T>
  119. void test(T t, const char* p)
  120. {
  121. test_bessel(t, p);
  122. }
  123. BOOST_AUTO_TEST_CASE(test_main)
  124. {
  125. using namespace boost::multiprecision;
  126. expected_results();
  127. //
  128. // Test at:
  129. // 18 decimal digits: tests 80-bit long double approximations
  130. // 30 decimal digits: tests 128-bit long double approximations
  131. // 35 decimal digits: tests arbitrary precision code
  132. //
  133. ALL_SMALL_TESTS
  134. }