diff options
Diffstat (limited to 'sysdeps/ieee754/ldbl-128/s_roundl.c')
| -rw-r--r-- | sysdeps/ieee754/ldbl-128/s_roundl.c | 52 |
1 files changed, 20 insertions, 32 deletions
diff --git a/sysdeps/ieee754/ldbl-128/s_roundl.c b/sysdeps/ieee754/ldbl-128/s_roundl.c index 98b448e4f3..22789cedf3 100644 --- a/sysdeps/ieee754/ldbl-128/s_roundl.c +++ b/sysdeps/ieee754/ldbl-128/s_roundl.c @@ -1,5 +1,5 @@ /* Round long double to integer away from zero. - Copyright (C) 1997-2016 Free Software Foundation, Inc. + Copyright (C) 1997-2018 Free Software Foundation, Inc. This file is part of the GNU C Library. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997 and Jakub Jelinek <jj@ultra.linux.cz>, 1999. @@ -21,16 +21,14 @@ #include <math.h> #include <math_private.h> +#include <libm-alias-ldouble.h> -static const long double huge = 1.0E4930L; - - -long double -__roundl (long double x) +_Float128 +__roundl (_Float128 x) { int32_t j0; - u_int64_t i1, i0; + uint64_t i1, i0; GET_LDOUBLE_WORDS64 (i0, i1, x); j0 = ((i0 >> 48) & 0x7fff) - 0x3fff; @@ -38,27 +36,21 @@ __roundl (long double x) { if (j0 < 0) { - if (huge + x > 0.0) - { - i0 &= 0x8000000000000000ULL; - if (j0 == -1) - i0 |= 0x3fff000000000000LL; - i1 = 0; - } + i0 &= 0x8000000000000000ULL; + if (j0 == -1) + i0 |= 0x3fff000000000000LL; + i1 = 0; } else { - u_int64_t i = 0x0000ffffffffffffLL >> j0; + uint64_t i = 0x0000ffffffffffffLL >> j0; if (((i0 & i) | i1) == 0) /* X is integral. */ return x; - if (huge + x > 0.0) - { - /* Raise inexact if x != 0. */ - i0 += 0x0000800000000000LL >> j0; - i0 &= ~i; - i1 = 0; - } + + i0 += 0x0000800000000000LL >> j0; + i0 &= ~i; + i1 = 0; } } else if (j0 > 111) @@ -71,23 +63,19 @@ __roundl (long double x) } else { - u_int64_t i = -1ULL >> (j0 - 48); + uint64_t i = -1ULL >> (j0 - 48); if ((i1 & i) == 0) /* X is integral. */ return x; - if (huge + x > 0.0) - { - /* Raise inexact if x != 0. */ - u_int64_t j = i1 + (1LL << (111 - j0)); - if (j < i1) - i0 += 1; - i1 = j; - } + uint64_t j = i1 + (1LL << (111 - j0)); + if (j < i1) + i0 += 1; + i1 = j; i1 &= ~i; } SET_LDOUBLE_WORDS64 (x, i0, i1); return x; } -weak_alias (__roundl, roundl) +libm_alias_ldouble (__round, round) |