X-Git-Url: http://erislabs.net/gitweb/?a=blobdiff_plain;f=lib%2Flogl.c;h=de46bd357b4a944dc701a640c9cdd3daf55c0cba;hb=620ece7758e453efebf7f52a1ab9fc1c49a5020c;hp=52adaecc17d08469521f4180efeb227c23601ac6;hpb=febae7ee2b8b61030ed32d2b5b53316b240434f2;p=gnulib.git diff --git a/lib/logl.c b/lib/logl.c index 52adaecc1..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. * @@ -185,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; @@ -215,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; @@ -239,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). */