1*05a0b428SJohn Marino /* @(#)w_lgamma.c 5.1 93/09/24 */
2*05a0b428SJohn Marino /*
3*05a0b428SJohn Marino  * ====================================================
4*05a0b428SJohn Marino  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
5*05a0b428SJohn Marino  *
6*05a0b428SJohn Marino  * Developed at SunPro, a Sun Microsystems, Inc. business.
7*05a0b428SJohn Marino  * Permission to use, copy, modify, and distribute this
8*05a0b428SJohn Marino  * software is freely granted, provided that this notice
9*05a0b428SJohn Marino  * is preserved.
10*05a0b428SJohn Marino  * ====================================================
11*05a0b428SJohn Marino  */
12*05a0b428SJohn Marino 
13*05a0b428SJohn Marino /* double lgamma(double x)
14*05a0b428SJohn Marino  * Return the logarithm of the Gamma function of x.
15*05a0b428SJohn Marino  *
16*05a0b428SJohn Marino  * Method: call lgamma_r
17*05a0b428SJohn Marino  */
18*05a0b428SJohn Marino 
19*05a0b428SJohn Marino #include <float.h>
20*05a0b428SJohn Marino #include <math.h>
21*05a0b428SJohn Marino 
22*05a0b428SJohn Marino #include "math_private.h"
23*05a0b428SJohn Marino 
24*05a0b428SJohn Marino double
lgamma(double x)25*05a0b428SJohn Marino lgamma(double x)
26*05a0b428SJohn Marino {
27*05a0b428SJohn Marino 	return lgamma_r(x,&signgam);
28*05a0b428SJohn Marino }
29*05a0b428SJohn Marino 
30*05a0b428SJohn Marino #if	LDBL_MANT_DIG == DBL_MANT_DIG
31*05a0b428SJohn Marino __strong_alias(lgammal, lgamma);
32*05a0b428SJohn Marino #endif	/* LDBL_MANT_DIG == DBL_MANT_DIG */
33