pio2_hi
return (pio2_hi + pio2_lo); /* if |x| < 2**-57 */
return (pio2_hi - (x - (pio2_lo - x * r)));
return (x * pio2_hi + x * pio2_lo);
t = pio2_hi - (2.0 * (s + s * w) - pio2_lo);