diff options
Diffstat (limited to 'sysdeps/ieee754/ldbl-96/e_j1l.c')
| -rw-r--r-- | sysdeps/ieee754/ldbl-96/e_j1l.c | 19 |
1 files changed, 10 insertions, 9 deletions
diff --git a/sysdeps/ieee754/ldbl-96/e_j1l.c b/sysdeps/ieee754/ldbl-96/e_j1l.c index e8a7349cf4..581615d563 100644 --- a/sysdeps/ieee754/ldbl-96/e_j1l.c +++ b/sysdeps/ieee754/ldbl-96/e_j1l.c @@ -75,6 +75,7 @@ #include <float.h> #include <math.h> #include <math_private.h> +#include <math-underflow.h> static long double pone (long double), qone (long double); @@ -112,7 +113,7 @@ __ieee754_j1l (long double x) { long double z, c, r, s, ss, cc, u, v, y; int32_t ix; - u_int32_t se; + uint32_t se; GET_LDOUBLE_EXP (se, x); ix = se & 0x7fff; @@ -137,12 +138,12 @@ __ieee754_j1l (long double x) * y1(x) = 1/sqrt(pi) * (P(1,x)*ss + Q(1,x)*cc) / sqrt(x) */ if (__glibc_unlikely (ix > 0x4080)) - z = (invsqrtpi * cc) / __ieee754_sqrtl (y); + z = (invsqrtpi * cc) / sqrtl (y); else { u = pone (y); v = qone (y); - z = invsqrtpi * (u * cc - v * ss) / __ieee754_sqrtl (y); + z = invsqrtpi * (u * cc - v * ss) / sqrtl (y); } if (se & 0x8000) return -z; @@ -195,7 +196,7 @@ __ieee754_y1l (long double x) { long double z, s, c, ss, cc, u, v; int32_t ix; - u_int32_t se, i0, i1; + uint32_t se, i0, i1; GET_LDOUBLE_WORDS (se, i0, i1, x); ix = se & 0x7fff; @@ -231,12 +232,12 @@ __ieee754_y1l (long double x) * to compute the worse one. */ if (__glibc_unlikely (ix > 0x4080)) - z = (invsqrtpi * ss) / __ieee754_sqrtl (x); + z = (invsqrtpi * ss) / sqrtl (x); else { u = pone (x); v = qone (x); - z = invsqrtpi * (u * ss + v * cc) / __ieee754_sqrtl (x); + z = invsqrtpi * (u * ss + v * cc) / sqrtl (x); } return z; } @@ -362,7 +363,7 @@ pone (long double x) const long double *p, *q; long double z, r, s; int32_t ix; - u_int32_t se, i0, i1; + uint32_t se, i0, i1; GET_LDOUBLE_WORDS (se, i0, i1, x); ix = se & 0x7fff; @@ -507,9 +508,9 @@ static long double qone (long double x) { const long double *p, *q; - static long double s, r, z; + long double s, r, z; int32_t ix; - u_int32_t se, i0, i1; + uint32_t se, i0, i1; GET_LDOUBLE_WORDS (se, i0, i1, x); ix = se & 0x7fff; |