13f708241SDavid Schultz 23a8617a8SJordan K. Hubbard /* 33a8617a8SJordan K. Hubbard * ==================================================== 43a8617a8SJordan K. Hubbard * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 53a8617a8SJordan K. Hubbard * 63f708241SDavid Schultz * Developed at SunSoft, a Sun Microsystems, Inc. business. 73a8617a8SJordan K. Hubbard * Permission to use, copy, modify, and distribute this 83a8617a8SJordan K. Hubbard * software is freely granted, provided that this notice 93a8617a8SJordan K. Hubbard * is preserved. 103a8617a8SJordan K. Hubbard * ==================================================== 113f708241SDavid Schultz * 123a8617a8SJordan K. Hubbard */ 133a8617a8SJordan K. Hubbard 14*99843eb8SSteve Kargl /* lgamma(x) 153a8617a8SJordan K. Hubbard * Return the logarithm of the Gamma function of x. 163a8617a8SJordan K. Hubbard * 17*99843eb8SSteve Kargl * Method: call lgamma_r 183a8617a8SJordan K. Hubbard */ 193a8617a8SJordan K. Hubbard 20f7efd14dSSteve Kargl #include <float.h> 21f7efd14dSSteve Kargl 223a8617a8SJordan K. Hubbard #include "math.h" 233a8617a8SJordan K. Hubbard #include "math_private.h" 243a8617a8SJordan K. Hubbard 253a8617a8SJordan K. Hubbard extern int signgam; 263a8617a8SJordan K. Hubbard 2759b19ff1SAlfred Perlstein double lgamma(double x)28*99843eb8SSteve Kargllgamma(double x) 293a8617a8SJordan K. Hubbard { 30*99843eb8SSteve Kargl return lgamma_r(x,&signgam); 313a8617a8SJordan K. Hubbard } 32f7efd14dSSteve Kargl 33f7efd14dSSteve Kargl #if (LDBL_MANT_DIG == 53) 34f7efd14dSSteve Kargl __weak_reference(lgamma, lgammal); 35f7efd14dSSteve Kargl #endif 36