DBL_FRACLBITS
r <<= DBL_FRACLBITS;
exp -= DBL_FRACLBITS;
for (; s > s0 + sigfigs - (DBL_FRACLBITS / 4) - 1 && s > s0; s--) {