Update.

1999-03-09  Andreas Schwab  <schwab@issan.cs.uni-dortmund.de>

	* stdlib/strtod.c (round_and_return): Fix rounding of denormals.
	(DENORM_EXP): New define.  Use it as the exponent for a denormal.
	* stdlib/strtold.c: Moved to ...
	* sysdeps/generic/strtold.c: ... here.
	* sysdeps/m68k/strtold.c: New file.

1999-03-09  Andreas Schwab  <schwab@issan.cs.uni-dortmund.de>

	* stdio-common/printf_fphex.c (__printf_fphex): Print no decimal
	point if the precision is zero, even if the mantissa is non-zero.

diff --git a/ChangeLog b/ChangeLog
index 5dbcd4d..0548fb3 100644 (file)
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,16 @@
+1999-03-09  Andreas Schwab  <schwab@issan.cs.uni-dortmund.de>
+
+       * stdlib/strtod.c (round_and_return): Fix rounding of denormals.
+       (DENORM_EXP): New define.  Use it as the exponent for a denormal.
+       * stdlib/strtold.c: Moved to ...
+       * sysdeps/generic/strtold.c: ... here.
+       * sysdeps/m68k/strtold.c: New file.
+
+1999-03-09  Andreas Schwab  <schwab@issan.cs.uni-dortmund.de>
+
+       * stdio-common/printf_fphex.c (__printf_fphex): Print no decimal
+       point if the precision is zero, even if the mantissa is non-zero.
+
 1999-03-09  Roland McGrath  <roland@baalperazim.frob.com>
 
        * sysdeps/mach/hurd/ioctl.c (__ioctl): Move the marshalling of the

diff --git a/stdlib/strtod.c b/stdlib/strtod.c
index 5b15859..0de96a5 100644 (file)
--- a/stdlib/strtod.c
+++ b/stdlib/strtod.c
@@ -226,7 +226,7 @@ round_and_return (mp_limb_t *retval, int exponent, int negative,
          int i;
 
          round_limb = retval[RETURN_LIMB_SIZE - 1];
-         round_bit = BITS_PER_MP_LIMB - 1;
+         round_bit = (MANT_DIG - 1) % BITS_PER_MP_LIMB;
          for (i = 0; i < RETURN_LIMB_SIZE; ++i)
            more_bits |= retval[i] != 0;
          MPN_ZERO (retval, RETURN_LIMB_SIZE);
@@ -254,7 +254,13 @@ round_and_return (mp_limb_t *retval, int exponent, int negative,
           round_bit = shift - 1;
          (void) __mpn_rshift (retval, retval, RETURN_LIMB_SIZE, shift);
        }
-      exponent = MIN_EXP - 2;
+      /* This is a hook for the m68k long double format, where the
+        exponent bias is the same for normalized and denormalized
+        numbers.  */
+#ifndef DENORM_EXP
+# define DENORM_EXP (MIN_EXP - 2)
+#endif
+      exponent = DENORM_EXP;
     }
 
   if ((round_limb & (((mp_limb_t) 1) << round_bit)) != 0
@@ -273,7 +279,7 @@ round_and_return (mp_limb_t *retval, int exponent, int negative,
          retval[RETURN_LIMB_SIZE - 1]
            |= ((mp_limb_t) 1) << ((MANT_DIG - 1) % BITS_PER_MP_LIMB);
        }
-      else if (exponent == MIN_EXP - 2
+      else if (exponent == DENORM_EXP
               && (retval[RETURN_LIMB_SIZE - 1]
                   & (((mp_limb_t) 1) << ((MANT_DIG - 1) % BITS_PER_MP_LIMB)))
               != 0)

diff --git a/stdlib/strtold.c b/sysdeps/generic/strtold.c
similarity index 100%
rename from stdlib/strtold.c
rename to sysdeps/generic/strtold.c

diff --git a/sysdeps/m68k/strtold.c b/sysdeps/m68k/strtold.c
new file mode 100644 (file)
index 0000000..dd7fbce
--- /dev/null
+++ b/sysdeps/m68k/strtold.c
@@ -0,0 +1,2 @@
+#define DENORM_EXP (MIN_EXP - 1)
+#include <sysdeps/generic/strtold.c>