184 static long double stirf(long double x)  in stirf()  argument
188 w = 1.0L/x;  in stirf()
190 if( x > 1024.0L )  in stirf()
200 y = expl(x);  in stirf()
201 if( x > MAXSTIR )  in stirf()
203 	v = powl( x, 0.5L * x - 0.25L );  in stirf()
208 	y = powl( x, x - 0.5L ) / y;  in stirf()
215 tgammal(long double x)  in tgammal()  argument
220 if( isnan(x) )  in tgammal()
221         return(x + x);  in tgammal()
222 if(x == (long double) INFINITY)  in tgammal()
224 if(x == -(long double) INFINITY)  in tgammal()
225 	return __math_invalidl(x);  in tgammal()
226 if( x == 0.0L )  in tgammal()
227 	return __math_divzerol(signbitl_inline(x));  in tgammal()
228 q = fabsl(x);  in tgammal()
234                 if (x < 0.0L)  in tgammal()
235                         return __math_invalidl(x);  in tgammal()
238 	if( x < 0.0L )  in tgammal()
242 			return __math_invalidl(x);  in tgammal()
262 		z = stirf(x);  in tgammal()
268 while( x >= 3.0L )  in tgammal()
270 	x -= 1.0L;  in tgammal()
271 	z *= x;  in tgammal()
274 while( x < -0.03125L )  in tgammal()
276 	z /= x;  in tgammal()
277 	x += 1.0L;  in tgammal()
280 if( x <= 0.03125L )  in tgammal()
283 while( x < 2.0L )  in tgammal()
285 	z /= x;  in tgammal()
286 	x += 1.0L;  in tgammal()
289 if( x == 2.0L )  in tgammal()
292 x -= 2.0L;  in tgammal()
293 p = __polevll( x, P, 7 );  in tgammal()
294 q = __polevll( x, Q, 8 );  in tgammal()
299 if( x == 0.0L )  in tgammal()
300         return __math_invalidl(x);  in tgammal()
303         q = check_oflowl(1.0L/x);  in tgammal()
306                 if( x < 0.0L )  in tgammal()
308                         x = -x;  in tgammal()
309                         q = z / __polevll( x, SN, 8 ) * (-q);  in tgammal()
312                         q = z / __polevll( x, S, 8 ) * q;  in tgammal()