summaryrefslogtreecommitdiff
path: root/libm/powerpc/w_scalb.c
diff options
context:
space:
mode:
authorEric Andersen <andersen@codepoet.org>2005-03-06 07:11:53 +0000
committerEric Andersen <andersen@codepoet.org>2005-03-06 07:11:53 +0000
commitc4e44e97f8562254d9da898f6ed7e6cb4d8a3ce4 (patch)
tree6c61f83ac5b94085222b3eda8d731309d61be99b /libm/powerpc/w_scalb.c
parentd4fad9c64ee518be51ecb40662af69b405a49556 (diff)
Trim off whitespace
Diffstat (limited to 'libm/powerpc/w_scalb.c')
-rw-r--r--libm/powerpc/w_scalb.c32
1 files changed, 16 insertions, 16 deletions
diff --git a/libm/powerpc/w_scalb.c b/libm/powerpc/w_scalb.c
index 98a2b7d3f..c93c74b68 100644
--- a/libm/powerpc/w_scalb.c
+++ b/libm/powerpc/w_scalb.c
@@ -1,26 +1,26 @@
/***********************************************************************
** File: scalb.c
-**
+**
** Contains: C source code for implementations of floating-point
** scalb functions defined in header <fp.h>. In
** particular, this file contains implementations of
** functions scalb and scalbl for double and long double
** formats on PowerPC platforms.
-**
+**
** Written by: Jon Okada, SANEitation Engineer, ext. 4-4838
-**
+**
** Copyright: © 1992 by Apple Computer, Inc., all rights reserved
-**
+**
** Change History ( most recent first ):
**
** 28 May 97 ali made an speed improvement for large n,
** removed scalbl.
** 12 Dec 92 JPO First created.
-**
+**
***********************************************************************/
-typedef union
- {
+typedef union
+ {
struct {
#if defined(__BIG_ENDIAN__)
unsigned long int hi;
@@ -41,35 +41,35 @@ static const double twoToM1022 = 2.225073858507201383e-308; // 0x1p-1022
double scalb( double x, long int n ) returns its argument x scaled
by the factor 2^m. NaNs, signed zeros, and infinities are propagated
by this function regardless of the value of n.
-
+
Exceptions: OVERFLOW/INEXACT or UNDERFLOW inexact may occur;
INVALID for signaling NaN inputs ( quiet NaN returned ).
-
+
Calls: none.
***********************************************************************/
double scalb ( double x, int n )
{
DblInHex xInHex;
-
+
xInHex.words.lo = 0UL; // init. low half of xInHex
-
- if ( n > 1023 )
+
+ if ( n > 1023 )
{ // large positive scaling
if ( n > 2097 ) // huge scaling
return ( ( x * twoTo1023 ) * twoTo1023 ) * twoTo1023;
- while ( n > 1023 )
+ while ( n > 1023 )
{ // scale reduction loop
x *= twoTo1023; // scale x by 2^1023
n -= 1023; // reduce n by 1023
}
}
-
- else if ( n < -1022 )
+
+ else if ( n < -1022 )
{ // large negative scaling
if ( n < -2098 ) // huge negative scaling
return ( ( x * twoToM1022 ) * twoToM1022 ) * twoToM1022;
- while ( n < -1022 )
+ while ( n < -1022 )
{ // scale reduction loop
x *= twoToM1022; // scale x by 2^( -1022 )
n += 1022; // incr n by 1022