Return to w_lgamma_r.c CVS log | Up to [cvs.NetBSD.org] / src / lib / libm / src |
File: [cvs.NetBSD.org] / src / lib / libm / src / w_lgamma_r.c (download)
Revision 1.8, Fri Jul 2 15:37:44 1999 UTC (24 years, 9 months ago) by simonb
More trailing white space. |
/* @(#)wr_lgamma.c 5.1 93/09/24 */ /* * ==================================================== * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. * * Developed at SunPro, a Sun Microsystems, Inc. business. * Permission to use, copy, modify, and distribute this * software is freely granted, provided that this notice * is preserved. * ==================================================== */ #include <sys/cdefs.h> #if defined(LIBM_SCCS) && !defined(lint) __RCSID("$NetBSD: w_lgamma_r.c,v 1.8 1999/07/02 15:37:44 simonb Exp $"); #endif /* * wrapper double lgamma_r(double x, int *signgamp) */ #include "math.h" #include "math_private.h" double lgamma_r __P((double, int *)); #ifdef __STDC__ double lgamma_r(double x, int *signgamp) /* wrapper lgamma_r */ #else double lgamma_r(x,signgamp) /* wrapper lgamma_r */ double x; int *signgamp; #endif { #ifdef _IEEE_LIBM return __ieee754_lgamma_r(x,signgamp); #else double y; y = __ieee754_lgamma_r(x,signgamp); if(_LIB_VERSION == _IEEE_) return y; if(!finite(y)&&finite(x)) { if(floor(x)==x&&x<=0.0) return __kernel_standard(x,x,15); /* lgamma pole */ else return __kernel_standard(x,x,14); /* lgamma overflow */ } else return y; #endif }