204   long double y, z;  in sin_pi()  local
213   y = -x;			/* x is assume negative */  in sin_pi()
219   z = floorl (y);  in sin_pi()
220   if (z != y)  in sin_pi()
222       y  *= 0.5;  in sin_pi()
223       y = 2.0*(y - floorl(y));		/* y = |x| mod 2.0 */  in sin_pi()
224       n = (int) (y*4.0);  in sin_pi()
230 	  y = zero; n = 0;		/* y must be even */  in sin_pi()
235 	  z = y + two63;	/* exact */  in sin_pi()
238 	y  = n;  in sin_pi()
246       y = sinl (pi * y);  in sin_pi()
250       y = cosl (pi * (half - y));  in sin_pi()
254       y = sinl (pi * (one - y));  in sin_pi()
258       y = -cosl (pi * (y - 1.5));  in sin_pi()
261       y = sinl (pi * (y - 2.0));  in sin_pi()
264   return -y;  in sin_pi()
271   long double t, y, z, nadj, p, p1, p2, q, r, w;  in lgammal()  local
326 	      y = x - one;  in lgammal()
331 	      y = x - (tc - one);  in lgammal()
337 	      y = x;  in lgammal()
347 	      y = x - 2.0;  in lgammal()
353 	      y = x - tc;  in lgammal()
359 	      y = x - one;  in lgammal()
366 	  p1 = a0 + y * (a1 + y * (a2 + y * (a3 + y * (a4 + y * a5))));  in lgammal()
367 	  p2 = b0 + y * (b1 + y * (b2 + y * (b3 + y * (b4 + y))));  in lgammal()
368 	  r += half * y + y * p1/p2;  in lgammal()
371     p1 = g0 + y * (g1 + y * (g2 + y * (g3 + y * (g4 + y * (g5 + y * g6)))));  in lgammal()
372     p2 = h0 + y * (h1 + y * (h2 + y * (h3 + y * (h4 + y * (h5 + y)))));  in lgammal()
373     p = tt + y * p1/p2;  in lgammal()
377  p1 = y * (u0 + y * (u1 + y * (u2 + y * (u3 + y * (u4 + y * (u5 + y * u6))))));  in lgammal()
378       p2 = v0 + y * (v1 + y * (v2 + y * (v3 + y * (v4 + y * (v5 + y)))));  in lgammal()
379 	  r += (-half * y + p1 / p2);  in lgammal()
387       y = x - (double) i;  in lgammal()
388   p = y *  in lgammal()
389      (s0 + y * (s1 + y * (s2 + y * (s3 + y * (s4 + y * (s5 + y * s6))))));  in lgammal()
390   q = r0 + y * (r1 + y * (r2 + y * (r3 + y * (r4 + y * (r5 + y * (r6 + y))))));  in lgammal()
391       r = half * y + p / q;  in lgammal()
396 	  z *= (y + 6.0);	/* FALLTHRU */  in lgammal()
398 	  z *= (y + 5.0);	/* FALLTHRU */  in lgammal()
400 	  z *= (y + 4.0);	/* FALLTHRU */  in lgammal()
402 	  z *= (y + 3.0);	/* FALLTHRU */  in lgammal()
404 	  z *= (y + 2.0);	/* FALLTHRU */  in lgammal()
414       y = z * z;  in lgammal()
416 	  + y * (w2 + y * (w3 + y * (w4 + y * (w5 + y * (w6 + y * w7))))));  in lgammal()