aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorPaul E. Murphy <murphyp@linux.vnet.ibm.com>2016-08-31 13:17:35 -0500
committerPaul E. Murphy <murphyp@linux.vnet.ibm.com>2016-08-31 17:17:03 -0500
commitde6b6d14e9b4d8e9a7333addca55163d9ebb7ff6 (patch)
treeb5af69384978b9de2e942e3192c2bd503be656f1
parent15089e046b6c71bbefe29687fe4c7e569c9e1c03 (diff)
downloadglibc-de6b6d14e9b4d8e9a7333addca55163d9ebb7ff6.zip
glibc-de6b6d14e9b4d8e9a7333addca55163d9ebb7ff6.tar.gz
glibc-de6b6d14e9b4d8e9a7333addca55163d9ebb7ff6.tar.bz2
ldbl-128: Cleanup e_gammal_r.c after _Float128 rename
-rw-r--r--ChangeLog6
-rw-r--r--sysdeps/ieee754/ldbl-128/e_gammal_r.c12
2 files changed, 12 insertions, 6 deletions
diff --git a/ChangeLog b/ChangeLog
index 8d9483b..4fd9342 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -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));