--- _ldtoa.c.orig 2004-06-03 15:17:18.000000000 -0700
+++ _ldtoa.c 2005-10-08 22:43:25.000000000 -0700
@@ -61,14 +61,34 @@
char *ret;
union IEEEl2bits u;
uint32_t bits[(LDBL_MANT_DIG + 31) / 32];
+#if defined(__ppc__) || defined(__ppc64__)
+ int type;
+#endif /* defined(__ppc__) || defined(__ppc64__) */
u.e = *ld;
+#if defined(__ppc__) || defined(__ppc64__)
+ /*
+ * Subnormal head-tail doubles don't seem to be converted correctly
+ * by gdtoa. So we multiply by 10^32 to make them normal then
+ * subtract 32 from the exponent later.
+ */
+ if ((type = __fpclassify(u.e)) == FP_NORMAL && __fpclassifyd(u.d[1]) == FP_SUBNORMAL)
+ type = FP_SUBNORMAL;
+ if (type == FP_SUBNORMAL)
+ u.e *= 1.0e32L;
+#endif /* defined(__ppc__) || defined(__ppc64__) */
*sign = u.bits.sign;
be = u.bits.exp - (LDBL_MAX_EXP - 1) - (LDBL_MANT_DIG - 1);
LDBL_TO_ARRAY32(u, bits);
+#if defined(__ppc__) || defined(__ppc64__)
+ switch (type) {
+ case FP_SUBNORMAL:
+#else /* !defined(__ppc__) && !defined(__ppc64__) */
switch (fpclassify(u.e)) {
+#endif /* defined(__ppc__) || defined(__ppc64__) */
case FP_NORMAL:
+ case FP_SUPERNORMAL:
kind = STRTOG_Normal;
#ifdef LDBL_IMPLICIT_NBIT
bits[LDBL_MANT_DIG / 32] |= 1 << ((LDBL_MANT_DIG - 1) % 32);
@@ -77,12 +97,12 @@
case FP_ZERO:
kind = STRTOG_Zero;
break;
+#if !defined(__ppc__) && !defined(__ppc64__)
case FP_SUBNORMAL:
kind = STRTOG_Denormal;
-#ifdef LDBL_IMPLICIT_NBIT
be++;
-#endif
break;
+#endif /* !defined(__ppc__) && !defined(__ppc64__) */
case FP_INFINITE:
kind = STRTOG_Infinite;
break;
@@ -96,5 +116,9 @@
ret = gdtoa(&fpi, be, (ULong *)bits, &kind, mode, ndigits, decpt, rve);
if (*decpt == -32768)
*decpt = INT_MAX;
+#if defined(__ppc__) || defined(__ppc64__)
+ else if (type == FP_SUBNORMAL)
+ *decpt -= 32;
+#endif /* defined(__ppc__) || defined(__ppc64__) */
return ret;
}