diff options
Diffstat (limited to 'sysdeps/ieee754/flt-32/lgamma_negf.c')
| -rw-r--r-- | sysdeps/ieee754/flt-32/lgamma_negf.c | 9 |
1 files changed, 5 insertions, 4 deletions
diff --git a/sysdeps/ieee754/flt-32/lgamma_negf.c b/sysdeps/ieee754/flt-32/lgamma_negf.c index 01edb0b8de..c2c00eb5b2 100644 --- a/sysdeps/ieee754/flt-32/lgamma_negf.c +++ b/sysdeps/ieee754/flt-32/lgamma_negf.c @@ -1,5 +1,5 @@ /* lgammaf 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 float lgamma_zeros[][2] = { @@ -196,7 +197,7 @@ __lgamma_negf (float x, int *signgamp) { /* Determine the half-integer region X lies in, handle exact integers and determine the sign of the result. */ - int i = __floorf (-2 * x); + int i = floorf (-2 * x); if ((i & 1) == 0 && i == -2 * x) return 1.0f / 0.0f; float xn = ((i & 1) == 0 ? -i / 2 : (-i - 1) / 2); @@ -213,7 +214,7 @@ __lgamma_negf (float x, int *signgamp) approximations to an adjusted version of the gamma function. */ if (i < 2) { - int j = __floorf (-8 * x) - 16; + int j = floorf (-8 * x) - 16; float xm = (-33 - 2 * j) * 0.0625f; float x_adj = x - xm; size_t deg = poly_deg[j]; |