This is the mail archive of the
gcc-bugs@gcc.gnu.org
mailing list for the GCC project.
[Bug libstdc++/83566] cyl_bessel_j returns wrong result for x>1000 for high orders.
- From: "emsr at gcc dot gnu.org" <gcc-bugzilla at gcc dot gnu dot org>
- To: gcc-bugs at gcc dot gnu dot org
- Date: Sun, 18 Nov 2018 18:32:58 +0000
- Subject: [Bug libstdc++/83566] cyl_bessel_j returns wrong result for x>1000 for high orders.
- Auto-submitted: auto-generated
- References: <bug-83566-4@http.gcc.gnu.org/bugzilla/>
https://gcc.gnu.org/bugzilla/show_bug.cgi?id=83566
--- Comment #3 from emsr at gcc dot gnu.org ---
Author: emsr
Date: Sun Nov 18 18:32:26 2018
New Revision: 266252
URL: https://gcc.gnu.org/viewcvs?rev=266252&root=gcc&view=rev
Log:
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.
Modified:
trunk/libstdc++-v3/ChangeLog
trunk/libstdc++-v3/include/tr1/bessel_function.tcc
trunk/libstdc++-v3/testsuite/special_functions/08_cyl_bessel_j/check_value.cc
trunk/libstdc++-v3/testsuite/special_functions/10_cyl_neumann/check_value.cc
trunk/libstdc++-v3/testsuite/tr1/5_numerical_facilities/special_functions/09_cyl_bessel_j/check_value.cc
trunk/libstdc++-v3/testsuite/tr1/5_numerical_facilities/special_functions/11_cyl_neumann/check_value.cc