From 0ccf5430f2abac4a7d4242dfb309c5fb27b4f79d Mon Sep 17 00:00:00 2001 From: Christopher Faylor Date: Wed, 11 Feb 2004 00:47:52 +0000 Subject: * libm/mathfp/er_gamma.c (gamma): Add new non-reentrant function. * libm/mathfp/er_lgamma.c (lgamma): Ditto. * libm/mathfp/erf_gamma.c (gammaf): Ditto. * libm/mathfp/erf_lgamma.c (lgammaf): Ditto. --- newlib/libm/mathfp/er_gamma.c | 5 +++++ newlib/libm/mathfp/er_lgamma.c | 6 ++++++ newlib/libm/mathfp/erf_gamma.c | 7 ++++++- newlib/libm/mathfp/erf_lgamma.c | 5 +++++ 4 files changed, 22 insertions(+), 1 deletion(-) (limited to 'newlib/libm') diff --git a/newlib/libm/mathfp/er_gamma.c b/newlib/libm/mathfp/er_gamma.c index 63758bb43..0060ed031 100644 --- a/newlib/libm/mathfp/er_gamma.c +++ b/newlib/libm/mathfp/er_gamma.c @@ -30,3 +30,8 @@ { return exp (lgamma_r(x,signgamp)); } + +double gamma(double x) +{ + return gamma_r(x, &(_REENT_SIGNGAM(_REENT))); +} diff --git a/newlib/libm/mathfp/er_lgamma.c b/newlib/libm/mathfp/er_lgamma.c index 65d60026e..fd49c252e 100644 --- a/newlib/libm/mathfp/er_lgamma.c +++ b/newlib/libm/mathfp/er_lgamma.c @@ -422,3 +422,9 @@ static double zero= 0.00000000000000000000e+00; if(hx<0) r = nadj - r; return r; } + +double +lgamma(double x) +{ + return lgamma_r(x, &(_REENT_SIGNGAM(_REENT))); +} diff --git a/newlib/libm/mathfp/erf_gamma.c b/newlib/libm/mathfp/erf_gamma.c index 16be0d040..246645d25 100644 --- a/newlib/libm/mathfp/erf_gamma.c +++ b/newlib/libm/mathfp/erf_gamma.c @@ -16,7 +16,7 @@ /* gammaf_r(x, signgamp) * Reentrant version of the logarithm of the Gamma function - * with user provide pointer for the sign of Gamma(x). + * with user provided pointer for the sign of Gamma(x). * * Method: See lgammaf_r */ @@ -32,3 +32,8 @@ { return expf (lgammaf_r(x,signgamp)); } + +float gammaf(float x) +{ + return gammaf_r(x, &(_REENT_SIGNGAM(_REENT))); +} diff --git a/newlib/libm/mathfp/erf_lgamma.c b/newlib/libm/mathfp/erf_lgamma.c index 64bec0aca..4f092147b 100644 --- a/newlib/libm/mathfp/erf_lgamma.c +++ b/newlib/libm/mathfp/erf_lgamma.c @@ -244,3 +244,8 @@ static float zero= 0.0000000000e+00; if(hx<0) r = nadj - r; return r; } + +float lgammaf(float x) +{ + return lgammaf_r(x, &(_REENT_SIGNGAM(_REENT))); +} -- cgit v1.2.3