diff options
author | Wilco Dijkstra <wdijkstr@arm.com> | 2018-03-15 18:05:03 +0000 |
---|---|---|
committer | Wilco Dijkstra <wdijkstr@arm.com> | 2018-03-15 19:21:36 +0000 |
commit | f67a8147b0ac274655c6c5c3b4c8bcb02be5745b (patch) | |
tree | 1fe21849a6ae923ddfd2c81e8e70ecd06fd96806 /sysdeps/ieee754/ldbl-128 | |
parent | 1294b1892e19d70e9e4dca0a2f3e39497f262a42 (diff) | |
download | glibc-f67a8147b0ac274655c6c5c3b4c8bcb02be5745b.tar.gz glibc-f67a8147b0ac274655c6c5c3b4c8bcb02be5745b.tar.xz glibc-f67a8147b0ac274655c6c5c3b4c8bcb02be5745b.zip |
Rename all __ieee754_sqrt(f/l) calls to sqrt(f/l)
Use sqrt(f/l) to enable inlining by GCC - if inlining doesn't happen, the asm redirect ensures we will still call __ieee754_sqrt(f/l). * sysdeps/ieee754/dbl-64/e_acosh.c (__ieee754_acosh): Use sqrt. * sysdeps/ieee754/dbl-64/e_gamma_r.c (gamma_positive): Likewise. * sysdeps/ieee754/dbl-64/e_hypot.c (__ieee754_hypot): Likewise. * sysdeps/ieee754/dbl-64/e_j0.c (__ieee754_j0): Likewise. * sysdeps/ieee754/dbl-64/e_j1.c (__ieee754_j1): Likewise. * sysdeps/ieee754/dbl-64/e_jn.c (__ieee754_jn): Likewise. * sysdeps/ieee754/dbl-64/s_asinh.c (__asinh): Likewise. * sysdeps/ieee754/dbl-64/wordsize-64/e_acosh.c (__ieee754_acosh): Likewise. * sysdeps/ieee754/flt-32/e_acosf.c (__ieee754_acosf): Likewise. * sysdeps/ieee754/flt-32/e_acoshf.c (__ieee754_acoshf): Likewise. * sysdeps/ieee754/flt-32/e_asinf.c (__ieee754_asinf): Likewise. * sysdeps/ieee754/flt-32/e_gammaf_r.c (gammaf_positive): Likewise. * sysdeps/ieee754/flt-32/e_hypotf.c (__ieee754_hypotf): Likewise. * sysdeps/ieee754/flt-32/e_j0f.c (__ieee754_j0f): Likewise. * sysdeps/ieee754/flt-32/e_j1f.c (__ieee754_j1f): Likewise. * sysdeps/ieee754/flt-32/e_powf.c (__ieee754_powf): Likewise. * sysdeps/ieee754/flt-32/s_asinhf.c (__asinhf): Likewise. * sysdeps/ieee754/ldbl-128/e_acoshl.c (__ieee754_acoshl): Use sqrtl. * sysdeps/ieee754/ldbl-128/e_acosl.c (__ieee754_acosl): Likewise. * sysdeps/ieee754/ldbl-128/e_asinl.c (__ieee754_asinl): Likewise. * sysdeps/ieee754/ldbl-128/e_gammal_r.c (gammal_positive): Likewise. * sysdeps/ieee754/ldbl-128/e_hypotl.c (__ieee754_hypotl): Likewise. * sysdeps/ieee754/ldbl-128/e_j0l.c (__ieee754_j0l): Likewise. * sysdeps/ieee754/ldbl-128/e_j1l.c (__ieee754_j1l): Likewise. * sysdeps/ieee754/ldbl-128/e_jnl.c (__ieee754_jnl): Likewise. * sysdeps/ieee754/ldbl-128/e_powl.c (__ieee754_powl): Likewise. * sysdeps/ieee754/ldbl-128/s_asinhl.c (__ieee754_asinhl): Likewise. * sysdeps/ieee754/ldbl-128ibm/e_acoshl.c (__ieee754_acoshl): Likewise. * sysdeps/ieee754/ldbl-128ibm/e_acosl.c (__ieee754_acosl): Likewise. * sysdeps/ieee754/ldbl-128ibm/e_asinl.c (__ieee754_asinl): Likewise. * sysdeps/ieee754/ldbl-128ibm/e_gammal_r.c (gammal_positive): Likewise. * sysdeps/ieee754/ldbl-128ibm/e_hypotl.c (__ieee754_hypotl): Likewise. * sysdeps/ieee754/ldbl-128ibm/e_j0l.c (__ieee754_j0l): Likewise. * sysdeps/ieee754/ldbl-128ibm/e_j1l.c (__ieee754_j1l): Likewise * sysdeps/ieee754/ldbl-128ibm/e_jnl.c (__ieee754_jnl): Likewise. * sysdeps/ieee754/ldbl-128ibm/e_powl.c (__ieee754_powl): Likewise. * sysdeps/ieee754/ldbl-128ibm/s_asinhl.c (__ieee754_asinhl): Likewise. * sysdeps/ieee754/ldbl-96/e_acoshl.c (__ieee754_acoshl): Use sqrtl. * sysdeps/ieee754/ldbl-96/e_asinl.c (__ieee754_asinl): Likewise. * sysdeps/ieee754/ldbl-96/e_gammal_r.c (gammal_positive): Likewise. * sysdeps/ieee754/ldbl-96/e_hypotl.c (__ieee754_hypotl): Likewise. * sysdeps/ieee754/ldbl-96/e_j0l.c (__ieee754_j0l): Likewise. * sysdeps/ieee754/ldbl-96/e_j1l.c (__ieee754_j1l): Likewise. * sysdeps/ieee754/ldbl-96/e_jnl.c (__ieee754_jnl): Likewise. * sysdeps/ieee754/ldbl-96/s_asinhl.c (__ieee754_asinhl): Likewise. * sysdeps/m68k/m680x0/fpu/e_pow.c (__ieee754_pow): Likewise. * sysdeps/powerpc/fpu/e_hypot.c (__ieee754_hypot): Likewise. * sysdeps/powerpc/fpu/e_hypotf.c (__ieee754_hypotf): Likewise.
Diffstat (limited to 'sysdeps/ieee754/ldbl-128')
-rw-r--r-- | sysdeps/ieee754/ldbl-128/e_acoshl.c | 4 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/e_acosl.c | 4 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/e_asinl.c | 2 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/e_gammal_r.c | 2 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/e_hypotl.c | 4 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/e_j0l.c | 8 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/e_j1l.c | 8 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/e_jnl.c | 4 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/e_powl.c | 2 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/s_asinhl.c | 4 |
10 files changed, 21 insertions, 21 deletions
diff --git a/sysdeps/ieee754/ldbl-128/e_acoshl.c b/sysdeps/ieee754/ldbl-128/e_acoshl.c index 10576c9a46..c292759ac6 100644 --- a/sysdeps/ieee754/ldbl-128/e_acoshl.c +++ b/sysdeps/ieee754/ldbl-128/e_acoshl.c @@ -52,10 +52,10 @@ __ieee754_acoshl(_Float128 x) return 0; /* acosh(1) = 0 */ } else if (hx > 0x4000000000000000LL) { /* 2**28 > x > 2 */ t=x*x; - return __ieee754_logl(2*x-one/(x+__ieee754_sqrtl(t-one))); + return __ieee754_logl(2*x-one/(x+sqrtl(t-one))); } else { /* 1<x<2 */ t = x-one; - return __log1pl(t+__sqrtl(2*t+t*t)); + return __log1pl(t+sqrtl(2*t+t*t)); } } strong_alias (__ieee754_acoshl, __acoshl_finite) diff --git a/sysdeps/ieee754/ldbl-128/e_acosl.c b/sysdeps/ieee754/ldbl-128/e_acosl.c index 342ea5f47d..0dc23dee2b 100644 --- a/sysdeps/ieee754/ldbl-128/e_acosl.c +++ b/sysdeps/ieee754/ldbl-128/e_acosl.c @@ -51,7 +51,7 @@ * if x is NaN, return x itself; * if |x|>1, return NaN with invalid signal. * - * Functions needed: __ieee754_sqrtl. + * Functions needed: sqrtl. */ #include <math.h> @@ -270,7 +270,7 @@ __ieee754_acosl (_Float128 x) else { /* |x| >= .625 */ z = (one - u.value) * 0.5; - s = __ieee754_sqrtl (z); + s = sqrtl (z); /* Compute an extended precision square root from the Newton iteration s -> 0.5 * (s + z / s). The change w from s to the improved value is diff --git a/sysdeps/ieee754/ldbl-128/e_asinl.c b/sysdeps/ieee754/ldbl-128/e_asinl.c index 1edf1c05a1..53c2f1fada 100644 --- a/sysdeps/ieee754/ldbl-128/e_asinl.c +++ b/sysdeps/ieee754/ldbl-128/e_asinl.c @@ -231,7 +231,7 @@ __ieee754_asinl (_Float128 x) return x + x * w; } - s = __ieee754_sqrtl (t); + s = sqrtl (t); if (ix >= 0x3ffef333) /* |x| > 0.975 */ { w = p / q; diff --git a/sysdeps/ieee754/ldbl-128/e_gammal_r.c b/sysdeps/ieee754/ldbl-128/e_gammal_r.c index 1d415d4400..95b9fe0677 100644 --- a/sysdeps/ieee754/ldbl-128/e_gammal_r.c +++ b/sysdeps/ieee754/ldbl-128/e_gammal_r.c @@ -107,7 +107,7 @@ gammal_positive (_Float128 x, int *exp2_adj) _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) + * sqrtl (2 * M_PIl / x_adj) / prod); exp_adj += x_eps * __ieee754_logl (x_adj); _Float128 bsum = gamma_coeff[NCOEFF - 1]; diff --git a/sysdeps/ieee754/ldbl-128/e_hypotl.c b/sysdeps/ieee754/ldbl-128/e_hypotl.c index 87160e4230..dd82b3aec8 100644 --- a/sysdeps/ieee754/ldbl-128/e_hypotl.c +++ b/sysdeps/ieee754/ldbl-128/e_hypotl.c @@ -116,7 +116,7 @@ __ieee754_hypotl(_Float128 x, _Float128 y) t1 = 0; SET_LDOUBLE_MSW64(t1,ha); t2 = a-t1; - w = __ieee754_sqrtl(t1*t1-(b*(-b)-t2*(a+t1))); + w = sqrtl(t1*t1-(b*(-b)-t2*(a+t1))); } else { a = a+a; y1 = 0; @@ -125,7 +125,7 @@ __ieee754_hypotl(_Float128 x, _Float128 y) t1 = 0; SET_LDOUBLE_MSW64(t1,ha+0x0001000000000000LL); t2 = a - t1; - w = __ieee754_sqrtl(t1*y1-(w*(-w)-(t1*y2+t2*b))); + w = sqrtl(t1*y1-(w*(-w)-(t1*y2+t2*b))); } if(k!=0) { uint64_t high; diff --git a/sysdeps/ieee754/ldbl-128/e_j0l.c b/sysdeps/ieee754/ldbl-128/e_j0l.c index 9190e27c2a..5c699956b6 100644 --- a/sysdeps/ieee754/ldbl-128/e_j0l.c +++ b/sysdeps/ieee754/ldbl-128/e_j0l.c @@ -723,7 +723,7 @@ __ieee754_j0l (_Float128 x) } if (xx > L(0x1p256)) - return ONEOSQPI * cc / __ieee754_sqrtl (xx); + return ONEOSQPI * cc / sqrtl (xx); xinv = 1 / xx; z = xinv * xinv; @@ -786,7 +786,7 @@ __ieee754_j0l (_Float128 x) p = 1 + z * p; q = z * xinv * q; q = q - L(0.125) * xinv; - z = ONEOSQPI * (p * cc - q * ss) / __ieee754_sqrtl (xx); + z = ONEOSQPI * (p * cc - q * ss) / sqrtl (xx); return z; } strong_alias (__ieee754_j0l, __j0l_finite) @@ -868,7 +868,7 @@ _Float128 } if (xx > L(0x1p256)) - return ONEOSQPI * ss / __ieee754_sqrtl (x); + return ONEOSQPI * ss / sqrtl (x); xinv = 1 / xx; z = xinv * xinv; @@ -931,7 +931,7 @@ _Float128 p = 1 + z * p; q = z * xinv * q; q = q - L(0.125) * xinv; - z = ONEOSQPI * (p * ss + q * cc) / __ieee754_sqrtl (x); + z = ONEOSQPI * (p * ss + q * cc) / sqrtl (x); return z; } strong_alias (__ieee754_y0l, __y0l_finite) diff --git a/sysdeps/ieee754/ldbl-128/e_j1l.c b/sysdeps/ieee754/ldbl-128/e_j1l.c index e7af40ef2a..bff3a5e88c 100644 --- a/sysdeps/ieee754/ldbl-128/e_j1l.c +++ b/sysdeps/ieee754/ldbl-128/e_j1l.c @@ -736,7 +736,7 @@ __ieee754_j1l (_Float128 x) if (xx > L(0x1p256)) { - z = ONEOSQPI * cc / __ieee754_sqrtl (xx); + z = ONEOSQPI * cc / sqrtl (xx); if (x < 0) z = -z; return z; @@ -803,7 +803,7 @@ __ieee754_j1l (_Float128 x) p = 1 + z * p; q = z * q; q = q * xinv + L(0.375) * xinv; - z = ONEOSQPI * (p * cc - q * ss) / __ieee754_sqrtl (xx); + z = ONEOSQPI * (p * cc - q * ss) / sqrtl (xx); if (x < 0) z = -z; return z; @@ -892,7 +892,7 @@ __ieee754_y1l (_Float128 x) } if (xx > L(0x1p256)) - return ONEOSQPI * ss / __ieee754_sqrtl (xx); + return ONEOSQPI * ss / sqrtl (xx); xinv = 1 / xx; z = xinv * xinv; @@ -955,7 +955,7 @@ __ieee754_y1l (_Float128 x) p = 1 + z * p; q = z * q; q = q * xinv + L(0.375) * xinv; - z = ONEOSQPI * (p * ss + q * cc) / __ieee754_sqrtl (xx); + z = ONEOSQPI * (p * ss + q * cc) / sqrtl (xx); return z; } strong_alias (__ieee754_y1l, __y1l_finite) diff --git a/sysdeps/ieee754/ldbl-128/e_jnl.c b/sysdeps/ieee754/ldbl-128/e_jnl.c index 7ae671bc3e..635b4e7e42 100644 --- a/sysdeps/ieee754/ldbl-128/e_jnl.c +++ b/sysdeps/ieee754/ldbl-128/e_jnl.c @@ -149,7 +149,7 @@ __ieee754_jnl (int n, _Float128 x) temp = c - s; break; } - b = invsqrtpi * temp / __ieee754_sqrtl (x); + b = invsqrtpi * temp / sqrtl (x); } else { @@ -385,7 +385,7 @@ __ieee754_ynl (int n, _Float128 x) temp = s + c; break; } - b = invsqrtpi * temp / __ieee754_sqrtl (x); + b = invsqrtpi * temp / sqrtl (x); } else { diff --git a/sysdeps/ieee754/ldbl-128/e_powl.c b/sysdeps/ieee754/ldbl-128/e_powl.c index 2aef955b46..4418ad7c38 100644 --- a/sysdeps/ieee754/ldbl-128/e_powl.c +++ b/sysdeps/ieee754/ldbl-128/e_powl.c @@ -233,7 +233,7 @@ __ieee754_powl (_Float128 x, _Float128 y) if (hy == 0x3ffe0000) { /* y is 0.5 */ if (hx >= 0) /* x >= +0 */ - return __ieee754_sqrtl (x); + return sqrtl (x); } } diff --git a/sysdeps/ieee754/ldbl-128/s_asinhl.c b/sysdeps/ieee754/ldbl-128/s_asinhl.c index 307ca0023c..92a482658d 100644 --- a/sysdeps/ieee754/ldbl-128/s_asinhl.c +++ b/sysdeps/ieee754/ldbl-128/s_asinhl.c @@ -65,12 +65,12 @@ __asinhl (_Float128 x) else if (ix >0x40000000) { /* 2^ 54 > |x| > 2.0 */ t = u.value; - w = __ieee754_logl (2.0 * t + one / (__ieee754_sqrtl (x * x + one) + t)); + w = __ieee754_logl (2.0 * t + one / (sqrtl (x * x + one) + t)); } else { /* 2.0 > |x| > 2 ^ -56 */ t = x * x; - w = __log1pl (u.value + t / (one + __ieee754_sqrtl (one + t))); + w = __log1pl (u.value + t / (one + sqrtl (one + t))); } if (sign & 0x80000000) return -w; |