1 /* w_gammaf.c -- float version of w_gamma.c.
2 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
6 * ====================================================
7 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
9 * Developed at SunPro, a Sun Microsystems, Inc. business.
10 * Permission to use, copy, modify, and distribute this
11 * software is freely granted, provided that this notice
13 * ====================================================
18 #include <math_private.h>
24 float y
= __ieee754_gammaf_r(x
,&local_signgam
);
26 if(__glibc_unlikely (!isfinite (y
) || y
== 0)
27 && (isfinite (x
) || isinf (x
) < 0)
28 && _LIB_VERSION
!= _IEEE_
) {
31 return __kernel_standard_f(x
, x
, 150);
32 else if(__floorf(x
)==x
&&x
<0.0f
)
34 return __kernel_standard_f(x
, x
, 141);
36 /* tgammaf underflow */
39 /* tgammaf overflow */
40 return __kernel_standard_f(x
, x
, 140);
42 return local_signgam
< 0 ? - y
: y
;
44 weak_alias (__tgammaf
, tgammaf
)