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 /* gamma_r(x, signgamp)
153a8617a8SJordan K. Hubbard * Reentrant version of the logarithm of the Gamma function
163a8617a8SJordan K. Hubbard * with user provide pointer for the sign of Gamma(x).
173a8617a8SJordan K. Hubbard *
18*99843eb8SSteve Kargl * Method: See lgamma_r
193a8617a8SJordan K. Hubbard */
203a8617a8SJordan K. Hubbard
213a8617a8SJordan K. Hubbard #include "math.h"
223a8617a8SJordan K. Hubbard #include "math_private.h"
233a8617a8SJordan K. Hubbard
24a82bbc73SAlfred Perlstein double
gamma_r(double x,int * signgamp)25*99843eb8SSteve Kargl gamma_r(double x, int *signgamp)
263a8617a8SJordan K. Hubbard {
27*99843eb8SSteve Kargl return lgamma_r(x,signgamp);
283a8617a8SJordan K. Hubbard }
29