mingw-wsl-patches

Fork

# HG changeset patch
# Parent c4d8313f54f1f8a436c4ca0fd692e675e91746d6
diff --git a/mingwrt/mingwex/math/powl.c b/mingwrt/mingwex/math/powl.c
--- a/mingwrt/mingwex/math/powl.c
+++ b/mingwrt/mingwex/math/powl.c
@@ -450,34 +450,41 @@ long double powl( long double x, long do
 {
   /* double F, Fa, Fb, G, Ga, Gb, H, Ha, Hb */
   int i, nflg, iyflg, yoddint;
   long e;
 
-  if( y == 0.0L )
+  if( (x == 1.0L) || (y == 0.0L) )
     return( 1.0L );
 
 #ifdef NANS
   if( isnanl(x) )
   {
-    _SET_ERRNO (EDOM);
+    //_SET_ERRNO (EDOM);
     return( x );
   }
   if( isnanl(y) )
   {
-    _SET_ERRNO (EDOM);
+    //_SET_ERRNO (EDOM);
     return( y );
   }
 #endif
 
   if( y == 1.0L )
     return( x );
 
-  if( isinfl(y) && (x == -1.0L || x == 1.0L) )
-    return( y );
+  if( isinfl(y) )
+  {
+    if( x == -1.0L )
+      return( 1.0L );
+    else if( (-1.0L < x) && (x < 1.0L) )
+      return signbitl(y) ? INFINITYL : 0.0L;
+    else
+      return signbitl(y) ? 0.0L : INFINITYL;
+  }
 
-  if( x == 1.0L )
-    return( 1.0L );
+  if( (x == 0.0L) && (y < 0.0L) )
+    _SET_ERRNO (ERANGE);
 
   if( y >= MAXNUML )
   {
     _SET_ERRNO (ERANGE);
 #ifdef INFINITIES