diff options
Diffstat (limited to 'sysdeps/ieee754/dbl-64/lgamma_neg.c')
| -rw-r--r-- | sysdeps/ieee754/dbl-64/lgamma_neg.c | 9 |
1 files changed, 5 insertions, 4 deletions
diff --git a/sysdeps/ieee754/dbl-64/lgamma_neg.c b/sysdeps/ieee754/dbl-64/lgamma_neg.c index 5bb2f10c71..84e7022af5 100644 --- a/sysdeps/ieee754/dbl-64/lgamma_neg.c +++ b/sysdeps/ieee754/dbl-64/lgamma_neg.c @@ -1,5 +1,5 @@ /* lgamma expanding around zeros. - Copyright (C) 2015-2018 Free Software Foundation, Inc. + Copyright (C) 2015-2019 Free Software Foundation, Inc. This file is part of the GNU C Library. The GNU C Library is free software; you can redistribute it and/or @@ -14,12 +14,13 @@ You should have received a copy of the GNU Lesser General Public License along with the GNU C Library; if not, see - <http://www.gnu.org/licenses/>. */ + <https://www.gnu.org/licenses/>. */ #include <float.h> #include <math.h> #include <math-narrow-eval.h> #include <math_private.h> +#include <fenv_private.h> static const double lgamma_zeros[][2] = { @@ -281,7 +282,7 @@ __lgamma_neg (double x, int *signgamp) { /* Determine the half-integer region X lies in, handle exact integers and determine the sign of the result. */ - int i = __floor (-2 * x); + int i = floor (-2 * x); if ((i & 1) == 0 && i == -2 * x) return 1.0 / 0.0; double xn = ((i & 1) == 0 ? -i / 2 : (-i - 1) / 2); @@ -298,7 +299,7 @@ __lgamma_neg (double x, int *signgamp) approximations to an adjusted version of the gamma function. */ if (i < 2) { - int j = __floor (-8 * x) - 16; + int j = floor (-8 * x) - 16; double xm = (-33 - 2 * j) * 0.0625; double x_adj = x - xm; size_t deg = poly_deg[j]; |