From f517e06ab71e07b24bf085063b0ba2d072da4a56 Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Wed, 9 Dec 2015 22:37:08 +0000 Subject: Fix ldbl-128ibm sinhl spurious overflows (bug 19350). The ldbl-128ibm implementation of sinhl uses a slightly too small overflow threshold (similar to bug 16407 for coshl). This patch fixes it to use a safe threshold (so that values whose high part is above the value compared with definitely result in an overflow in all rounding modes). Tested for powerpc. [BZ #19350] * sysdeps/ieee754/ldbl-128ibm/e_sinhl.c (__ieee754_sinhl): Increase overflow threshold. --- ChangeLog | 4 ++++ sysdeps/ieee754/ldbl-128ibm/e_sinhl.c | 2 +- 2 files changed, 5 insertions(+), 1 deletion(-) diff --git a/ChangeLog b/ChangeLog index b76c54b..1c7b198 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,5 +1,9 @@ 2015-12-09 Joseph Myers + [BZ #19350] + * sysdeps/ieee754/ldbl-128ibm/e_sinhl.c (__ieee754_sinhl): + Increase overflow threshold. + [BZ #19349] * sysdeps/ieee754/ldbl-128ibm/s_tanhl.c (__tanhl): Return argument when small. diff --git a/sysdeps/ieee754/ldbl-128ibm/e_sinhl.c b/sysdeps/ieee754/ldbl-128ibm/e_sinhl.c index 29700ad..67d9d24 100644 --- a/sysdeps/ieee754/ldbl-128ibm/e_sinhl.c +++ b/sysdeps/ieee754/ldbl-128ibm/e_sinhl.c @@ -67,7 +67,7 @@ __ieee754_sinhl(long double x) if (ix < 0x40862e42fefa39efLL) return h*__ieee754_expl(fabsl(x)); /* |x| in [log(maxdouble), overflowthresold] */ - if (ix <= 0x408633ce8fb9f87dLL) { + if (ix <= 0x408633ce8fb9f87eLL) { w = __ieee754_expl(0.5*fabsl(x)); t = h*w; return t*w; -- cgit v1.1