diff options
Diffstat (limited to 'libquadmath/math/llrintq.c')
-rw-r--r-- | libquadmath/math/llrintq.c | 43 |
1 files changed, 38 insertions, 5 deletions
diff --git a/libquadmath/math/llrintq.c b/libquadmath/math/llrintq.c index eef31d823b6..a6a0ae64bd3 100644 --- a/libquadmath/math/llrintq.c +++ b/libquadmath/math/llrintq.c @@ -1,9 +1,9 @@ /* Round argument to nearest integral value according to current rounding direction. - Copyright (C) 1997, 1999, 2006 Free Software Foundation, Inc. + Copyright (C) 1997-2017 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. + Jakub Jelinek <jj@ultra.linux.cz>, 1999. The GNU C Library is free software; you can redistribute it and/or modify it under the terms of the GNU Lesser General Public @@ -33,7 +33,7 @@ llrintq (__float128 x) { int32_t j0; uint64_t i0,i1; - volatile __float128 w; + __float128 w; __float128 t; long long int result; int sx; @@ -46,8 +46,23 @@ llrintq (__float128 x) if (j0 < (int32_t) (8 * sizeof (long long int)) - 1) { - w = two112[sx] + x; - t = w - two112[sx]; +#if defined FE_INVALID || defined FE_INEXACT + /* X < LLONG_MAX + 1 implied by J0 < 63. */ + if (x > (__float128) LLONG_MAX) + { + /* In the event of overflow we must raise the "invalid" + exception, but not "inexact". */ + t = nearbyintq (x); +#ifdef USE_FENV_H + feraiseexcept (t == LLONG_MAX ? FE_INEXACT : FE_INVALID); +#endif + } + else +#endif + { + w = two112[sx] + x; + t = w - two112[sx]; + } GET_FLT128_WORDS64 (i0, i1, t); j0 = ((i0 >> 48) & 0x7fff) - 0x3fff; i0 &= 0x0000ffffffffffffLL; @@ -62,6 +77,24 @@ llrintq (__float128 x) } else { + /* The number is too large. Unless it rounds to LLONG_MIN, + FE_INVALID must be raised and the return value is + unspecified. */ +#if defined FE_INVALID || defined FE_INEXACT + if (x < (__float128) LLONG_MIN + && x > (__float128) LLONG_MIN - 1.0Q) + { + /* If truncation produces LLONG_MIN, the cast will not raise + the exception, but may raise "inexact". */ + t = nearbyintq (x); +#ifdef USE_FENV_H + feraiseexcept (t == LLONG_MIN ? FE_INEXACT : FE_INVALID); +#endif + return LLONG_MIN; + } + +#endif + /* The number is too large. It is left implementation defined what happens. */ return (long long int) x; |