pio2_1
pio2_1 = 1.570796326734125614166, /* 2^ 0 * 1.921FB54400000 */
t = fabs(x) - pio2_1;
r = t - fn * pio2_1;
y = (y + pio2_1) + pio2_t;
y = (y - n * pio2_1) - n * pio2_t;
y = (y - pio2_1) - pio2_t;
y = (y - pio2_1) - pio2_t;
y = (y + pio2_1) + pio2_t;
y = (y - n * pio2_1) - n * pio2_t;
y = (y + pio2_1) + pio2_t;
y = (y - n * pio2_1) - n * pio2_t;
y = (y - pio2_1) - pio2_t;
y = (y - pio2_1) - pio2_t;
y = (y + pio2_1) + pio2_t;
y = (y - n * pio2_1) - n * pio2_t;
a1 = x1 - fn1 * pio2_1;
a0 = x0 - fn0 * pio2_1;
a0 = x0 - fn0 * pio2_1;
a1 = x1 - fn1 * pio2_1;
a2 = x2 - fn2 * pio2_1;
pio2_1 = 1.570796326734125614166, /* first 33 bits of pi/2 */
a0 = x0 - fn0 * pio2_1;
a1 = x1 - fn1 * pio2_1;
a2 = x2 - fn2 * pio2_1;
pio2_1 = 1.570796326734125614166,
a1 = x1 - fn1 * pio2_1;
a0 = x0 - fn0 * pio2_1;
y##N = (y##N - z##N * pio2_1) - z##N * pio2_t; \
y0 = (y0 - z0 * pio2_1) - z0 * pio2_t;
y1 = (y1 - z1 * pio2_1) - z1 * pio2_t;
y2 = (y2 - z2 * pio2_1) - z2 * pio2_t;
y3 = (y3 - z3 * pio2_1) - z3 * pio2_t;
a1 = x1 - fn1 * pio2_1;
a0 = x0 - fn0 * pio2_1;
pio2_1 = 1.570796326734125614166,
a0 = x0 - fn0 * pio2_1;
a1 = x1 - fn1 * pio2_1;
a2 = x2 - fn2 * pio2_1;
a0 = x0 - fn0 * pio2_1;
a1 = x1 - fn1 * pio2_1;
a2 = x2 - fn2 * pio2_1;
pio2_1 = 1.570796326734125614166,
a1 = x1 - fn1 * pio2_1;
a0 = x0 - fn0 * pio2_1;
a1 = x1 - fn1 * pio2_1;
a0 = x0 - fn0 * pio2_1;
pio2_1 = 1.570796326734125614166, /* first 33 bits of pi/2 */
a0 = x0 - fn0 * pio2_1;
a1 = x1 - fn1 * pio2_1;
a2 = x2 - fn2 * pio2_1;
y##N = (y##N - z##N * pio2_1) - z##N * pio2_t; \
y0 = (y0 - z0 * pio2_1) - z0 * pio2_t;
y1 = (y1 - z1 * pio2_1) - z1 * pio2_t;
y2 = (y2 - z2 * pio2_1) - z2 * pio2_t;
y3 = (y3 - z3 * pio2_1) - z3 * pio2_t;
y##N = (y##N - z##N * pio2_1) - z##N * pio2_t; \
y0 = (y0 - z0 * pio2_1) - z0 * pio2_t;
y1 = (y1 - z1 * pio2_1) - z1 * pio2_t;
y2 = (y2 - z2 * pio2_1) - z2 * pio2_t;
y3 = (y3 - z3 * pio2_1) - z3 * pio2_t;