From owner-svn-src-all@FreeBSD.ORG  Thu Dec  4 15:57:59 2014
Return-Path: <owner-svn-src-all@FreeBSD.ORG>
Delivered-To: svn-src-all@freebsd.org
Received: from mx1.freebsd.org (mx1.freebsd.org
 [IPv6:2001:1900:2254:206a::19:1])
 (using TLSv1.2 with cipher AECDH-AES256-SHA (256/256 bits))
 (No client certificate requested)
 by hub.freebsd.org (Postfix) with ESMTPS id 58280F4C;
 Thu,  4 Dec 2014 15:57:59 +0000 (UTC)
Received: from svn.freebsd.org (svn.freebsd.org
 [IPv6:2001:1900:2254:2068::e6a:0])
 (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits))
 (Client did not present a certificate)
 by mx1.freebsd.org (Postfix) with ESMTPS id 4215FA3E;
 Thu,  4 Dec 2014 15:57:59 +0000 (UTC)
Received: from svn.freebsd.org ([127.0.1.70])
 by svn.freebsd.org (8.14.9/8.14.9) with ESMTP id sB4FvxCL093591;
 Thu, 4 Dec 2014 15:57:59 GMT (envelope-from kargl@FreeBSD.org)
Received: (from kargl@localhost)
 by svn.freebsd.org (8.14.9/8.14.9/Submit) id sB4FvxH8093590;
 Thu, 4 Dec 2014 15:57:59 GMT (envelope-from kargl@FreeBSD.org)
Message-Id: <201412041557.sB4FvxH8093590@svn.freebsd.org>
X-Authentication-Warning: svn.freebsd.org: kargl set sender to
 kargl@FreeBSD.org using -f
From: Steve Kargl <kargl@FreeBSD.org>
Date: Thu, 4 Dec 2014 15:57:59 +0000 (UTC)
To: src-committers@freebsd.org, svn-src-all@freebsd.org,
 svn-src-head@freebsd.org
Subject: svn commit: r275476 - head/lib/msun/src
X-SVN-Group: head
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
X-BeenThere: svn-src-all@freebsd.org
X-Mailman-Version: 2.1.18-1
Precedence: list
List-Id: "SVN commit messages for the entire src tree \(except for &quot;
 user&quot; and &quot; projects&quot; \)" <svn-src-all.freebsd.org>
List-Unsubscribe: <http://lists.freebsd.org/mailman/options/svn-src-all>,
 <mailto:svn-src-all-request@freebsd.org?subject=unsubscribe>
List-Archive: <http://lists.freebsd.org/pipermail/svn-src-all/>
List-Post: <mailto:svn-src-all@freebsd.org>
List-Help: <mailto:svn-src-all-request@freebsd.org?subject=help>
List-Subscribe: <http://lists.freebsd.org/mailman/listinfo/svn-src-all>,
 <mailto:svn-src-all-request@freebsd.org?subject=subscribe>
X-List-Received-Date: Thu, 04 Dec 2014 15:57:59 -0000

Author: kargl
Date: Thu Dec  4 15:57:58 2014
New Revision: 275476
URL: https://svnweb.freebsd.org/changeset/base/275476

Log:
  Fix a 20+ year bug by using an appropriate constant for
  the transition from one asymptotic approximation to another
  for the zeroth order Bessel and Neumann functions.
  
  Reviewed by:	bde

Modified:
  head/lib/msun/src/e_j0f.c

Modified: head/lib/msun/src/e_j0f.c
==============================================================================
--- head/lib/msun/src/e_j0f.c	Thu Dec  4 15:54:45 2014	(r275475)
+++ head/lib/msun/src/e_j0f.c	Thu Dec  4 15:57:58 2014	(r275476)
@@ -62,7 +62,7 @@ __ieee754_j0f(float x)
 	 * j0(x) = 1/sqrt(pi) * (P(0,x)*cc - Q(0,x)*ss) / sqrt(x)
 	 * y0(x) = 1/sqrt(pi) * (P(0,x)*ss + Q(0,x)*cc) / sqrt(x)
 	 */
-		if(ix>0x80000000) z = (invsqrtpi*cc)/sqrtf(x);
+		if(ix>0x54000000) z = (invsqrtpi*cc)/sqrtf(x);
 		else {
 		    u = pzerof(x); v = qzerof(x);
 		    z = invsqrtpi*(u*cc-v*ss)/sqrtf(x);
@@ -136,7 +136,7 @@ __ieee754_y0f(float x)
                     if ((s*c)<zero) cc = z/ss;
                     else            ss = z/cc;
                 }
-                if(ix>0x80000000) z = (invsqrtpi*ss)/sqrtf(x);
+                if(ix>0x54800000) z = (invsqrtpi*ss)/sqrtf(x);
                 else {
                     u = pzerof(x); v = qzerof(x);
                     z = invsqrtpi*(u*ss+v*cc)/sqrtf(x);