Print this page

        

*** 492,502 **** if (n > 78 || y == 0.0L) { #else if (n > 122 || y == 0.0L) { #endif ! XFSCALE(x, 0x3fff - (ix >> 16)); i = ((ix & 0xffff) + 0x100) >> 9; /* 7.5 bits of x */ zk = 1.0L + ((long double) i) * 0.0078125L; z = x - zk; dk = (double)z; --- 492,502 ---- if (n > 78 || y == 0.0L) { #else if (n > 122 || y == 0.0L) { #endif ! XFSCALE(x, (0x3fff - (ix >> 16))); i = ((ix & 0xffff) + 0x100) >> 9; /* 7.5 bits of x */ zk = 1.0L + ((long double) i) * 0.0078125L; z = x - zk; dk = (double)z;
*** 515,526 **** w = k_log_NKzl(n, k, zh, er); } else { /* * compute z = x*x + y*y */ ! XFSCALE(x, 0x3fff - (ix >> 16)); ! XFSCALE(y, 0x3fff - n - (iy >> 16)); ix = (ix & 0xffff) | 0x3fff0000; iy = (iy & 0xffff) | (0x3fff0000 - (n << 16)); nx -= 0x3fff; t1 = x * x; t2 = y * y; wh = x; --- 515,526 ---- w = k_log_NKzl(n, k, zh, er); } else { /* * compute z = x*x + y*y */ ! XFSCALE(x, (0x3fff - (ix >> 16))); ! XFSCALE(y, (0x3fff - n - (iy >> 16))); ix = (ix & 0xffff) | 0x3fff0000; iy = (iy & 0xffff) | (0x3fff0000 - (n << 16)); nx -= 0x3fff; t1 = x * x; t2 = y * y; wh = x;