1/* w_gammaf.c -- float version of w_gamma.c.
2 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
3 */
4
5/*
6 * ====================================================
7 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8 *
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
12 * is preserved.
13 * ====================================================
14 */
15
16#include <sys/cdefs.h>
17#if defined(LIBM_SCCS) && !defined(lint)
18__RCSID("$NetBSD: w_gammaf.c,v 1.7 2001/01/06 00:15:00 christos Exp $");
19#endif
20
21#include "math.h"
22#include "math_private.h"
23
24float
25gammaf(float x)
26{
27#ifdef _IEEE_LIBM
28	return __ieee754_lgammaf_r(x,&signgam);
29#else
30        float y;
31        y = __ieee754_lgammaf_r(x,&signgam);
32        if(_LIB_VERSION == _IEEE_) return y;
33        if(!finitef(y)&&finitef(x)) {
34            if(floorf(x)==x&&x<=(float)0.0)
35	        /* gammaf pole */
36                return (float)__kernel_standard((double)x,(double)x,141);
37            else
38	        /* gammaf overflow */
39                return (float)__kernel_standard((double)x,(double)x,140);
40        } else
41            return y;
42#endif
43}
44