twom53
*pe = (nrx * twom53 * *x < half - rmx)? 1 : 2;
ay *= twom53;
ay *= twom53;
ay *= twom53;
if (onep1u * onep1u != onep1u + twom53) {
twom53 = 1.11022302462515654042e-16, /* 0x3ca00000 0 = 2**-53 */
ay *= twom53;
ay *= twom53;
} else if (onep1u * onep1u != onep1u + twom53) {
twom53 = 1.11022302462515654042e-16, /* 0x3ca00000 0 = 2**-53 */
ay *= twom53;