Back to index

glibc  2.9
w_lgamma_r.c
Go to the documentation of this file.
00001 /* @(#)wr_lgamma.c 5.1 93/09/24 */
00002 /*
00003  * ====================================================
00004  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
00005  *
00006  * Developed at SunPro, a Sun Microsystems, Inc. business.
00007  * Permission to use, copy, modify, and distribute this
00008  * software is freely granted, provided that this notice
00009  * is preserved.
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  * wrapper double lgamma_r(double x, int *signgamp)
00019  */
00020 
00021 #include <math.h>
00022 #include "math_private.h"
00023 
00024 
00025 #ifdef __STDC__
00026        double __lgamma_r(double x, int *signgamp) /* wrapper lgamma_r */
00027 #else
00028        double __lgamma_r(x,signgamp)              /* wrapper lgamma_r */
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); /* lgamma pole */
00041             else
00042                 return __kernel_standard(x,x,14); /* lgamma overflow */
00043         } else
00044             return y;
00045 #endif
00046 }
00047 weak_alias (__lgamma_r, lgamma_r)
00048 #ifdef NO_LONG_DOUBLE
00049 strong_alias (__lgamma_r, __lgammal_r)
00050 weak_alias (__lgamma_r, lgammal_r)
00051 #endif