From owner-freebsd-bugs@FreeBSD.ORG Fri Feb 26 04:20:03 2010 Return-Path: Delivered-To: freebsd-bugs@hub.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:4f8:fff6::34]) by hub.freebsd.org (Postfix) with ESMTP id 60BE8106564A for ; Fri, 26 Feb 2010 04:20:03 +0000 (UTC) (envelope-from gnats@FreeBSD.org) Received: from freefall.freebsd.org (freefall.freebsd.org [IPv6:2001:4f8:fff6::28]) by mx1.freebsd.org (Postfix) with ESMTP id 4FD738FC20 for ; Fri, 26 Feb 2010 04:20:03 +0000 (UTC) Received: from freefall.freebsd.org (localhost [127.0.0.1]) by freefall.freebsd.org (8.14.3/8.14.3) with ESMTP id o1Q4K37A044887 for ; Fri, 26 Feb 2010 04:20:03 GMT (envelope-from gnats@freefall.freebsd.org) Received: (from gnats@localhost) by freefall.freebsd.org (8.14.3/8.14.3/Submit) id o1Q4K3SG044886; Fri, 26 Feb 2010 04:20:03 GMT (envelope-from gnats) Date: Fri, 26 Feb 2010 04:20:03 GMT Message-Id: <201002260420.o1Q4K3SG044886@freefall.freebsd.org> To: freebsd-bugs@FreeBSD.org From: "Steven G. Kargl" Cc: Subject: Re: bin/144306: Nasty bug in jn(3) X-BeenThere: freebsd-bugs@freebsd.org X-Mailman-Version: 2.1.5 Precedence: list Reply-To: "Steven G. Kargl" List-Id: Bug reports List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Fri, 26 Feb 2010 04:20:03 -0000 The following reply was made to PR bin/144306; it has been noted by GNATS. From: "Steven G. Kargl" To: FreeBSD-gnats-submit@FreeBSD.org, freebsd-bugs@FreeBSD.org Cc: Subject: Re: bin/144306: Nasty bug in jn(3) Date: Thu, 25 Feb 2010 20:15:33 -0800 (PST) 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/