Remove trailing whitespace.

4 files changed, 54 insertions, 54 deletions

diff --git a/lib/libm/common/atan2.c b/lib/libm/common/atan2.c
index 958a154..b847a1d 100644
--- a/lib/libm/common/atan2.c
+++ b/lib/libm/common/atan2.c
@@ -38,21 +38,21 @@ static char sccsid[] = "@(#)atan2.c	8.1 (Berkeley) 6/4/93";
 /* ATAN2(Y,X)
  * RETURN ARG (X+iY)
  * DOUBLE PRECISION (VAX D format 56 bits, IEEE DOUBLE 53 BITS)
- * CODED IN C BY K.C. NG, 1/8/85; 
+ * CODED IN C BY K.C. NG, 1/8/85;
  * REVISED BY K.C. NG on 2/7/85, 2/13/85, 3/7/85, 3/30/85, 6/29/85.
  *
  * Required system supported functions :
  *	copysign(x,y)
  *	scalb(x,y)
  *	logb(x)
- *	
+ *
  * Method :
  *	1. Reduce y to positive by atan2(y,x)=-atan2(-y,x).
- *	2. Reduce x to positive by (if x and y are unexceptional): 
+ *	2. Reduce x to positive by (if x and y are unexceptional):
  *		ARG (x+iy) = arctan(y/x)   	   ... if x > 0,
  *		ARG (x+iy) = pi - arctan[y/(-x)]   ... if x < 0,
- *	3. According to the integer k=4t+0.25 truncated , t=y/x, the argument 
- *	   is further reduced to one of the following intervals and the 
+ *	3. According to the integer k=4t+0.25 truncated , t=y/x, the argument
+ *	   is further reduced to one of the following intervals and the
  *	   arctangent of y/x is evaluated by the corresponding formula:
  *
  *         [0,7/16]	   atan(y/x) = t - t^3*(a1+t^2*(a2+...(a10+t^2*a11)...)
@@ -76,32 +76,32 @@ static char sccsid[] = "@(#)atan2.c	8.1 (Berkeley) 6/4/93";
  *	ARG( (anything but,0,NaN, and INF),+-INF ) is +-PI/2;
  *
  * Accuracy:
- *	atan2(y,x) returns (PI/pi) * the exact ARG (x+iy) nearly rounded, 
+ *	atan2(y,x) returns (PI/pi) * the exact ARG (x+iy) nearly rounded,
  *	where
  *
  *	in decimal:
- *		pi = 3.141592653589793 23846264338327 ..... 
+ *		pi = 3.141592653589793 23846264338327 .....
  *    53 bits   PI = 3.141592653589793 115997963 ..... ,
- *    56 bits   PI = 3.141592653589793 227020265 ..... ,  
+ *    56 bits   PI = 3.141592653589793 227020265 ..... ,
  *
  *	in hexadecimal:
  *		pi = 3.243F6A8885A308D313198A2E....
  *    53 bits   PI = 3.243F6A8885A30  =  2 * 1.921FB54442D18	error=.276ulps
  *    56 bits   PI = 3.243F6A8885A308 =  4 * .C90FDAA22168C2    error=.206ulps
- *	
+ *
  *	In a test run with 356,000 random argument on [-1,1] * [-1,1] on a
  *	VAX, the maximum observed error was 1.41 ulps (units of the last place)
  *	compared with (PI/pi)*(the exact ARG(x+iy)).
  *
  * Note:
  *	We use machine PI (the true pi rounded) in place of the actual
- *	value of pi for all the trig and inverse trig functions. In general, 
- *	if trig is one of sin, cos, tan, then computed trig(y) returns the 
- *	exact trig(y*pi/PI) nearly rounded; correspondingly, computed arctrig 
- *	returns the exact arctrig(y)*PI/pi nearly rounded. These guarantee the 
+ *	value of pi for all the trig and inverse trig functions. In general,
+ *	if trig is one of sin, cos, tan, then computed trig(y) returns the
+ *	exact trig(y*pi/PI) nearly rounded; correspondingly, computed arctrig
+ *	returns the exact arctrig(y)*PI/pi nearly rounded. These guarantee the
  *	trig functions have period PI, and trig(arctrig(x)) returns x for
  *	all critical values x.
- *	
+ *
  * Constants:
  * The hexadecimal values are the intended ones for the following constants.
  * The decimal values may be used, provided that the compiler will convert
@@ -174,7 +174,7 @@ ic(a11,    1.6438029044759730479E-2  ,  -6,  1.0D52174A1BB54)
 
 double atan2(y,x)
 double  y,x;
-{  
+{
 	static const double zero=0, one=1, small=1.0E-9, big=1.0E18;
 	double t,z,signy,signx,hi,lo;
 	int k,m;
@@ -185,8 +185,8 @@ double  y,x;
 #endif	/* !defined(vax)&&!defined(tahoe) */
 
     /* copy down the sign of y and x */
-	signy = copysign(one,y) ;  
-	signx = copysign(one,x) ;  
+	signy = copysign(one,y) ;
+	signx = copysign(one,x) ;
 
     /* if x is 1.0, goto begin */
 	if(x==1) { y=copysign(y,one); t=y; if(finite(t)) goto begin;}
@@ -196,10 +196,10 @@ double  y,x;
 
     /* when x = 0 */
 	if(x==zero) return(copysign(PIo2,signy));
-	    
+
     /* when x is INF */
 	if(!finite(x))
-	    if(!finite(y)) 
+	    if(!finite(y))
 		return(copysign((signx==one)?PIo4:3*PIo4,signy));
 	    else
 		return(copysign((signx==one)?zero:PI,signy));
@@ -208,43 +208,43 @@ double  y,x;
 	if(!finite(y)) return(copysign(PIo2,signy));
 
     /* compute y/x */
-	x=copysign(x,one); 
-	y=copysign(y,one); 
-	if((m=(k=logb(y))-logb(x)) > 60) t=big+big; 
+	x=copysign(x,one);
+	y=copysign(y,one);
+	if((m=(k=logb(y))-logb(x)) > 60) t=big+big;
 	    else if(m < -80 ) t=y/x;
 	    else { t = y/x ; y = scalb(y,-k); x=scalb(x,-k); }
 
     /* begin argument reduction */
 begin:
-	if (t < 2.4375) {		 
+	if (t < 2.4375) {
 
 	/* truncate 4(t+1/16) to integer for branching */
 	    k = 4 * (t+0.0625);
 	    switch (k) {
 
 	    /* t is in [0,7/16] */
-	    case 0:                    
+	    case 0:
 	    case 1:
-		if (t < small) 
+		if (t < small)
 		    { big + small ;  /* raise inexact flag */
 		      return (copysign((signx>zero)?t:PI-t,signy)); }
 
 		hi = zero;  lo = zero;  break;
 
 	    /* t is in [7/16,11/16] */
-	    case 2:                    
+	    case 2:
 		hi = athfhi; lo = athflo;
 		z = x+x;
 		t = ( (y+y) - x ) / ( z +  y ); break;
 
 	    /* t is in [11/16,19/16] */
-	    case 3:                    
+	    case 3:
 	    case 4:
 		hi = PIo4; lo = zero;
 		t = ( y - x ) / ( x + y ); break;
 
 	    /* t is in [19/16,39/16] */
-	    default:                   
+	    default:
 		hi = at1fhi; lo = at1flo;
 		z = y-x; y=y+y+y; t = x+x;
 		t = ( (z+z)-x ) / ( t + y ); break;
@@ -252,7 +252,7 @@ begin:
 	}
 	/* end of if (t < 2.4375) */
 
-	else                           
+	else
 	{
 	    hi = PIo2; lo = zero;
 
@@ -260,7 +260,7 @@ begin:
 	    if (t <= big)  t = - x / y;
 
 	    /* t is in [big, INF] */
-	    else          
+	    else
 	      { big+small;	/* raise inexact flag */
 		t = zero; }
 	}
diff --git a/lib/libm/common/sincos.c b/lib/libm/common/sincos.c
index ab88560..fc35618 100644
--- a/lib/libm/common/sincos.c
+++ b/lib/libm/common/sincos.c
@@ -67,7 +67,7 @@ double x;
 }
 
 double
-cos(x) 
+cos(x)
 double x;
 {
 	double a,c,z,s = 1.0;
@@ -83,7 +83,7 @@ double x;
 		}
 		else {			/* ... in [PI/4,3PI/4] */
 			a = PIo2-a;
-			return a+a*sin__S(a*a);	/* rtn. S(PI/2-|x|) */ 
+			return a+a*sin__S(a*a);	/* rtn. S(PI/2-|x|) */
 		}
 	}
 	if (a < small) {
diff --git a/lib/libm/common/tan.c b/lib/libm/common/tan.c
index 61ed5c5..7b49bce 100644
--- a/lib/libm/common/tan.c
+++ b/lib/libm/common/tan.c
@@ -37,7 +37,7 @@ static char sccsid[] = "@(#)tan.c	8.1 (Berkeley) 6/4/93";
 
 #include "trig.h"
 double
-tan(x) 
+tan(x)
 double x;
 {
 	double a,z,ss,cc,c;
diff --git a/lib/libm/common/trig.h b/lib/libm/common/trig.h
index 9e05b0e..e31fb4c 100644
--- a/lib/libm/common/trig.h
+++ b/lib/libm/common/trig.h
@@ -67,7 +67,7 @@ static const double
 	zero = 0,
 	one = 1,
 	negone = -1,
-	half = 1.0/2.0, 
+	half = 1.0/2.0,
 	small = 1E-10,	/* 1+small**2 == 1; better values for small:
 			 *		small	= 1.5E-9 for VAX D
 			 *			= 1.2E-8 for IEEE Double
@@ -77,27 +77,27 @@ static const double
 
 /* sin__S(x*x) ... re-implemented as a macro
  * DOUBLE PRECISION (VAX D format 56 bits, IEEE DOUBLE 53 BITS)
- * STATIC KERNEL FUNCTION OF SIN(X), COS(X), AND TAN(X) 
- * CODED IN C BY K.C. NG, 1/21/85; 
+ * STATIC KERNEL FUNCTION OF SIN(X), COS(X), AND TAN(X)
+ * CODED IN C BY K.C. NG, 1/21/85;
  * REVISED BY K.C. NG on 8/13/85.
  *
  *	    sin(x*k) - x
  * RETURN  --------------- on [-PI/4,PI/4] , where k=pi/PI, PI is the rounded
- *	            x	
+ *	            x
  * value of pi in machine precision:
  *
  *	Decimal:
- *		pi = 3.141592653589793 23846264338327 ..... 
+ *		pi = 3.141592653589793 23846264338327 .....
  *    53 bits   PI = 3.141592653589793 115997963 ..... ,
- *    56 bits   PI = 3.141592653589793 227020265 ..... ,  
+ *    56 bits   PI = 3.141592653589793 227020265 ..... ,
  *
  *	Hexadecimal:
  *		pi = 3.243F6A8885A308D313198A2E....
  *    53 bits   PI = 3.243F6A8885A30  =  2 * 1.921FB54442D18
- *    56 bits   PI = 3.243F6A8885A308 =  4 * .C90FDAA22168C2    
+ *    56 bits   PI = 3.243F6A8885A308 =  4 * .C90FDAA22168C2
  *
  * Method:
- *	1. Let z=x*x. Create a polynomial approximation to 
+ *	1. Let z=x*x. Create a polynomial approximation to
  *	    (sin(k*x)-x)/x  =  z*(S0 + S1*z^1 + ... + S5*z^5).
  *	Then
  *      sin__S(x*x) = z*(S0 + S1*z^1 + ... + S5*z^5)
@@ -105,8 +105,8 @@ static const double
  *	The coefficient S's are obtained by a special Remez algorithm.
  *
  * Accuracy:
- *	In the absence of rounding error, the approximation has absolute error 
- *	less than 2**(-61.11) for VAX D FORMAT, 2**(-57.45) for IEEE DOUBLE. 
+ *	In the absence of rounding error, the approximation has absolute error
+ *	less than 2**(-61.11) for VAX D FORMAT, 2**(-57.45) for IEEE DOUBLE.
  *
  * Constants:
  * The hexadecimal values are the intended ones for the following constants.
@@ -149,28 +149,28 @@ ic(S5,  1.5868926979889205164E-10 , -33,  1.5CF61DF672B13)
 
 /* cos__C(x*x) ... re-implemented as a macro
  * DOUBLE PRECISION (VAX D FORMAT 56 BITS, IEEE DOUBLE 53 BITS)
- * STATIC KERNEL FUNCTION OF SIN(X), COS(X), AND TAN(X) 
- * CODED IN C BY K.C. NG, 1/21/85; 
+ * STATIC KERNEL FUNCTION OF SIN(X), COS(X), AND TAN(X)
+ * CODED IN C BY K.C. NG, 1/21/85;
  * REVISED BY K.C. NG on 8/13/85.
  *
- *	   		    x*x	
+ *	   		    x*x
  * RETURN   cos(k*x) - 1 + ----- on [-PI/4,PI/4],  where k = pi/PI,
- *	  		     2	
+ *	  		     2
  * PI is the rounded value of pi in machine precision :
  *
  *	Decimal:
- *		pi = 3.141592653589793 23846264338327 ..... 
+ *		pi = 3.141592653589793 23846264338327 .....
  *    53 bits   PI = 3.141592653589793 115997963 ..... ,
- *    56 bits   PI = 3.141592653589793 227020265 ..... ,  
+ *    56 bits   PI = 3.141592653589793 227020265 ..... ,
  *
  *	Hexadecimal:
  *		pi = 3.243F6A8885A308D313198A2E....
  *    53 bits   PI = 3.243F6A8885A30  =  2 * 1.921FB54442D18
- *    56 bits   PI = 3.243F6A8885A308 =  4 * .C90FDAA22168C2    
+ *    56 bits   PI = 3.243F6A8885A308 =  4 * .C90FDAA22168C2
  *
  *
  * Method:
- *	1. Let z=x*x. Create a polynomial approximation to 
+ *	1. Let z=x*x. Create a polynomial approximation to
  *	    cos(k*x)-1+z/2  =  z*z*(C0 + C1*z^1 + ... + C5*z^5)
  *	then
  *      cos__C(z) =  z*z*(C0 + C1*z^1 + ... + C5*z^5)
@@ -178,9 +178,9 @@ ic(S5,  1.5868926979889205164E-10 , -33,  1.5CF61DF672B13)
  *	The coefficient C's are obtained by a special Remez algorithm.
  *
  * Accuracy:
- *	In the absence of rounding error, the approximation has absolute error 
- *	less than 2**(-64) for VAX D FORMAT, 2**(-58.3) for IEEE DOUBLE. 
- *	
+ *	In the absence of rounding error, the approximation has absolute error
+ *	less than 2**(-64) for VAX D FORMAT, 2**(-58.3) for IEEE DOUBLE.
+ *
  *
  * Constants:
  * The hexadecimal values are the intended ones for the following constants.