118 	double z, r;  in GT1()  local
120 	z = y * y;  in GT1()
121 	r = TZ1 * y + z * ((P10 + y * P11 + z * P12) + (z * y) * (P13 + y *  in GT1()
122 		P14 + z * P15));  in GT1()
129 	double z;  in GT2()  local
131 	z = y * y;  in GT2()
132 	return (GZ2 + z * ((P20 + y * P21 + z * P22) + (z * y) * (P23 + y *  in GT2()
133 		P24 + z * P25)));  in GT2()
139 double z, r;  in GT3()  local
141 	z = y * y;  in GT3()
142 	r = TZ3 * y + z * ((P30 + y * P31 + z * P32) + (z * y) * (P33 + y *  in GT3()
309 	double ss, zz, z, t1, t2, w, y, u;  in large_gam()  local
319 	__HI(z) = (ix & 0xffffc000) | 0x2000;	/* z[j]=1+j/64+1/128 */  in large_gam()
320 	__LO(z) = 0;  in large_gam()
322 	t1 = y + z;  in large_gam()
323 	t2 = y - z;  in large_gam()
331 	z = one / x;  in large_gam()
332 	zz = z * z;  in large_gam()
333 	w = ((x - half) * ss + hln2pi) + z * (GP0 + zz * GP1 + (zz * zz) * GP2);  in large_gam()
339 	z = w - (double) k *ln2_32;  in large_gam()
340 	zz = S[j] * (one + z + (z * z) * (Et1 + z * Et2));  in large_gam()
360 	double z;  in kpsin()  local
362 	z = x * x;  in kpsin()
363 	return (x + (x * z) * ((ks[0] + z * ks[1]) + (z * z) * (ks[2] + z *  in kpsin()
383 	double z;  in kpcos()  local
385 	z = x * x;  in kpcos()
386 	return (kc[0] + z * (kc[1] + z * kc[2] + (z * z) * kc[3]));  in kpcos()
455 	double x, y, z;  in tgammaf()  local
512 			z = -tiny;  in tgammaf()
514 			z = tiny;  in tgammaf()
515 		return ((float)z);  in tgammaf()
527 	z = y - (double) j;  in tgammaf()
528 	if (z > 0.3183098861837906715377675)  in tgammaf()
529 		if (z > 0.6816901138162093284622325)  in tgammaf()
530 			ss = kpsin(one - z);  in tgammaf()
532 			ss = kpcos(0.5 - z);  in tgammaf()
534 		ss = kpsin(z);  in tgammaf()
540 		ww = gam_n(j + 1, z);  in tgammaf()