diff options
author | Andreas Jaeger <aj@suse.de> | 2001-11-26 08:05:53 +0000 |
---|---|---|
committer | Andreas Jaeger <aj@suse.de> | 2001-11-26 08:05:53 +0000 |
commit | 60a06b7ce15c3fd911c59693ee24fda12404f3a7 (patch) | |
tree | e1a7c69715a316fc852702877aadc135e445bd47 | |
parent | 9992bc08aa26c1cf00f6dad46c5be8e56265aa21 (diff) | |
download | glibc-60a06b7ce15c3fd911c59693ee24fda12404f3a7.tar glibc-60a06b7ce15c3fd911c59693ee24fda12404f3a7.tar.gz glibc-60a06b7ce15c3fd911c59693ee24fda12404f3a7.tar.bz2 glibc-60a06b7ce15c3fd911c59693ee24fda12404f3a7.zip |
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.
-rw-r--r-- | ChangeLog | 19 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/e_coshl.c | 3 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/e_j0l.c | 4 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/e_j1l.c | 11 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/e_lgammal_r.c | 4 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/s_erfl.c | 10 |
6 files changed, 33 insertions, 18 deletions
@@ -1,3 +1,22 @@ +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. + 2001-11-22 Andreas Jaeger <aj@suse.de> * sysdeps/ieee754/ldbl-96/s_ilogbl.c (__ilogbl): Add brace to 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; |