diff options
Diffstat (limited to 'math/w_lgamma_r_compat.c')
| -rw-r--r-- | math/w_lgamma_r_compat.c | 38 |
1 files changed, 38 insertions, 0 deletions
diff --git a/math/w_lgamma_r_compat.c b/math/w_lgamma_r_compat.c new file mode 100644 index 0000000000..6578b2c76c --- /dev/null +++ b/math/w_lgamma_r_compat.c @@ -0,0 +1,38 @@ +/* @(#)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. + * ==================================================== + */ + +/* + * wrapper double lgamma_r(double x, int *signgamp) + */ + +#include <math.h> +#include <math_private.h> + + +double +__lgamma_r(double x, int *signgamp) +{ + double y = __ieee754_lgamma_r(x,signgamp); + if(__builtin_expect(!isfinite(y), 0) + && isfinite(x) && _LIB_VERSION != _IEEE_) + return __kernel_standard(x, x, + __floor(x)==x&&x<=0.0 + ? 15 /* lgamma pole */ + : 14); /* lgamma overflow */ + + return y; +} +weak_alias (__lgamma_r, lgamma_r) +#ifdef NO_LONG_DOUBLE +strong_alias (__lgamma_r, __lgammal_r) +weak_alias (__lgamma_r, lgammal_r) +#endif |