freebsd-skq/lib
Ulrich Spörlein e61ffaea2a 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
..
bind Prep for the 9.6-ESV-R2 update 2010-10-31 04:45:25 +00:00
clang When the make target is 'install', don't descend into the clang 2010-10-30 16:53:42 +00:00
csu Remove two .endp's without matching .proc in lib/csu/ia64/crtn.S. 2010-10-15 21:40:20 +00:00
libalias Remove redundant WARNS?=6 overrides and inherit the WARNS setting from 2010-03-02 18:44:08 +00:00
libarchive If the Zip reader doesn't see a PK signature block 2010-11-07 03:40:37 +00:00
libauditd Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libautofs Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libbegemot Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libblocksruntime Import libcompiler_rt into HEAD and add Makefiles. 2010-11-11 15:13:11 +00:00
libbluetooth Add the following Linux BlueZ compatibility macros: htobs(), htobl(), 2010-09-22 23:41:02 +00:00
libbsm Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libbsnmp Always assign WARNS using ?= 2010-03-02 16:58:04 +00:00
libbz2 Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libc Sync with OpenBSD, primarily better signal and terminal handling. 2010-11-13 10:38:06 +00:00
libc_r MFtbemd: 2010-08-23 22:24:11 +00:00
libcalendar Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libcam camlib.c: update one overlooked comment 2010-10-11 21:34:35 +00:00
libcom_err
libcompat Small style(9) fix: use tabs instead of spaces. 2010-08-08 08:19:23 +00:00
libcompiler_rt Revert to libgcc for sparc64. 2010-11-12 15:52:27 +00:00
libcrypt Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libdevinfo Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libdevstat Spelling fixes. 2010-08-03 17:40:09 +00:00
libdisk mdoc: drop redundant .Pp and .LP calls 2010-10-08 12:40:16 +00:00
libdwarf Removed redundant -I. from CFLAGS and "yes" from WITHOUT_MAN. 2010-02-25 22:16:30 +00:00
libedit libedit: Try to map <Delete> to ed-delete-next-char. 2010-09-05 16:12:10 +00:00
libefi Spelling fixes. 2010-08-03 17:40:09 +00:00
libelf mdoc: drop redundant .Pp and .LP calls 2010-10-08 12:40:16 +00:00
libexpat Fix buildworld -DNO_CLEAN when using with Perforce, which marks files as 2010-08-12 20:46:49 +00:00
libfetch Move variable declarations into the conditional block where they are 2010-10-24 01:05:10 +00:00
libftpio mdoc: consistently spell our email addresses <foo@FreeBSD.org> 2010-05-19 08:57:53 +00:00
libgeom Remove code duplication by introducing static gctl_param_add() function which 2010-10-21 10:38:14 +00:00
libgpib Fix typos, spelling, formatting and mdoc mistakes found by Nobuyuki while 2010-08-16 15:18:30 +00:00
libgssapi mdoc: drop redundant .Pp and .LP calls 2010-10-08 12:40:16 +00:00
libipsec Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libipx Fix incorrect usage of 'assure' and 'insure'. 2010-08-28 16:32:01 +00:00
libjail Find a jail's type as part of jailparam_init rather than waiting until 2010-10-27 21:01:53 +00:00
libkiconv Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libkse Merge from tbemd: 2010-06-13 01:27:29 +00:00
libkvm amd64: introduce minidump version 2 2010-11-11 18:35:28 +00:00
liblzma Update xz to release 5.0.0 2010-11-12 16:53:17 +00:00
libmagic Use CFLAGS and LDFLAGS when compiling mkmagic. This ensures that all the 2010-08-01 12:10:32 +00:00
libmd Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libmemstat Spelling fixes. 2010-08-03 17:40:09 +00:00
libmilter Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libmp Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libncp Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libnetgraph bring the NGM_IPFW_COOKIE back into ng_ipfw.h, libnetgraph expects 2009-12-28 12:29:13 +00:00
libngatm Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libopie Let libopie use utmpx instead of utmp. 2010-01-13 18:02:30 +00:00
libpam Bump .Dd date. 2010-05-03 09:49:42 +00:00
libpcap Regenerate our lib/libpcap/config.h for libpcap 1.1.1. 2010-10-29 22:57:14 +00:00
libpkg - Take libinstall.a out of pkg_install and make it a proper shared library. 2010-04-23 11:07:43 +00:00
libpmc Fix manpage markup. 2010-11-06 10:54:33 +00:00
libproc Ignore EINTR when calling waitpid. 2010-09-18 23:38:21 +00:00
libradius Add links for libradius(3) functions. 2010-10-18 12:35:10 +00:00
librpcsec_gss Fix typos, spelling, formatting and mdoc mistakes found by Nobuyuki while 2010-08-16 15:18:30 +00:00
librpcsvc Use ANSI function declarations in librpcsvc. 2009-02-26 20:32:11 +00:00
librt To support stack unwinding for cancellation points, add -fexceptions flag 2010-09-25 01:57:47 +00:00
librtld_db Spelling fixes. 2010-08-03 17:40:09 +00:00
libsbuf Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libsdp Fix typos, spelling, formatting and mdoc mistakes found by Nobuyuki while 2010-08-16 15:18:30 +00:00
libsm Enable the use of nanosleep() instead of using pause() and signals. 2010-03-04 05:53:06 +00:00
libsmb Bump the version of all non-symbol-versioned shared libraries in 2009-07-19 17:25:24 +00:00
libsmdb Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libsmutil Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libstand mdoc: drop redundant .Pp and .LP calls 2010-10-08 12:40:16 +00:00
libtacplus mdoc: drop redundant .Pp and .LP calls 2010-10-08 12:40:16 +00:00
libtelnet Rename all symbols in libmp(3) to mp_*, just like Solaris. 2009-02-26 21:43:15 +00:00
libthr Use sysctl kern.sched.cpusetsize to retrieve size of kernel cpuset. 2010-11-02 02:13:13 +00:00
libthread_db Repair some build breakage introduced in r211725 and garbage collect some 2010-08-28 15:03:11 +00:00
libufs Search beyond the first 1/8th of inodes. 2010-05-01 18:56:45 +00:00
libugidfw mdoc: order prologue macros consistently by Dd/Dt/Os 2010-04-14 19:08:06 +00:00
libulog Fix reference to nonexistent manpage getuid(3). 2010-10-11 20:30:57 +00:00
libusb - Add missing LibUSB API functions: 2010-10-14 20:50:33 +00:00
libusbhid - Add support for libusbhid in 32-bit compatibility mode. 2010-10-16 11:20:53 +00:00
libutil mdoc: drop redundant .Pp and .LP calls 2010-10-08 12:40:16 +00:00
libvgl Change libvgl's set4pixels() and set2lines() functions from plain 2010-10-04 18:16:38 +00:00
libwrap Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
liby Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libypclnt Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
libz Revert 212517 to restore pristine state of this file 2010-10-13 16:34:08 +00:00
msun Fix bug in jn(3) and jnf(3) that led to -inf results 2010-11-13 10:54:10 +00:00
ncurses Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00
Makefile Replace libgcc.a by libcompiler_rt.a. 2010-11-11 15:48:27 +00:00
Makefile.inc Build lib/ with WARNS=6 by default. 2010-01-02 09:58:07 +00:00