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 (!__finitef (y) || y == 0)
27 && (__finitef (x) || __isinff (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)