161 	struct Double t, u, v;  local
172 	t.a = v.a*u.a;			/* t = (x-.5)*(log(x)-1) */
173 	t.b = v.b*u.a + x*u.b;
175 	t.b += lns2pi_lo; t.b += p;
176 	u.a = lns2pi_hi + t.b; u.a += t.a;
177 	u.b = t.a - u.a;
178 	u.b += lns2pi_hi; u.b += t.b;
189 	double y, ym1, t;  local
204 		t = r.a*yy.a;
206 		r.a = t;
208 		r.b += (t - r.a);
223 	double t, d;  local
226 		t = x, TRUNC(t);
227 		d = (t+x)*(x-t);
228 		t *= t;
229 		xx.a = (t + x), TRUNC(xx.a);
230 		xx.b = x - xx.a; xx.b += t; xx.b += d;
231 		t = (one-x0); t += x;
232 		d = (one-x0); d -= t; d += x;
237 		t = x - x0;
238 		d = (-x0 -t); d += x;
240 	r = ratfun_gam(t, d);
253 	struct Double r, t;  local
260 	t.a = z, TRUNC(t.a);		/* t ~= z + c */
261 	t.b = (z - t.a) + c;
262 	t.b *= (t.a + z);
263 	q = (t.a *= t.a);		/* t = (z+c)^2 */
264 	TRUNC(t.a);
265 	t.b += (q - t.a);
268 	t.b = t.b*p + t.a*r.b + a0_lo;
269 	t.a *= r.a;			/* t = (z+c)^2*(P/Q) */
270 	r.a = t.a + a0_hi, TRUNC(r.a);
271 	r.b = ((a0_hi-r.a) + t.a) + t.b;