From owner-svn-src-stable-8@FreeBSD.ORG Tue Nov 23 18:48:57 2010 Return-Path: Delivered-To: svn-src-stable-8@freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:4f8:fff6::34]) by hub.freebsd.org (Postfix) with ESMTP id 9CE6D1065672; Tue, 23 Nov 2010 18:48:57 +0000 (UTC) (envelope-from uqs@FreeBSD.org) Received: from svn.freebsd.org (svn.freebsd.org [IPv6:2001:4f8:fff6::2c]) by mx1.freebsd.org (Postfix) with ESMTP id 8A36F8FC0C; Tue, 23 Nov 2010 18:48:57 +0000 (UTC) Received: from svn.freebsd.org (localhost [127.0.0.1]) by svn.freebsd.org (8.14.3/8.14.3) with ESMTP id oANImveW039482; Tue, 23 Nov 2010 18:48:57 GMT (envelope-from uqs@svn.freebsd.org) Received: (from uqs@localhost) by svn.freebsd.org (8.14.3/8.14.3/Submit) id oANImvgG039479; Tue, 23 Nov 2010 18:48:57 GMT (envelope-from uqs@svn.freebsd.org) Message-Id: <201011231848.oANImvgG039479@svn.freebsd.org> From: Ulrich Spoerlein Date: Tue, 23 Nov 2010 18:48:57 +0000 (UTC) To: src-committers@freebsd.org, svn-src-all@freebsd.org, svn-src-stable@freebsd.org, svn-src-stable-8@freebsd.org X-SVN-Group: stable-8 MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Cc: Subject: svn commit: r215757 - in stable: 7/lib/msun/src 8/lib/msun/src X-BeenThere: svn-src-stable-8@freebsd.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: SVN commit messages for only the 8-stable src tree List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 23 Nov 2010 18:48:57 -0000 Author: uqs Date: Tue Nov 23 18:48:57 2010 New Revision: 215757 URL: http://svn.freebsd.org/changeset/base/215757 Log: MFC r215237: Fix bug in jn(3) and jnf(3) that led to -inf results Modified: stable/8/lib/msun/src/e_jn.c stable/8/lib/msun/src/e_jnf.c Directory Properties: stable/8/lib/msun/ (props changed) Changes in other areas also in this revision: Modified: stable/7/lib/msun/src/e_jn.c stable/7/lib/msun/src/e_jnf.c Directory Properties: stable/7/lib/msun/ (props changed) Modified: stable/8/lib/msun/src/e_jn.c ============================================================================== --- stable/8/lib/msun/src/e_jn.c Tue Nov 23 18:47:43 2010 (r215756) +++ stable/8/lib/msun/src/e_jn.c Tue Nov 23 18:48:57 2010 (r215757) @@ -200,7 +200,12 @@ __ieee754_jn(int n, double x) } } } - 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; Modified: stable/8/lib/msun/src/e_jnf.c ============================================================================== --- stable/8/lib/msun/src/e_jnf.c Tue Nov 23 18:47:43 2010 (r215756) +++ stable/8/lib/msun/src/e_jnf.c Tue Nov 23 18:48:57 2010 (r215757) @@ -152,7 +152,12 @@ __ieee754_jnf(int n, float x) } } } - b = (t*__ieee754_j0f(x)/b); + z = __ieee754_j0f(x); + w = __ieee754_j1f(x); + if (fabsf(z) >= fabsf(w)) + b = (t*z/b); + else + b = (t*w/a); } } if(sgn==1) return -b; else return b;