diff options
author | Tulio Magno Quites Machado Filho <tuliom@linux.ibm.com> | 2018-05-29 17:52:24 -0300 |
---|---|---|
committer | Paul E. Murphy <murphyp@linux.vnet.ibm.com> | 2020-02-14 08:24:56 -0600 |
commit | c624d23260d5c136bebdd3a0734876527cedc505 (patch) | |
tree | c3b9039f6ef2e62892420a9b147668fc8924d759 /math/e_scalbf.c | |
parent | ad180676b83dc1782d407dbff57dabbaab0c1f71 (diff) | |
download | glibc-c624d23260d5c136bebdd3a0734876527cedc505.tar.gz glibc-c624d23260d5c136bebdd3a0734876527cedc505.tar.xz glibc-c624d23260d5c136bebdd3a0734876527cedc505.zip |
Add a generic scalb implementation
This is a preparatory patch to enable building a _Float128 variant to ease reuse when building a _Float128 variant to alias this long double only symbol. Notably, stubs are added where missing to the native _Float128 sysdep dir to prevent building these newly templated variants created inside the build directories. Also noteworthy are the changes around LIBM_SVID_COMPAT. These changes are not intuitive. The templated version is only enabled when !LIBM_SVID_COMPAT, and the compat version is predicated entirely on LIBM_SVID_COMPAT. Thus, exactly one is stubbed out entirely when building. The nldbl scalb compat files are updated to account for this. Likewise, fixup the reuse of m68k's e_scalb{f,l}.c to include it's override of e_scalb.c. Otherwise, the search path finds the templated copy in the build directory. This could be futher simplified by providing an overridden template, but I lack the hardware to verify.
Diffstat (limited to 'math/e_scalbf.c')
-rw-r--r-- | math/e_scalbf.c | 54 |
1 files changed, 0 insertions, 54 deletions
diff --git a/math/e_scalbf.c b/math/e_scalbf.c deleted file mode 100644 index 944cfbefc9..0000000000 --- a/math/e_scalbf.c +++ /dev/null @@ -1,54 +0,0 @@ -/* Copyright (C) 2011-2020 Free Software Foundation, Inc. - This file is part of the GNU C Library. - Contributed by Ulrich Drepper <drepper@gmail.com>, 2011. - - The GNU C Library is free software; you can redistribute it and/or - modify it under the terms of the GNU Lesser General Public - License as published by the Free Software Foundation; either - version 2.1 of the License, or (at your option) any later version. - - The GNU C Library is distributed in the hope that it will be useful, - but WITHOUT ANY WARRANTY; without even the implied warranty of - MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU - Lesser General Public License for more details. - - You should have received a copy of the GNU Lesser General Public - License along with the GNU C Library; if not, see - <https://www.gnu.org/licenses/>. */ - -#include <math.h> -#include <math_private.h> -#include <libm-alias-finite.h> - -static float -__attribute__ ((noinline)) -invalid_fn (float x, float fn) -{ - if (rintf (fn) != fn) - return (fn - fn) / (fn - fn); - else if (fn > 65000.0f) - return __scalbnf (x, 65000); - else - return __scalbnf (x,-65000); -} - - -float -__ieee754_scalbf (float x, float fn) -{ - if (__glibc_unlikely (isnan (x))) - return x * fn; - if (__glibc_unlikely (!isfinite (fn))) - { - if (isnan (fn) || fn > 0.0f) - return x * fn; - if (x == 0.0f) - return x; - return x / -fn; - } - if (__glibc_unlikely (fabsf (fn) >= 0x1p31f || (float) (int) fn != fn)) - return invalid_fn (x, fn); - - return __scalbnf (x, (int) fn); -} -libm_alias_finite (__ieee754_scalbf, __scalbf) |