aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--ChangeLog10
-rw-r--r--NEWS4
-rw-r--r--aclocal.m44
-rw-r--r--sysdeps/ieee754/ldbl-128/e_coshl.c103
-rw-r--r--sysdeps/powerpc/fpu/libm-test-ulps6
-rw-r--r--sysdeps/s390/fpu/libm-test-ulps6
6 files changed, 129 insertions, 4 deletions
diff --git a/ChangeLog b/ChangeLog
index 5bc9902cff..5a1923e63e 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,13 @@
+2001-06-11 Andreas Jaeger <aj@suse.de>
+
+ * sysdeps/ieee754/ldbl-128/e_coshl.c: New file. Contributed by
+ Stephen L Moshier <moshier@mediaone.net>.
+
+ * aclocal.m4: Quote AC_FD_MSG and AC_FD_CC.
+
+ * sysdeps/powerpc/fpu/libm-test-ulps: Add ULPs for tanh (-0.7).
+ * sysdeps/s390/fpu/libm-test-ulps: Likewise.
+
2001-06-11 Michael Deutschmann <michael@talamasca.ocis.net>
* rt/tst-aio4.c (do_test): Test whether rt signals are supported.
diff --git a/NEWS b/NEWS
index e4deae51b7..d0e7f899d7 100644
--- a/NEWS
+++ b/NEWS
@@ -8,8 +8,8 @@ Please send GNU C library bug reports using the `glibcbug' script to
Version 2.2.4
-* Stephen Moshier implemented expm1, log1p, acos, sinh for the 128-bit
- long double format.
+* Stephen Moshier implemented cosh, expm1, log1p, acos, sinh for the
+ 128-bit long double format.
Version 2.2.3
diff --git a/aclocal.m4 b/aclocal.m4
index 5affe451c4..ba9ca30e74 100644
--- a/aclocal.m4
+++ b/aclocal.m4
@@ -3,8 +3,8 @@ dnl which appears in configure.in before the sysdep configure scripts are run.
dnl Each sysdep configure.in does GLIBC_PROVIDES first, to avoid any
dnl AC_REQUIREs or AC_BEFOREs duplicating their code.
dnl
-define(AC_FD_MSG,6)dnl Autoconf lossage.
-define(AC_FD_CC,5)dnl Autoconf lossage.
+define([AC_FD_MSG],6)dnl Autoconf lossage.
+define([AC_FD_CC],5)dnl Autoconf lossage.
AC_DEFUN([GLIBC_PROVIDES], [dnl
AC_PROVIDE([AC_PROG_INSTALL])dnl
AC_PROVIDE([AC_PROG_RANLIB])dnl
diff --git a/sysdeps/ieee754/ldbl-128/e_coshl.c b/sysdeps/ieee754/ldbl-128/e_coshl.c
new file mode 100644
index 0000000000..7ffc4d5159
--- /dev/null
+++ b/sysdeps/ieee754/ldbl-128/e_coshl.c
@@ -0,0 +1,103 @@
+/*
+ * ====================================================
+ * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
+ *
+ * Developed at SunPro, a Sun Microsystems, Inc. business.
+ * Permission to use, copy, modify, and distribute this
+ * software is freely granted, provided that this notice
+ * is preserved.
+ * ====================================================
+ */
+
+/* Changes for 128-bit long double contributed by
+ Stephen L. Moshier <moshier@na-net.ornl.gov> */
+
+/* __ieee754_coshl(x)
+ * Method :
+ * mathematically coshl(x) if defined to be (exp(x)+exp(-x))/2
+ * 1. Replace x by |x| (coshl(x) = coshl(-x)).
+ * 2.
+ * [ exp(x) - 1 ]^2
+ * 0 <= x <= ln2/2 : coshl(x) := 1 + -------------------
+ * 2*exp(x)
+ *
+ * exp(x) + 1/exp(x)
+ * ln2/2 <= x <= 22 : coshl(x) := -------------------
+ * 2
+ * 22 <= x <= lnovft : coshl(x) := expl(x)/2
+ * lnovft <= x <= ln2ovft: coshl(x) := expl(x/2)/2 * expl(x/2)
+ * ln2ovft < x : coshl(x) := huge*huge (overflow)
+ *
+ * Special cases:
+ * coshl(x) is |x| if x is +INF, -INF, or NaN.
+ * only coshl(0)=1 is exact for finite x.
+ */
+
+#include "math.h"
+#include "math_private.h"
+
+#ifdef __STDC__
+static const long double one = 1.0, half = 0.5, huge = 1.0e4900L,
+ovf_thresh = 1.1357216553474703894801348310092223067821E4L;
+#else
+static long double one = 1.0, half = 0.5, huge = 1.0e4900L,
+ovf_thresh = 1.1357216553474703894801348310092223067821E4L;
+#endif
+
+#ifdef __STDC__
+long double
+__ieee754_coshl (long double x)
+#else
+long double
+__ieee754_coshl (x)
+ long double x;
+#endif
+{
+ long double t, w;
+ int32_t ex;
+ u_int32_t mx, lx;
+ ieee854_long_double_shape_type u;
+
+ u.value = x;
+ ex = u.parts32.w0 & 0x7fffffff;
+
+ /* Absolute value of x. */
+ u.parts32.w0 = ex;
+
+ /* x is INF or NaN */
+ if (ex >= 0x7fff0000)
+ return x * x;
+
+ /* |x| in [0,0.5*ln2], return 1+expm1l(|x|)^2/(2*expl(|x|)) */
+ if (ex < 0x3ffd62e4) /* 0.3465728759765625 */
+ {
+ t = __expm1l (u.value);
+ w = one + t;
+ if (ex < 0x3fc60000) /* |x| < 2^-57 */
+ return w; /* cosh(tiny) = 1 */
+
+ return one + (t * t) / (w + w);
+ }
+
+ /* |x| in [0.5*ln2,40], return (exp(|x|)+1/exp(|x|)/2; */
+ if (ex < 0x40044000)
+ {
+ t = __ieee754_expl (u.value);
+ return half * t + half / t;
+ }
+
+ /* |x| in [22, ln(maxdouble)] return half*exp(|x|) */
+ if (ex <= 0x400c62e3) /* 11356.375 */
+ return half * __ieee754_expl (u.value);
+
+ /* |x| in [log(maxdouble), overflowthresold] */
+ if (u.value <= ovf_thresh)
+ {
+ w = __ieee754_expl (half * u.value);
+ t = half * w;
+ return t * w;
+ }
+
+ /* |x| > overflowthresold, cosh(x) overflow */
+ return huge * huge;
+}
diff --git a/sysdeps/powerpc/fpu/libm-test-ulps b/sysdeps/powerpc/fpu/libm-test-ulps
index e985f3be2c..b1d43323e4 100644
--- a/sysdeps/powerpc/fpu/libm-test-ulps
+++ b/sysdeps/powerpc/fpu/libm-test-ulps
@@ -597,6 +597,12 @@ float: 1
idouble: 1
ifloat: 1
+Test "tanh (-0.7) == -0.60436777711716349631":
+double: 1
+float: 1
+idouble: 1
+ifloat: 1
+
# tgamma
Test "tgamma (-0.5) == -2 sqrt (pi)":
double: 1
diff --git a/sysdeps/s390/fpu/libm-test-ulps b/sysdeps/s390/fpu/libm-test-ulps
index c1463c92d8..675ceeb2a7 100644
--- a/sysdeps/s390/fpu/libm-test-ulps
+++ b/sysdeps/s390/fpu/libm-test-ulps
@@ -640,6 +640,12 @@ float: 1
idouble: 1
ifloat: 1
+Test "tanh (-0.7) == -0.60436777711716349631":
+double: 1
+float: 1
+idouble: 1
+ifloat: 1
+
# tgamma
Test "tgamma (-0.5) == -2 sqrt (pi)":
double: 1