aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/ieee754/ldbl-128ibm/s_sincosl.c
diff options
context:
space:
mode:
authorJoseph Myers <joseph@codesourcery.com>2015-02-11 23:17:25 +0000
committerJoseph Myers <joseph@codesourcery.com>2015-02-11 23:17:25 +0000
commitd435569cd626bccb9c27361202a279d0fd7fce88 (patch)
treed307669a2d90be6d52f58a9409f8f5f21575bbe1 /sysdeps/ieee754/ldbl-128ibm/s_sincosl.c
parentef0ffcfd13fea7544c2d83d222963b0da4963734 (diff)
downloadglibc-d435569cd626bccb9c27361202a279d0fd7fce88.zip
glibc-d435569cd626bccb9c27361202a279d0fd7fce88.tar.gz
glibc-d435569cd626bccb9c27361202a279d0fd7fce88.tar.bz2
Fix sincos errno setting (bug 15467).
This patch makes sincos set errno to EDOM when passed an infinity, similarly to sin and cos. Tested for x86_64, x86, powerpc and mips64. I don't know if the architecture-specific implementations for ia64 and m68k might need corresponding fixes. 2015-02-11 Joseph Myers <joseph@codesourcery.com> [BZ #15467] * sysdeps/ieee754/dbl-64/s_sincos.c: Include <errno.h>. (__sincos): Set errno to EDOM for infinite argument. * sysdeps/ieee754/flt-32/s_sincosf.c: Include <errno.h>. (SINCOSF_FUNC): Set errno to EDOM for infinite argument. * sysdeps/ieee754/ldbl-128/s_sincosl.c: Include <errno.h>. (__sincosl): Set errno to EDOM for infinite argument. * sysdeps/ieee754/ldbl-128ibm/s_sincosl.c: Include <errno.h>. (__sincosl): Set errno to EDOM for infinite argument. * sysdeps/ieee754/ldbl-96/s_sincosl.c: Include <errno.h>. (__sincosl): Set errno to EDOM for infinite argument. * math/libm-test.inc (sincos_test_data): Test errno setting.
Diffstat (limited to 'sysdeps/ieee754/ldbl-128ibm/s_sincosl.c')
-rw-r--r--sysdeps/ieee754/ldbl-128ibm/s_sincosl.c3
1 files changed, 3 insertions, 0 deletions
diff --git a/sysdeps/ieee754/ldbl-128ibm/s_sincosl.c b/sysdeps/ieee754/ldbl-128ibm/s_sincosl.c
index f52e733..8a49c54 100644
--- a/sysdeps/ieee754/ldbl-128ibm/s_sincosl.c
+++ b/sysdeps/ieee754/ldbl-128ibm/s_sincosl.c
@@ -18,6 +18,7 @@
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>
@@ -41,6 +42,8 @@ __sincosl (long double x, long double *sinx, long double *cosx)
{
/* sin(Inf or NaN) is NaN */
*sinx = *cosx = x - x;
+ if (__isinf_nsl (x))
+ __set_errno (EDOM);
}
else
{