120 	double p, z, thi, tlo, xhi, xlo;  in large_gam()  local
138 	tlo = xlo * u.a + x * u.b;  in large_gam()
141 	tlo += ln2pi_lo;  in large_gam()
142 	tlo += p;  in large_gam()
143 	u.a = ln2pi_hi + tlo;  in large_gam()
147 	u.b += tlo;  in large_gam()
182 	double p, q, thi, tlo;  in ratfun_gam()  local
192 	tlo = (z - thi) + c;  in ratfun_gam()
193 	tlo *= (thi + z);  in ratfun_gam()
199 	tlo += (q - thi);  in ratfun_gam()
205 	tlo = tlo * p + thi * r.b + a0_lo;  in ratfun_gam()
208 	r.b = ((a0_hi - r.a) + thi) + tlo;  in ratfun_gam()