gamma3
Log gamma functions, gamma function
Part of lgamma(3).
Usage
#include <math.h>
extern int
signgam;
double
lgamma(double x);
double
lgamma_r(double x, int *signgamp);
float
lgammaf(float x);
float
lgammaf_r(float x, int *signgamp);
long double
lgammal(long double x);
long double
lgammal_r(long double x, int *signgamp);
double
gamma(double x);
double
gamma_r(double x, int *signgamp);
float
gammaf(float x);
float
gammaf_r(float x, int *signgamp);
long double
tgamma(double x);
float
tgammaf(float x);
long double
tgammal(long double x);
Library
Math Library (libm, -lm)
Overview
lgamma(x), lgammaf(x), and lgammal(x) return ln|Γ(x)|. The external integer signgam returns the sign of Γ(x).
lgamma_r(x, signgamp), lgammaf_r(x, signgamp), and lgammal_r(x, signgamp) provide the same functionality as lgamma(x), lgammaf(x), and lgammal(x), but the caller must provide an integer to store the sign of Γ(x).
The tgamma(x), tgammaf(x), and tgammal(x) functions return Γ(x), with no effect on signgam.
gamma(), gammaf(), gamma_r(), and gammaf_r() are deprecated aliases for lgamma(), lgammaf(), lgamma_r(), and lgammaf_r(), respectively.
Idiosyncrasies
Do not use the expression “signgam∗exp(lgamma(x))” to compute g := Γ(x). Instead use a program like this (in C):
lg = lgamma(x); g = signgam∗exp(lg);
Only after lgamma() or lgammaf() has returned can signgam be correct.
For arguments in its range, tgamma() is preferred, as for positive arguments it is accurate to within one unit in the last place. Exponentiation of lgamma() will lose up to 10 significant bits.
Return values
gamma(), gammaf(), gammal(), gamma_r(), gammaf_r(), gammal_r(), lgamma(), lgammaf(), lgammal(), lgamma_r(), lgammaf_r(), and lgammal_r() return appropriate values unless an argument is out of range. Overflow will occur for sufficiently large positive values, and non-positive integers. For large non-integer negative values, tgamma() will underflow.
Releases
87 releases since FreeBSD 2.0.5-RELEASE; the text changed 49 times.