Back to index

glibc  2.9
w_lgammal_r.c
Go to the documentation of this file.
00001 /* w_lgammal_r.c -- long double version of w_lgamma_r.c.
00002  * Conversion to long double by Ulrich Drepper,
00003  * Cygnus Support, drepper@cygnus.com.
00004  */
00005 
00006 /*
00007  * ====================================================
00008  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
00009  *
00010  * Developed at SunPro, a Sun Microsystems, Inc. business.
00011  * Permission to use, copy, modify, and distribute this
00012  * software is freely granted, provided that this notice
00013  * is preserved.
00014  * ====================================================
00015  */
00016 
00017 #if defined(LIBM_SCCS) && !defined(lint)
00018 static char rcsid[] = "$NetBSD: $";
00019 #endif
00020 
00021 /*
00022  * wrapper long double lgammal_r(long double x, int *signgamp)
00023  */
00024 
00025 #include <math.h>
00026 #include "math_private.h"
00027 
00028 
00029 #ifdef __STDC__
00030        long double __lgammal_r(long double x, int *signgamp)
00031               /* wrapper lgamma_r */
00032 #else
00033        long double __lgammal_r(x,signgamp)             /* wrapper lgamma_r */
00034         long double x; int *signgamp;
00035 #endif
00036 {
00037 #ifdef _IEEE_LIBM
00038        return __ieee754_lgammal_r(x,signgamp);
00039 #else
00040         long double y;
00041         y = __ieee754_lgammal_r(x,signgamp);
00042         if(_LIB_VERSION == _IEEE_) return y;
00043         if(!__finitel(y)&&__finitel(x)) {
00044             if(__floorl(x)==x&&x<=0.0)
00045                 return __kernel_standard(x,x,215); /* lgamma pole */
00046             else
00047                 return __kernel_standard(x,x,214); /* lgamma overflow */
00048         } else
00049             return y;
00050 #endif
00051 }
00052 weak_alias (__lgammal_r, lgammal_r)