summaryrefslogtreecommitdiff
path: root/libc/sysdeps/ieee754/ldbl-128ibm/w_expl.c
diff options
context:
space:
mode:
authorjoseph <joseph@7b3dc134-2b1b-0410-93df-9e9f96275f8d>2013-04-28 15:21:50 +0000
committerjoseph <joseph@7b3dc134-2b1b-0410-93df-9e9f96275f8d>2013-04-28 15:21:50 +0000
commitc19459a92177fc8864b461e84d017a21c4006dc5 (patch)
treefcf20ad0222de01ad3a955edb3b65b06724a3a15 /libc/sysdeps/ieee754/ldbl-128ibm/w_expl.c
parent8751114637bcc3caaf16a4216da0afb84456558a (diff)
downloadeglibc2-c19459a92177fc8864b461e84d017a21c4006dc5.tar.gz
Merge changes between r22663 and r22954 from /fsf/trunk.
git-svn-id: svn://svn.eglibc.org/trunk@22955 7b3dc134-2b1b-0410-93df-9e9f96275f8d
Diffstat (limited to 'libc/sysdeps/ieee754/ldbl-128ibm/w_expl.c')
-rw-r--r--libc/sysdeps/ieee754/ldbl-128ibm/w_expl.c26
1 files changed, 22 insertions, 4 deletions
diff --git a/libc/sysdeps/ieee754/ldbl-128ibm/w_expl.c b/libc/sysdeps/ieee754/ldbl-128ibm/w_expl.c
index a5e72b217..70fe5f693 100644
--- a/libc/sysdeps/ieee754/ldbl-128ibm/w_expl.c
+++ b/libc/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);