85 	float y,z;  in sin_pif()  local
88 	y = -x;  in sin_pif()
90 	vz = y+0x1p23F;			/* depend on 0 <= y < 0x1p23 */  in sin_pif()
92 	if (z == y)  in sin_pif()
95 	vz = y+0x1p21F;  in sin_pif()
98 	if (z > y) {  in sin_pif()
103 	y = y - z + n * 0.25F;		/* y mod 2 */  in sin_pif()
106 	    case 0:   y =  __kernel_sindf(pi*y); break;  in sin_pif()
108 	    case 2:   y =  __kernel_cosdf(pi*((float)0.5-y)); break;  in sin_pif()
110 	    case 4:   y =  __kernel_sindf(pi*(one-y)); break;  in sin_pif()
112 	    case 6:   y = -__kernel_cosdf(pi*(y-(float)1.5)); break;  in sin_pif()
113 	    default:  y =  __kernel_sindf(pi*(y-(float)2.0)); break;  in sin_pif()
115 	return -y;  in sin_pif()
122 	float nadj,p,p1,p2,q,r,t,w,y,z;  in lgammaf_r()  local
159 		if(ix>=0x3f3b4a20) {y = one-x; i= 0;}  in lgammaf_r()
160 		else if(ix>=0x3e6d3308) {y= x-(tc-one); i=1;}  in lgammaf_r()
161 	  	else {y = x; i=2;}  in lgammaf_r()
164 	        if(ix>=0x3fdda618) {y=2-x;i=0;} /* [1.7316,2] */  in lgammaf_r()
165 	        else if(ix>=0x3F9da620) {y=x-tc;i=1;} /* [1.23,1.73] */  in lgammaf_r()
166 		else {y=x-one;i=2;}  in lgammaf_r()
170 		z = y*y;  in lgammaf_r()
173 		p  = y*p1+p2;  in lgammaf_r()
174 		r  += p-y/2; break;  in lgammaf_r()
176 		p = t0+y*t1+y*y*(t2+y*(t3+y*(t4+y*(t5+y*(t6+y*t7)))));  in lgammaf_r()
179 		p1 = y*(u0+y*(u1+y*u2));  in lgammaf_r()
180 		p2 = one+y*(v1+y*(v2+y*v3));  in lgammaf_r()
181 		r += p1/p2-y/2;  in lgammaf_r()
187 	    y = x-i;  in lgammaf_r()
188 	    p = y*(s0+y*(s1+y*(s2+y*s3)));  in lgammaf_r()
189 	    q = one+y*(r1+y*(r2+y*r3));  in lgammaf_r()
190 	    r = y/2+p/q;  in lgammaf_r()
193 	    case 7: z *= (y+6);		/* FALLTHRU */  in lgammaf_r()
194 	    case 6: z *= (y+5);		/* FALLTHRU */  in lgammaf_r()
195 	    case 5: z *= (y+4);		/* FALLTHRU */  in lgammaf_r()
196 	    case 4: z *= (y+3);		/* FALLTHRU */  in lgammaf_r()
197 	    case 3: z *= (y+2);		/* FALLTHRU */  in lgammaf_r()
204 	    y = z*z;  in lgammaf_r()
205 	    w = w0+z*(w1+y*w2);  in lgammaf_r()