161     __float64 y, z;  in sin_pi()  local
169     y = -x; /* x is assume negative */  in sin_pi()
175     z = floor64(y);  in sin_pi()
176     if (z != y) { /* inexact anyway */  in sin_pi()
177         y *= _F_64(0.5);  in sin_pi()
178         y = _F_64(2.0) * (y - floor64(y)); /* y = |x| mod 2.0 */  in sin_pi()
179         n = (__int32_t)(y * _F_64(4.0));  in sin_pi()
182             y = zero;  in sin_pi()
186                 z = y + two52; /* exact */  in sin_pi()
189             y = n;  in sin_pi()
195         y = __kernel_sin(pi * y, zero, 0);  in sin_pi()
199         y = __kernel_cos(pi * (_F_64(0.5) - y), zero);  in sin_pi()
203         y = __kernel_sin(pi * (one - y), zero, 0);  in sin_pi()
207         y = -__kernel_cos(pi * (y - _F_64(1.5)), zero);  in sin_pi()
210         y = __kernel_sin(pi * (y - _F_64(2.0)), zero, 0);  in sin_pi()
213     return -y;  in sin_pi()
219     __float64 t, y, z, nadj = 0.0, p, p1, p2, p3, q, r, w;  in __math_lgamma_r()  local
266                 y = one - x;  in __math_lgamma_r()
269                 y = x - (tc - one);  in __math_lgamma_r()
272                 y = x;  in __math_lgamma_r()
278                 y = 2.0 - x;  in __math_lgamma_r()
282                 y = x - tc;  in __math_lgamma_r()
286                 y = x - one;  in __math_lgamma_r()
292             z = y * y;  in __math_lgamma_r()
295             p = y * p1 + p2;  in __math_lgamma_r()
296             r += (p - 0.5 * y);  in __math_lgamma_r()
299             z = y * y;  in __math_lgamma_r()
300             w = z * y;  in __math_lgamma_r()
305             p = z * p1 - (tt - w * (p2 + y * p3));  in __math_lgamma_r()
309             p1 = y * (u0 + y * (u1 + y * (u2 + y * (u3 + y * (u4 + y * u5)))));  in __math_lgamma_r()
310             p2 = one + y * (v1 + y * (v2 + y * (v3 + y * (v4 + y * v5))));  in __math_lgamma_r()
311             r += (-0.5 * y + p1 / p2);  in __math_lgamma_r()
315         y = x - (__float64)i;  in __math_lgamma_r()
316         p = y * (s0 +  in __math_lgamma_r()
317                  y * (s1 + y * (s2 + y * (s3 + y * (s4 + y * (s5 + y * s6))))));  in __math_lgamma_r()
318         q = one + y * (r1 + y * (r2 + y * (r3 + y * (r4 + y * (r5 + y * r6)))));  in __math_lgamma_r()
319         r = half * y + p / q;  in __math_lgamma_r()
323             z *= (y + 6.0);  in __math_lgamma_r()
326             z *= (y + 5.0);  in __math_lgamma_r()
329             z *= (y + 4.0);  in __math_lgamma_r()
332             z *= (y + 3.0);  in __math_lgamma_r()
335             z *= (y + 2.0);  in __math_lgamma_r()
343         y = z * z;  in __math_lgamma_r()
344         w = w0 + z * (w1 + y * (w2 + y * (w3 + y * (w4 + y * (w5 + y * w6)))));  in __math_lgamma_r()