AuroraMiddleware/glibc
1
0
Fork 0
mirror of https://sourceware.org/git/glibc.git synced 2024-11-21 20:40:05 +00:00
Code Issues Packages Projects Releases Wiki Activity

ldbl-128: Cleanup e_gammal_r.c after _Float128 rename

Browse Source
This commit is contained in:
Paul E. Murphy 2016-08-31 13:17:35 -05:00
parent 15089e046b
commit de6b6d14e9
2 changed files with 12 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
ChangeLog
View File

@ -1,3 +1,9 @@
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:

12
sysdeps/ieee754/ldbl-128/e_gammal_r.c
View File

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