Lines Matching refs:x_adj
66 double x_adj = x - n; in gamma_positive() local
68 double prod = __gamma_product (x_adj, 0, n, &eps); in gamma_positive()
69 return (__ieee754_exp (__ieee754_lgamma_r (x_adj, &local_signgam)) in gamma_positive()
76 double x_adj = x; in gamma_positive() local
83 x_adj = math_narrow_eval (x + n); in gamma_positive()
84 x_eps = (x - (x_adj - n)); in gamma_positive()
85 prod = __gamma_product (x_adj - n, x_eps, n, &eps); in gamma_positive()
91 double x_adj_int = round (x_adj); in gamma_positive()
92 double x_adj_frac = x_adj - x_adj_int; in gamma_positive()
94 double x_adj_mant = __frexp (x_adj, &x_adj_log2); in gamma_positive()
102 mul_split (&h1, &l1, __ieee754_pow (x_adj_mant, x_adj), in gamma_positive()
104 mul_split (&h2, &l2, __ieee754_exp (-x_adj), sqrt (2 * M_PI / x_adj)); in gamma_positive()
108 double exp_adj = x_eps * __ieee754_log (x_adj); in gamma_positive()
110 double x_adj2 = x_adj * x_adj; in gamma_positive()
113 exp_adj += bsum / x_adj; in gamma_positive()