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>
19 #include <math-svid-compat.h>
25 float y
= __ieee754_gammaf_r(x
,&local_signgam
);
27 if(__glibc_unlikely (!isfinite (y
) || y
== 0)
28 && (isfinite (x
) || (isinf (x
) && x
< 0.0))
29 && _LIB_VERSION
!= _IEEE_
) {
32 return __kernel_standard_f(x
, x
, 150);
33 else if(__floorf(x
)==x
&&x
<0.0f
)
35 return __kernel_standard_f(x
, x
, 141);
37 /* tgammaf underflow */
40 /* tgammaf overflow */
41 return __kernel_standard_f(x
, x
, 140);
43 return local_signgam
< 0 ? - y
: y
;
45 weak_alias (__tgammaf
, tgammaf
)