]> git.ipfire.org Git - thirdparty/glibc.git/commitdiff
PowerPC floating point little-endian [1 of 15]
authorAlan Modra <amodra@gmail.com>
Sat, 17 Aug 2013 08:51:58 +0000 (18:21 +0930)
committerAlan Modra <amodra@gmail.com>
Fri, 4 Oct 2013 01:01:41 +0000 (10:31 +0930)
http://sourceware.org/ml/libc-alpha/2013-08/msg00081.html

This is the first of a series of patches to ban ieee854_long_double
and the ieee854_long_double macros when using IBM long double.  union
ieee854_long_double just isn't correct for IBM long double, especially
when little-endian, and pretending it is OK has allowed a number of
bugs to remain undetected in sysdeps/ieee754/ldbl-128ibm/.

This changes the few places in generic code that use it.

* stdio-common/printf_size.c (__printf_size): Don't use
union ieee854_long_double in fpnum union.
* stdio-common/printf_fphex.c (__printf_fphex): Likewise.  Use
signbit macro to retrieve sign from long double.
* stdio-common/printf_fp.c (___printf_fp): Use signbit macro to
retrieve sign from long double.
* sysdeps/ieee754/ldbl-128ibm/printf_fphex.c: Adjust for fpnum change.
* sysdeps/ieee754/ldbl-128/printf_fphex.c: Likewise.
* sysdeps/ieee754/ldbl-96/printf_fphex.c: Likewise.
* sysdeps/x86_64/fpu/printf_fphex.c: Likewise.
* math/test-misc.c (main): Don't use union ieee854_long_double.
ports/
* sysdeps/ia64/fpu/printf_fphex.c: Adjust for fpnum change.

ChangeLog
math/test-misc.c
ports/ChangeLog
ports/sysdeps/ia64/fpu/printf_fphex.c
stdio-common/printf_fp.c
stdio-common/printf_fphex.c
stdio-common/printf_size.c
sysdeps/ieee754/ldbl-128/printf_fphex.c
sysdeps/ieee754/ldbl-128ibm/printf_fphex.c
sysdeps/ieee754/ldbl-96/printf_fphex.c
sysdeps/x86_64/fpu/printf_fphex.c

index d6bbe3455f1a7f524462cae2455745048ef1390b..c4825619bb3bbb9402d49a7dab7dc05ee501718a 100644 (file)
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,17 @@
+2013-10-04  Alan Modra  <amodra@gmail.com>
+
+       * stdio-common/printf_size.c (__printf_size): Don't use
+       union ieee854_long_double in fpnum union.
+       * stdio-common/printf_fphex.c (__printf_fphex): Likewise.  Use
+       signbit macro to retrieve sign from long double.
+       * stdio-common/printf_fp.c (___printf_fp): Use signbit macro to
+       retrieve sign from long double.
+       * sysdeps/ieee754/ldbl-128ibm/printf_fphex.c: Adjust for fpnum change.
+       * sysdeps/ieee754/ldbl-128/printf_fphex.c: Likewise.
+       * sysdeps/ieee754/ldbl-96/printf_fphex.c: Likewise.
+       * sysdeps/x86_64/fpu/printf_fphex.c: Likewise.
+       * math/test-misc.c (main): Don't use union ieee854_long_double.
+
 2013-10-04  Alan Modra  <amodra@gmail.com>
 
        [BZ #15680]
index 27d673b823c2bc284b2d62b9c710ab925733c77f..f5276ebff95609606e40c0b7b6b110d97388774e 100644 (file)
@@ -722,300 +722,161 @@ main (void)
 
 #ifndef NO_LONG_DOUBLE
   {
-    union ieee854_long_double v1;
-    union ieee854_long_double v2;
-    long double ld;
+    long double v1, v2;
 
-    v1.d = ld = LDBL_MIN;
-    if (fpclassify (ld) != FP_NORMAL)
+    v1 = LDBL_MIN;
+    if (fpclassify (v1) != FP_NORMAL)
       {
-       printf ("fpclassify (LDBL_MIN) failed: %d\n", fpclassify (ld));
+       printf ("fpclassify (LDBL_MIN) failed: %d (%La)\n",
+               fpclassify (v1), v1);
        result = 1;
       }
-    ld = nextafterl (ld, LDBL_MIN / 2.0);
-    if (fpclassify (ld) != FP_SUBNORMAL)
+    v2 = nextafterl (v1, LDBL_MIN / 2.0);
+    if (fpclassify (v2) != FP_SUBNORMAL)
       {
        printf ("fpclassify (LDBL_MIN-epsilon) failed: %d (%La)\n",
-               fpclassify (ld), ld);
+               fpclassify (v2), v2);
        result = 1;
       }
-    v2.d = ld = nextafterl (ld, LDBL_MIN);
-    if (fpclassify (ld) != FP_NORMAL)
+    v2 = nextafterl (v2, LDBL_MIN);
+    if (fpclassify (v2) != FP_NORMAL)
       {
        printf ("fpclassify (LDBL_MIN-epsilon+epsilon) failed: %d (%La)\n",
-               fpclassify (ld), ld);
+               fpclassify (v2), v2);
        result = 1;
       }
 
-    if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
+    if (v1 != v2)
       {
-       printf ("LDBL_MIN: mantissa0 differs: %8x vs %8x\n",
-               v1.ieee.mantissa0, v2.ieee.mantissa0);
-       result = 1;
-      }
-    if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
-      {
-       printf ("LDBL_MIN: mantissa1 differs: %8x vs %8x\n",
-               v1.ieee.mantissa1, v2.ieee.mantissa1);
-       result = 1;
-      }
-    if (v1.ieee.exponent != v2.ieee.exponent)
-      {
-       printf ("LDBL_MIN: exponent differs: %4x vs %4x\n",
-               v1.ieee.exponent, v2.ieee.exponent);
-       result = 1;
-      }
-    if (v1.ieee.negative != v2.ieee.negative)
-      {
-       printf ("LDBL_MIN: negative differs: %d vs %d\n",
-               v1.ieee.negative, v2.ieee.negative);
+       printf ("LDBL_MIN-epsilon+epsilon != LDBL_MIN: %La vs %La\n", v2, v1);
        result = 1;
       }
 
-    v1.d = ld = -LDBL_MIN;
-    if (fpclassify (ld) != FP_NORMAL)
+    v1 = -LDBL_MIN;
+    if (fpclassify (v1) != FP_NORMAL)
       {
-       printf ("fpclassify (-LDBL_MIN) failed: %d\n", fpclassify (ld));
+       printf ("fpclassify (-LDBL_MIN) failed: %d (%La)\n",
+               fpclassify (v1), v1);
        result = 1;
       }
-    ld = nextafterl (ld, -LDBL_MIN / 2.0);
-    if (fpclassify (ld) != FP_SUBNORMAL)
+    v2 = nextafterl (v1, -LDBL_MIN / 2.0);
+    if (fpclassify (v2) != FP_SUBNORMAL)
       {
        printf ("fpclassify (-LDBL_MIN-epsilon) failed: %d (%La)\n",
-               fpclassify (ld), ld);
+               fpclassify (v2), v2);
        result = 1;
       }
-    v2.d = ld = nextafterl (ld, -LDBL_MIN);
-    if (fpclassify (ld) != FP_NORMAL)
+    v2 = nextafterl (v2, -LDBL_MIN);
+    if (fpclassify (v2) != FP_NORMAL)
       {
        printf ("fpclassify (-LDBL_MIN-epsilon+epsilon) failed: %d (%La)\n",
-               fpclassify (ld), ld);
+               fpclassify (v2), v2);
        result = 1;
       }
 
-    if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
-      {
-       printf ("-LDBL_MIN: mantissa0 differs: %8x vs %8x\n",
-               v1.ieee.mantissa0, v2.ieee.mantissa0);
-       result = 1;
-      }
-    if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
-      {
-       printf ("-LDBL_MIN: mantissa1 differs: %8x vs %8x\n",
-               v1.ieee.mantissa1, v2.ieee.mantissa1);
-       result = 1;
-      }
-    if (v1.ieee.exponent != v2.ieee.exponent)
+    if (v1 != v2)
       {
-       printf ("-LDBL_MIN: exponent differs: %4x vs %4x\n",
-               v1.ieee.exponent, v2.ieee.exponent);
-       result = 1;
-      }
-    if (v1.ieee.negative != v2.ieee.negative)
-      {
-       printf ("-LDBL_MIN: negative differs: %d vs %d\n",
-               v1.ieee.negative, v2.ieee.negative);
+       printf ("-LDBL_MIN-epsilon+epsilon != -LDBL_MIN: %La vs %La\n", v2, v1);
        result = 1;
       }
 
-    ld = LDBL_MAX;
-    if (fpclassify (ld) != FP_NORMAL)
+    v1 = LDBL_MAX;
+    if (fpclassify (v1) != FP_NORMAL)
       {
-       printf ("fpclassify (LDBL_MAX) failed: %d\n", fpclassify (ld));
+       printf ("fpclassify (LDBL_MAX) failed: %d (%La)\n",
+               fpclassify (v1), v1);
        result = 1;
       }
-    ld = nextafterl (ld, INFINITY);
-    if (fpclassify (ld) != FP_INFINITE)
+    v2 = nextafterl (v1, INFINITY);
+    if (fpclassify (v2) != FP_INFINITE)
       {
-       printf ("fpclassify (LDBL_MAX+epsilon) failed: %d\n", fpclassify (ld));
+       printf ("fpclassify (LDBL_MAX+epsilon) failed: %d (%La)\n",
+               fpclassify (v2), v2);
        result = 1;
       }
 
-    ld = -LDBL_MAX;
-    if (fpclassify (ld) != FP_NORMAL)
+    v1 = -LDBL_MAX;
+    if (fpclassify (v1) != FP_NORMAL)
       {
-       printf ("fpclassify (-LDBL_MAX) failed: %d\n", fpclassify (ld));
+       printf ("fpclassify (-LDBL_MAX) failed: %d (%La)\n",
+               fpclassify (v1), v1);
        result = 1;
       }
-    ld = nextafterl (ld, -INFINITY);
-    if (fpclassify (ld) != FP_INFINITE)
+    v2 = nextafterl (v1, -INFINITY);
+    if (fpclassify (v2) != FP_INFINITE)
       {
-       printf ("fpclassify (-LDBL_MAX-epsilon) failed: %d\n",
-               fpclassify (ld));
+       printf ("fpclassify (-LDBL_MAX-epsilon) failed: %d (%La)\n",
+               fpclassify (v2), v2);
        result = 1;
       }
 
-    v1.d = ld = 0.0625;
-    ld = nextafterl (ld, 0.0);
-    v2.d = ld = nextafterl (ld, 1.0);
+    v1 = 0.0625;
+    v2 = nextafterl (v1, 0.0);
+    v2 = nextafterl (v2, 1.0);
 
-    if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
-      {
-       printf ("0.0625L down: mantissa0 differs: %8x vs %8x\n",
-               v1.ieee.mantissa0, v2.ieee.mantissa0);
-       result = 1;
-      }
-    if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
-      {
-       printf ("0.0625L down: mantissa1 differs: %8x vs %8x\n",
-               v1.ieee.mantissa1, v2.ieee.mantissa1);
-       result = 1;
-      }
-    if (v1.ieee.exponent != v2.ieee.exponent)
-      {
-       printf ("0.0625L down: exponent differs: %4x vs %4x\n",
-               v1.ieee.exponent, v2.ieee.exponent);
-       result = 1;
-      }
-    if (v1.ieee.negative != v2.ieee.negative)
+    if (v1 != v2)
       {
-       printf ("0.0625L down: negative differs: %d vs %d\n",
-               v1.ieee.negative, v2.ieee.negative);
+       printf ("0.0625L-epsilon+epsilon != 0.0625L: %La vs %La\n", v2, v1);
        result = 1;
       }
 
-    v1.d = ld = 0.0625;
-    ld = nextafterl (ld, 1.0);
-    v2.d = ld = nextafterl (ld, 0.0);
+    v1 = 0.0625;
+    v2 = nextafterl (v1, 1.0);
+    v2 = nextafterl (v2, 0.0);
 
-    if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
-      {
-       printf ("0.0625L up: mantissa0 differs: %8x vs %8x\n",
-               v1.ieee.mantissa0, v2.ieee.mantissa0);
-       result = 1;
-      }
-    if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
-      {
-       printf ("0.0625L up: mantissa1 differs: %8x vs %8x\n",
-               v1.ieee.mantissa1, v2.ieee.mantissa1);
-       result = 1;
-      }
-    if (v1.ieee.exponent != v2.ieee.exponent)
-      {
-       printf ("0.0625L up: exponent differs: %4x vs %4x\n",
-               v1.ieee.exponent, v2.ieee.exponent);
-       result = 1;
-      }
-    if (v1.ieee.negative != v2.ieee.negative)
+    if (v1 != v2)
       {
-       printf ("0.0625L up: negative differs: %d vs %d\n",
-               v1.ieee.negative, v2.ieee.negative);
+       printf ("0.0625L+epsilon-epsilon != 0.0625L: %La vs %La\n", v2, v1);
        result = 1;
       }
 
-    v1.d = ld = -0.0625;
-    ld = nextafterl (ld, 0.0);
-    v2.d = ld = nextafterl (ld, -1.0);
+    v1 = -0.0625;
+    v2 = nextafterl (v1, 0.0);
+    v2 = nextafterl (v2, -1.0);
 
-    if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
+    if (v1 != v2)
       {
-       printf ("-0.0625L up: mantissa0 differs: %8x vs %8x\n",
-               v1.ieee.mantissa0, v2.ieee.mantissa0);
-       result = 1;
-      }
-    if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
-      {
-       printf ("-0.0625L up: mantissa1 differs: %8x vs %8x\n",
-               v1.ieee.mantissa1, v2.ieee.mantissa1);
-       result = 1;
-      }
-    if (v1.ieee.exponent != v2.ieee.exponent)
-      {
-       printf ("-0.0625L up: exponent differs: %4x vs %4x\n",
-               v1.ieee.exponent, v2.ieee.exponent);
-       result = 1;
-      }
-    if (v1.ieee.negative != v2.ieee.negative)
-      {
-       printf ("-0.0625L up: negative differs: %d vs %d\n",
-               v1.ieee.negative, v2.ieee.negative);
+       printf ("-0.0625L+epsilon-epsilon != -0.0625L: %La vs %La\n", v2, v1);
        result = 1;
       }
 
-    v1.d = ld = -0.0625;
-    ld = nextafterl (ld, -1.0);
-    v2.d = ld = nextafterl (ld, 0.0);
+    v1 = -0.0625;
+    v2 = nextafterl (v1, -1.0);
+    v2 = nextafterl (v2, 0.0);
 
-    if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
-      {
-       printf ("-0.0625L down: mantissa0 differs: %8x vs %8x\n",
-               v1.ieee.mantissa0, v2.ieee.mantissa0);
-       result = 1;
-      }
-    if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
-      {
-       printf ("-0.0625L down: mantissa1 differs: %8x vs %8x\n",
-               v1.ieee.mantissa1, v2.ieee.mantissa1);
-       result = 1;
-      }
-    if (v1.ieee.exponent != v2.ieee.exponent)
+    if (v1 != v2)
       {
-       printf ("-0.0625L down: exponent differs: %4x vs %4x\n",
-               v1.ieee.exponent, v2.ieee.exponent);
-       result = 1;
-      }
-    if (v1.ieee.negative != v2.ieee.negative)
-      {
-       printf ("-0.0625L down: negative differs: %d vs %d\n",
-               v1.ieee.negative, v2.ieee.negative);
+       printf ("-0.0625L-epsilon+epsilon != -0.0625L: %La vs %La\n", v2, v1);
        result = 1;
       }
 
-    v1.d = ld = 0.0;
-    ld = nextafterl (ld, 1.0);
-    v2.d = nextafterl (ld, -1.0);
+    v1 = 0.0;
+    v2 = nextafterl (v1, 1.0);
+    v2 = nextafterl (v2, -1.0);
 
-    if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
-      {
-       printf ("0.0L up: mantissa0 differs: %8x vs %8x\n",
-               v1.ieee.mantissa0, v2.ieee.mantissa0);
-       result = 1;
-      }
-    if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
-      {
-       printf ("0.0L up: mantissa1 differs: %8x vs %8x\n",
-               v1.ieee.mantissa1, v2.ieee.mantissa1);
-       result = 1;
-      }
-    if (v1.ieee.exponent != v2.ieee.exponent)
+    if (v1 != v2)
       {
-       printf ("0.0L up: exponent differs: %4x vs %4x\n",
-               v1.ieee.exponent, v2.ieee.exponent);
+       printf ("0.0+epsilon-epsilon != 0.0L: %La vs %La\n", v2, v1);
        result = 1;
       }
-    if (0 != v2.ieee.negative)
+    if (signbit (v2))
       {
-       printf ("0.0L up: negative differs: 0 vs %d\n",
-               v2.ieee.negative);
+       printf ("0.0+epsilon-epsilon is negative\n");
        result = 1;
       }
 
-    v1.d = ld = 0.0;
-    ld = nextafterl (ld, -1.0);
-    v2.d = nextafterl (ld, 1.0);
+    v1 = 0.0;
+    v2 = nextafterl (v1, -1.0);
+    v2 = nextafterl (v2, 1.0);
 
-    if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
-      {
-       printf ("0.0L down: mantissa0 differs: %8x vs %8x\n",
-               v1.ieee.mantissa0, v2.ieee.mantissa0);
-       result = 1;
-      }
-    if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
+    if (v1 != v2)
       {
-       printf ("0.0L down: mantissa1 differs: %8x vs %8x\n",
-               v1.ieee.mantissa1, v2.ieee.mantissa1);
-       result = 1;
-      }
-    if (v1.ieee.exponent != v2.ieee.exponent)
-      {
-       printf ("0.0L down: exponent differs: %4x vs %4x\n",
-               v1.ieee.exponent, v2.ieee.exponent);
+       printf ("0.0-epsilon+epsilon != 0.0L: %La vs %La\n", v2, v1);
        result = 1;
       }
-    if (1 != v2.ieee.negative)
+    if (!signbit (v2))
       {
-       printf ("0.0L down: negative differs: 1 vs %d\n",
-               v2.ieee.negative);
+       printf ("0.0-epsilon+epsilon is positive\n");
        result = 1;
       }
 
index 9b638017f84170e013167460b916db9b9996c350..d1875f4fd99b74e590b271d56e6b19cef6dcdeae 100644 (file)
@@ -1,3 +1,7 @@
+2013-19-04  Alan Modra  <amodra@gmail.com>
+
+       * sysdeps/ia64/fpu/printf_fphex.c: Adjust for fpnum change.
+
 2013-09-02  Mike Frysinger  <vapier@gentoo.org>
 
        * sysdeps/unix/sysv/linux/hppa/syscalls.list (fanotify_mark): New
index ca021427eca4e054196816016d761d23836fafc5..0698cdad10159f871a85912371338470081cdb28 100644 (file)
@@ -25,9 +25,11 @@ do {                                                                       \
       /* The "strange" 80 bit format on ia64 has an explicit                 \
         leading digit in the 64 bit mantissa.  */                            \
       unsigned long long int num;                                            \
+      union ieee854_long_double u;                                           \
+      u.d = fpnum.ldbl;                                                              \
                                                                              \
-      num = (((unsigned long long int) fpnum.ldbl.ieee.mantissa0) << 32              \
-            | fpnum.ldbl.ieee.mantissa1);                                    \
+      num = (((unsigned long long int) u.ieee.mantissa0) << 32               \
+            | u.ieee.mantissa1);                                             \
                                                                              \
       zero_mantissa = num == 0;                                                      \
                                                                              \
@@ -49,8 +51,8 @@ do {                                                                        \
                                                                              \
       /* We have 3 bits from the mantissa in the leading nibble.             \
         Therefore we are here using `IEEE854_LONG_DOUBLE_BIAS + 3'.  */      \
-      exponent = fpnum.ldbl.ieee.exponent;                                   \
-                                                                            \
+      exponent = u.ieee.exponent;                                            \
+                                                                             \
       if (exponent == 0)                                                     \
        {                                                                     \
          if (zero_mantissa)                                                  \
index e20eab6a345c8a2202140c78d207a92f137b8b6f..2b93e6c57ac135a0bd5460d5dd75578ef22bc494 100644 (file)
@@ -332,8 +332,7 @@ ___printf_fp (FILE *fp,
       int res;
       if (__isnanl (fpnum.ldbl))
        {
-         union ieee854_long_double u = { .d = fpnum.ldbl };
-         is_neg = u.ieee.negative != 0;
+         is_neg = signbit (fpnum.ldbl);
          if (isupper (info->spec))
            {
              special = "NAN";
index 3da2eec40cdc4a1399e68555a69117adda8c91de..50b6fbfbc7a4179890775574530734d84ff717b1 100644 (file)
@@ -93,7 +93,7 @@ __printf_fphex (FILE *fp,
   union
     {
       union ieee754_double dbl;
-      union ieee854_long_double ldbl;
+      long double ldbl;
     }
   fpnum;
 
@@ -162,12 +162,11 @@ __printf_fphex (FILE *fp,
 #ifndef __NO_LONG_DOUBLE_MATH
   if (info->is_long_double && sizeof (long double) > sizeof (double))
     {
-      fpnum.ldbl.d = *(const long double *) args[0];
+      fpnum.ldbl = *(const long double *) args[0];
 
       /* Check for special values: not a number or infinity.  */
-      if (__isnanl (fpnum.ldbl.d))
+      if (__isnanl (fpnum.ldbl))
        {
-         negative = fpnum.ldbl.ieee.negative != 0;
          if (isupper (info->spec))
            {
              special = "NAN";
@@ -181,8 +180,7 @@ __printf_fphex (FILE *fp,
        }
       else
        {
-         int res = __isinfl (fpnum.ldbl.d);
-         if (res)
+         if (__isinfl (fpnum.ldbl))
            {
              if (isupper (info->spec))
                {
@@ -194,11 +192,9 @@ __printf_fphex (FILE *fp,
                  special = "inf";
                  wspecial = L"inf";
                }
-             negative = res < 0;
            }
-         else
-           negative = signbit (fpnum.ldbl.d);
        }
+      negative = signbit (fpnum.ldbl);
     }
   else
 #endif /* no long double */
index 2c496e5b2f0f271e388983a1140e54d65e03db59..dfb3a539b92f70ffc6ce5aafed627e4a1f7ca052 100644 (file)
@@ -103,7 +103,7 @@ __printf_size (FILE *fp, const struct printf_info *info,
   union
     {
       union ieee754_double dbl;
-      union ieee854_long_double ldbl;
+      long double ldbl;
     }
   fpnum;
   const void *ptr = &fpnum;
@@ -123,25 +123,25 @@ __printf_size (FILE *fp, const struct printf_info *info,
 #ifndef __NO_LONG_DOUBLE_MATH
   if (info->is_long_double && sizeof (long double) > sizeof (double))
     {
-      fpnum.ldbl.d = *(const long double *) args[0];
+      fpnum.ldbl = *(const long double *) args[0];
 
       /* Check for special values: not a number or infinity.  */
-      if (__isnanl (fpnum.ldbl.d))
+      if (__isnanl (fpnum.ldbl))
        {
          special = "nan";
          wspecial = L"nan";
          // fpnum_sign = 0;    Already zero
        }
-      else if ((res = __isinfl (fpnum.ldbl.d)))
+      else if ((res = __isinfl (fpnum.ldbl)))
        {
          fpnum_sign = res;
          special = "inf";
          wspecial = L"inf";
        }
       else
-       while (fpnum.ldbl.d >= divisor && tag[1] != '\0')
+       while (fpnum.ldbl >= divisor && tag[1] != '\0')
          {
-           fpnum.ldbl.d /= divisor;
+           fpnum.ldbl /= divisor;
            ++tag;
          }
     }
index c9e09a4b7421b7983321461f639e46eac0d89e91..e82228a53e850e3c74faee798077420c62e47fd5 100644 (file)
@@ -24,13 +24,15 @@ do {                                                                              \
         digits we use only the implicit digits for the number before         \
         the decimal point.  */                                               \
       unsigned long long int num0, num1;                                     \
+      union ieee854_long_double u;                                           \
+      u.d = fpnum.ldbl;                                                              \
                                                                              \
       assert (sizeof (long double) == 16);                                   \
                                                                              \
-      num0 = (((unsigned long long int) fpnum.ldbl.ieee.mantissa0) << 32      \
-            | fpnum.ldbl.ieee.mantissa1);                                    \
-      num1 = (((unsigned long long int) fpnum.ldbl.ieee.mantissa2) << 32      \
-            | fpnum.ldbl.ieee.mantissa3);                                    \
+      num0 = (((unsigned long long int) u.ieee.mantissa0) << 32                      \
+            | u.ieee.mantissa1);                                             \
+      num1 = (((unsigned long long int) u.ieee.mantissa2) << 32                      \
+            | u.ieee.mantissa3);                                             \
                                                                              \
       zero_mantissa = (num0|num1) == 0;                                              \
                                                                              \
@@ -75,9 +77,9 @@ do {                                                                        \
          *--wnumstr = L'0';                                                  \
        }                                                                     \
                                                                              \
-      leading = fpnum.ldbl.ieee.exponent == 0 ? '0' : '1';                   \
+      leading = u.ieee.exponent == 0 ? '0' : '1';                            \
                                                                              \
-      exponent = fpnum.ldbl.ieee.exponent;                                   \
+      exponent = u.ieee.exponent;                                            \
                                                                              \
       if (exponent == 0)                                                     \
        {                                                                     \
index 453c2be7b83894be00aa9dcba28cd7531ab6a5f9..e0ec422b01ff9b9894c1f2c94b67d82ac69a9401 100644 (file)
@@ -27,14 +27,14 @@ do {                                                                              \
       unsigned long long hi, lo;                                             \
       int ediff;                                                             \
       union ibm_extended_long_double u;                                              \
-      u.ld = fpnum.ldbl.d;                                                   \
+      u.ld = fpnum.ldbl;                                                     \
                                                                              \
       assert (sizeof (long double) == 16);                                   \
                                                                              \
       lo = ((long long)u.d[1].ieee.mantissa0 << 32) | u.d[1].ieee.mantissa1;  \
       hi = ((long long)u.d[0].ieee.mantissa0 << 32) | u.d[0].ieee.mantissa1;  \
       lo <<= 7; /* pre-shift lo to match ieee854.  */                        \
-      /* If the lower double is not a denomal or zero then set the hidden     \
+      /* If the lower double is not a denormal or zero then set the hidden    \
         53rd bit.  */                                                        \
       if (u.d[1].ieee.exponent != 0)                                         \
        lo |= (1ULL << (52 + 7));                                             \
index f356a4843e3689fb0b1146fe13039d47f5f5c82b..715c93b500c9259d797ee1c9ee357c4e27836756 100644 (file)
@@ -25,11 +25,13 @@ do {                                                                              \
       /* The "strange" 80 bit format on ix86 and m68k has an explicit        \
         leading digit in the 64 bit mantissa.  */                            \
       unsigned long long int num;                                            \
+      union ieee854_long_double u;                                           \
+      u.d = fpnum.ldbl;                                                              \
                                                                              \
       assert (sizeof (long double) == 12);                                   \
                                                                              \
-      num = (((unsigned long long int) fpnum.ldbl.ieee.mantissa0) << 32              \
-            | fpnum.ldbl.ieee.mantissa1);                                    \
+      num = (((unsigned long long int) u.ieee.mantissa0) << 32               \
+            | u.ieee.mantissa1);                                             \
                                                                              \
       zero_mantissa = num == 0;                                                      \
                                                                              \
@@ -62,7 +64,7 @@ do {                                                                        \
                                                                              \
       /* We have 3 bits from the mantissa in the leading nibble.             \
         Therefore we are here using `IEEE854_LONG_DOUBLE_BIAS + 3'.  */      \
-      exponent = fpnum.ldbl.ieee.exponent;                                   \
+      exponent = u.ieee.exponent;                                            \
                                                                              \
       if (exponent == 0)                                                     \
        {                                                                     \
index c85d1f79fbd42fc1f41bc22649d179bd9a12beee..be55f9cf6bb6a48bc48845c0228fc3c4e1ed566d 100644 (file)
@@ -25,10 +25,11 @@ do {                                                                              \
       /* The "strange" 80 bit format on ix86 and m68k has an explicit        \
         leading digit in the 64 bit mantissa.  */                            \
       unsigned long long int num;                                            \
+      union ieee854_long_double u;                                           \
+      u.d = fpnum.ldbl;                                                              \
                                                                              \
-                                                                             \
-      num = (((unsigned long long int) fpnum.ldbl.ieee.mantissa0) << 32              \
-            | fpnum.ldbl.ieee.mantissa1);                                    \
+      num = (((unsigned long long int) u.ieee.mantissa0) << 32               \
+            | u.ieee.mantissa1);                                             \
                                                                              \
       zero_mantissa = num == 0;                                                      \
                                                                              \
@@ -61,7 +62,7 @@ do {                                                                        \
                                                                              \
       /* We have 3 bits from the mantissa in the leading nibble.             \
         Therefore we are here using `IEEE854_LONG_DOUBLE_BIAS + 3'.  */      \
-      exponent = fpnum.ldbl.ieee.exponent;                                   \
+      exponent = u.ieee.exponent;                                            \
                                                                              \
       if (exponent == 0)                                                     \
        {                                                                     \