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.
This commit is contained in:
Andreas Jaeger 2001-11-26 08:05:53 +00:00
parent 9992bc08aa
commit 60a06b7ce1
6 changed files with 33 additions and 18 deletions

View File

@ -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

View File

@ -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);

View File

@ -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;

View File

@ -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;

View File

@ -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;

View File

@ -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;