summary refs log tree commit diff
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