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;
_______________________________________________
svn-src-all@freebsd.org mailing list
http://lists.freebsd.org/mailman/listinfo/svn-src-all
To unsubscribe, send any mail to "svn-src-all-unsubscr...@freebsd.org"
