1 files changed, 28 insertions, 46 deletions

diff --git a/newlib/libm/mathfp/s_pow.c b/newlib/libm/mathfp/s_pow.c
index 90d9d0b5f..7c0a38a20 100644
--- a/newlib/libm/mathfp/s_pow.c
+++ b/newlib/libm/mathfp/s_pow.c
@@ -24,7 +24,7 @@ TRAD_SYNOPSIS
         float <[x]>, <[y]>;
 
 DESCRIPTION
-        <<pow>> and <<powf>> calculate <[x]> raised to the exponent <[y]>.
+        <<pow>> and <<powf>> calculate <[x]> raised to the exp1.0nt <[y]>.
         @tex
         (That is, $x^y$.)
         @end tex
@@ -52,68 +52,57 @@ PORTABILITY
 
 double pow (double x, double y)
 {
-  double d, k, t, r = 1.0;
-  int n, sign, exponent_is_even_int = 0;
+  double d, t, r = 1.0;
+  int n, k, sign = 0;
   __uint32_t px;
 
   GET_HIGH_WORD (px, x);
 
   k = modf (y, &d);
-
-  if (k == 0.0)
+  if (k == 0.0) 
     {
-      /* Exponent y is an integer. */
       if (modf (ldexp (y, -1), &t))
-        {
-          /* y is odd. */
-          exponent_is_even_int = 0;
-        }
+        sign = 0;
       else
-        {
-          /* y is even. */
-          exponent_is_even_int = 1;
-        }
+        sign = 1; 
     }
 
-  if (x == 0.0)
-    {
-      if (y <= 0.0)
-        errno = EDOM;
-    }
+  if (x == 0.0 && y <= 0.0) 
+    errno = EDOM;
+
   else if ((t = y * log (fabs (x))) >= BIGX) 
     {
       errno = ERANGE;
       if (px & 0x80000000) 
         {
-          /* x is negative. */
-          if (k) 
+          if (!k) 
             {
-              /* y is not an integer. */
               errno = EDOM;
               x = 0.0;
             }
-          else if (exponent_is_even_int)
-            x = z_infinity.d;
-          else
+          else if (sign)
             x = -z_infinity.d;
+          else
+            x =  z_infinity.d;
         }
-      else
-        {
-          x = z_infinity.d;
-        }
-    }
+
+    else 
+      x = z_infinity.d;
+  }
+
   else if (t < SMALLX)
     {
       errno = ERANGE;
       x = 0.0;
     }
+
   else 
     {
-      if ( !k && fabs(d) <= 32767 ) 
+      if ( k && fabs(d) <= 32767 ) 
         {
           n = (int) d;
 
-          if ((sign = (n < 0)))
+          if (sign = (n < 0))
             n = -n;
 
           while ( n > 0 ) 
@@ -129,14 +118,13 @@ double pow (double x, double y)
 
           return r;
         }
+		
       else 
         {
           if ( px & 0x80000000 ) 
             {
-              /* x is negative. */
-              if ( k ) 
+              if ( !k ) 
                 {
-                  /* y is not an integer. */
                   errno = EDOM;
                   return 0.0;
                 }
@@ -144,19 +132,13 @@ double pow (double x, double y)
 
           x = exp (t);
 
-          if (!exponent_is_even_int)
-            {
-              if (px & 0x80000000)
-                {
-                  /* y is an odd integer, and x is negative,
-                     so the result is negative. */
-                  GET_HIGH_WORD (px, x);
-                  px |= 0x80000000;
-                  SET_HIGH_WORD (x, px);
-                }
+          if ( sign ) 
+            { 
+              px ^= 0x80000000;
+              SET_HIGH_WORD (x, px);
             }
         }
-    }
+      }
 
   return x;
 }