Update.

2001-11-26  Stephen L Moshier <moshier@mediaone.net>

	* sysdeps/ieee754/ldbl-128/e_log10l.c (L102B): Fix typo in
	initialization.

	* sysdeps/ieee754/ldbl-128/e_coshl.c (__ieee754_coshl): Remove
	unused variables.  Fix threshold for returning 1.0.

	* sysdeps/ieee754/ldbl-128/e_lgammal_r.c (neval, deval): Constify
	argument.

	* sysdeps/ieee754/ldbl-128/s_erfl.c (neval, deval): Likewise, and
	remove unused variables.

	* sysdeps/ieee754/ldbl-128/e_j0l.c (neval, deval): Likewise.

	* sysdeps/ieee754/ldbl-128/e_j1l.c (neval, deval): Likewise, and
	include local math header files.

5 files changed, 14 insertions, 18 deletions

diff --git a/sysdeps/ieee754/ldbl-128/e_coshl.c b/sysdeps/ieee754/ldbl-128/e_coshl.c
index 7ffc4d5159..90c5f4fd6a 100644
--- a/sysdeps/ieee754/ldbl-128/e_coshl.c
+++ b/sysdeps/ieee754/ldbl-128/e_coshl.c
@@ -55,7 +55,6 @@ __ieee754_coshl (x)
 {
   long double t, w;
   int32_t ex;
-  u_int32_t mx, lx;
   ieee854_long_double_shape_type u;
 
   u.value = x;
@@ -73,7 +72,7 @@ __ieee754_coshl (x)
     {
       t = __expm1l (u.value);
       w = one + t;
-      if (ex < 0x3fc60000) /* |x| < 2^-57 */
+      if (ex < 0x3fb80000) /* |x| < 2^-116 */
 	return w;		/* cosh(tiny) = 1 */
 
       return one + (t * t) / (w + w);
diff --git a/sysdeps/ieee754/ldbl-128/e_j0l.c b/sysdeps/ieee754/ldbl-128/e_j0l.c
index a804f0a276..dbe414d2b1 100644
--- a/sysdeps/ieee754/ldbl-128/e_j0l.c
+++ b/sysdeps/ieee754/ldbl-128/e_j0l.c
@@ -625,7 +625,7 @@ static const long double Q2_2r3D[NQ2_2r3D + 1] = {
 /* Evaluate P[n] x^n  +  P[n-1] x^(n-1)  +  ...  +  P[0] */
 
 static long double
-neval (long double x, long double *p, int n)
+neval (long double x, const long double *p, int n)
 {
   long double y;
 
@@ -643,7 +643,7 @@ neval (long double x, long double *p, int n)
 /* Evaluate x^n+1  +  P[n] x^(n)  +  P[n-1] x^(n-1)  +  ...  +  P[0] */
 
 static long double
-deval (long double x, long double *p, int n)
+deval (long double x, const long double *p, int n)
 {
   long double y;
 
diff --git a/sysdeps/ieee754/ldbl-128/e_j1l.c b/sysdeps/ieee754/ldbl-128/e_j1l.c
index ad899be3de..052f2c18eb 100644
--- a/sysdeps/ieee754/ldbl-128/e_j1l.c
+++ b/sysdeps/ieee754/ldbl-128/e_j1l.c
@@ -81,11 +81,8 @@
 
 /* Copyright 2001 by Stephen L. Moshier (moshier@na-net.onrl.gov).  */
 
-extern long double fabsl (long double);
-extern long double cosl (long double);
-extern long double sinl (long double);
-extern long double sqrtl (long double);
-extern long double logl (long double);
+#include "math.h"
+#include "math_private.h"
 
 /* 1 / sqrt(pi) */
 static const long double ONEOSQPI = 5.6418958354775628694807945156077258584405E-1L;
@@ -634,7 +631,7 @@ static const long double Q2_2r3D[NQ2_2r3D + 1] = {
 /* Evaluate P[n] x^n  +  P[n-1] x^(n-1)  +  ...  +  P[0] */
 
 static long double
-neval (long double x, long double *p, int n)
+neval (long double x, const long double *p, int n)
 {
   long double y;
 
@@ -652,7 +649,7 @@ neval (long double x, long double *p, int n)
 /* Evaluate x^n+1  +  P[n] x^(n)  +  P[n-1] x^(n-1)  +  ...  +  P[0] */
 
 static long double
-deval (long double x, long double *p, int n)
+deval (long double x, const long double *p, int n)
 {
   long double y;
 
diff --git a/sysdeps/ieee754/ldbl-128/e_lgammal_r.c b/sysdeps/ieee754/ldbl-128/e_lgammal_r.c
index 9537f12b5c..b13ca91c3f 100644
--- a/sysdeps/ieee754/ldbl-128/e_lgammal_r.c
+++ b/sysdeps/ieee754/ldbl-128/e_lgammal_r.c
@@ -707,7 +707,7 @@ static const long double RDr9[NRDr9 + 1] =
 /* Evaluate P[n] x^n  +  P[n-1] x^(n-1)  +  ...  +  P[0] */
 
 static long double
-neval (long double x, long double *p, int n)
+neval (long double x, const long double *p, int n)
 {
   long double y;
 
@@ -725,7 +725,7 @@ neval (long double x, long double *p, int n)
 /* Evaluate x^n+1  +  P[n] x^(n)  +  P[n-1] x^(n-1)  +  ...  +  P[0] */
 
 static long double
-deval (long double x, long double *p, int n)
+deval (long double x, const long double *p, int n)
 {
   long double y;
 
diff --git a/sysdeps/ieee754/ldbl-128/s_erfl.c b/sysdeps/ieee754/ldbl-128/s_erfl.c
index efa0e8fcb3..ef8eea6875 100644
--- a/sysdeps/ieee754/ldbl-128/s_erfl.c
+++ b/sysdeps/ieee754/ldbl-128/s_erfl.c
@@ -84,7 +84,7 @@
 /* Evaluate P[n] x^n  +  P[n-1] x^(n-1)  +  ...  +  P[0] */
 
 static long double
-neval (long double x, long double *p, int n)
+neval (long double x, const long double *p, int n)
 {
   long double y;
 
@@ -102,7 +102,7 @@ neval (long double x, long double *p, int n)
 /* Evaluate x^n+1  +  P[n] x^(n)  +  P[n-1] x^(n-1)  +  ...  +  P[0] */
 
 static long double
-deval (long double x, long double *p, int n)
+deval (long double x, const long double *p, int n)
 {
   long double y;
 
@@ -746,7 +746,7 @@ __erfl (x)
 #endif
 {
   long double a, y, z;
-  int32_t i, ix, sign, flag;
+  int32_t i, ix, sign;
   ieee854_long_double_shape_type u;
 
   u.value = x;
@@ -804,8 +804,8 @@ weak_alias (__erf, erfl)
        x;
 #endif
 {
-  long double P, Q, s, y, z, p, r;
-  int32_t i, ix, sign, flag;
+  long double y, z, p, r;
+  int32_t i, ix, sign;
   ieee854_long_double_shape_type u;
 
   u.value = x;