diff options
author | Denis Vlasenko <vda.linux@googlemail.com> | 2008-12-29 16:00:21 +0000 |
---|---|---|
committer | Denis Vlasenko <vda.linux@googlemail.com> | 2008-12-29 16:00:21 +0000 |
commit | d4f94ef391e602b6c2df453283ab2aeca1905eb2 (patch) | |
tree | 9d5924cfddc21406325bbe7aa769c5c7c2cfcedc /libm/w_gamma.c | |
parent | da3fc969d3d4e34ede64c347ac2399684c7cbcfe (diff) |
libm/*: delete many incarnations of logarithmic gamma function,
aliasing them instead af appropriate. Also alias drem to remainder.
text data bss dec hex filename
- 43568 176 4 43748 aae4 lib/libm.so
+ 43444 176 4 43624 aa68 lib/libm.so
Diffstat (limited to 'libm/w_gamma.c')
-rw-r--r-- | libm/w_gamma.c | 36 |
1 files changed, 0 insertions, 36 deletions
diff --git a/libm/w_gamma.c b/libm/w_gamma.c deleted file mode 100644 index 246c8f764..000000000 --- a/libm/w_gamma.c +++ /dev/null @@ -1,36 +0,0 @@ -/* - * ==================================================== - * 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. - * ==================================================== - */ - -/* double gamma(double x) - * Return the logarithm of the Gamma function of x. - * - * Method: call gamma_r - */ - -#include <math.h> -#include "math_private.h" - -double gamma(double x) -{ -#ifdef _IEEE_LIBM - return __ieee754_lgamma_r(x, &signgam); -#else - double y = __ieee754_lgamma_r(x, &signgam); - if (_LIB_VERSION == _IEEE_) - return y; - if (!isfinite(y) && isfinite(x)) { - if (floor(x) == x && x <= 0.0) - return __kernel_standard(x, x, 41); /* gamma pole */ - return __kernel_standard(x, x, 40); /* gamma overflow */ - } - return y; -#endif -} |