--- _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; }