pio2hi
y = pio2hi; p = pio2lo;
y = -pio2hi; p = -pio2lo;
return (hx >= 0 ? pio2hi - y : y - pio2hi);
pio2hi = 1.570796326794896619231321691639751398740e+0000L,
return (pio2hi - y);
return (y - pio2hi);
pio2hi = 1.570796371e+0000F,
x = pio2hi - y;
x = y - pio2hi;
pio2hi = 1.570796326794896558e+00,
t = pio2hi - y;
*err = -(y - (pio2hi - t));
t = y - pio2hi;
*err = y - (pio2hi + t);
s = pio2hi - v;
t -= (v - (pio2hi - s));
s = v - pio2hi;
t = -(t - (v - (s + pio2hi)));