1151 	double t3, t4, y, z;  in GT1()  local
1155 	z = y * y;  in GT1()
1156 	t3 = (z * (P10 + y * ((P11 + y * P12) + z * (P13 + y * P14)))) /  in GT1()
1157 		(Q10 + y * ((Q11 + y * Q12) + z * ((Q13 + Q14 * y) + z * Q15)));  in GT1()
1170 	double t3, y, z;  in GT2()  local
1174 	z = y * y;  in GT2()
1175 	t3 = (z * (P20 + y * P21 + z * (P22 + y * P23))) /  in GT2()
1176 		(Q20 + (y * ((Q21 + Q22 * y) + z * Q23) +  in GT2()
1177 		(z * z) * ((Q24 + Q25 * y) + z * Q26))) + GZ2_l;  in GT2()
1187 	double t3, t4, y, z;  in GT3()  local
1191 	z = y * y;  in GT3()
1192 	t3 = (z * (P30 + y * ((P31 + y * P32) + z * (P33 + y * P34)))) /  in GT3()
1193 		(Q30 + y * ((Q31 + y * Q32) + z * ((Q33 + Q34 * y) + z * Q35)));  in GT3()
1211 	double z, t1, t2, t3, z2, t5, w, y, u, r, z4, v, t24 = 16777216.0,  in large_gam()  local
1250 	__HI(z) = (ix & 0xffffc000) | 0x2000;	/* z[j]=1+j/64+1/128 */  in large_gam()
1251 	__LO(z) = 0;  in large_gam()
1253 	t1 = y + z;  in large_gam()
1254 	t2 = y - z;  in large_gam()
1268 		t2 = t5 + r * ((two * t2 - v * t1) - v * (y - (t1 - z)));  in large_gam()
1278 	z = one / x;  in large_gam()
1282 	z2 = z * z;  in large_gam()
1289 	w_l = z * (GP0 + t1) + w;  in large_gam()
1305 	z = w_h - w_l;  in large_gam()
1306 	z2 = z * z;  in large_gam()
1309 	t3 = w_h - (w_l - (t1 + z * t2));  in large_gam()
1335 	double z, t1, t2, t3, t4;  in kpsin()  local
1338 	z = x * x;  in kpsin()
1340 	t1 = z * x;  in kpsin()
1341 	t2 = z * z;  in kpsin()
1343 	t3 = (t1 * z) * ((ks[1] + z * ks[2] + t2 * ks[3]) + (z * t2) *  in kpsin()
1344 		(ks[4] + z * ks[5] + t2 * ks[6]));  in kpsin()
1373 	double z, t1, t2, t3, t4, x4, x8;  in kpcos()  local
1376 	z = x * x;  in kpcos()
1379 	x4 = z * z;  in kpcos()
1380 	t2 = npi_2_l * z + npi_2_h * (x + t1) * (x - t1);  in kpcos()
1381 	t3 = one_pi_l + x4 * ((kc[1] + z * kc[2]) + x4 * (kc[3] + z *  in kpcos()
1408 	double r1, r2, t2, z, xh, xl, yh, yl, zh, z1, z2, zl, x5, wh, wl;  in gam_n()  local
1444 		z = x + one;	/* may not be exact */  in gam_n()
1445 		zh = (double) ((float) z);  in gam_n()
1447 		rr.l = z * yy.l + (x - (zh - one)) * yy.h;  in gam_n()
1452 		z = z1 * z2;  in gam_n()
1453 		xh = (double) ((float) z);  in gam_n()
1457 		rr.l = z * yy.l + xl * yy.h;  in gam_n()
1467 		z = z1 * z2;  in gam_n()
1468 		xh = (double) ((float) z);  in gam_n()
1471 		rr.l = z * yy.l + xl * yy.h;  in gam_n()
1476 		z = z1 * z2;  in gam_n()
1478 		yh = (double) ((float) z);  in gam_n()
1480 		z2 = z - 2.0;  in gam_n()
1481 		z *= z2;  in gam_n()
1482 		xh = (double) ((float) z);  in gam_n()
1485 		rr.l = z * yy.l + xl * yy.h;  in gam_n()
1490 		z = z1 * z2;  in gam_n()
1492 		yh = (double) ((float) z);  in gam_n()
1495 		z2 = z - 2.0;  in gam_n()
1497 		z *= z2;  in gam_n()
1498 		xh = (double) ((float) z);  in gam_n()
1506 		rr.l = z * wl + xl * wh;  in gam_n()
1511 		z = z2 * z1;  in gam_n()
1513 		yh = (double) ((float) z);	/* yh+yl = (x+3)(x+4) */  in gam_n()
1516 		z2 = z - 2.0;	/* z2 = (x+2)*(x+5) */  in gam_n()
1517 		z *= z2;  in gam_n()
1518 		xh = (double) ((float) z);  in gam_n()
1526 		rr.l = z * wl + xl * wh;  in gam_n()
1534 	double t, t1, t2, t3, t4, t5, w, y, z, z1, z2, z3, z5;  in tgamma()  local
1550 		z = x / tiny;  in tgamma()
1551 		return (z * z);  in tgamma()
1605 		z = tiny / x;  in tgamma()
1607 			z = -z;  in tgamma()
1608 		return (z * tiny);  in tgamma()
1619 	z = y - (double) j;  in tgamma()
1620 	if (z > 0.3183098861837906715377675)  in tgamma()
1621 		if (z > 0.6816901138162093284622325)  in tgamma()
1622 			ss = kpsin(one - z);  in tgamma()
1624 			ss = kpcos(0.5 - z);  in tgamma()
1626 		ss = kpsin(z);  in tgamma()
1635 		ww = gam_n(j + 1, z);  in tgamma()
1646 					ww = gam_n(j, z);  in tgamma()
1674 	z = t2 - z2;  in tgamma()
1678 		hx = __HI(z);  in tgamma()
1684 			__HI(z) = ix ^ i;  in tgamma()
1688 				z = -tiny * tiny;  in tgamma()
1690 				z = tiny * tiny;  in tgamma()
1697 			__HI(z) = ix ^ i;  in tgamma()
1698 			z *= t;  in tgamma()
1701 	return (z);  in tgamma()