diff options
| author | Paul E. Murphy <murphyp@linux.vnet.ibm.com> | 2016-08-31 13:17:35 -0500 |
|---|---|---|
| committer | Paul E. Murphy <murphyp@linux.vnet.ibm.com> | 2016-08-31 17:17:03 -0500 |
| commit | de6b6d14e9b4d8e9a7333addca55163d9ebb7ff6 (patch) | |
| tree | b5af69384978b9de2e942e3192c2bd503be656f1 | |
| parent | 15089e046b6c71bbefe29687fe4c7e569c9e1c03 (diff) | |
| download | glibc-de6b6d14e9b4d8e9a7333addca55163d9ebb7ff6.zip glibc-de6b6d14e9b4d8e9a7333addca55163d9ebb7ff6.tar.gz glibc-de6b6d14e9b4d8e9a7333addca55163d9ebb7ff6.tar.bz2 | |
ldbl-128: Cleanup e_gammal_r.c after _Float128 rename
| -rw-r--r-- | ChangeLog | 6 | ||||
| -rw-r--r-- | sysdeps/ieee754/ldbl-128/e_gammal_r.c | 12 |
2 files changed, 12 insertions, 6 deletions
@@ -1,5 +1,11 @@ 2016-08-31 Paul E. Murphy <murphyp@linux.vnet.ibm.com> + * sysdeps/ieee754/ldbl-128/e_gammal_r.c: + (gammal_positive): Fix spacing. + (__ieee754_gammal_r): Likewise. + +2016-08-31 Paul E. Murphy <murphyp@linux.vnet.ibm.com> + * sysdeps/ieee754/ldbl-128/math_ldbl.h: (_Float128): Define as long double, always. (L): Apply long double type suffix, always. diff --git a/sysdeps/ieee754/ldbl-128/e_gammal_r.c b/sysdeps/ieee754/ldbl-128/e_gammal_r.c index 43282fe..ac68b31 100644 --- a/sysdeps/ieee754/ldbl-128/e_gammal_r.c +++ b/sysdeps/ieee754/ldbl-128/e_gammal_r.c @@ -105,10 +105,10 @@ gammal_positive (_Float128 x, int *exp2_adj) } *exp2_adj = x_adj_log2 * (int) x_adj_int; _Float128 ret = (__ieee754_powl (x_adj_mant, x_adj) - * __ieee754_exp2l (x_adj_log2 * x_adj_frac) - * __ieee754_expl (-x_adj) - * __ieee754_sqrtl (2 * M_PIl / x_adj) - / prod); + * __ieee754_exp2l (x_adj_log2 * x_adj_frac) + * __ieee754_expl (-x_adj) + * __ieee754_sqrtl (2 * M_PIl / x_adj) + / prod); exp_adj += x_eps * __ieee754_logl (x_adj); _Float128 bsum = gamma_coeff[NCOEFF - 1]; _Float128 x_adj2 = x_adj * x_adj; @@ -188,8 +188,8 @@ __ieee754_gammal_r (_Float128 x, int *signgamp) if (frac > 0.5L) frac = 1.0L - frac; _Float128 sinpix = (frac <= 0.25L - ? __sinl (M_PIl * frac) - : __cosl (M_PIl * (0.5L - frac))); + ? __sinl (M_PIl * frac) + : __cosl (M_PIl * (0.5L - frac))); int exp2_adj; ret = M_PIl / (-x * sinpix * gammal_positive (-x, &exp2_adj)); |