diff options
| -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)); |