*** 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;