1 /* @(#)w_gamma.c 5.1 93/09/24 */
3 * ====================================================
4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
6 * Developed at SunPro, a Sun Microsystems, Inc. business.
7 * Permission to use, copy, modify, and distribute this
8 * software is freely granted, provided that this notice
10 * ====================================================
13 #if defined(LIBM_SCCS) && !defined(lint)
14 static char rcsid
[] = "$NetBSD: w_gamma.c,v 1.7 1995/11/20 22:06:43 jtc Exp $";
17 /* double gamma(double x)
18 * Return the logarithm of the Gamma function of x or the Gamma function of x,
19 * depending on the library mode.
23 #include "math_private.h"
26 double __gamma(double x
)
34 if (_LIB_VERSION
== _SVID_
)
36 y
= __ieee754_lgamma_r(x
,&signgam
);
38 if(!__finite(y
)&&__finite(x
)) {
39 if(__floor(x
)==x
&&x
<=(double)0.0)
41 return __kernel_standard(x
,x
,15);
44 return __kernel_standard(x
,x
,14);
51 y
= __ieee754_gamma_r(x
,&local_signgam
);
52 if (local_signgam
< 0) y
= -y
;
56 if(_LIB_VERSION
== _IEEE_
) return y
;
58 if(!__finite(y
)&&__finite(x
)) {
59 if(__floor(x
)==x
&&x
<=0.0)
60 return __kernel_standard(x
,x
,41); /* gamma pole */
62 return __kernel_standard(x
,x
,40); /* gamma overflow */
68 weak_alias (__gamma
, gamma
)
70 strong_alias (__gamma
, __gammal
)
71 weak_alias (__gamma
, gammal
)