PR libstdc++/83566 - cyl_bessel_j returns wrong result for x>1000
authorMichele Pezzutti <mpezz@tiscali.it>
Sun, 18 Nov 2018 18:32:26 +0000 (19:32 +0100)
committerEdward Smith-Rowland <emsr@gcc.gnu.org>
Sun, 18 Nov 2018 18:32:26 +0000 (18:32 +0000)
commitbee39274cb9ffb5415ab50fa102673b0394d4f39
tree3a8f7894243809aefb24e4e78105559254d4f539
parentcb40e8071e0fa1a86183f6b97ebeb3355e5e1722
PR libstdc++/83566 - cyl_bessel_j returns wrong result for x>1000

2018-11-16  Michele Pezzutti <mpezz@tiscali.it>
    Edward Smith-Rowland  <3dw4rd@verizon.net>

PR libstdc++/83566 - cyl_bessel_j returns wrong result for x>1000
for high orders.
* include/tr1/bessel_function.tcc: Perform no fewer than nu/2 iterations
of the asymptotic series (nu is the Bessel order).
* testsuite/tr1/5_numerical_facilities/special_functions/
09_cyl_bessel_j/check_value.cc: Add tests at nu=100, 1000<=x<=2000.
* testsuite/tr1/5_numerical_facilities/special_functions/
11_cyl_neumann/check_value.cc: Ditto.
* testsuite/special_functions/08_cyl_bessel_j/check_value.cc: Ditto.
* testsuite/special_functions/10_cyl_neumann/check_value.cc: Ditto.

Co-Authored-By: Edward Smith-Rowland <3dw4rd@verizon.net>
From-SVN: r266252
libstdc++-v3/ChangeLog
libstdc++-v3/include/tr1/bessel_function.tcc
libstdc++-v3/testsuite/special_functions/08_cyl_bessel_j/check_value.cc
libstdc++-v3/testsuite/special_functions/10_cyl_neumann/check_value.cc
libstdc++-v3/testsuite/tr1/5_numerical_facilities/special_functions/09_cyl_bessel_j/check_value.cc
libstdc++-v3/testsuite/tr1/5_numerical_facilities/special_functions/11_cyl_neumann/check_value.cc