X-Git-Url: http://erislabs.net/gitweb/?a=blobdiff_plain;f=lib%2Fexpl.c;h=407da5fb883a68bc04287320038895876502ec32;hb=41ffca9515f41dc46cc2e0296b67db431cb075fc;hp=202d35952ec8465227ed876442e437b825deffa5;hpb=57fdfd3f8ec62b105c53bcdf6f127c35c7fe7391;p=gnulib.git diff --git a/lib/expl.c b/lib/expl.c index 202d35952..407da5fb8 100644 --- a/lib/expl.c +++ b/lib/expl.c @@ -1,7 +1,7 @@ /* Emulation for expl. Contributed by Paolo Bonzini - Copyright 2002, 2003, 2007 Free Software Foundation, Inc. + Copyright 2002, 2003, 2007, 2009, 2010 Free Software Foundation, Inc. This file is part of gnulib. @@ -72,16 +72,16 @@ expl (long double x) x -= exponent / 8.0L; if (x > 0.0625) - { - exponent++; - x -= 0.125L; - } + { + exponent++; + x -= 0.125L; + } if (exponent < 0) { t = 0.8824969025845954028648921432290507362220L; /* e^-0.25 */ - exponent = -exponent; - } + exponent = -exponent; + } else t = 1.1331484530668263168290072278117938725655L; /* e^0.25 */ @@ -97,8 +97,8 @@ expl (long double x) } /* Approximate (e^x - 1)/x, using a seventh-degree polynomial, - with maximum error in [-2^-16-2^-53,2^-16+2^-53] - less than 4.8e-39. */ + with maximum error in [-2^-16-2^-53,2^-16+2^-53] + less than 4.8e-39. */ x22 = x + x*x*(P1+x*(P2+x*(P3+x*(P4+x*(P5+x*P6))))); return result + result * x22; @@ -107,11 +107,11 @@ expl (long double x) else if (x < himark) { if (x + x == x) - /* e^-inf == 0, with no error. */ - return 0; + /* e^-inf == 0, with no error. */ + return 0; else - /* Underflow */ - return TINY * TINY; + /* Underflow */ + return TINY * TINY; } else /* Return x, if x is a NaN or Inf; or overflow, otherwise. */ @@ -122,14 +122,14 @@ expl (long double x) int main (void) { - printf ("%.16Lg\n", expl(1.0L)); - printf ("%.16Lg\n", expl(-1.0L)); - printf ("%.16Lg\n", expl(2.0L)); - printf ("%.16Lg\n", expl(4.0L)); - printf ("%.16Lg\n", expl(-2.0L)); - printf ("%.16Lg\n", expl(-4.0L)); - printf ("%.16Lg\n", expl(0.0625L)); - printf ("%.16Lg\n", expl(0.3L)); - printf ("%.16Lg\n", expl(0.6L)); + printf ("%.16Lg\n", expl (1.0L)); + printf ("%.16Lg\n", expl (-1.0L)); + printf ("%.16Lg\n", expl (2.0L)); + printf ("%.16Lg\n", expl (4.0L)); + printf ("%.16Lg\n", expl (-2.0L)); + printf ("%.16Lg\n", expl (-4.0L)); + printf ("%.16Lg\n", expl (0.0625L)); + printf ("%.16Lg\n", expl (0.3L)); + printf ("%.16Lg\n", expl (0.6L)); } #endif