qiurui
/
glibc
mirror of git://sourceware.org/git/glibc.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

PowerPC floating point little-endian [2 of 15] 

http://sourceware.org/ml/libc-alpha/2013-08/msg00082.html

This patch replaces occurrences of GET_LDOUBLE_* and SET_LDOUBLE_*
macros, and union ieee854_long_double_shape_type in ldbl-128ibm/,
and a stray one in the 32-bit fpu support.  These files have no
significant changes apart from rewriting the long double bit access.

	* sysdeps/ieee754/ldbl-128ibm/math_ldbl.h (ldbl_high): Define.
	* sysdeps/ieee754/ldbl-128ibm/e_acoshl.c (__ieee754_acoshl): Rewrite
	all uses of ieee854 long double macros and unions.
	* 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_atanhl.c (__ieee754_atanhl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/e_coshl.c (__ieee754_coshl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/e_log2l.c (__ieee754_log2l): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/e_rem_pio2l.c (__ieee754_rem_pio2l):
	Likewise.
	* sysdeps/ieee754/ldbl-128ibm/e_sinhl.c (__ieee754_sinhl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/k_cosl.c (__kernel_cosl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/k_sincosl.c (__kernel_sincosl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/k_sinl.c (__kernel_sinl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_asinhl.c (__asinhl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_atanl.c (__atanl): Likewise.
	Simplify sign and nan test too.
	* sysdeps/ieee754/ldbl-128ibm/s_cosl.c (__cosl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_fabsl.c (__fabsl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_finitel.c (___finitel): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_fpclassifyl.c (___fpclassifyl):
	Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_isnanl.c (___isnanl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_issignalingl.c (__issignalingl):
	Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_logbl.c (__logbl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_signbitl.c (___signbitl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_sincosl.c (__sincosl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_sinl.c (__sinl): Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_tanl.c (__tanl): Likewise.
	* sysdeps/powerpc/powerpc32/power7/fpu/s_logbl.c (__logbl): Likewise.
This commit is contained in:
Alan Modra 2013-08-17 18:24:05 +09:30
parent 1b6adf888d
commit 4ebd120cd9
27 changed files with 176 additions and 89 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

34
ChangeLog
View File

@ -1,3 +1,37 @@
2013-10-04 Alan Modra <amodra@gmail.com>
* sysdeps/ieee754/ldbl-128ibm/math_ldbl.h (ldbl_high): Define.
* sysdeps/ieee754/ldbl-128ibm/e_acoshl.c (__ieee754_acoshl): Rewrite
all uses of ieee854 long double macros and unions.
* 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_atanhl.c (__ieee754_atanhl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/e_coshl.c (__ieee754_coshl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/e_log2l.c (__ieee754_log2l): Likewise.
* sysdeps/ieee754/ldbl-128ibm/e_rem_pio2l.c (__ieee754_rem_pio2l):
Likewise.
* sysdeps/ieee754/ldbl-128ibm/e_sinhl.c (__ieee754_sinhl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/k_cosl.c (__kernel_cosl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/k_sincosl.c (__kernel_sincosl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/k_sinl.c (__kernel_sinl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_asinhl.c (__asinhl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_atanl.c (__atanl): Likewise.
Simplify sign and nan test too.
* sysdeps/ieee754/ldbl-128ibm/s_cosl.c (__cosl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_fabsl.c (__fabsl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_finitel.c (___finitel): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_fpclassifyl.c (___fpclassifyl):
Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_isnanl.c (___isnanl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_issignalingl.c (__issignalingl):
Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_logbl.c (__logbl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_signbitl.c (___signbitl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_sincosl.c (__sincosl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_sinl.c (__sinl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_tanl.c (__tanl): Likewise.
* sysdeps/powerpc/powerpc32/power7/fpu/s_logbl.c (__logbl): Likewise.
2013-10-04 Alan Modra <amodra@gmail.com> 2013-10-04 Alan Modra <amodra@gmail.com>
* stdio-common/printf_size.c (__printf_size): Don't use * stdio-common/printf_size.c (__printf_size): Don't use

8
sysdeps/ieee754/ldbl-128ibm/e_acoshl.c
View File

@ -36,8 +36,12 @@ __ieee754_acoshl(long double x)
{ {
long double t; long double t;
int64_t hx; int64_t hx;
u_int64_t lx; uint64_t lx;
GET_LDOUBLE_WORDS64(hx,lx,x); double xhi, xlo;
ldbl_unpack (x, &xhi, &xlo);
EXTRACT_WORDS64 (hx, xhi);
EXTRACT_WORDS64 (lx, xlo);
if(hx<0x3ff0000000000000LL) { /* x < 1 */ if(hx<0x3ff0000000000000LL) { /* x < 1 */
return (x-x)/(x-x); return (x-x)/(x-x);
} else if(hx >=0x41b0000000000000LL) { /* x > 2**28 */ } else if(hx >=0x41b0000000000000LL) { /* x > 2**28 */

36
sysdeps/ieee754/ldbl-128ibm/e_acosl.c
View File

@ -151,26 +151,25 @@ static const long double
long double long double
__ieee754_acosl (long double x) __ieee754_acosl (long double x)
{ {
long double z, r, w, p, q, s, t, f2; long double a, z, r, w, p, q, s, t, f2;
ieee854_long_double_shape_type u;
u.value = __builtin_fabsl (x); a = __builtin_fabsl (x);
if (u.value == 1.0L) if (a == 1.0L)
{ {
if (x > 0.0L) if (x > 0.0L)
return 0.0; /* acos(1) = 0 */ return 0.0; /* acos(1) = 0 */
else else
return (2.0 * pio2_hi) + (2.0 * pio2_lo); /* acos(-1)= pi */ return (2.0 * pio2_hi) + (2.0 * pio2_lo); /* acos(-1)= pi */
} }
else if (u.value > 1.0L) else if (a > 1.0L)
{ {
return (x - x) / (x - x); /* acos(|x| > 1) is NaN */ return (x - x) / (x - x); /* acos(|x| > 1) is NaN */
} }
if (u.value < 0.5L) if (a < 0.5L)
{ {
if (u.value < 6.938893903907228e-18L) /* |x| < 2**-57 */ if (a < 6.938893903907228e-18L) /* |x| < 2**-57 */
return pio2_hi + pio2_lo; return pio2_hi + pio2_lo;
if (u.value < 0.4375L) if (a < 0.4375L)
{ {
/* Arcsine of x. */ /* Arcsine of x. */
z = x * x; z = x * x;
@ -199,7 +198,7 @@ __ieee754_acosl (long double x)
return z; return z;
} }
/* .4375 <= |x| < .5 */ /* .4375 <= |x| < .5 */
t = u.value - 0.4375L; t = a - 0.4375L;
p = ((((((((((P10 * t p = ((((((((((P10 * t
+ P9) * t + P9) * t
+ P8) * t + P8) * t
@ -230,9 +229,9 @@ __ieee754_acosl (long double x)
r = acosr4375 + r; r = acosr4375 + r;
return r; return r;
} }
else if (u.value < 0.625L) else if (a < 0.625L)
{ {
t = u.value - 0.5625L; t = a - 0.5625L;
p = ((((((((((rS10 * t p = ((((((((((rS10 * t
+ rS9) * t + rS9) * t
+ rS8) * t + rS8) * t
@ -264,7 +263,9 @@ __ieee754_acosl (long double x)
} }
else else
{ /* |x| >= .625 */ { /* |x| >= .625 */
z = (one - u.value) * 0.5; double shi, slo;
z = (one - a) * 0.5;
s = __ieee754_sqrtl (z); s = __ieee754_sqrtl (z);
/* Compute an extended precision square root from /* Compute an extended precision square root from
the Newton iteration s -> 0.5 * (s + z / s). the Newton iteration s -> 0.5 * (s + z / s).
@ -273,12 +274,11 @@ __ieee754_acosl (long double x)
Express s = f1 + f2 where f1 * f1 is exactly representable. Express s = f1 + f2 where f1 * f1 is exactly representable.
w = (z - s^2)/2s = (z - f1^2 - 2 f1 f2 - f2^2)/2s . w = (z - s^2)/2s = (z - f1^2 - 2 f1 f2 - f2^2)/2s .
s + w has extended precision. */ s + w has extended precision. */
u.value = s; ldbl_unpack (s, &shi, &slo);
u.parts32.w2 = 0; a = shi;
u.parts32.w3 = 0; f2 = slo;
f2 = s - u.value; w = z - a * a;
w = z - u.value * u.value; w = w - 2.0 * a * f2;
w = w - 2.0 * u.value * f2;
w = w - f2 * f2; w = w - f2 * f2;
w = w / (2.0 * s); w = w / (2.0 * s);
/* Arcsine of s. */ /* Arcsine of s. */

26
sysdeps/ieee754/ldbl-128ibm/e_asinl.c
View File

@ -131,19 +131,18 @@ static const long double
long double long double
__ieee754_asinl (long double x) __ieee754_asinl (long double x)
{ {
long double t, w, p, q, c, r, s; long double a, t, w, p, q, c, r, s;
int flag; int flag;
ieee854_long_double_shape_type u;
flag = 0; flag = 0;
u.value = __builtin_fabsl (x); a = __builtin_fabsl (x);
if (u.value == 1.0L) /* |x|>= 1 */ if (a == 1.0L) /* |x|>= 1 */
return x * pio2_hi + x * pio2_lo; /* asin(1)=+-pi/2 with inexact */ return x * pio2_hi + x * pio2_lo; /* asin(1)=+-pi/2 with inexact */
else if (u.value >= 1.0L) else if (a >= 1.0L)
return (x - x) / (x - x); /* asin(|x|>1) is NaN */ return (x - x) / (x - x); /* asin(|x|>1) is NaN */
else if (u.value < 0.5L) else if (a < 0.5L)
{ {
if (u.value < 6.938893903907228e-18L) /* |x| < 2**-57 */ if (a < 6.938893903907228e-18L) /* |x| < 2**-57 */
{ {
if (huge + x > one) if (huge + x > one)
return x; /* return x with inexact if x!=0 */ return x; /* return x with inexact if x!=0 */
@ -155,9 +154,9 @@ __ieee754_asinl (long double x)
flag = 1; flag = 1;
} }
} }
else if (u.value < 0.625L) else if (a < 0.625L)
{ {
t = u.value - 0.5625; t = a - 0.5625;
p = ((((((((((rS10 * t p = ((((((((((rS10 * t
+ rS9) * t + rS9) * t
+ rS8) * t + rS8) * t
@ -190,7 +189,7 @@ __ieee754_asinl (long double x)
else else
{ {
/* 1 > |x| >= 0.625 */ /* 1 > |x| >= 0.625 */
w = one - u.value; w = one - a;
t = w * 0.5; t = w * 0.5;
} }
@ -223,17 +222,14 @@ __ieee754_asinl (long double x)
} }
s = __ieee754_sqrtl (t); s = __ieee754_sqrtl (t);
if (u.value > 0.975L) if (a > 0.975L)
{ {
w = p / q; w = p / q;
t = pio2_hi - (2.0 * (s + s * w) - pio2_lo); t = pio2_hi - (2.0 * (s + s * w) - pio2_lo);
} }
else else
{ {
u.value = s; w = ldbl_high (s);
u.parts32.w3 = 0;
u.parts32.w2 = 0;
w = u.value;
c = (t - w * w) / (s + w); c = (t - w * w) / (s + w);
r = p / q; r = p / q;
p = 2.0 * s * r - (pio2_lo - 2.0 * c); p = 2.0 * s * r - (pio2_lo - 2.0 * c);

6
sysdeps/ieee754/ldbl-128ibm/e_atanhl.c
View File

@ -40,8 +40,10 @@ __ieee754_atanhl(long double x)
{ {
long double t; long double t;
int64_t hx,ix; int64_t hx,ix;
u_int64_t lx __attribute__ ((unused)); double xhi;
GET_LDOUBLE_WORDS64(hx,lx,x);
xhi = ldbl_high (x);
EXTRACT_WORDS64 (hx, xhi);
ix = hx&0x7fffffffffffffffLL; ix = hx&0x7fffffffffffffffLL;
if (ix >= 0x3ff0000000000000LL) { /* |x|>=1 */ if (ix >= 0x3ff0000000000000LL) { /* |x|>=1 */
if (ix > 0x3ff0000000000000LL) if (ix > 0x3ff0000000000000LL)

4
sysdeps/ieee754/ldbl-128ibm/e_coshl.c
View File

@ -41,9 +41,11 @@ __ieee754_coshl (long double x)
{ {
long double t,w; long double t,w;
int64_t ix; int64_t ix;
double xhi;
/* High word of |x|. */ /* High word of |x|. */
GET_LDOUBLE_MSW64(ix,x); xhi = ldbl_high (x);
EXTRACT_WORDS64 (ix, xhi);
ix &= 0x7fffffffffffffffLL; ix &= 0x7fffffffffffffffLL;
/* x is INF or NaN */ /* x is INF or NaN */

8
sysdeps/ieee754/ldbl-128ibm/e_log2l.c
View File

@ -177,11 +177,13 @@ __ieee754_log2l (x)
long double z; long double z;
long double y; long double y;
int e; int e;
int64_t hx, lx; int64_t hx;
double xhi;
/* Test for domain */ /* Test for domain */
GET_LDOUBLE_WORDS64 (hx, lx, x); xhi = ldbl_high (x);
if (((hx & 0x7fffffffffffffffLL) | (lx & 0x7fffffffffffffffLL)) == 0) EXTRACT_WORDS64 (hx, xhi);
if ((hx & 0x7fffffffffffffffLL) == 0)
return (-1.0L / (x - x)); return (-1.0L / (x - x));
if (hx < 0) if (hx < 0)
return (x - x) / (x - x); return (x - x) / (x - x);

5
sysdeps/ieee754/ldbl-128ibm/e_rem_pio2l.c
View File

@ -200,10 +200,11 @@ int32_t __ieee754_rem_pio2l(long double x, long double *y)
double tx[8]; double tx[8];
int exp; int exp;
int64_t n, ix, hx, ixd; int64_t n, ix, hx, ixd;
u_int64_t lx __attribute__ ((unused));
u_int64_t lxd; u_int64_t lxd;
double xhi;
GET_LDOUBLE_WORDS64 (hx, lx, x); xhi = ldbl_high (x);
EXTRACT_WORDS64 (hx, xhi);
ix = hx & 0x7fffffffffffffffLL; ix = hx & 0x7fffffffffffffffLL;
if (ix <= 0x3fe921fb54442d10LL) /* x in <-pi/4, pi/4> */ if (ix <= 0x3fe921fb54442d10LL) /* x in <-pi/4, pi/4> */
{ {

4
sysdeps/ieee754/ldbl-128ibm/e_sinhl.c
View File

@ -38,9 +38,11 @@ __ieee754_sinhl(long double x)
{ {
long double t,w,h; long double t,w,h;
int64_t ix,jx; int64_t ix,jx;
double xhi;
/* High word of |x|. */ /* High word of |x|. */
GET_LDOUBLE_MSW64(jx,x); xhi = ldbl_high (x);
EXTRACT_WORDS64 (jx, xhi);
ix = jx&0x7fffffffffffffffLL; ix = jx&0x7fffffffffffffffLL;
/* x is INF or NaN */ /* x is INF or NaN */

10
sysdeps/ieee754/ldbl-128ibm/k_cosl.c
View File

@ -81,8 +81,11 @@ __kernel_cosl(long double x, long double y)
{ {
long double h, l, z, sin_l, cos_l_m1; long double h, l, z, sin_l, cos_l_m1;
int64_t ix; int64_t ix;
u_int32_t tix, hix, index; uint32_t tix, hix, index;
GET_LDOUBLE_MSW64 (ix, x); double xhi, hhi;
xhi = ldbl_high (x);
EXTRACT_WORDS64 (ix, xhi);
tix = ((u_int64_t)ix) >> 32; tix = ((u_int64_t)ix) >> 32;
tix &= ~0x80000000; /* tix = |x|'s high 32 bits */ tix &= ~0x80000000; /* tix = |x|'s high 32 bits */
if (tix < 0x3fc30000) /* |x| < 0.1484375 */ if (tix < 0x3fc30000) /* |x| < 0.1484375 */
@ -136,7 +139,8 @@ __kernel_cosl(long double x, long double y)
case 2: index = (hix - 0x3fc30000) >> 14; break; case 2: index = (hix - 0x3fc30000) >> 14; break;
} }
*/ */
SET_LDOUBLE_WORDS64(h, ((u_int64_t)hix) << 32, 0); INSERT_WORDS64 (hhi, ((uint64_t)hix) << 32);
h = hhi;
l = y - (h - x); l = y - (h - x);
z = l * l; z = l * l;
sin_l = l*(ONE+z*(SSIN1+z*(SSIN2+z*(SSIN3+z*(SSIN4+z*SSIN5))))); sin_l = l*(ONE+z*(SSIN1+z*(SSIN2+z*(SSIN3+z*(SSIN4+z*SSIN5)))));

12
sysdeps/ieee754/ldbl-128ibm/k_sincosl.c
View File

@ -100,9 +100,12 @@ __kernel_sincosl(long double x, long double y, long double *sinx, long double *c
{ {
long double h, l, z, sin_l, cos_l_m1; long double h, l, z, sin_l, cos_l_m1;
int64_t ix; int64_t ix;
u_int32_t tix, hix, index; uint32_t tix, hix, index;
GET_LDOUBLE_MSW64 (ix, x); double xhi, hhi;
tix = ((u_int64_t)ix) >> 32;
xhi = ldbl_high (x);
EXTRACT_WORDS64 (ix, xhi);
tix = ((uint64_t)ix) >> 32;
tix &= ~0x80000000; /* tix = |x|'s high 32 bits */ tix &= ~0x80000000; /* tix = |x|'s high 32 bits */
if (tix < 0x3fc30000) /* |x| < 0.1484375 */ if (tix < 0x3fc30000) /* |x| < 0.1484375 */
{ {
@ -164,7 +167,8 @@ __kernel_sincosl(long double x, long double y, long double *sinx, long double *c
case 2: index = (hix - 0x3fc30000) >> 14; break; case 2: index = (hix - 0x3fc30000) >> 14; break;
} }
*/ */
SET_LDOUBLE_WORDS64(h, ((u_int64_t)hix) << 32, 0); INSERT_WORDS64 (hhi, ((uint64_t)hix) << 32);
h = hhi;
if (iy) if (iy)
l = y - (h - x); l = y - (h - x);
else else

8
sysdeps/ieee754/ldbl-128ibm/k_sinl.c
View File

@ -82,7 +82,10 @@ __kernel_sinl(long double x, long double y, int iy)
long double h, l, z, sin_l, cos_l_m1; long double h, l, z, sin_l, cos_l_m1;
int64_t ix; int64_t ix;
u_int32_t tix, hix, index; u_int32_t tix, hix, index;
GET_LDOUBLE_MSW64 (ix, x); double xhi, hhi;
xhi = ldbl_high (x);
EXTRACT_WORDS64 (ix, xhi);
tix = ((u_int64_t)ix) >> 32; tix = ((u_int64_t)ix) >> 32;
tix &= ~0x80000000; /* tix = |x|'s high 32 bits */ tix &= ~0x80000000; /* tix = |x|'s high 32 bits */
if (tix < 0x3fc30000) /* |x| < 0.1484375 */ if (tix < 0x3fc30000) /* |x| < 0.1484375 */
@ -132,7 +135,8 @@ __kernel_sinl(long double x, long double y, int iy)
case 2: index = (hix - 0x3fc30000) >> 14; break; case 2: index = (hix - 0x3fc30000) >> 14; break;
} }
*/ */
SET_LDOUBLE_WORDS64(h, ((u_int64_t)hix) << 32, 0); INSERT_WORDS64 (hhi, ((uint64_t)hix) << 32);
h = hhi;
if (iy) if (iy)
l = (ix < 0 ? -y : y) - (h - x); l = (ix < 0 ? -y : y) - (h - x);
else else

3
sysdeps/ieee754/ldbl-128ibm/math_ldbl.h
View File

@ -192,6 +192,9 @@ default_ldbl_unpack (long double l, double *a, double *aa)
# define ldbl_unpack default_ldbl_unpack # define ldbl_unpack default_ldbl_unpack
#endif #endif
/* Extract high double. */
#define ldbl_high(x) ((double) x)
/* Convert a finite long double to canonical form. /* Convert a finite long double to canonical form.
Does not handle +/-Inf properly. */ Does not handle +/-Inf properly. */
static inline void static inline void

5
sysdeps/ieee754/ldbl-128ibm/s_asinhl.c
View File

@ -38,7 +38,10 @@ long double __asinhl(long double x)
{ {
long double t,w; long double t,w;
int64_t hx,ix; int64_t hx,ix;
GET_LDOUBLE_MSW64(hx,x); double xhi;
xhi = ldbl_high (x);
EXTRACT_WORDS64 (hx, xhi);
ix = hx&0x7fffffffffffffffLL; ix = hx&0x7fffffffffffffffLL;
if(ix>=0x7ff0000000000000LL) return x+x; /* x is inf or NaN */ if(ix>=0x7ff0000000000000LL) return x+x; /* x is inf or NaN */
if(ix< 0x3e20000000000000LL) { /* |x|<2**-29 */ if(ix< 0x3e20000000000000LL) { /* |x|<2**-29 */

15
sysdeps/ieee754/ldbl-128ibm/s_atanl.c
View File

@ -173,23 +173,20 @@ static const long double
long double long double
__atanl (long double x) __atanl (long double x)
{ {
int k, sign; int32_t k, sign, lx;
long double t, u, p, q; long double t, u, p, q;
ieee854_long_double_shape_type s; double xhi;
s.value = x; xhi = ldbl_high (x);
k = s.parts32.w0; EXTRACT_WORDS (k, lx, xhi);
if (k & 0x80000000) sign = k & 0x80000000;
sign = 1;
else
sign = 0;
/* Check for IEEE special cases. */ /* Check for IEEE special cases. */
k &= 0x7fffffff; k &= 0x7fffffff;
if (k >= 0x7ff00000) if (k >= 0x7ff00000)
{ {
/* NaN. */ /* NaN. */
if ((k & 0xfffff) | s.parts32.w1 ) if (((k - 0x7ff00000) | lx) != 0)
return (x + x); return (x + x);
/* Infinity. */ /* Infinity. */

4
sysdeps/ieee754/ldbl-128ibm/s_cosl.c
View File

@ -53,9 +53,11 @@ long double __cosl(long double x)
{ {
long double y[2],z=0.0L; long double y[2],z=0.0L;
int64_t n, ix; int64_t n, ix;
double xhi;
/* High word of x. */ /* High word of x. */
GET_LDOUBLE_MSW64(ix,x); xhi = ldbl_high (x);
EXTRACT_WORDS64 (ix, xhi);
/* |x| ~< pi/4 */ /* |x| ~< pi/4 */
ix &= 0x7fffffffffffffffLL; ix &= 0x7fffffffffffffffLL;

10
sysdeps/ieee754/ldbl-128ibm/s_fabsl.c
View File

@ -29,10 +29,16 @@ static char rcsid[] = "$NetBSD: $";
long double __fabsl(long double x) long double __fabsl(long double x)
{ {
u_int64_t hx, lx; u_int64_t hx, lx;
GET_LDOUBLE_WORDS64(hx,lx,x); double xhi, xlo;
ldbl_unpack (x, &xhi, &xlo);
EXTRACT_WORDS64 (hx, xhi);
EXTRACT_WORDS64 (lx, xlo);
lx = lx ^ ( hx & 0x8000000000000000LL ); lx = lx ^ ( hx & 0x8000000000000000LL );
hx = hx & 0x7fffffffffffffffLL; hx = hx & 0x7fffffffffffffffLL;
SET_LDOUBLE_WORDS64(x,hx,lx); INSERT_WORDS64 (xhi, hx);
INSERT_WORDS64 (xlo, lx);
x = ldbl_pack (xhi, xlo);
return x; return x;
} }
long_double_symbol (libm, __fabsl, fabsl); long_double_symbol (libm, __fabsl, fabsl);

12
sysdeps/ieee754/ldbl-128ibm/s_finitel.c
View File

@ -29,10 +29,14 @@ static char rcsid[] = "$NetBSD: $";
int int
___finitel (long double x) ___finitel (long double x)
{ {
int64_t hx; uint64_t hx;
GET_LDOUBLE_MSW64(hx,x); double xhi;
return (int)((u_int64_t)((hx&0x7fffffffffffffffLL)
-0x7ff0000000000000LL)>>63); xhi = ldbl_high (x);
EXTRACT_WORDS64 (hx, xhi);
hx &= 0x7fffffffffffffffLL;
hx -= 0x7ff0000000000000LL;
return hx >> 63;
} }
hidden_ver (___finitel, __finitel) hidden_ver (___finitel, __finitel)
weak_alias (___finitel, ____finitel) weak_alias (___finitel, ____finitel)

5
sysdeps/ieee754/ldbl-128ibm/s_fpclassifyl.c
View File

@ -46,8 +46,10 @@ ___fpclassifyl (long double x)
{ {
u_int64_t hx, lx; u_int64_t hx, lx;
int retval = FP_NORMAL; int retval = FP_NORMAL;
double xhi, xlo;
GET_LDOUBLE_WORDS64 (hx, lx, x); ldbl_unpack (x, &xhi, &xlo);
EXTRACT_WORDS64 (hx, xhi);
if ((hx & 0x7ff0000000000000ULL) == 0x7ff0000000000000ULL) { if ((hx & 0x7ff0000000000000ULL) == 0x7ff0000000000000ULL) {
/* +/-NaN or +/-Inf */ /* +/-NaN or +/-Inf */
if (hx & 0x000fffffffffffffULL) { if (hx & 0x000fffffffffffffULL) {
@ -65,6 +67,7 @@ ___fpclassifyl (long double x)
retval = FP_NORMAL; retval = FP_NORMAL;
} else { } else {
if ((hx & 0x7ff0000000000000ULL) == 0x0360000000000000ULL) { if ((hx & 0x7ff0000000000000ULL) == 0x0360000000000000ULL) {
EXTRACT_WORDS64 (lx, xlo);
if ((lx & 0x7fffffffffffffff) /* lower is non-zero */ if ((lx & 0x7fffffffffffffff) /* lower is non-zero */
&& ((lx^hx) & 0x8000000000000000ULL)) { /* and sign differs */ && ((lx^hx) & 0x8000000000000000ULL)) { /* and sign differs */
/* +/- denormal */ /* +/- denormal */

14
sysdeps/ieee754/ldbl-128ibm/s_isnanl.c
View File

@ -29,12 +29,14 @@ static char rcsid[] = "$NetBSD: $";
int int
___isnanl (long double x) ___isnanl (long double x)
{ {
int64_t hx; uint64_t hx;
int64_t lx __attribute__ ((unused)); double xhi;
GET_LDOUBLE_WORDS64(hx,lx,x);
hx &= 0x7fffffffffffffffLL; xhi = ldbl_high (x);
hx = 0x7ff0000000000000LL - hx; EXTRACT_WORDS64 (hx, xhi);
return (int)((u_int64_t)hx>>63); hx &= 0x7fffffffffffffffLL;
hx = 0x7ff0000000000000LL - hx;
return (int) (hx >> 63);
} }
hidden_ver (___isnanl, __isnanl) hidden_ver (___isnanl, __isnanl)
#ifndef IS_IN_libm #ifndef IS_IN_libm

7
sysdeps/ieee754/ldbl-128ibm/s_issignalingl.c
View File

@ -22,10 +22,13 @@
int int
__issignalingl (long double x) __issignalingl (long double x)
{ {
u_int64_t xi; uint64_t xi;
/* For inspecting NaN status, we only have to look at the first of the pair /* For inspecting NaN status, we only have to look at the first of the pair
of IEEE 754 64-bit precision numbers. */ of IEEE 754 64-bit precision numbers. */
GET_LDOUBLE_MSW64 (xi, x); double xhi;
xhi = ldbl_high (x);
EXTRACT_WORDS64 (xi, xhi);
#ifdef HIGH_ORDER_BIT_IS_SET_FOR_SNAN #ifdef HIGH_ORDER_BIT_IS_SET_FOR_SNAN
# error untested # error untested
/* We only have to care about the high-order bit of x's significand, because /* We only have to care about the high-order bit of x's significand, because

5
sysdeps/ieee754/ldbl-128ibm/s_logbl.c
View File

@ -27,9 +27,10 @@ long double
__logbl (long double x) __logbl (long double x)
{ {
int64_t hx, rhx; int64_t hx, rhx;
int64_t lx __attribute__ ((unused)); double xhi;
GET_LDOUBLE_WORDS64 (hx, lx, x); xhi = ldbl_high (x);
EXTRACT_WORDS64 (hx, xhi);
hx &= 0x7fffffffffffffffLL; /* high |x| */ hx &= 0x7fffffffffffffffLL; /* high |x| */
if (hx == 0) if (hx == 0)
return -1.0 / fabs (x); return -1.0 / fabs (x);

4
sysdeps/ieee754/ldbl-128ibm/s_signbitl.c
View File

@ -25,8 +25,10 @@ int
___signbitl (long double x) ___signbitl (long double x)
{ {
int64_t e; int64_t e;
double xhi;
GET_LDOUBLE_MSW64 (e, x); xhi = ldbl_high (x);
EXTRACT_WORDS64 (e, xhi);
return e < 0; return e < 0;
} }
#ifdef IS_IN_libm #ifdef IS_IN_libm

4
sysdeps/ieee754/ldbl-128ibm/s_sincosl.c
View File

@ -27,9 +27,11 @@ void
__sincosl (long double x, long double *sinx, long double *cosx) __sincosl (long double x, long double *sinx, long double *cosx)
{ {
int64_t ix; int64_t ix;
double xhi;
/* High word of x. */ /* High word of x. */
GET_LDOUBLE_MSW64 (ix, x); xhi = ldbl_high (x);
EXTRACT_WORDS64 (ix, xhi);
/* |x| ~< pi/4 */ /* |x| ~< pi/4 */
ix &= 0x7fffffffffffffffLL; ix &= 0x7fffffffffffffffLL;

4
sysdeps/ieee754/ldbl-128ibm/s_sinl.c
View File

@ -53,9 +53,11 @@ long double __sinl(long double x)
{ {
long double y[2],z=0.0L; long double y[2],z=0.0L;
int64_t n, ix; int64_t n, ix;
double xhi;
/* High word of x. */ /* High word of x. */
GET_LDOUBLE_MSW64(ix,x); xhi = ldbl_high (x);
EXTRACT_WORDS64 (ix, xhi);
/* |x| ~< pi/4 */ /* |x| ~< pi/4 */
ix &= 0x7fffffffffffffffLL; ix &= 0x7fffffffffffffffLL;

4
sysdeps/ieee754/ldbl-128ibm/s_tanl.c
View File

@ -53,9 +53,11 @@ long double __tanl(long double x)
{ {
long double y[2],z=0.0L; long double y[2],z=0.0L;
int64_t n, ix; int64_t n, ix;
double xhi;
/* High word of x. */ /* High word of x. */
GET_LDOUBLE_MSW64(ix,x); xhi = ldbl_high (x);
EXTRACT_WORDS64 (ix, xhi);
/* |x| ~< pi/4 */ /* |x| ~< pi/4 */
ix &= 0x7fffffffffffffffLL; ix &= 0x7fffffffffffffffLL;

8
sysdeps/powerpc/powerpc32/power7/fpu/s_logbl.c
View File

@ -35,14 +35,14 @@ static const union {
long double long double
__logbl (long double x) __logbl (long double x)
{ {
double xh, xl; double xh;
double ret; double ret;
if (__builtin_expect (x == 0.0L, 0)) if (__builtin_expect (x == 0.0L, 0))
/* Raise FE_DIVBYZERO and return -HUGE_VAL[LF]. */ /* Raise FE_DIVBYZERO and return -HUGE_VAL[LF]. */
return -1.0L / __builtin_fabsl (x); return -1.0L / __builtin_fabsl (x);
ldbl_unpack (x, &xh, &xl); xh = ldbl_high (x);
/* ret = x & 0x7ff0000000000000; */ /* ret = x & 0x7ff0000000000000; */
asm ( asm (
"xxland %x0,%x1,%x2\n" "xxland %x0,%x1,%x2\n"
@ -58,9 +58,9 @@ __logbl (long double x)
{ {
/* POSIX specifies that denormal number is treated as /* POSIX specifies that denormal number is treated as
though it were normalized. */ though it were normalized. */
int64_t lx, hx; int64_t hx;
GET_LDOUBLE_WORDS64 (hx, lx, x); EXTRACT_WORDS64 (hx, xh);
return (long double) (-1023 - (__builtin_clzll (hx) - 12)); return (long double) (-1023 - (__builtin_clzll (hx) - 12));
} }
/* Test to avoid logb_downward (0.0) == -0.0. */ /* Test to avoid logb_downward (0.0) == -0.0. */