Lines Matching refs:x
30 static double sinpi(double x) in sinpi() argument
36 x = x * 0.5; in sinpi()
37 x = 2 * (x - floor(x)); in sinpi()
40 n = 4 * x; in sinpi()
42 x -= n * 0.5; in sinpi()
44 x *= pi; in sinpi()
48 return __sin(x, 0, 0); in sinpi()
50 return __cos(x, 0); in sinpi()
52 return __sin(-x, 0, 0); in sinpi()
54 return -__cos(x, 0); in sinpi()
89 static double S(double x) in S() argument
95 if (x < 8) in S()
97 num = num * x + Snum[i]; in S()
98 den = den * x + Sden[i]; in S()
102 num = num / x + Snum[i]; in S()
103 den = den / x + Sden[i]; in S()
108 double tgamma(double x) in tgamma() argument
110 union {double f; uint64_t i;} u = {x}; in tgamma()
119 return x + INFINITY; in tgamma()
122 return 1/x; in tgamma()
126 if (x == floor(x)) { in tgamma()
129 if (x <= sizeof fact/sizeof *fact) in tgamma()
130 return fact[(int)x - 1]; in tgamma()
137 FORCE_EVAL((float)(0x1p-126/x)); in tgamma()
138 if (floor(x) * 0.5 == floor(x * 0.5)) in tgamma()
142 x *= 0x1p1023; in tgamma()
143 return x; in tgamma()
146 absx = sign ? -x : x; in tgamma()
160 if (x < 0) { in tgamma()
174 double __lgamma_r(double x, int *sign) in __lgamma_r() argument
181 if (!isfinite(x)) in __lgamma_r()
183 return x*x; in __lgamma_r()
186 if (x == floor(x) && x <= 2) { in __lgamma_r()
189 if (x <= 0) in __lgamma_r()
194 absx = fabs(x); in __lgamma_r()
198 *sign = 1 - 2*!!signbit(x); in __lgamma_r()
204 x = tgamma(x); in __lgamma_r()
205 *sign = 1 - 2*!!signbit(x); in __lgamma_r()
206 return log(fabs(x)); in __lgamma_r()
212 if (x < 0) { in __lgamma_r()
214 x = sinpi(absx); in __lgamma_r()
215 *sign = 2*!!signbit(x) - 1; in __lgamma_r()
216 r = log(pi/(fabs(x)*absx)) - r; in __lgamma_r()