This is the mail archive of the
libc-alpha@sourceware.org
mailing list for the glibc project.
[PATCH 2.19] BZ16427: Fix ldbl-128 exp overflows
- From: Andreas Krebbel <krebbel at linux dot vnet dot ibm dot com>
- To: libc-alpha at sourceware dot org
- Date: Tue, 14 Jan 2014 09:29:48 +0100
- Subject: [PATCH 2.19] BZ16427: Fix ldbl-128 exp overflows
- Authentication-results: sourceware.org; auth=none
Hi,
the attached patch fixes the exp failures with long double 128 on
S/390.
Ok to apply?
Bye,
-Andreas-
2014-01-14 Andreas Krebbel <Andreas.Krebbel@de.ibm.com>
[BZ #16427]
* sysdeps/ieee754/ldbl-128/w_expl.c (__expl): Invoke the non-IEEE
handling only for numbers special also in the IEEE case.
commit 1e74b6daed6c95628a4710f64d140880457c9c44
Author: Andreas Krebbel <krebbel@linux.vnet.ibm.com>
Date: Mon Jan 13 14:41:47 2014 +0000
S/390: Fix ldbl-128 exp overflows.
diff --git a/sysdeps/ieee754/ldbl-128/w_expl.c b/sysdeps/ieee754/ldbl-128/w_expl.c
index 10193be..af37aba 100644
--- a/sysdeps/ieee754/ldbl-128/w_expl.c
+++ b/sysdeps/ieee754/ldbl-128/w_expl.c
@@ -25,24 +25,16 @@ static char rcsid[] = "$NetBSD: $";
#include <math.h>
#include <math_private.h>
-static const long double
-o_threshold= 1.1356523406294143949491931077970763428449E4L,
-u_threshold= -1.1433462743336297878837243843452621503410E4;
-
long double __expl(long double x) /* wrapper exp */
{
#ifdef _IEEE_LIBM
return __ieee754_expl(x);
#else
- long double z;
- z = __ieee754_expl(x);
- if(_LIB_VERSION == _IEEE_) return z;
- if(__finitel(x)) {
- if(x>o_threshold)
- return __kernel_standard_l(x,x,206); /* exp overflow */
- else if(x<u_threshold)
- return __kernel_standard_l(x,x,207); /* exp underflow */
- }
+ long double z = __ieee754_expl (x);
+ if (__builtin_expect (!__finitel (z) || z == 0, 0)
+ && __finitel (x) && _LIB_VERSION != _IEEE_)
+ return __kernel_standard_l (x, x, 206 + !!__signbitl (x));
+
return z;
#endif
}