[Bug libstdc++/83566] cyl_bessel_j returns wrong result for x>1000 for high orders.
mpezz at tiscali dot it
gcc-bugzilla@gcc.gnu.org
Mon Jan 1 02:40:00 GMT 2018
https://gcc.gnu.org/bugzilla/show_bug.cgi?id=83566
Michele Pezzutti <mpezz at tiscali dot it> changed:
What |Removed |Added
----------------------------------------------------------------------------
CC| |mpezz at tiscali dot it
--- Comment #2 from Michele Pezzutti <mpezz at tiscali dot it> ---
Series expansion in __cyl_bessel_jn_asymp() shall not be truncated at the first
terms. At least nu/2 terms shall be added, in order to have a guaranteed error
bound.
Patch proposed to mailing list.
More information about the Gcc-bugs
mailing list