int main () { test_float (); test_double (); { DECL_LONG_DOUBLE_ROUNDING BEGIN_LONG_DOUBLE_ROUNDING (); test_long_double (); END_LONG_DOUBLE_ROUNDING (); } return 0; }
long double ldexpl (long double x, int exp) { long double factor; int bit; DECL_LONG_DOUBLE_ROUNDING BEGIN_LONG_DOUBLE_ROUNDING (); /* Check for zero, nan and infinity. */ if (!(isnanl (x) || x + x == x)) { if (exp < 0) { exp = -exp; factor = 0.5L; } else 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; } } END_LONG_DOUBLE_ROUNDING (); return x; }