svn commit: r270932 - head/lib/msun/src

[Steve Kargl]

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;