00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
#if defined(LIBM_SCCS) && !defined(lint)
00014
static char rcsid[] =
"$NetBSD: w_lgamma_r.c,v 1.6 1995/05/10 20:49:27 jtc Exp $";
00015
#endif
00016
00017
00018
00019
00020
00021
#include "math.h"
00022
#include "mathP.h"
00023
00024
00025
#ifdef __STDC__
00026
double lgamma_r(
double x,
int *signgamp)
00027 #
else
00028 double lgamma_r(x,signgamp)
00029 double x;
int *signgamp;
00030 #endif
00031 {
00032
#ifdef _IEEE_LIBM
00033
return __ieee754_lgamma_r(x,signgamp);
00034
#else
00035
double y;
00036 y =
__ieee754_lgamma_r(x,signgamp);
00037
if(
_LIB_VERSION == _IEEE_)
return y;
00038
if(!
finite(y)&&
finite(x)) {
00039
if(
floor(x)==x&&x<=0.0)
00040
return __kernel_standard(x,x,15);
00041
else
00042
return __kernel_standard(x,x,14);
00043 }
else
00044
return y;
00045
#endif
00046
}