163 	double y,z;  in sin_pi()  local
166 	y = -x;  in sin_pi()
168 	vz = y+0x1p52;			/* depend on 0 <= y < 0x1p52 */  in sin_pi()
169 	z = vz-0x1p52;			/* rint(y) for the above range */  in sin_pi()
170 	if (z == y)  in sin_pi()
173 	vz = y+0x1p50;  in sin_pi()
174 	GET_LOW_WORD(n,vz);		/* bits for rounded y (units 0.25) */  in sin_pi()
175 	z = vz-0x1p50;			/* y rounded to a multiple of 0.25 */  in sin_pi()
176 	if (z > y) {  in sin_pi()
180 	n &= 7;				/* octant of y mod 2 */  in sin_pi()
181 	y = y - z + n * 0.25;		/* y mod 2 */  in sin_pi()
184 	    case 0:   y =  __kernel_sin(pi*y,zero,0); break;  in sin_pi()
186 	    case 2:   y =  __kernel_cos(pi*(0.5-y),zero); break;  in sin_pi()
188 	    case 4:   y =  __kernel_sin(pi*(one-y),zero,0); break;  in sin_pi()
190 	    case 6:   y = -__kernel_cos(pi*(y-1.5),zero); break;  in sin_pi()
191 	    default:  y =  __kernel_sin(pi*(y-2.0),zero,0); break;  in sin_pi()
193 	return -y;  in sin_pi()
200 	double nadj,p,p1,p2,p3,q,r,t,w,y,z;  in lgamma_r()  local
237 		if(ix>=0x3FE76944) {y = one-x; i= 0;}  in lgamma_r()
238 		else if(ix>=0x3FCDA661) {y= x-(tc-one); i=1;}  in lgamma_r()
239 	  	else {y = x; i=2;}  in lgamma_r()
242 	        if(ix>=0x3FFBB4C3) {y=2.0-x;i=0;} /* [1.7316,2] */  in lgamma_r()
243 	        else if(ix>=0x3FF3B4C4) {y=x-tc;i=1;} /* [1.23,1.73] */  in lgamma_r()
244 		else {y=x-one;i=2;}  in lgamma_r()
248 		z = y*y;  in lgamma_r()
251 		p  = y*p1+p2;  in lgamma_r()
252 		r  += p-y/2; break;  in lgamma_r()
254 		z = y*y;  in lgamma_r()
255 		w = z*y;  in lgamma_r()
259 		p  = z*p1-(tt-w*(p2+y*p3));  in lgamma_r()
262 		p1 = y*(u0+y*(u1+y*(u2+y*(u3+y*(u4+y*u5)))));  in lgamma_r()
263 		p2 = one+y*(v1+y*(v2+y*(v3+y*(v4+y*v5))));  in lgamma_r()
264 		r += p1/p2-y/2;  in lgamma_r()
270 	    y = x-i;  in lgamma_r()
271 	    p = y*(s0+y*(s1+y*(s2+y*(s3+y*(s4+y*(s5+y*s6))))));  in lgamma_r()
272 	    q = one+y*(r1+y*(r2+y*(r3+y*(r4+y*(r5+y*r6)))));  in lgamma_r()
273 	    r = y/2+p/q;  in lgamma_r()
276 	    case 7: z *= (y+6);		/* FALLTHRU */  in lgamma_r()
277 	    case 6: z *= (y+5);		/* FALLTHRU */  in lgamma_r()
278 	    case 5: z *= (y+4);		/* FALLTHRU */  in lgamma_r()
279 	    case 4: z *= (y+3);		/* FALLTHRU */  in lgamma_r()
280 	    case 3: z *= (y+2);		/* FALLTHRU */  in lgamma_r()
287 	    y = z*z;  in lgamma_r()
288 	    w = w0+z*(w1+y*(w2+y*(w3+y*(w4+y*(w5+y*w6)))));  in lgamma_r()