Date: Thu, 25 Feb 2010 20:15:33 -0800 (PST) From: "Steven G. Kargl" <kargl@troutmask.apl.washington.edu> To: FreeBSD-gnats-submit@FreeBSD.org, freebsd-bugs@FreeBSD.org Subject: Re: bin/144306: Nasty bug in jn(3) Message-ID: <201002260415.o1Q4FXC6029992@troutmask.apl.washington.edu> In-Reply-To: <201002260020.o1Q0K40O036894@freefall.freebsd.org>
next in thread | previous in thread | raw e-mail | index | archive | help
This patch fixes the problem. Index: e_jn.c =================================================================== --- e_jn.c (revision 204219) +++ e_jn.c (working copy) @@ -200,7 +200,12 @@ } } } - b = (t*__ieee754_j0(x)/b); + z = __ieee754_j0(x); + w = __ieee754_j1(x); + if (fabs(z) >= fabs(w)) + b = (t*z/b); + else + b = (t*w/a); } } if(sgn==1) return -b; else return b; -- Steve http://troutmask.apl.washington.edu/~kargl/
Want to link to this message? Use this URL: <https://mail-archive.FreeBSD.org/cgi/mid.cgi?201002260415.o1Q4FXC6029992>