about summary refs log tree commit diff
path: root/sysdeps/ieee754/ldbl-128ibm/w_expl.c
diff options
context:
space:
mode:
authorAdhemerval Zanella <azanella@linux.vnet.ibm.com>2013-03-22 12:39:10 -0300
committerAdhemerval Zanella <azanella@linux.vnet.ibm.com>2013-03-22 12:39:10 -0300
commite42a38dd9dd4bbeb0fbd6e99f35d796ba37b6879 (patch)
tree6d86cd75766421301163df58b8dd3284acb441ad /sysdeps/ieee754/ldbl-128ibm/w_expl.c
parent2e0fb52187504fad6657c9462ce650a540d5e387 (diff)
downloadglibc-e42a38dd9dd4bbeb0fbd6e99f35d796ba37b6879.tar.gz
glibc-e42a38dd9dd4bbeb0fbd6e99f35d796ba37b6879.tar.xz
glibc-e42a38dd9dd4bbeb0fbd6e99f35d796ba37b6879.zip
BZ#13889: expl (709.75) wrongly overflows for ldbl-128ibm
The patch increase the high value to check if expl overflows. Current
high mark value is not really correct, the algorithm accepts high values.
It also adds a correct wrapper function to check for overflow and underflow.
Diffstat (limited to 'sysdeps/ieee754/ldbl-128ibm/w_expl.c')
-rw-r--r--sysdeps/ieee754/ldbl-128ibm/w_expl.c26
1 files changed, 22 insertions, 4 deletions
diff --git a/sysdeps/ieee754/ldbl-128ibm/w_expl.c b/sysdeps/ieee754/ldbl-128ibm/w_expl.c
index a5e72b2170..70fe5f693e 100644
--- a/sysdeps/ieee754/ldbl-128ibm/w_expl.c
+++ b/sysdeps/ieee754/ldbl-128ibm/w_expl.c
@@ -1,6 +1,24 @@
-/* Looks like we can use ieee854 w_expl.c as is for IBM extended format. */
+#include <math.h>
+#include <math_private.h>
 #include <math_ldbl_opt.h>
-#undef weak_alias
-#define weak_alias(n,a)
-#include <sysdeps/ieee754/ldbl-128/w_expl.c>
+
+static const long double o_thres = 709.78271289338399678773454114191496482L;
+static const long double u_thres = -744.44007192138126231410729844608163411L;
+
+long double __expl(long double x)	/* wrapper exp  */
+{
+  long double z;
+  z = __ieee754_expl(x);
+  if (_LIB_VERSION == _IEEE_)
+    return z;
+  if (__finitel(x))
+    {
+      if (x >= o_thres)
+	return __kernel_standard_l(x,x,206); /* exp overflow  */
+      else if (x <= u_thres)
+	return __kernel_standard_l(x,x,207); /* exp underflow  */
+    }
+  return z;
+}
+hidden_def (__expl)
 long_double_symbol (libm, __expl, expl);