8795b4a443
This patch fixes the default mode of scalb to set errno (bugs 6803 and 6804). Previously, the _LIB_VERSION == _SVID_ mode would set errno but only in some relevant cases, and with various peculiarities (such as errno setting when an exact infinity or zero result arises with an argument to scalb being an infinity). This patch leaves this mode bug-compatible, while making the default mode set errno in accordance with normal practice (so an exact infinity from an infinite argument is not an error, and nor is an exact zero result). gen-libm-test.pl is taught new notation such as ERRNO_PLUS_OFLOW to facilitate writing the tests of errno setting for underflow / overflow in libm-test.inc. Note that bug 6803 also covers scalbn and scalbln, but this patch only addresses the scalb parts of that bug (along with the whole of bug 6804). Tested x86_64 and x86. [BZ #6803] [BZ #6804] * math/w_scalb.c (__scalb): For non-SVID mode, check result and set errno as appropriate. * math/w_scalbf.c (__scalbf): Likewise. * math/w_scalbl.c (__scalbl): Likewise. * math/gen-libm-test.pl (parse_args): Handle ERRNO_PLUS_OFLOW, ERRNO_MINUS_OFLOW, ERRNO_PLUS_UFLOW and ERRNO_MINUS_UFLOW. * math/libm-test.inc (scalb_test_data): Add errno expectations. Add more NaN tests.
77 lines
2.0 KiB
C
77 lines
2.0 KiB
C
/* Copyright (C) 2011-2014 Free Software Foundation, Inc.
|
|
This file is part of the GNU C Library.
|
|
Contributed by Ulrich Drepper <drepper@gmail.com>, 2011.
|
|
|
|
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 <errno.h>
|
|
#include <math.h>
|
|
#include <math_private.h>
|
|
|
|
|
|
static long double
|
|
__attribute__ ((noinline))
|
|
sysv_scalbl (long double x, long double fn)
|
|
{
|
|
long double z = __ieee754_scalbl (x, fn);
|
|
|
|
if (__glibc_unlikely (__isinfl (z)))
|
|
{
|
|
if (__finitel (x))
|
|
return __kernel_standard_l (x, fn, 232); /* scalb overflow */
|
|
else
|
|
__set_errno (ERANGE);
|
|
}
|
|
else if (__builtin_expect (z == 0.0L, 0) && z != x)
|
|
return __kernel_standard_l (x, fn, 233); /* scalb underflow */
|
|
|
|
return z;
|
|
}
|
|
|
|
|
|
/* Wrapper scalbl */
|
|
long double
|
|
__scalbl (long double x, long double fn)
|
|
{
|
|
if (__glibc_unlikely (_LIB_VERSION == _SVID_))
|
|
return sysv_scalbl (x, fn);
|
|
else
|
|
{
|
|
long double z = __ieee754_scalbl (x, fn);
|
|
|
|
if (__glibc_unlikely (!__finitel (z) || z == 0.0L))
|
|
{
|
|
if (__isnanl (z))
|
|
{
|
|
if (!__isnanl (x) && !__isnanl (fn))
|
|
__set_errno (EDOM);
|
|
}
|
|
else if (__isinf_nsl (z))
|
|
{
|
|
if (!__isinf_nsl (x) && !__isinf_nsl (fn))
|
|
__set_errno (ERANGE);
|
|
}
|
|
else
|
|
{
|
|
/* z == 0. */
|
|
if (x != 0.0L && !__isinf_nsl (fn))
|
|
__set_errno (ERANGE);
|
|
}
|
|
}
|
|
return z;
|
|
}
|
|
}
|
|
weak_alias (__scalbl, scalbl)
|