From owner-svn-src-head@FreeBSD.ORG Mon Sep 1 16:24:26 2014
Return-Path: <owner-svn-src-head@FreeBSD.ORG>
Delivered-To: svn-src-head@freebsd.org
Received: from mx1.freebsd.org (mx1.freebsd.org [8.8.178.115])
(using TLSv1 with cipher ADH-AES256-SHA (256/256 bits))
(No client certificate requested)
by hub.freebsd.org (Postfix) with ESMTPS id 25B41947;
Mon, 1 Sep 2014 16:24:26 +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 EB8F31FBE;
Mon, 1 Sep 2014 16:24:25 +0000 (UTC)
Received: from svn.freebsd.org ([127.0.1.70])
by svn.freebsd.org (8.14.9/8.14.9) with ESMTP id s81GOPRt019856;
Mon, 1 Sep 2014 16:24:25 GMT (envelope-from kargl@FreeBSD.org)
Received: (from kargl@localhost)
by svn.freebsd.org (8.14.9/8.14.9/Submit) id s81GOPdM019854;
Mon, 1 Sep 2014 16:24:25 GMT (envelope-from kargl@FreeBSD.org)
Message-Id: <201409011624.s81GOPdM019854@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: Mon, 1 Sep 2014 16:24:25 +0000 (UTC)
To: src-committers@freebsd.org, svn-src-all@freebsd.org,
svn-src-head@freebsd.org
Subject: svn commit: r270932 - 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-head@freebsd.org
X-Mailman-Version: 2.1.18-1
Precedence: list
List-Id: SVN commit messages for the src tree for head/-current
<svn-src-head.freebsd.org>
List-Unsubscribe: <http://lists.freebsd.org/mailman/options/svn-src-head>,
<mailto:svn-src-head-request@freebsd.org?subject=unsubscribe>
List-Archive: <http://lists.freebsd.org/pipermail/svn-src-head/>
List-Post: <mailto:svn-src-head@freebsd.org>
List-Help: <mailto:svn-src-head-request@freebsd.org?subject=help>
List-Subscribe: <http://lists.freebsd.org/mailman/listinfo/svn-src-head>,
<mailto:svn-src-head-request@freebsd.org?subject=subscribe>
X-List-Received-Date: Mon, 01 Sep 2014 16:24:26 -0000
Author: kargl
Date: Mon Sep 1 16:24:25 2014
New Revision: 270932
URL: http://svnweb.freebsd.org/changeset/base/270932
Log:
Fix a tab that somehow became 8 spaces.
Remove parentheses in a return statement to be consistent with the
rest of the file.
Rename sin_pi() in the float version to sin_pif().
Remove large comment that precedes sin_pif(). The comment
duplicates a comment in e_lgamma_r.c where the algorithm
is documented.
Requested by: bde
Modified:
head/lib/msun/src/e_lgamma_r.c
head/lib/msun/src/e_lgammaf_r.c
Modified: head/lib/msun/src/e_lgamma_r.c
==============================================================================
--- head/lib/msun/src/e_lgamma_r.c Mon Sep 1 16:12:29 2014 (r270931)
+++ head/lib/msun/src/e_lgamma_r.c Mon Sep 1 16:24:25 2014 (r270932)
@@ -171,9 +171,9 @@ sin_pi(double x)
y = -x;
vz = y+0x1p52; /* depend on 0 <= y < 0x1p52 */
- z = vz-0x1p52; /* rint(y) for the above range */
+ z = vz-0x1p52; /* rint(y) for the above range */
if (z == y)
- return (zero);
+ return zero;
vz = y+0x1p50;
GET_LOW_WORD(n,vz); /* bits for rounded y (units 0.25) */
Modified: head/lib/msun/src/e_lgammaf_r.c
==============================================================================
--- head/lib/msun/src/e_lgammaf_r.c Mon Sep 1 16:12:29 2014 (r270931)
+++ head/lib/msun/src/e_lgammaf_r.c Mon Sep 1 16:24:25 2014 (r270932)
@@ -89,13 +89,8 @@ w6 = -1.6309292987e-03; /* 0xbad5c4e8 *
static const float zero= 0.0000000000e+00;
-/*
- * Compute sin(pi*x) without actually doing the pi*x multiplication.
- * sin_pi(x) is only called for x < 0 and |x| < 2**(p-1) where p is
- * the precision of x.
- */
static float
-sin_pi(float x)
+sin_pif(float x)
{
volatile float vz;
float y,z;
@@ -106,7 +101,7 @@ sin_pi(float x)
vz = y+0x1p23F; /* depend on 0 <= y < 0x1p23 */
z = vz-0x1p23F; /* rintf(y) for the above range */
if (z == y)
- return (zero);
+ return zero;
vz = y+0x1p21F;
GET_FLOAT_WORD(n,vz); /* bits for rounded y (units 0.25) */
@@ -155,7 +150,7 @@ __ieee754_lgammaf_r(float x, int *signga
if(hx<0) {
if(ix>=0x4b000000) /* |x|>=2**23, must be -integer */
return one/zero;
- t = sin_pi(x);
+ t = sin_pif(x);
if(t==zero) return one/zero; /* -integer */
nadj = __ieee754_logf(pi/fabsf(t*x));
if(t<zero) *signgamp = -1;