X-Git-Url: http://erislabs.net/gitweb/?a=blobdiff_plain;f=lib%2Flogl.c;h=de46bd357b4a944dc701a640c9cdd3daf55c0cba;hb=bf18c9b4d8ab7f7e8db272594b5d2c8a8dd06b64;hp=4694645813a14faf21ab167a5b3337670ca0b0f3;hpb=6c38d72875ec56c4b2c250717d458e70b2cb97ef;p=gnulib.git diff --git a/lib/logl.c b/lib/logl.c index 469464581..de46bd357 100644 --- a/lib/logl.c +++ b/lib/logl.c @@ -18,7 +18,7 @@ /* Specification. */ #include -/* 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,7 +185,7 @@ 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 t; @@ -217,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; @@ -241,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). */