mirror of
https://sourceware.org/git/glibc.git
synced 2024-11-16 01:50:11 +00:00
078d1cf8ac
C99 and C11 allow but do not require ceil, floor, round and trunc to raise the "inexact" exception for noninteger arguments. TS 18661-1 requires that this exception not be raised by these functions. This aligns them with general IEEE semantics, where "inexact" is only raised if the final step of rounding the infinite-precision result to the result type is inexact; for these functions, the infinite-precision integer result is always representable in the result type, so "inexact" should never be raised. The generic implementations of ceil, floor and round functions contain code to force "inexact" to be raised. This patch removes it for round functions to align them with TS 18661-1 in this regard. The tests *are* updated by this patch; there are fewer architecture-specific versions than for ceil and floor, and I fixed the powerpc ones some time ago. If any others still have the issue, as shown by tests for round failing with spurious exceptions, they can be fixed separately by architecture maintainers or others. Tested for x86_64, x86 and mips64. [BZ #15479] * sysdeps/ieee754/dbl-64/s_round.c (huge): Remove variable. (__round): Do not force "inexact" exception. * sysdeps/ieee754/dbl-64/wordsize-64/s_round.c (huge): Remove variable. (__round): Do not force "inexact" exception. * sysdeps/ieee754/flt-32/s_roundf.c (huge): Remove variable. (__roundf): Do not force "inexact" exception. * sysdeps/ieee754/ldbl-128/s_roundl.c (huge): Remove variable. (__roundl): Do not force "inexact" exception. * sysdeps/ieee754/ldbl-96/s_roundl.c (huge): Remove variable. (__roundl): Do not force "inexact" exception. * math/libm-test.inc (round_test_data): Do not allow spurious "inexact" exceptions.
93 lines
2.0 KiB
C
93 lines
2.0 KiB
C
/* Round long double to integer away from zero.
|
|
Copyright (C) 1997-2016 Free Software Foundation, Inc.
|
|
This file is part of the GNU C Library.
|
|
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
|
|
|
|
The GNU C Library is free software; you can redistribute it and/or
|
|
modify it under the terms of the GNU Lesser General Public
|
|
License as published by the Free Software Foundation; either
|
|
version 2.1 of the License, or (at your option) any later version.
|
|
|
|
The GNU C Library is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
Lesser General Public License for more details.
|
|
|
|
You should have received a copy of the GNU Lesser General Public
|
|
License along with the GNU C Library; if not, see
|
|
<http://www.gnu.org/licenses/>. */
|
|
|
|
#include <math.h>
|
|
|
|
#include <math_private.h>
|
|
|
|
|
|
long double
|
|
__roundl (long double x)
|
|
{
|
|
int32_t j0;
|
|
u_int32_t se, i1, i0;
|
|
|
|
GET_LDOUBLE_WORDS (se, i0, i1, x);
|
|
j0 = (se & 0x7fff) - 0x3fff;
|
|
if (j0 < 31)
|
|
{
|
|
if (j0 < 0)
|
|
{
|
|
se &= 0x8000;
|
|
i0 = i1 = 0;
|
|
if (j0 == -1)
|
|
{
|
|
se |= 0x3fff;
|
|
i0 = 0x80000000;
|
|
}
|
|
}
|
|
else
|
|
{
|
|
u_int32_t i = 0x7fffffff >> j0;
|
|
if (((i0 & i) | i1) == 0)
|
|
/* X is integral. */
|
|
return x;
|
|
|
|
u_int32_t j = i0 + (0x40000000 >> j0);
|
|
if (j < i0)
|
|
se += 1;
|
|
i0 = (j & ~i) | 0x80000000;
|
|
i1 = 0;
|
|
}
|
|
}
|
|
else if (j0 > 62)
|
|
{
|
|
if (j0 == 0x4000)
|
|
/* Inf or NaN. */
|
|
return x + x;
|
|
else
|
|
return x;
|
|
}
|
|
else
|
|
{
|
|
u_int32_t i = 0xffffffff >> (j0 - 31);
|
|
if ((i1 & i) == 0)
|
|
/* X is integral. */
|
|
return x;
|
|
|
|
u_int32_t j = i1 + (1 << (62 - j0));
|
|
if (j < i1)
|
|
{
|
|
u_int32_t k = i0 + 1;
|
|
if (k < i0)
|
|
{
|
|
se += 1;
|
|
k |= 0x80000000;
|
|
}
|
|
i0 = k;
|
|
}
|
|
i1 = j;
|
|
i1 &= ~i;
|
|
}
|
|
|
|
SET_LDOUBLE_WORDS (x, se, i0, i1);
|
|
return x;
|
|
}
|
|
weak_alias (__roundl, roundl)
|