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;