X-Git-Url: http://erislabs.net/gitweb/?a=blobdiff_plain;f=lib%2Fldexpl.c;h=4389f761b2442169d409ccb04bfee067fa4e4e47;hb=b4c2b236d0007a7f9226e6ccef82e9e37bdb2b3a;hp=60474136bbf6a2350dd708bac28d9759018630fc;hpb=57fdfd3f8ec62b105c53bcdf6f127c35c7fe7391;p=gnulib.git diff --git a/lib/ldexpl.c b/lib/ldexpl.c index 60474136b..4389f761b 100644 --- a/lib/ldexpl.c +++ b/lib/ldexpl.c @@ -1,7 +1,7 @@ /* Emulation for ldexpl. Contributed by Paolo Bonzini - Copyright 2002, 2003, 2007 Free Software Foundation, Inc. + Copyright 2002-2003, 2007-2011 Free Software Foundation, Inc. This file is part of gnulib. @@ -25,10 +25,9 @@ #include #include "fpucw.h" -#include "isnanl.h" long double -ldexpl(long double x, int exp) +ldexpl (long double x, int exp) { long double factor; int bit; @@ -40,25 +39,25 @@ ldexpl(long double x, int exp) if (!(isnanl (x) || x + x == x)) { if (exp < 0) - { - exp = -exp; - factor = 0.5L; - } + { + exp = -exp; + factor = 0.5L; + } else - factor = 2.0L; + factor = 2.0L; if (exp > 0) - for (bit = 1;;) - { - /* Invariant: Here bit = 2^i, factor = 2^-2^i or = 2^2^i, - and bit <= exp. */ - if (exp & bit) - x *= factor; - bit <<= 1; - if (bit > exp) - break; - factor = factor * factor; - } + for (bit = 1;;) + { + /* Invariant: Here bit = 2^i, factor = 2^-2^i or = 2^2^i, + and bit <= exp. */ + if (exp & bit) + x *= factor; + bit <<= 1; + if (bit > exp) + break; + factor = factor * factor; + } } END_LONG_DOUBLE_ROUNDING (); @@ -73,6 +72,6 @@ main (void) long double x; int y; for (y = 0; y < 29; y++) - printf ("%5d %.16Lg %.16Lg\n", y, ldexpl(0.8L, y), ldexpl(0.8L, -y) * ldexpl(0.8L, y)); + printf ("%5d %.16Lg %.16Lg\n", y, ldexpl (0.8L, y), ldexpl (0.8L, -y) * ldexpl (0.8L, y)); } #endif