Correctly determine whether pow is available in libc on AIX 7 with xlc.

[gnulib.git] / lib / logl.c

diff --git a/lib/logl.c b/lib/logl.c
index ffa6a95..de46bd3 100644 (file)
--- a/lib/logl.c
+++ b/lib/logl.c
@@ -18,7 +18,7 @@
 /* Specification.  */
 #include <math.h>
 
-/*                                                     logll.c
+/*                                                      logll.c
  *
  * Natural logarithm for 128-bit long double precision.
  *
@@ -62,8 +62,6 @@
  *
  */
 
-#include "isnanl.h"
-
 /* log(1+x) = x - .5 x^2 + x^3 l(x)
    -.0078125 <= x <= +.0078125
    peak relative error 1.2e-37 */
@@ -187,11 +185,10 @@ static const long double
   ln2b = 1.4286068203094172321214581765680755001344E-6L;
 
 long double
-logl(long double x)
+logl (long double x)
 {
   long double z, y, w;
-  long double u, t;
-  unsigned int m;
+  long double t;
   int k, e;
 
   /* Check for IEEE special cases.  */
@@ -218,7 +215,7 @@ logl(long double x)
     }
 
   /* Extract exponent and reduce domain to 0.703125 <= u < 1.40625  */
-  x = frexpl(x, &e);
+  x = frexpl (x, &e);
   if (x < 0.703125L)
     {
       x += x;
@@ -242,17 +239,17 @@ logl(long double x)
   /* Series expansion of log(1+z).  */
   w = z * z;
   y = ((((((((((((l15 * z
-                 + l14) * z
-                + l13) * z
-               + l12) * z
-              + l11) * z
-             + l10) * z
-            + l9) * z
-           + l8) * z
-          + l7) * z
-         + l6) * z
-        + l5) * z
-       + l4) * z
+                  + l14) * z
+                 + l13) * z
+                + l12) * z
+               + l11) * z
+              + l10) * z
+             + l9) * z
+            + l8) * z
+           + l7) * z
+          + l6) * z
+         + l5) * z
+        + l4) * z
        + l3) * z * w;
   y -= 0.5 * w;
   y += e * ln2b;  /* Base 2 exponent offset times ln(2).  */