202   long double y, z;  in sin_pi()  local
211   y = -x;			/* x is assume negative */  in sin_pi()
217   z = floorl (y);  in sin_pi()
218   if (z != y)  in sin_pi()
220       y  *= 0.5L;  in sin_pi()
221       y = 2.0L*(y - floorl(y));		/* y = |x| mod 2.0 */  in sin_pi()
222       n = (int) (y*4.0L);  in sin_pi()
228 	  y = zero; n = 0;		/* y must be even */  in sin_pi()
233 	  z = y + two63;	/* exact */  in sin_pi()
236 	y  = n;  in sin_pi()
244       y = sinl (pi * y);  in sin_pi()
248       y = cosl (pi * (half - y));  in sin_pi()
252       y = sinl (pi * (one - y));  in sin_pi()
256       y = -cosl (pi * (y - 1.5L));  in sin_pi()
259       y = sinl (pi * (y - 2.0L));  in sin_pi()
262   return -y;  in sin_pi()
269   long double t, y, z, nadj = 0, p, p1, p2, q, r, w;  in lgammal_r()  local
324 	      y = x - one;  in lgammal_r()
329 	      y = x - (tc - one);  in lgammal_r()
335 	      y = x;  in lgammal_r()
345 	      y = x - 2.0L;  in lgammal_r()
351 	      y = x - tc;  in lgammal_r()
357 	      y = x - one;  in lgammal_r()
364 	  p1 = a0 + y * (a1 + y * (a2 + y * (a3 + y * (a4 + y * a5))));  in lgammal_r()
365 	  p2 = b0 + y * (b1 + y * (b2 + y * (b3 + y * (b4 + y))));  in lgammal_r()
366 	  r += half * y + y * p1/p2;  in lgammal_r()
369     p1 = g0 + y * (g1 + y * (g2 + y * (g3 + y * (g4 + y * (g5 + y * g6)))));  in lgammal_r()
370     p2 = h0 + y * (h1 + y * (h2 + y * (h3 + y * (h4 + y * (h5 + y)))));  in lgammal_r()
371     p = tt + y * p1/p2;  in lgammal_r()
375  p1 = y * (u0 + y * (u1 + y * (u2 + y * (u3 + y * (u4 + y * (u5 + y * u6))))));  in lgammal_r()
376       p2 = v0 + y * (v1 + y * (v2 + y * (v3 + y * (v4 + y * (v5 + y)))));  in lgammal_r()
377 	  r += (-half * y + p1 / p2);  in lgammal_r()
385       y = x - (long double) i;  in lgammal_r()
386   p = y *  in lgammal_r()
387      (s0 + y * (s1 + y * (s2 + y * (s3 + y * (s4 + y * (s5 + y * s6))))));  in lgammal_r()
388   q = r0 + y * (r1 + y * (r2 + y * (r3 + y * (r4 + y * (r5 + y * (r6 + y))))));  in lgammal_r()
389       r = half * y + p / q;  in lgammal_r()
394 	  z *= (y + 6.0L);  in lgammal_r()
397 	  z *= (y + 5.0L);  in lgammal_r()
400 	  z *= (y + 4.0L);  in lgammal_r()
403 	  z *= (y + 3.0L);  in lgammal_r()
406 	  z *= (y + 2.0L);  in lgammal_r()
416       y = z * z;  in lgammal_r()
418 	  + y * (w2 + y * (w3 + y * (w4 + y * (w5 + y * (w6 + y * w7))))));  in lgammal_r()