freebsd-skq/lib/msun
uqs 34589aa4cb Fix bug in jn(3) and jnf(3) that led to -inf results
Explanation by Steve:
jn[f](n,x) for certain ranges of x uses downward recursion to compute
the value of the function.  The recursion sequence that is generated is
proportional to the actual desired value, so a normalization step is
taken.  This normalization is j0[f](x) divided by the zeroth sequence
member.  As Bruce notes, near the zeros of j0[f](x) the computed value
can have giga-ULP inaccuracy. I found for the 1st zero of j0f(x) only
the leading decimal digit is correct.  The solution to the issue is
fairly straight forward.  The zeros of j0(x) and j1(x) never coincide,
so as j0(x) approaches a zero, the normalization constant switches to
j1[f](x) divided by the 2nd sequence member.  The expectation is that
j1[f](x) is a more accurately computed value.

PR:		bin/144306
Submitted by:	Steven G. Kargl <kargl@troutmask.apl.washington.edu>
Reviewed by:	bde
MFC after:	7 days
2010-11-13 10:54:10 +00:00
..
amd64 Use __FBSDID() instead of RCSID() in most .S files under lib/msun/i386, 2010-10-01 20:14:36 +00:00
arm
bsdsrc Eliminate some warnings. 2008-02-22 02:26:51 +00:00
i387 Use __FBSDID() instead of RCSID() in most .S files under lib/msun/i386, 2010-10-01 20:14:36 +00:00
ia64 Use, in uncovered part, the END() macro in order to improve debugging. 2009-05-25 14:37:10 +00:00
ld80 On i386, gcc truncates long double constants to double precision 2008-08-02 03:56:22 +00:00
ld128 Add implementations of acosl(), asinl(), atanl(), atan2l(), 2008-07-31 22:41:26 +00:00
man Fix typos, spelling, formatting and mdoc mistakes found by Nobuyuki while 2010-08-16 15:18:30 +00:00
mips Symbol.map is handled by cpp, so use C-style comments 2008-05-03 21:16:08 +00:00
powerpc Let PowerPC world optionally build with -msoft-float. For FPU-less PowerPC 2008-02-24 19:22:53 +00:00
sparc64
src Fix bug in jn(3) and jnf(3) that led to -inf results 2010-11-13 10:54:10 +00:00
Makefile This is exactly the same as the .else, so remove it. 2010-09-13 04:23:23 +00:00
Symbol.map Introduce __isnanf() as an alias for isnanf(), and make the isnan() 2010-06-12 17:32:05 +00:00