invpio2
invpio2 = 0.636619772367581343075535, /* 2^ -1 * 1.45F306DC9C883 */
n = (int)(t * invpio2 + half);
w = y * invpio2;
w = y * invpio2;
w = y * invpio2;
w = y * invpio2;
n1 = (int) (x1 * invpio2 + half[xsb1]);
n0 = (int) (x0 * invpio2 + half[xsb0]);
n0 = (int) (x0 * invpio2 + half[xsb0]);
n1 = (int) (x1 * invpio2 + half[xsb1]);
n2 = (int) (x2 * invpio2 + half[xsb2]);
invpio2 = 0.636619772367581343075535, /* 53 bits of pi/2 */
n0 = (int) (x0 * invpio2 + half[xsb0]);
n1 = (int) (x1 * invpio2 + half[xsb1]);
n2 = (int) (x2 * invpio2 + half[xsb2]);
invpio2 = 0.636619772367581343075535,
n1 = (int) (x1 * invpio2 + half[xsb1]);
n0 = (int) (x0 * invpio2 + half[xsb0]);
z##N = y##N * invpio2 + c3two51; \
z0 = y0 * invpio2 + c3two51;
z1 = y1 * invpio2 + c3two51;
z2 = y2 * invpio2 + c3two51;
z3 = y3 * invpio2 + c3two51;
n1 = (int) (x1 * invpio2 + half[xsb1]);
n0 = (int) (x0 * invpio2 + half[xsb0]);
invpio2 = 0.636619772367581343075535,
n0 = (int) (x0 * invpio2 + half[xsb0]);
n1 = (int) (x1 * invpio2 + half[xsb1]);
n2 = (int) (x2 * invpio2 + half[xsb2]);
n0 = (int) (x0 * invpio2 + half[xsb0]);
n1 = (int) (x1 * invpio2 + half[xsb1]);
n2 = (int) (x2 * invpio2 + half[xsb2]);
invpio2 = 0.636619772367581343075535,
n1 = (int) (x1 * invpio2 + half[xsb1]);
n0 = (int) (x0 * invpio2 + half[xsb0]);
n1 = (int) (x1 * invpio2 + half[xsb1]);
n0 = (int) (x0 * invpio2 + half[xsb0]);
invpio2 = 0.636619772367581343075535, /* 53 bits of pi/2 */
n0 = (int) (x0 * invpio2 + half[xsb0]);
n1 = (int) (x1 * invpio2 + half[xsb1]);
n2 = (int) (x2 * invpio2 + half[xsb2]);
z##N = y##N * invpio2 + c3two51; \
z0 = y0 * invpio2 + c3two51;
z1 = y1 * invpio2 + c3two51;
z2 = y2 * invpio2 + c3two51;
z3 = y3 * invpio2 + c3two51;
z##N = y##N * invpio2 + c3two51; \
z0 = y0 * invpio2 + c3two51;
z1 = y1 * invpio2 + c3two51;
z2 = y2 * invpio2 + c3two51;
z3 = y3 * invpio2 + c3two51;