From 60a06b7ce15c3fd911c59693ee24fda12404f3a7 Mon Sep 17 00:00:00 2001 From: Andreas Jaeger Date: Mon, 26 Nov 2001 08:05:53 +0000 Subject: Update. 2001-11-26 Stephen L Moshier * 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. --- sysdeps/ieee754/ldbl-128/e_coshl.c | 3 +-- sysdeps/ieee754/ldbl-128/e_j0l.c | 4 ++-- sysdeps/ieee754/ldbl-128/e_j1l.c | 11 ++++------- sysdeps/ieee754/ldbl-128/e_lgammal_r.c | 4 ++-- sysdeps/ieee754/ldbl-128/s_erfl.c | 10 +++++----- 5 files changed, 14 insertions(+), 18 deletions(-) (limited to 'sysdeps/ieee754') 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; -- cgit 1.4.1