DBL_MANT_DIG
res_ptr[N - 1] |= 1L << (DBL_MANT_DIG - 1 - ((N - 1) * BITS_PER_MP_LIMB));
- (DBL_MANT_DIG - ((N - 1) * BITS_PER_MP_LIMB)))
<< (DBL_MANT_DIG - 32)) - 1);
<< (DBL_MANT_DIG - 32)) - 1);
to_shift = 1 + p.fracsize * BITS_PER_MP_LIMB - DBL_MANT_DIG;