1 /* Double-precision x^y function.
56 /* Compute y+TAIL = log(x) where the rounded result is y and TAIL has about
63   double_t z, r, y, invc, logc, logctail, kd, hi, t1, t2, lo, lo1, lo2, p;  in log_inline()  local
124   y = hi + lo;  in log_inline()
125   *tail = hi - y + lo;  in log_inline()
126   return y;  in log_inline()
153   double_t scale, y;  in specialcase()  local
160       y = 0x1p1009 * (scale + scale * tmp);  in specialcase()
161       return check_oflow (y);  in specialcase()
167   y = scale + scale * tmp;  in specialcase()
171   if (fabs (y) < 1.0)  in specialcase()
174       /* Round y to the right precision before scaling it into the subnormal  in specialcase()
179       if (y < 0.0)  in specialcase()
181       lo = scale - y + scale * tmp;  in specialcase()
182       hi = one + y;  in specialcase()
183       lo = one - hi + y + lo;  in specialcase()
184       y = eval_as_double (hi + lo) - one;  in specialcase()
186       if (y == 0.0)  in specialcase()
187 	y = asfloat64 (sbits & 0x8000000000000000);  in specialcase()
191   y = 0x1p-1022 * y;  in specialcase()
192   return check_uflow (y);  in specialcase()
300 pow (double x, double y)  in pow()  argument
307   iy = asuint64 (y);  in pow()
309   topy = top12 (y);  in pow()
313       /* Note: if |y| > 1075 * ln2 * 2^53 ~= 0x1.749p62 then pow(x,y) = inf/0  in pow()
314 	 and if |y| < 2^-54 / 1075 ~= 0x1.e7b6p-65 then pow(x,y) = +-1.  */  in pow()
316 	 (|y| < 0x1p-65 or |y| >= 0x1p63 or nan).  */  in pow()
320 	    return issignaling64_inline (x) ? x + y : 1.0;  in pow()
322 	    return issignaling64_inline (y) ? x + y : 1.0;  in pow()
325 	    return x + y;  in pow()
329 	    return 0.0; /* |x|<1 && y==inf or |x|>1 && y==-inf.  */  in pow()
330 	  return y * y;  in pow()
346       /* Here x and y are non-zero finite.  */  in pow()
360 	  /* Note: sign_bias == 0 here because y is not odd.  */  in pow()
365 	      /* |y| < 2^-65, x^y ~= 1 + y*log(x).  */  in pow()
367 		return ix > asuint64 (1.0) ? 1.0 + y : 1.0 - y;  in pow()
387   ehi = y * hi;  in pow()
388   elo = y * lo + fma (y, hi, -ehi);  in pow()
391   double_t ylo = y - yhi;  in pow()
395   elo = ylo * lhi + y * llo; /* |elo| < |ehi| * 2^-25.  */  in pow()