]> sourceware.org Git - glibc.git/commitdiff
PowerPC floating point little-endian [2 of 15]
authorAlan Modra <amodra@gmail.com>
Sat, 17 Aug 2013 08:54:05 +0000 (18:24 +0930)
committerAlan Modra <amodra@gmail.com>
Fri, 4 Oct 2013 01:02:19 +0000 (10:32 +0930)
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.

27 files changed:
ChangeLog
sysdeps/ieee754/ldbl-128ibm/e_acoshl.c
sysdeps/ieee754/ldbl-128ibm/e_acosl.c
sysdeps/ieee754/ldbl-128ibm/e_asinl.c
sysdeps/ieee754/ldbl-128ibm/e_atanhl.c
sysdeps/ieee754/ldbl-128ibm/e_coshl.c
sysdeps/ieee754/ldbl-128ibm/e_log2l.c
sysdeps/ieee754/ldbl-128ibm/e_rem_pio2l.c
sysdeps/ieee754/ldbl-128ibm/e_sinhl.c
sysdeps/ieee754/ldbl-128ibm/k_cosl.c
sysdeps/ieee754/ldbl-128ibm/k_sincosl.c
sysdeps/ieee754/ldbl-128ibm/k_sinl.c
sysdeps/ieee754/ldbl-128ibm/math_ldbl.h
sysdeps/ieee754/ldbl-128ibm/s_asinhl.c
sysdeps/ieee754/ldbl-128ibm/s_atanl.c
sysdeps/ieee754/ldbl-128ibm/s_cosl.c
sysdeps/ieee754/ldbl-128ibm/s_fabsl.c
sysdeps/ieee754/ldbl-128ibm/s_finitel.c
sysdeps/ieee754/ldbl-128ibm/s_fpclassifyl.c
sysdeps/ieee754/ldbl-128ibm/s_isnanl.c
sysdeps/ieee754/ldbl-128ibm/s_issignalingl.c
sysdeps/ieee754/ldbl-128ibm/s_logbl.c
sysdeps/ieee754/ldbl-128ibm/s_signbitl.c
sysdeps/ieee754/ldbl-128ibm/s_sincosl.c
sysdeps/ieee754/ldbl-128ibm/s_sinl.c
sysdeps/ieee754/ldbl-128ibm/s_tanl.c
sysdeps/powerpc/powerpc32/power7/fpu/s_logbl.c

index c4825619bb3bbb9402d49a7dab7dc05ee501718a..4daa80f53b27cbe7f91e615b7649264b2fd85533 100644 (file)
--- a/ChangeLog
+++ b/ChangeLog
@@ -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>
 
        * stdio-common/printf_size.c (__printf_size): Don't use
index abc78a35bd0ac5c1c0a67548a4f8929eb9f1ea73..8a4a5bb7b9796a3ed57481110f9fb4545052484f 100644 (file)
@@ -36,8 +36,12 @@ __ieee754_acoshl(long double x)
 {
        long double t;
        int64_t hx;
-       u_int64_t lx;
-       GET_LDOUBLE_WORDS64(hx,lx,x);
+       uint64_t lx;
+       double xhi, xlo;
+
+       ldbl_unpack (x, &xhi, &xlo);
+       EXTRACT_WORDS64 (hx, xhi);
+       EXTRACT_WORDS64 (lx, xlo);
        if(hx<0x3ff0000000000000LL) {           /* x < 1 */
            return (x-x)/(x-x);
        } else if(hx >=0x41b0000000000000LL) {  /* x > 2**28 */
index 5d2af30346b446434473c6ad4af97725802ae9ed..86639930c96b0f59a7444264e4b141df96d956ad 100644 (file)
@@ -151,26 +151,25 @@ static const long double
 long double
 __ieee754_acosl (long double x)
 {
-  long double z, r, w, p, q, s, t, f2;
-  ieee854_long_double_shape_type u;
+  long double a, z, r, w, p, q, s, t, f2;
 
-  u.value = __builtin_fabsl (x);
-  if (u.value == 1.0L)
+  a = __builtin_fabsl (x);
+  if (a == 1.0L)
     {
       if (x > 0.0L)
        return 0.0;             /* acos(1) = 0  */
       else
        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 */
     }
-  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;
-      if (u.value < 0.4375L)
+      if (a < 0.4375L)
        {
          /* Arcsine of x.  */
          z = x * x;
@@ -199,7 +198,7 @@ __ieee754_acosl (long double x)
          return z;
        }
       /* .4375 <= |x| < .5 */
-      t = u.value - 0.4375L;
+      t = a - 0.4375L;
       p = ((((((((((P10 * t
                    + P9) * t
                   + P8) * t
@@ -230,9 +229,9 @@ __ieee754_acosl (long double x)
        r = acosr4375 + 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
                    + rS9) * t
                   + rS8) * t
@@ -264,7 +263,9 @@ __ieee754_acosl (long double x)
     }
   else
     {                          /* |x| >= .625 */
-      z = (one - u.value) * 0.5;
+      double shi, slo;
+
+      z = (one - a) * 0.5;
       s = __ieee754_sqrtl (z);
       /* Compute an extended precision square root from
         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.
          w = (z - s^2)/2s = (z - f1^2 - 2 f1 f2 - f2^2)/2s .
          s + w has extended precision.  */
-      u.value = s;
-      u.parts32.w2 = 0;
-      u.parts32.w3 = 0;
-      f2 = s - u.value;
-      w = z - u.value * u.value;
-      w = w - 2.0 * u.value * f2;
+      ldbl_unpack (s, &shi, &slo);
+      a = shi;
+      f2 = slo;
+      w = z - a * a;
+      w = w - 2.0 * a * f2;
       w = w - f2 * f2;
       w = w / (2.0 * s);
       /* Arcsine of s.  */
index b395439495e0c769fb3b55ee284c100ea894b2ae..99a5b85fa024185afca69a544fa8bf3a6291e510 100644 (file)
@@ -131,19 +131,18 @@ static const long double
 long double
 __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;
-  ieee854_long_double_shape_type u;
 
   flag = 0;
-  u.value = __builtin_fabsl (x);
-  if (u.value == 1.0L) /* |x|>= 1 */
+  a = __builtin_fabsl (x);
+  if (a == 1.0L)       /* |x|>= 1 */
     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 */
-  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)
            return x;           /* return x with inexact if x!=0 */
@@ -155,9 +154,9 @@ __ieee754_asinl (long double x)
          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
                    + rS9) * t
                   + rS8) * t
@@ -190,7 +189,7 @@ __ieee754_asinl (long double x)
   else
     {
       /* 1 > |x| >= 0.625 */
-      w = one - u.value;
+      w = one - a;
       t = w * 0.5;
     }
 
@@ -223,17 +222,14 @@ __ieee754_asinl (long double x)
     }
 
   s = __ieee754_sqrtl (t);
-  if (u.value > 0.975L)
+  if (a > 0.975L)
     {
       w = p / q;
       t = pio2_hi - (2.0 * (s + s * w) - pio2_lo);
     }
   else
     {
-      u.value = s;
-      u.parts32.w3 = 0;
-      u.parts32.w2 = 0;
-      w = u.value;
+      w = ldbl_high (s);
       c = (t - w * w) / (s + w);
       r = p / q;
       p = 2.0 * s * r - (pio2_lo - 2.0 * c);
index f35182f03e5819be643bd8d7f219ca6a81e16f9e..29f2e92072128713361e309cc0fd34703cb4b9e6 100644 (file)
@@ -40,8 +40,10 @@ __ieee754_atanhl(long double x)
 {
        long double t;
        int64_t hx,ix;
-       u_int64_t lx __attribute__ ((unused));
-       GET_LDOUBLE_WORDS64(hx,lx,x);
+       double xhi;
+
+       xhi = ldbl_high (x);
+       EXTRACT_WORDS64 (hx, xhi);
        ix = hx&0x7fffffffffffffffLL;
        if (ix >= 0x3ff0000000000000LL) { /* |x|>=1 */
            if (ix > 0x3ff0000000000000LL)
index 3e8e1875c6cd576c500a11744456aaade02e403a..05683bc02feaec95eea92d2129c4b24e499fa60c 100644 (file)
@@ -41,9 +41,11 @@ __ieee754_coshl (long double x)
 {
        long double t,w;
        int64_t ix;
+       double xhi;
 
     /* High word of |x|. */
-       GET_LDOUBLE_MSW64(ix,x);
+       xhi = ldbl_high (x);
+       EXTRACT_WORDS64 (ix, xhi);
        ix &= 0x7fffffffffffffffLL;
 
     /* x is INF or NaN */
index f0098f6c73b67a22b2af46114a8ba55eaac301a0..323ded0c0fa5e2e151fa4837ed9fa05249757eba 100644 (file)
@@ -177,11 +177,13 @@ __ieee754_log2l (x)
   long double z;
   long double y;
   int e;
-  int64_t hx, lx;
+  int64_t hx;
+  double xhi;
 
 /* Test for domain */
-  GET_LDOUBLE_WORDS64 (hx, lx, x);
-  if (((hx & 0x7fffffffffffffffLL) | (lx & 0x7fffffffffffffffLL)) == 0)
+  xhi = ldbl_high (x);
+  EXTRACT_WORDS64 (hx, xhi);
+  if ((hx & 0x7fffffffffffffffLL) == 0)
     return (-1.0L / (x - x));
   if (hx < 0)
     return (x - x) / (x - x);
index 8885def4e150669c74931cb57bb89b8a960bbd44..36bc03226b9c836cbf1878ecd64409b61d89ea44 100644 (file)
@@ -200,10 +200,11 @@ int32_t __ieee754_rem_pio2l(long double x, long double *y)
   double tx[8];
   int exp;
   int64_t n, ix, hx, ixd;
-  u_int64_t lx __attribute__ ((unused));
   u_int64_t lxd;
+  double xhi;
 
-  GET_LDOUBLE_WORDS64 (hx, lx, x);
+  xhi = ldbl_high (x);
+  EXTRACT_WORDS64 (hx, xhi);
   ix = hx & 0x7fffffffffffffffLL;
   if (ix <= 0x3fe921fb54442d10LL)      /* x in <-pi/4, pi/4> */
     {
index 4e8481c4199c4aec261e29468b2ed41b3ab283d3..1790bef87eafe63bfade9228de266312731f086e 100644 (file)
@@ -38,9 +38,11 @@ __ieee754_sinhl(long double x)
 {
        long double t,w,h;
        int64_t ix,jx;
+       double xhi;
 
     /* High word of |x|. */
-       GET_LDOUBLE_MSW64(jx,x);
+       xhi = ldbl_high (x);
+       EXTRACT_WORDS64 (jx, xhi);
        ix = jx&0x7fffffffffffffffLL;
 
     /* x is INF or NaN */
index 0b81782fdb96ac1f6bb5c14254398f8861361e1b..046f3b573c9dc18e5188b23d4d7d344707692035 100644 (file)
@@ -81,8 +81,11 @@ __kernel_cosl(long double x, long double y)
 {
   long double h, l, z, sin_l, cos_l_m1;
   int64_t ix;
-  u_int32_t tix, hix, index;
-  GET_LDOUBLE_MSW64 (ix, x);
+  uint32_t tix, hix, index;
+  double xhi, hhi;
+
+  xhi = ldbl_high (x);
+  EXTRACT_WORDS64 (ix, xhi);
   tix = ((u_int64_t)ix) >> 32;
   tix &= ~0x80000000;                  /* tix = |x|'s high 32 bits */
   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;
        }
 */
-      SET_LDOUBLE_WORDS64(h, ((u_int64_t)hix) << 32, 0);
+      INSERT_WORDS64 (hhi, ((uint64_t)hix) << 32);
+      h = hhi;
       l = y - (h - x);
       z = l * l;
       sin_l = l*(ONE+z*(SSIN1+z*(SSIN2+z*(SSIN3+z*(SSIN4+z*SSIN5)))));
index fc1ead65979825e326ff49d414cef1af7da3f962..3ba9d7e907657a7cd4729332056c549945bbd242 100644 (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;
   int64_t ix;
-  u_int32_t tix, hix, index;
-  GET_LDOUBLE_MSW64 (ix, x);
-  tix = ((u_int64_t)ix) >> 32;
+  uint32_t tix, hix, index;
+  double xhi, hhi;
+
+  xhi = ldbl_high (x);
+  EXTRACT_WORDS64 (ix, xhi);
+  tix = ((uint64_t)ix) >> 32;
   tix &= ~0x80000000;                  /* tix = |x|'s high 32 bits */
   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;
        }
 */
-      SET_LDOUBLE_WORDS64(h, ((u_int64_t)hix) << 32, 0);
+      INSERT_WORDS64 (hhi, ((uint64_t)hix) << 32);
+      h = hhi;
       if (iy)
        l = y - (h - x);
       else
index f17c0ae5db9ce1234ea684b08b085a741fedabdf..b12ea134d58fd17366f1e5bc7c76212fede4da42 100644 (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;
   int64_t ix;
   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 &= ~0x80000000;                  /* tix = |x|'s high 32 bits */
   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;
        }
 */
-      SET_LDOUBLE_WORDS64(h, ((u_int64_t)hix) << 32, 0);
+      INSERT_WORDS64 (hhi, ((uint64_t)hix) << 32);
+      h = hhi;
       if (iy)
        l = (ix < 0 ? -y : y) - (h - x);
       else
index 4bb49c8dd0847aa9b61f61ad9f4c20f58e6474de..8adb081c38d5319909b764189357acaa22f117f1 100644 (file)
@@ -192,6 +192,9 @@ default_ldbl_unpack (long double l, double *a, double *aa)
 # define ldbl_unpack default_ldbl_unpack
 #endif
 
+/* Extract high double.  */
+#define ldbl_high(x) ((double) x)
+
 /* Convert a finite long double to canonical form.
    Does not handle +/-Inf properly.  */
 static inline void
index a833457eab3a114d493e6315cbb5495a71f21af7..63c6edbb1a05d501fdf619d08b7d45cc6184633c 100644 (file)
@@ -38,7 +38,10 @@ long double __asinhl(long double x)
 {
        long double t,w;
        int64_t hx,ix;
-       GET_LDOUBLE_MSW64(hx,x);
+       double xhi;
+
+       xhi = ldbl_high (x);
+       EXTRACT_WORDS64 (hx, xhi);
        ix = hx&0x7fffffffffffffffLL;
        if(ix>=0x7ff0000000000000LL) return x+x;        /* x is inf or NaN */
        if(ix< 0x3e20000000000000LL) {  /* |x|<2**-29 */
index 2a36d16bb410ab0a27e63c65ed6f28b2ff647694..41dde23998db56f421dc5a25bae4828cbecc1c35 100644 (file)
@@ -173,23 +173,20 @@ static const long double
 long double
 __atanl (long double x)
 {
-  int k, sign;
+  int32_t k, sign, lx;
   long double t, u, p, q;
-  ieee854_long_double_shape_type s;
+  double xhi;
 
-  s.value = x;
-  k = s.parts32.w0;
-  if (k & 0x80000000)
-    sign = 1;
-  else
-    sign = 0;
+  xhi = ldbl_high (x);
+  EXTRACT_WORDS (k, lx, xhi);
+  sign = k & 0x80000000;
 
   /* Check for IEEE special cases.  */
   k &= 0x7fffffff;
   if (k >= 0x7ff00000)
     {
       /* NaN. */
-      if ((k & 0xfffff) | s.parts32.w1 )
+      if (((k - 0x7ff00000) | lx) != 0)
        return (x + x);
 
       /* Infinity. */
index 23148392f12f282ec26b2f948210447f58d3793f..54c6cc77d227119cc1f5532b8f61ea344690e2a9 100644 (file)
@@ -53,9 +53,11 @@ long double __cosl(long double x)
 {
        long double y[2],z=0.0L;
        int64_t n, ix;
+       double xhi;
 
     /* High word of x. */
-       GET_LDOUBLE_MSW64(ix,x);
+       xhi = ldbl_high (x);
+       EXTRACT_WORDS64 (ix, xhi);
 
     /* |x| ~< pi/4 */
        ix &= 0x7fffffffffffffffLL;
index 99146d8021470a99e82ee17d887326bebacd15bc..c801c97065251a007e1f90ec840954155a72defb 100644 (file)
@@ -29,10 +29,16 @@ static char rcsid[] = "$NetBSD: $";
 long double __fabsl(long double x)
 {
        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 );
        hx = hx & 0x7fffffffffffffffLL;
-       SET_LDOUBLE_WORDS64(x,hx,lx);
+       INSERT_WORDS64 (xhi, hx);
+       INSERT_WORDS64 (xlo, lx);
+       x = ldbl_pack (xhi, xlo);
        return x;
 }
 long_double_symbol (libm, __fabsl, fabsl);
index 8edb34154d743445765e22c3d57af3f62390773d..7b4655fadb4207baf6c7d3e86cc36aa485a40bed 100644 (file)
@@ -29,10 +29,14 @@ static char rcsid[] = "$NetBSD: $";
 int
 ___finitel (long double x)
 {
-       int64_t hx;
-       GET_LDOUBLE_MSW64(hx,x);
-       return (int)((u_int64_t)((hx&0x7fffffffffffffffLL)
-                                -0x7ff0000000000000LL)>>63);
+  uint64_t hx;
+  double xhi;
+
+  xhi = ldbl_high (x);
+  EXTRACT_WORDS64 (hx, xhi);
+  hx &= 0x7fffffffffffffffLL;
+  hx -= 0x7ff0000000000000LL;
+  return hx >> 63;
 }
 hidden_ver (___finitel, __finitel)
 weak_alias (___finitel, ____finitel)
index f4a90b08c740e9c7467951667fea5dd4406e348c..90586e822e7574159bbcb55520ffafb922a30495 100644 (file)
@@ -46,8 +46,10 @@ ___fpclassifyl (long double x)
 {
   u_int64_t hx, lx;
   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) {
       /* +/-NaN or +/-Inf */
       if (hx & 0x000fffffffffffffULL) {
@@ -65,6 +67,7 @@ ___fpclassifyl (long double x)
              retval = FP_NORMAL;
          } else {
              if ((hx & 0x7ff0000000000000ULL) == 0x0360000000000000ULL) {
+                 EXTRACT_WORDS64 (lx, xlo);
                  if ((lx & 0x7fffffffffffffff) /* lower is non-zero */
                  && ((lx^hx) & 0x8000000000000000ULL)) { /* and sign differs */
                      /* +/- denormal */
index 264dec745ec7bb7a4c266648ae2c8af1134defc8..d12f1d3bf55a58cc593bb9e746d72681ff7246ed 100644 (file)
@@ -29,12 +29,14 @@ static char rcsid[] = "$NetBSD: $";
 int
 ___isnanl (long double x)
 {
-       int64_t hx;
-       int64_t lx __attribute__ ((unused));
-       GET_LDOUBLE_WORDS64(hx,lx,x);
-       hx &= 0x7fffffffffffffffLL;
-       hx = 0x7ff0000000000000LL - hx;
-       return (int)((u_int64_t)hx>>63);
+  uint64_t hx;
+  double xhi;
+
+  xhi = ldbl_high (x);
+  EXTRACT_WORDS64 (hx, xhi);
+  hx &= 0x7fffffffffffffffLL;
+  hx = 0x7ff0000000000000LL - hx;
+  return (int) (hx >> 63);
 }
 hidden_ver (___isnanl, __isnanl)
 #ifndef IS_IN_libm
index 96fab1aff1bdf7da2747da2e57fcd8102d9872c1..bdd58f8f252eb6247331b846a990f2b746fb0c5f 100644 (file)
 int
 __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
      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
 # error untested
   /* We only have to care about the high-order bit of x's significand, because
index 6cbfcfa1cc02a59cb2e4922ab61c9e24fa9b4552..e1402886121b772c94a3e039f53c3a7544e885a4 100644 (file)
@@ -27,9 +27,10 @@ long double
 __logbl (long double x)
 {
   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| */
   if (hx == 0)
     return -1.0 / fabs (x);
index ee4aea6cf45b6df5b945252f00a9e87d827e093a..aecb1fd792ede2075e1932f9390c322949a54cdf 100644 (file)
@@ -25,8 +25,10 @@ int
 ___signbitl (long double x)
 {
   int64_t e;
+  double xhi;
 
-  GET_LDOUBLE_MSW64 (e, x);
+  xhi = ldbl_high (x);
+  EXTRACT_WORDS64 (e, xhi);
   return e < 0;
 }
 #ifdef IS_IN_libm
index 3b1e547bdc8f15c28463e45a099a1e398a15084e..a9e2f3d19a5d2dcb0dd40d321b608e9e91be6703 100644 (file)
@@ -27,9 +27,11 @@ void
 __sincosl (long double x, long double *sinx, long double *cosx)
 {
   int64_t ix;
+  double xhi;
 
   /* High word of x. */
-  GET_LDOUBLE_MSW64 (ix, x);
+  xhi = ldbl_high (x);
+  EXTRACT_WORDS64 (ix, xhi);
 
   /* |x| ~< pi/4 */
   ix &= 0x7fffffffffffffffLL;
index 6fec16f85183a8d54756de8420c35225007fd65d..087921a9132dd31d4c81ee1869fc67fd7000280a 100644 (file)
@@ -53,9 +53,11 @@ long double __sinl(long double x)
 {
        long double y[2],z=0.0L;
        int64_t n, ix;
+       double xhi;
 
     /* High word of x. */
-       GET_LDOUBLE_MSW64(ix,x);
+       xhi = ldbl_high (x);
+       EXTRACT_WORDS64 (ix, xhi);
 
     /* |x| ~< pi/4 */
        ix &= 0x7fffffffffffffffLL;
index 9967d0c2063d8dd8ac62a4613f884dbb1abeedd3..66b8a0621eb2fed0a6028f3d94e41534c1b7a333 100644 (file)
@@ -53,9 +53,11 @@ long double __tanl(long double x)
 {
        long double y[2],z=0.0L;
        int64_t n, ix;
+       double xhi;
 
     /* High word of x. */
-       GET_LDOUBLE_MSW64(ix,x);
+       xhi = ldbl_high (x);
+       EXTRACT_WORDS64 (ix, xhi);
 
     /* |x| ~< pi/4 */
        ix &= 0x7fffffffffffffffLL;
index e008ed0c3daf984f64f78176b115325a1056ee14..1c82577f5d5eddfe8d2dd133f931bf87503ebde4 100644 (file)
@@ -35,14 +35,14 @@ static const union {
 long double
 __logbl (long double x)
 {
-  double xh, xl;
+  double xh;
   double ret;
 
   if (__builtin_expect (x == 0.0L, 0))
     /* Raise FE_DIVBYZERO and return -HUGE_VAL[LF].  */
     return -1.0L / __builtin_fabsl (x);
 
-  ldbl_unpack (x, &xh, &xl);
+  xh = ldbl_high (x);
   /* ret = x & 0x7ff0000000000000;  */
   asm (
     "xxland %x0,%x1,%x2\n"
@@ -58,9 +58,9 @@ __logbl (long double x)
     {
       /* POSIX specifies that denormal number is treated as
          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));
     }
   /* Test to avoid logb_downward (0.0) == -0.0.  */
This page took 0.147718 seconds and 5 git commands to generate.