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 }