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 Marinolgamma(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