Fix ldbl-128ibm sinhl spurious overflows (bug 19350)
Joseph Myers
joseph@codesourcery.com
Wed Dec 9 22:38:00 GMT 2015
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.
2015-12-09 Joseph Myers <joseph@codesourcery.com>
[BZ #19350]
* sysdeps/ieee754/ldbl-128ibm/e_sinhl.c (__ieee754_sinhl):
Increase overflow threshold.
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;
--
Joseph S. Myers
joseph@codesourcery.com
More information about the Libc-alpha
mailing list