]>
git.ipfire.org Git - thirdparty/glibc.git/blob - sysdeps/libm-ieee754/w_gammaf.c
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 * ====================================================
16 #if defined(LIBM_SCCS) && !defined(lint)
17 static char rcsid
[] = "$NetBSD: w_gammaf.c,v 1.4 1995/11/20 22:06:48 jtc Exp $";
21 #include "math_private.h"
24 float __gammaf(float x
)
32 if (_LIB_VERSION
== _SVID_
)
33 y
= __ieee754_lgammaf_r(x
,&signgam
);
36 y
= __ieee754_gammaf_r(x
,&signgam
);
37 if (signgam
< 0) y
= -y
;
41 if(_LIB_VERSION
== _IEEE_
) return y
;
43 if(!__finitef(y
)&&__finitef(x
)) {
44 if(__floorf(x
)==x
&&x
<=(float)0.0)
46 return (float)__kernel_standard((double)x
,(double)x
,141);
49 return (float)__kernel_standard((double)x
,(double)x
,140);
54 weak_alias (__gammaf
, gammaf
)