Lines Matching refs:z
178 double z, p, x1; in large_lgam() local
193 z = x1*x1; in large_lgam()
194 p = pb0+z*(pb1+z*(pb2+z*(pb3+z*(pb4+z*(pb5+z*(pb6+z*pb7)))))); in large_lgam()
212 double y, z, t, r = 0, p, q, hi, lo; in small_lgam() local
223 z = t - x0_lo; in small_lgam()
224 p = r0+z*(r1+z*(r2+z*(r3+z*(r4+z*r5)))); in small_lgam()
225 q = s0+z*(s1+z*(s2+z*(s3+z*s4))); in small_lgam()
226 r = t*(z*(p/q) - x0_lo); in small_lgam()
228 z = 1.0; in small_lgam()
230 case 6: z = (y + 5); in small_lgam()
231 case 5: z *= (y + 4); in small_lgam()
232 case 4: z *= (y + 3); in small_lgam()
233 case 3: z *= (y + 2); in small_lgam()
234 rr = __log__D(z); in small_lgam()
248 z = y-t; in small_lgam()
251 r = lo*y + z*hi; /* q + r = y*(a0+p/q) */ in small_lgam()
253 z = 1.0; in small_lgam()
255 case 6: z = (y + 5); in small_lgam()
256 case 5: z *= (y + 4); in small_lgam()
257 case 4: z *= (y + 3); in small_lgam()
258 case 3: z *= (y + 2); in small_lgam()
259 rr = __log__D(z); in small_lgam()
279 double y, z, one = 1.0, zero = 0.0; in neg_lgam() local
295 z = floor(x + .5); in neg_lgam()
296 if (z == x) { /* convention: G(-(integer)) -> +Inf */ in neg_lgam()
306 z = fabs(x + z); /* 0 < z <= .5 */ in neg_lgam()
307 if (z < .25) in neg_lgam()
308 z = sin(M_PI*z); in neg_lgam()
310 z = cos(M_PI*(0.5-z)); in neg_lgam()
311 z = log(M_PI/(z*x)); in neg_lgam()
313 return (z - y); in neg_lgam()