161 	double y,z;  in sin_pi()  local
168 	y = -x;		/* x is assume negative */  in sin_pi()
174 	z = floor(y);  in sin_pi()
175 	if(z!=y) {				/* inexact anyway */  in sin_pi()
176 	    y  *= 0.5;  in sin_pi()
177 	    y   = 2.0*(y - floor(y));		/* y = |x| mod 2.0 */  in sin_pi()
178 	    n   = (int) (y*4.0);  in sin_pi()
181                 y = zero; n = 0;                 /* y must be even */  in sin_pi()
183                 if(ix<0x43300000) z = y+two52;	/* exact */  in sin_pi()
186                 y  = n;  in sin_pi()
191 	    case 0:   y =  __kernel_sin(pi*y,zero,0); break;  in sin_pi()
193 	    case 2:   y =  __kernel_cos(pi*(0.5-y),zero); break;  in sin_pi()
195 	    case 4:   y =  __kernel_sin(pi*(one-y),zero,0); break;  in sin_pi()
197 	    case 6:   y = -__kernel_cos(pi*(y-1.5),zero); break;  in sin_pi()
198 	    default:  y =  __kernel_sin(pi*(y-2.0),zero,0); break;  in sin_pi()
200 	return -y;  in sin_pi()
207 	double t,y,z,nadj,p,p1,p2,p3,q,r,w;  in __ieee754_lgamma_r()  local
239 		if(ix>=0x3FE76944) {y = one-x; i= 0;}  in __ieee754_lgamma_r()
240 		else if(ix>=0x3FCDA661) {y= x-(tc-one); i=1;}  in __ieee754_lgamma_r()
241 	  	else {y = x; i=2;}  in __ieee754_lgamma_r()
244 	        if(ix>=0x3FFBB4C3) {y=2.0-x;i=0;} /* [1.7316,2] */  in __ieee754_lgamma_r()
245 	        else if(ix>=0x3FF3B4C4) {y=x-tc;i=1;} /* [1.23,1.73] */  in __ieee754_lgamma_r()
246 		else {y=x-one;i=2;}  in __ieee754_lgamma_r()
250 		z = y*y;  in __ieee754_lgamma_r()
253 		p  = y*p1+p2;  in __ieee754_lgamma_r()
254 		r  += (p-0.5*y); break;  in __ieee754_lgamma_r()
256 		z = y*y;  in __ieee754_lgamma_r()
257 		w = z*y;  in __ieee754_lgamma_r()
261 		p  = z*p1-(tt-w*(p2+y*p3));  in __ieee754_lgamma_r()
264 		p1 = y*(u0+y*(u1+y*(u2+y*(u3+y*(u4+y*u5)))));  in __ieee754_lgamma_r()
265 		p2 = one+y*(v1+y*(v2+y*(v3+y*(v4+y*v5))));  in __ieee754_lgamma_r()
266 		r += (-0.5*y + p1/p2);  in __ieee754_lgamma_r()
272 	    y = x-(double)i;  in __ieee754_lgamma_r()
273 	    p = y*(s0+y*(s1+y*(s2+y*(s3+y*(s4+y*(s5+y*s6))))));  in __ieee754_lgamma_r()
274 	    q = one+y*(r1+y*(r2+y*(r3+y*(r4+y*(r5+y*r6)))));  in __ieee754_lgamma_r()
275 	    r = half*y+p/q;  in __ieee754_lgamma_r()
278 	    case 7: z *= (y+6.0);	/* FALLTHRU */  in __ieee754_lgamma_r()
279 	    case 6: z *= (y+5.0);	/* FALLTHRU */  in __ieee754_lgamma_r()
280 	    case 5: z *= (y+4.0);	/* FALLTHRU */  in __ieee754_lgamma_r()
281 	    case 4: z *= (y+3.0);	/* FALLTHRU */  in __ieee754_lgamma_r()
282 	    case 3: z *= (y+2.0);	/* FALLTHRU */  in __ieee754_lgamma_r()
289 	    y = z*z;  in __ieee754_lgamma_r()
290 	    w = w0+z*(w1+y*(w2+y*(w3+y*(w4+y*(w5+y*w6)))));  in __ieee754_lgamma_r()