Functions
mirkl_error_t *	mirml_lngamma_f (double x, double res, double err, double *sgn)

double	mirml_lngamma (double x)

mirkl_error_t *	mirml_gamma_f (double x, double res, double err)

mirkl_error_t *	mirml_lnfact_f (unsigned int n, double res, double err)

double	mirml_lngamma_stirling (double x)

Detailed Description

Macro Definition Documentation

#define MIRML_FACT_NMAX 170

Maximum n for which n! does not overflow

#define MIRML_GAMMA_XMAX 171.0

The maximum x such that gamma(x) does not overflow.

double mirml_lngamma_stirling ( double x )

Return natural log of gamma function.

Parameters

[in] x Argument, must be > 0.

Note: Computes natural logarithm of Gamma(x) using an asymptotic Stirling's expansion, see Abramowitz & Stegun, Handbook of Mathematical Functions, 1964 [Eq. 6.1.41]