picnoir
/
glibc
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Add narrowing add functions. 

This patch adds the narrowing add functions from TS 18661-1 to glibc's
libm: fadd, faddl, daddl, f32addf64, f32addf32x, f32xaddf64 for all
configurations; f32addf64x, f32addf128, f64addf64x, f64addf128,
f32xaddf64x, f32xaddf128, f64xaddf128 for configurations with
_Float64x and _Float128; __nldbl_daddl for ldbl-opt.  As discussed for
the build infrastructure patch, tgmath.h support is deliberately
deferred, and FP_FAST_* macros are not applicable without optimized
function implementations.

Function implementations are added for all relevant pairs of formats
(including certain cases of a format and itself where more than one
type has that format).  The main implementations use round-to-odd, or
a trivial computation in the case where both formats are the same or
where the wider format is IBM long double (in which case we don't
attempt to be correctly rounding).  The sysdeps/ieee754/soft-fp
implementations use soft-fp, and are used automatically for
configurations without exceptions and rounding modes by virtue of
existing Implies files.  As previously discussed, optimized versions
for particular architectures are possible, but not included.

i386 gets a special version of f32xaddf64 to avoid problems with
double rounding (similar to the existing fdim version), since this
function must round just once without an intermediate rounding to long
double.  (No such special version is needed for any other function,
because the nontrivial functions use round-to-odd, which does the
intermediate computation with the rounding mode set to round-to-zero,
and double rounding is OK except in round-to-nearest mode, so is OK
for that intermediate round-to-zero computation.)  mul and div will
need slightly different special versions for i386 (using round-to-odd
on long double instead of precision control) because of the
possibility of inexact intermediate results in the subnormal range for
double.

To reduce duplication among the different function implementations,
math-narrow.h gets macros CHECK_NARROW_ADD, NARROW_ADD_ROUND_TO_ODD
and NARROW_ADD_TRIVIAL.

In the trivial cases and for any architecture-specific optimized
implementations, the overhead of the errno setting might be
significant, but I think that's best handled through compiler built-in
functions rather than providing separate no-errno versions in glibc
(and likewise there are no __*_finite entry points for these function
provided, __*_finite effectively being no-errno versions at present in
most cases).

Tested for x86_64 and x86, with both GCC 6 and GCC 7.  Tested for
mips64 (all three ABIs, both hard and soft float) and powerpc with GCC
7.  Tested with build-many-glibcs.py with both GCC 6 and GCC 7.

	* math/Makefile (libm-narrow-fns): Add add.
	(libm-test-funcs-narrow): Likewise.
	* math/Versions (GLIBC_2.28): Add narrowing add functions.
	* math/bits/mathcalls-narrow.h (add): Use __MATHCALL_NARROW .
	* math/gen-auto-libm-tests.c (test_functions): Add add.
	* math/math-narrow.h (CHECK_NARROW_ADD): New macro.
	(NARROW_ADD_ROUND_TO_ODD): Likewise.
	(NARROW_ADD_TRIVIAL): Likewise.
	* sysdeps/ieee754/float128/float128_private.h (__faddl): New
	macro.
	(__daddl): Likewise.
	* sysdeps/ieee754/ldbl-opt/Makefile (libnldbl-calls): Add fadd and
	dadd.
	(CFLAGS-nldbl-dadd.c): New variable.
	(CFLAGS-nldbl-fadd.c): Likewise.
	* sysdeps/ieee754/ldbl-opt/Versions (GLIBC_2.28): Add
	__nldbl_daddl.
	* sysdeps/ieee754/ldbl-opt/nldbl-compat.h (__nldbl_daddl): New
	prototype.
	* manual/arith.texi (Misc FP Arithmetic): Document fadd, faddl,
	daddl, fMaddfN, fMaddfNx, fMxaddfN and fMxaddfNx.
	* math/auto-libm-test-in: Add tests of add.
	* math/auto-libm-test-out-narrow-add: New generated file.
	* math/libm-test-narrow-add.inc: New file.
	* sysdeps/i386/fpu/s_f32xaddf64.c: Likewise.
	* sysdeps/ieee754/dbl-64/s_f32xaddf64.c: Likewise.
	* sysdeps/ieee754/dbl-64/s_fadd.c: Likewise.
	* sysdeps/ieee754/float128/s_f32addf128.c: Likewise.
	* sysdeps/ieee754/float128/s_f64addf128.c: Likewise.
	* sysdeps/ieee754/float128/s_f64xaddf128.c: Likewise.
	* sysdeps/ieee754/ldbl-128/s_daddl.c: Likewise.
	* sysdeps/ieee754/ldbl-128/s_f64xaddf128.c: Likewise.
	* sysdeps/ieee754/ldbl-128/s_faddl.c: Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_daddl.c: Likewise.
	* sysdeps/ieee754/ldbl-128ibm/s_faddl.c: Likewise.
	* sysdeps/ieee754/ldbl-96/s_daddl.c: Likewise.
	* sysdeps/ieee754/ldbl-96/s_faddl.c: Likewise.
	* sysdeps/ieee754/ldbl-opt/nldbl-dadd.c: Likewise.
	* sysdeps/ieee754/ldbl-opt/nldbl-fadd.c: Likewise.
	* sysdeps/ieee754/soft-fp/s_daddl.c: Likewise.
	* sysdeps/ieee754/soft-fp/s_fadd.c: Likewise.
	* sysdeps/ieee754/soft-fp/s_faddl.c: Likewise.
	* sysdeps/powerpc/fpu/libm-test-ulps: Update.
	* sysdeps/mach/hurd/i386/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/aarch64/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/alpha/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/arm/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/hppa/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/i386/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/ia64/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/m68k/coldfire/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/m68k/m680x0/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/microblaze/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/mips/mips32/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/mips/mips64/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/nios2/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/powerpc/powerpc32/fpu/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/powerpc/powerpc32/nofpu/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/powerpc/powerpc64/libm-le.abilist: Likewise.
	* sysdeps/unix/sysv/linux/powerpc/powerpc64/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/riscv/rv64/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/s390/s390-32/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/s390/s390-64/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/sh/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/sparc/sparc32/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/sparc/sparc64/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/tile/tilegx32/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/tile/tilegx64/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/x86_64/64/libm.abilist: Likewise.
	* sysdeps/unix/sysv/linux/x86_64/x32/libm.abilist: Likewise.
This commit is contained in:
Joseph Myers 2018-02-10 02:08:43 +00:00
parent a85b70dbd2
commit d8742dd82f
61 changed files with 24818 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

71
ChangeLog
View File

@ -1,5 +1,76 @@
2018-02-10 Joseph Myers <joseph@codesourcery.com>
* math/Makefile (libm-narrow-fns): Add add.
(libm-test-funcs-narrow): Likewise.
* math/Versions (GLIBC_2.28): Add narrowing add functions.
* math/bits/mathcalls-narrow.h (add): Use __MATHCALL_NARROW .
* math/gen-auto-libm-tests.c (test_functions): Add add.
* math/math-narrow.h (CHECK_NARROW_ADD): New macro.
(NARROW_ADD_ROUND_TO_ODD): Likewise.
(NARROW_ADD_TRIVIAL): Likewise.
* sysdeps/ieee754/float128/float128_private.h (__faddl): New
macro.
(__daddl): Likewise.
* sysdeps/ieee754/ldbl-opt/Makefile (libnldbl-calls): Add fadd and
dadd.
(CFLAGS-nldbl-dadd.c): New variable.
(CFLAGS-nldbl-fadd.c): Likewise.
* sysdeps/ieee754/ldbl-opt/Versions (GLIBC_2.28): Add
__nldbl_daddl.
* sysdeps/ieee754/ldbl-opt/nldbl-compat.h (__nldbl_daddl): New
prototype.
* manual/arith.texi (Misc FP Arithmetic): Document fadd, faddl,
daddl, fMaddfN, fMaddfNx, fMxaddfN and fMxaddfNx.
* math/auto-libm-test-in: Add tests of add.
* math/auto-libm-test-out-narrow-add: New generated file.
* math/libm-test-narrow-add.inc: New file.
* sysdeps/i386/fpu/s_f32xaddf64.c: Likewise.
* sysdeps/ieee754/dbl-64/s_f32xaddf64.c: Likewise.
* sysdeps/ieee754/dbl-64/s_fadd.c: Likewise.
* sysdeps/ieee754/float128/s_f32addf128.c: Likewise.
* sysdeps/ieee754/float128/s_f64addf128.c: Likewise.
* sysdeps/ieee754/float128/s_f64xaddf128.c: Likewise.
* sysdeps/ieee754/ldbl-128/s_daddl.c: Likewise.
* sysdeps/ieee754/ldbl-128/s_f64xaddf128.c: Likewise.
* sysdeps/ieee754/ldbl-128/s_faddl.c: Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_daddl.c: Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_faddl.c: Likewise.
* sysdeps/ieee754/ldbl-96/s_daddl.c: Likewise.
* sysdeps/ieee754/ldbl-96/s_faddl.c: Likewise.
* sysdeps/ieee754/ldbl-opt/nldbl-dadd.c: Likewise.
* sysdeps/ieee754/ldbl-opt/nldbl-fadd.c: Likewise.
* sysdeps/ieee754/soft-fp/s_daddl.c: Likewise.
* sysdeps/ieee754/soft-fp/s_fadd.c: Likewise.
* sysdeps/ieee754/soft-fp/s_faddl.c: Likewise.
* sysdeps/powerpc/fpu/libm-test-ulps: Update.
* sysdeps/mach/hurd/i386/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/aarch64/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/alpha/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/arm/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/hppa/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/i386/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/ia64/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/m68k/coldfire/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/m68k/m680x0/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/microblaze/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/mips/mips32/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/mips/mips64/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/nios2/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/powerpc/powerpc32/fpu/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/powerpc/powerpc32/nofpu/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/powerpc/powerpc64/libm-le.abilist: Likewise.
* sysdeps/unix/sysv/linux/powerpc/powerpc64/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/riscv/rv64/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/s390/s390-32/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/s390/s390-64/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/sh/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/sparc/sparc32/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/sparc/sparc64/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/tile/tilegx32/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/tile/tilegx64/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/x86_64/64/libm.abilist: Likewise.
* sysdeps/unix/sysv/linux/x86_64/x32/libm.abilist: Likewise.
* sysdeps/powerpc/powerpc64le/Makefile [$(subdir) = math]
(f128-pairs): New variable.
[$(subdir) = math] ($(foreach suf,$(all-object-suffixes),$(foreach

6
NEWS
View File

@ -9,7 +9,11 @@ Version 2.28
Major new features:
[Add new features here]
* <math.h> functions that round their results to a narrower type are added
from TS 18661-1:2014 and TS 18661-3:2015:
- fadd, faddl, daddl and corresponding fMaddfN, fMaddfNx, fMxaddfN and
fMxaddfNx functions.
Deprecated and removed features, and other changes affecting compatibility:

19
manual/arith.texi
View File

@ -2191,6 +2191,25 @@ In @theglibc{}, this always means the operation is implemented in
hardware.
@end deftypefun
@deftypefun float fadd (double @var{x}, double @var{y})
@deftypefunx float faddl (long double @var{x}, long double @var{y})
@deftypefunx double daddl (long double @var{x}, long double @var{y})
@deftypefunx _FloatM fMaddfN (_Float@var{N} @var{x}, _Float@var{N} @var{y})
@deftypefunx _FloatM fMaddfNx (_Float@var{N}x @var{x}, _Float@var{N}x @var{y})
@deftypefunx _FloatMx fMxaddfN (_Float@var{N} @var{x}, _Float@var{N} @var{y})
@deftypefunx _FloatMx fMxaddfNx (_Float@var{N}x @var{x}, _Float@var{N}x @var{y})
@standards{TS 18661-1:2014, math.h}
@standardsx{fMaddfN, TS 18661-3:2015, math.h}
@standardsx{fMaddfNx, TS 18661-3:2015, math.h}
@standardsx{fMxaddfN, TS 18661-3:2015, math.h}
@standardsx{fMxaddfNx, TS 18661-3:2015, math.h}
@safety{@prelim{}@mtsafe{}@assafe{}@acsafe{}}
These functions, from TS 18661-1:2014 and TS 18661-3:2015, return
@math{@var{x} + @var{y}}, rounded once to the return type of the
function without any intermediate rounding to the type of the
arguments.
@end deftypefun
@node Complex Numbers
@section Complex Numbers
@pindex complex.h

4
math/Makefile
View File

@ -89,7 +89,7 @@ libm-compat-calls = \
w_lgammaF_r_compat w_lgammaF_compat2 w_expF_compat \
w_lgamma_compatF k_standardF
libm-narrow-fns =
libm-narrow-fns = add
libm-narrow-types-basic = s_fF s_f32xFf64
libm-narrow-types-ldouble-yes = s_fFl s_dFl
libm-narrow-types-float128-yes = s_f32Ff128 s_f64Ff128 s_f64xFf128
@ -284,7 +284,7 @@ libm-test-funcs-noauto = canonicalize ceil cimag conj copysign cproj creal \
scalbln scalbn setpayload setpayloadsig signbit \
significand totalorder totalordermag trunc ufromfp \
ufromfpx
libm-test-funcs-narrow =
libm-test-funcs-narrow = add
libm-test-funcs-all = $(libm-test-funcs-auto) $(libm-test-funcs-noauto)
libm-test-c-auto = $(foreach f,$(libm-test-funcs-auto),libm-test-$(f).c)
libm-test-c-noauto = $(foreach f,$(libm-test-funcs-noauto),libm-test-$(f).c)

9
math/Versions
View File

@ -551,4 +551,13 @@ libm {
y1f32; y1f64; y1f32x;
ynf32; ynf64; ynf32x;
}
GLIBC_2.28 {
# Functions not involving _Float64x or _Float128, for all configurations.
fadd; faddl; daddl;
f32addf32x; f32addf64; f32xaddf64;
# Functions involving _Float64x or _Float128, for some configurations.
f32addf64x; f32addf128;
f32xaddf64x; f32xaddf128; f64addf64x; f64addf128;
f64xaddf128;
}
}

95
math/auto-libm-test-in
View File

@ -142,6 +142,101 @@ acosh 0x1.1808eep+0
acosh 0x1.1052c4p+0
acosh max no-test-inline
add 0 0
add 0 -0
add -0 0
add -0 -0
add max max
add max -max
add -max max
add -max -max
add min min missing-underflow:arg-ibm128
add min -min missing-underflow:arg-ibm128
add -min min missing-underflow:arg-ibm128
add -min -min missing-underflow:arg-ibm128
add min_subnorm min_subnorm missing-underflow:arg-ibm128
add min_subnorm -min_subnorm missing-underflow:arg-ibm128
add -min_subnorm min_subnorm missing-underflow:arg-ibm128
add -min_subnorm -min_subnorm missing-underflow:arg-ibm128
add 1 2
add 1 -2
add -1 2
add -1 -2
add 100.5 0.75
add 100.5 -0.75
add -100.5 0.75
add -100.5 -0.75
add 1 0x1p-23
add 1 0x1.7fp-23
add 1 0x1.8p-23
add 1 0x1.81p-23
add 1 0x1p-24
add 1 0x1.1p-24
add 1 0x0.fp-24
add 1 min
add 1 -min
add 1 min_subnorm
add 1 -min_subnorm
add -1 min
add -1 -min
add -1 min_subnorm
add -1 -min_subnorm
# Cases where larger argument is half way between two values of a
# floating-point format, so that double rounding would sometimes yield
# the wrong result.
add 0x1.000001p0 min
add 0x1.000001p0 -min
add 0x1.000001p0 min_subnorm
add 0x1.000001p0 -min_subnorm
add -0x1.000001p0 min
add -0x1.000001p0 -min
add -0x1.000001p0 min_subnorm
add -0x1.000001p0 -min_subnorm
add 0x1.00000000000008p0 min
add 0x1.00000000000008p0 -min
add 0x1.00000000000008p0 min_subnorm
add 0x1.00000000000008p0 -min_subnorm
add -0x1.00000000000008p0 min
add -0x1.00000000000008p0 -min
add -0x1.00000000000008p0 min_subnorm
add -0x1.00000000000008p0 -min_subnorm
add 0x1.0000000000000001p0 min
add 0x1.0000000000000001p0 -min
add 0x1.0000000000000001p0 min_subnorm
add 0x1.0000000000000001p0 -min_subnorm
add -0x1.0000000000000001p0 min
add -0x1.0000000000000001p0 -min
add -0x1.0000000000000001p0 min_subnorm
add -0x1.0000000000000001p0 -min_subnorm
# Two values representable in the same format with sum very close to
# half way between two representable values. In particular, verify
# that there is no intermediate rounding to 64 bits when adding two
# double values.
add 1 0x1.000002p-24
add 1 -0x1.000002p-24
add 1 0x0.ffffffp-24
add 1 -0x0.ffffffp-24
add 0x1.000002p0 0x1.000002p-24
add 0x1.000002p0 -0x1.000002p-24
add 0x1.000002p0 0x0.ffffffp-24
add 0x1.000002p0 -0x0.ffffffp-24
add 1 0x1.0000000000001p-53
add 1 -0x1.0000000000001p-53
add 1 0x0.fffffffffffff8p-53
add 1 -0x0.fffffffffffff8p-53
add 0x1.0000000000001p0 0x1.0000000000001p-53
add 0x1.0000000000001p0 -0x1.0000000000001p-53
add 0x1.0000000000001p0 0x0.fffffffffffff8p-53
add 0x1.0000000000001p0 -0x0.fffffffffffff8p-53
add 1 0x1.0000000000000002p-64
add 1 -0x1.0000000000000002p-64
add 1 0x0.ffffffffffffffffp-64
add 1 -0x0.ffffffffffffffffp-64
add 0x1.0000000000000002p0 0x1.0000000000000002p-64
add 0x1.0000000000000002p0 -0x1.0000000000000002p-64
add 0x1.0000000000000002p0 0x0.ffffffffffffffffp-64
add 0x1.0000000000000002p0 -0x0.ffffffffffffffffp-64
asin 0
asin -0
asin 0.5

23487
math/auto-libm-test-out-narrow-add Normal file
View File

File diff suppressed because it is too large Load Diff

3
math/bits/mathcalls-narrow.h
View File

@ -19,3 +19,6 @@
#ifndef _MATH_H
# error "Never include <bits/mathcalls-narrow.h> directly; include <math.h> instead."
#endif
/* Add. */
__MATHCALL_NARROW (__MATHCALL_NAME (add), __MATHCALL_REDIR_NAME (add), 2);

1
math/gen-auto-libm-tests.c
View File

@ -531,6 +531,7 @@ static test_function test_functions[] =
{
FUNC_mpfr_f_f ("acos", mpfr_acos, false),
FUNC_mpfr_f_f ("acosh", mpfr_acosh, false),
FUNC_mpfr_ff_f ("add", mpfr_add, true),
FUNC_mpfr_f_f ("asin", mpfr_asin, false),
FUNC_mpfr_f_f ("asinh", mpfr_asinh, false),
FUNC_mpfr_f_f ("atan", mpfr_atan, false),

179
math/libm-test-narrow-add.inc Normal file
View File

@ -0,0 +1,179 @@
/* Test narrowing add.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#include "libm-test-driver.c"
static const struct test_aa_f_data add_test_data[] =
{
TEST_aa_f (add, arg_qnan_value, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_qnan_value, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_qnan_value, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_qnan_value, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_qnan_value, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_qnan_value, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_qnan_value, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_qnan_value, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_snan_value, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_snan_value, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_snan_value, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_snan_value, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_qnan_value, arg_plus_zero, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_qnan_value, arg_minus_zero, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_qnan_value, arg_plus_infty, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_qnan_value, arg_minus_infty, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_qnan_value, arg_min_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_qnan_value, -arg_min_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_qnan_value, arg_min_subnorm_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_qnan_value, -arg_min_subnorm_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_qnan_value, arg_max_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_qnan_value, -arg_max_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_plus_zero, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_minus_zero, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_plus_infty, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_minus_infty, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_min_value, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_min_value, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_min_subnorm_value, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_min_subnorm_value, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_max_value, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_max_value, arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_qnan_value, arg_plus_zero, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_qnan_value, arg_minus_zero, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_qnan_value, arg_plus_infty, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_qnan_value, arg_minus_infty, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_qnan_value, arg_min_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_qnan_value, -arg_min_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_qnan_value, arg_min_subnorm_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_qnan_value, -arg_min_subnorm_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_qnan_value, arg_max_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_qnan_value, -arg_max_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_plus_zero, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_minus_zero, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_plus_infty, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_minus_infty, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_min_value, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_min_value, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_min_subnorm_value, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_min_subnorm_value, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_max_value, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_max_value, -arg_qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_snan_value, arg_plus_zero, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_snan_value, arg_minus_zero, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_snan_value, arg_plus_infty, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_snan_value, arg_minus_infty, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_snan_value, arg_min_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_snan_value, -arg_min_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_snan_value, arg_min_subnorm_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_snan_value, -arg_min_subnorm_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_snan_value, arg_max_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_snan_value, -arg_max_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_plus_zero, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_minus_zero, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_plus_infty, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_minus_infty, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_min_value, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_min_value, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_min_subnorm_value, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_min_subnorm_value, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_max_value, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_max_value, arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, arg_plus_zero, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, arg_minus_zero, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, arg_plus_infty, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, arg_minus_infty, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, arg_min_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, -arg_min_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, arg_min_subnorm_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, -arg_min_subnorm_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, arg_max_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_snan_value, -arg_max_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_plus_zero, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_minus_zero, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_plus_infty, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_minus_infty, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_min_value, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_min_value, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_min_subnorm_value, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_min_subnorm_value, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_max_value, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, -arg_max_value, -arg_snan_value, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION),
TEST_aa_f (add, arg_plus_infty, arg_plus_infty, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_plus_infty, arg_minus_infty, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION|ERRNO_EDOM),
TEST_aa_f (add, arg_minus_infty, arg_plus_infty, qnan_value, NO_INEXACT_EXCEPTION|INVALID_EXCEPTION|ERRNO_EDOM),
TEST_aa_f (add, arg_minus_infty, arg_minus_infty, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_plus_infty, arg_plus_zero, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_plus_infty, arg_minus_zero, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_plus_infty, arg_min_value, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_plus_infty, -arg_min_value, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_plus_infty, arg_min_subnorm_value, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_plus_infty, -arg_min_subnorm_value, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_plus_infty, arg_max_value, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED|XFAIL_ROUNDING_IBM128_LIBGCC),
TEST_aa_f (add, arg_plus_infty, -arg_max_value, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED|XFAIL_ROUNDING_IBM128_LIBGCC),
TEST_aa_f (add, arg_plus_zero, arg_plus_infty, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_minus_zero, arg_plus_infty, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_min_value, arg_plus_infty, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_min_value, arg_plus_infty, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_min_subnorm_value, arg_plus_infty, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_min_subnorm_value, arg_plus_infty, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_max_value, arg_plus_infty, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_max_value, arg_plus_infty, plus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_minus_infty, arg_plus_zero, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_minus_infty, arg_minus_zero, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_minus_infty, arg_min_value, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_minus_infty, -arg_min_value, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_minus_infty, arg_min_subnorm_value, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_minus_infty, -arg_min_subnorm_value, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_minus_infty, arg_max_value, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED|XFAIL_ROUNDING_IBM128_LIBGCC),
TEST_aa_f (add, arg_minus_infty, -arg_max_value, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED|XFAIL_ROUNDING_IBM128_LIBGCC),
TEST_aa_f (add, arg_plus_zero, arg_minus_infty, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_minus_zero, arg_minus_infty, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_min_value, arg_minus_infty, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_min_value, arg_minus_infty, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_min_subnorm_value, arg_minus_infty, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_min_subnorm_value, arg_minus_infty, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, arg_max_value, arg_minus_infty, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
TEST_aa_f (add, -arg_max_value, arg_minus_infty, minus_infty, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED),
AUTO_TESTS_aa_f (add),
};
static void
add_test (void)
{
ALL_RM_TEST (add, 1, add_test_data, RUN_TEST_LOOP_aa_f, END);
}
static void
do_test (void)
{
add_test ();
}
/*
* Local Variables:
* mode:c
* End:
*/

58
math/math-narrow.h
View File

@ -49,6 +49,64 @@
u.d; \
})
/* Check for error conditions from a narrowing add function returning
RET with arguments X and Y and set errno as needed. Overflow and
underflow can occur for finite arguments and a domain error for
infinite ones. */
#define CHECK_NARROW_ADD(RET, X, Y) \
do \
{ \
if (!isfinite (RET)) \
{ \
if (isnan (RET)) \
{ \
if (!isnan (X) && !isnan (Y)) \
__set_errno (EDOM); \
} \
else if (isfinite (X) && isfinite (Y)) \
__set_errno (ERANGE); \
} \
else if ((RET) == 0 && (X) != -(Y)) \
__set_errno (ERANGE); \
} \
while (0)
/* Implement narrowing add using round-to-odd. The arguments are X
and Y, the return type is TYPE and UNION, MANTISSA and SUFFIX are
as for ROUND_TO_ODD. */
#define NARROW_ADD_ROUND_TO_ODD(X, Y, TYPE, UNION, SUFFIX, MANTISSA) \
do \
{ \
TYPE ret; \
\
/* Ensure a zero result is computed in the original rounding \
mode. */ \
if ((X) == -(Y)) \
ret = (TYPE) ((X) + (Y)); \
else \
ret = (TYPE) ROUND_TO_ODD (math_opt_barrier (X) + (Y), \
UNION, SUFFIX, MANTISSA); \
\
CHECK_NARROW_ADD (ret, (X), (Y)); \
return ret; \
} \
while (0)
/* Implement a narrowing add function that is not actually narrowing
or where no attempt is made to be correctly rounding (the latter
only applies to IBM long double). The arguments are X and Y and
the return type is TYPE. */
#define NARROW_ADD_TRIVIAL(X, Y, TYPE) \
do \
{ \
TYPE ret; \
\
ret = (TYPE) ((X) + (Y)); \
CHECK_NARROW_ADD (ret, (X), (Y)); \
return ret; \
} \
while (0)
/* The following macros declare aliases for a narrowing function. The
sole argument is the base name of a family of functions, such as
"add". If any platform changes long double format after the

41
sysdeps/i386/fpu/s_f32xaddf64.c Normal file
View File

@ -0,0 +1,41 @@
/* Add _Float64 values, converting the result to _Float32x. i386 version.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#include <math.h>
#include <fpu_control.h>
#include <math_private.h>
#include <math-narrow.h>
_Float32x
__f32xaddf64 (_Float64 x, _Float64 y)
{
/* To avoid double rounding, set double precision for the addition.
math_narrow_eval is still needed to eliminate excess range in the
case of overflow. If the result of the addition is in the
subnormal range for double, it is exact, so no issues of double
rounding for subnormals arise. */
fpu_control_t cw, cw_double;
_FPU_GETCW (cw);
cw_double = (cw & ~_FPU_EXTENDED) | _FPU_DOUBLE;
_FPU_SETCW (cw_double);
_Float32x ret = math_narrow_eval (x + y);
_FPU_SETCW (cw);
CHECK_NARROW_ADD (ret, x, y);
return ret;
}
libm_alias_float32x_float64 (add)

30
sysdeps/ieee754/dbl-64/s_f32xaddf64.c Normal file
View File

@ -0,0 +1,30 @@
/* Add _Float64 values, converting the result to _Float32x.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#define daddl __hide_daddl
#include <math.h>
#undef daddl
#include <math-narrow.h>
_Float32x
__f32xaddf64 (_Float64 x, _Float64 y)
{
NARROW_ADD_TRIVIAL (x, y, _Float32x);
}
libm_alias_float32x_float64 (add)

34
sysdeps/ieee754/dbl-64/s_fadd.c Normal file
View File

@ -0,0 +1,34 @@
/* Add double values, narrowing the result to float.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#define f32addf64 __hide_f32addf64
#define f32addf32x __hide_f32addf32x
#define faddl __hide_faddl
#include <math.h>
#undef f32addf64
#undef f32addf32x
#undef faddl
#include <math-narrow.h>
float
__fadd (double x, double y)
{
NARROW_ADD_ROUND_TO_ODD (x, y, float, union ieee754_double, , mantissa1);
}
libm_alias_float_double (add)

3
sysdeps/ieee754/float128/float128_private.h
View File

@ -256,6 +256,9 @@
#define __truncl __truncf128
#define __x2y2m1l __x2y2m1f128
#define __faddl __f32addf128
#define __daddl __f64addf128
/* __nexttowardf128 is not _Float128 API. */
#define __nexttowardl __nexttowardf128_do_not_use
#define nexttowardl nexttowardf128_do_not_use

6
sysdeps/ieee754/float128/s_f32addf128.c Normal file
View File

@ -0,0 +1,6 @@
#define f32addf64x __hide_f32addf64x
#define f32addf128 __hide_f32addf128
#include <float128_private.h>
#undef f32addf64x
#undef f32addf128
#include "../ldbl-128/s_faddl.c"

10
sysdeps/ieee754/float128/s_f64addf128.c Normal file
View File

@ -0,0 +1,10 @@
#define f32xaddf64x __hide_f32xaddf64x
#define f32xaddf128 __hide_f32xaddf128
#define f64addf64x __hide_f64addf64x
#define f64addf128 __hide_f64addf128
#include <float128_private.h>
#undef f32xaddf64x
#undef f32xaddf128
#undef f64addf64x
#undef f64addf128
#include "../ldbl-128/s_daddl.c"

2
sysdeps/ieee754/float128/s_f64xaddf128.c Normal file
View File

@ -0,0 +1,2 @@
#include <float128_private.h>
#include "../ldbl-128/s_f64xaddf128.c"

37
sysdeps/ieee754/ldbl-128/s_daddl.c Normal file
View File

@ -0,0 +1,37 @@
/* Add long double (ldbl-128) values, narrowing the result to double.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#define f32xaddf64x __hide_f32xaddf64x
#define f32xaddf128 __hide_f32xaddf128
#define f64addf64x __hide_f64addf64x
#define f64addf128 __hide_f64addf128
#include <math.h>
#undef f32xaddf64x
#undef f32xaddf128
#undef f64addf64x
#undef f64addf128
#include <math-narrow.h>
double
__daddl (_Float128 x, _Float128 y)
{
NARROW_ADD_ROUND_TO_ODD (x, y, double, union ieee854_long_double, l,
mantissa3);
}
libm_alias_double_ldouble (add)

38
sysdeps/ieee754/ldbl-128/s_f64xaddf128.c Normal file
View File

@ -0,0 +1,38 @@
/* Add _Float128 values, converting the result to _Float64x.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#include <math.h>
#include <math-narrow.h>
/* math_ldbl.h defines _Float128 to long double for this directory,
but when they are different, this function must be defined with
_Float128 arguments to avoid defining an alias with an incompatible
type. */
#undef _Float128
_Float64x
__f64xaddf128 (_Float128 x, _Float128 y)
{
#if __HAVE_FLOAT64X_LONG_DOUBLE && __HAVE_DISTINCT_FLOAT128
NARROW_ADD_ROUND_TO_ODD (x, y, _Float64x, union ieee854_long_double, l,
mantissa3);
#else
NARROW_ADD_TRIVIAL (x, y, _Float64x);
#endif
}
libm_alias_float64x_float128 (add)

33
sysdeps/ieee754/ldbl-128/s_faddl.c Normal file
View File

@ -0,0 +1,33 @@
/* Add long double (ldbl-128) values, narrowing the result to float.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#define f32addf64x __hide_f32addf64x
#define f32addf128 __hide_f32addf128
#include <math.h>
#undef f32addf64x
#undef f32addf128
#include <math-narrow.h>
float
__faddl (_Float128 x, _Float128 y)
{
NARROW_ADD_ROUND_TO_ODD (x, y, float, union ieee854_long_double, l,
mantissa3);
}
libm_alias_float_ldouble (add)

27
sysdeps/ieee754/ldbl-128ibm/s_daddl.c Normal file
View File

@ -0,0 +1,27 @@
/* Add long double (ldbl-128ibm) values, narrowing the result to double.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#include <math.h>
#include <math-narrow.h>
double
__daddl (long double x, long double y)
{
NARROW_ADD_TRIVIAL (x, y, double);
}
libm_alias_double_ldouble (add)

27
sysdeps/ieee754/ldbl-128ibm/s_faddl.c Normal file
View File

@ -0,0 +1,27 @@
/* Add long double (ldbl-128ibm) values, narrowing the result to float.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#include <math.h>
#include <math-narrow.h>
float
__faddl (long double x, long double y)
{
NARROW_ADD_TRIVIAL (x, y, float);
}
libm_alias_float_ldouble (add)

33
sysdeps/ieee754/ldbl-96/s_daddl.c Normal file
View File

@ -0,0 +1,33 @@
/* Add long double (ldbl-96) values, narrowing the result to double.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#define f32xaddf64x __hide_f32xaddf64x
#define f64addf64x __hide_f64addf64x
#include <math.h>
#undef f32xaddf64x
#undef f64addf64x
#include <math-narrow.h>
double
__daddl (long double x, long double y)
{
NARROW_ADD_ROUND_TO_ODD (x, y, double, union ieee854_long_double, l,
mantissa1);
}
libm_alias_double_ldouble (add)

31
sysdeps/ieee754/ldbl-96/s_faddl.c Normal file
View File

@ -0,0 +1,31 @@
/* Add long double (ldbl-96) values, narrowing the result to float.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#define f32addf64x __hide_f32addf64x
#include <math.h>
#undef f32addf64x
#include <math-narrow.h>
float
__faddl (long double x, long double y)
{
NARROW_ADD_ROUND_TO_ODD (x, y, float, union ieee854_long_double, l,
mantissa1);
}
libm_alias_float_ldouble (add)

4
sysdeps/ieee754/ldbl-opt/Makefile
View File

@ -44,7 +44,7 @@ libnldbl-calls = asprintf dprintf fprintf fscanf fwprintf fwscanf iovfscanf \
isoc99_vwscanf isoc99_vfwscanf isoc99_vswscanf \
nextup nextdown totalorder totalordermag getpayload \
canonicalize setpayload setpayloadsig llogb fmaxmag fminmag \
roundeven fromfp ufromfp fromfpx ufromfpx
roundeven fromfp ufromfp fromfpx ufromfpx fadd dadd
libnldbl-routines = $(libnldbl-calls:%=nldbl-%)
libnldbl-inhibit-o = $(object-suffixes)
libnldbl-static-only-routines = $(libnldbl-routines)
@ -86,6 +86,7 @@ CFLAGS-nldbl-csinh.c = -fno-builtin-csinhl
CFLAGS-nldbl-csqrt.c = -fno-builtin-csqrtl
CFLAGS-nldbl-ctan.c = -fno-builtin-ctanl
CFLAGS-nldbl-ctanh.c = -fno-builtin-ctanhl
CFLAGS-nldbl-dadd.c = -fno-builtin-daddl
CFLAGS-nldbl-erf.c = -fno-builtin-erfl
CFLAGS-nldbl-erfc.c = -fno-builtin-erfcl
CFLAGS-nldbl-exp.c = -fno-builtin-expl
@ -93,6 +94,7 @@ CFLAGS-nldbl-exp10.c = -fno-builtin-exp10l
CFLAGS-nldbl-exp2.c = -fno-builtin-exp2l
CFLAGS-nldbl-expm1.c = -fno-builtin-expm1l
CFLAGS-nldbl-fabs.c = -fno-builtin-fabsl
CFLAGS-nldbl-fadd.c = -fno-builtin-faddl
CFLAGS-nldbl-fdim.c = -fno-builtin-fdiml
CFLAGS-nldbl-finite.c = -fno-builtin-finitel
CFLAGS-nldbl-floor.c = -fno-builtin-floorl

6
sysdeps/ieee754/ldbl-opt/Versions
View File

@ -97,4 +97,10 @@ libm {
sincosl; tanl; tanhl; truncl; expl; __finitel; __signbitl;
__fpclassifyl; nexttowardf; nexttoward; __nldbl_nexttowardf;
}
GLIBC_2.28 {
# Functions taking long double = double argument and rounding
# result to double (same as f32x*f64 functions, but those names
# are not reserved in TS 18661-1).
__nldbl_daddl;
}
}

3
sysdeps/ieee754/ldbl-opt/nldbl-compat.h
View File

@ -101,5 +101,8 @@ extern int __nldbl___obstack_vprintf_chk (struct obstack *, int, const char *,
__gnuc_va_list) __THROW;
extern void __nldbl___vsyslog_chk (int, int, const char *, va_list);
/* The original declarations of these were hidden by the including
file. */
extern double __nldbl_daddl (double, double) __THROW;
#endif /* __NLDBL_COMPAT_H */

28
sysdeps/ieee754/ldbl-opt/nldbl-dadd.c Normal file
View File

@ -0,0 +1,28 @@
/* Compatibility routine for IEEE double as long double for dadd.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#define daddl __hide_daddl
#include "nldbl-compat.h"
#undef daddl
double
attribute_hidden
daddl (double x, double y)
{
return __nldbl_daddl (x, y);
}

28
sysdeps/ieee754/ldbl-opt/nldbl-fadd.c Normal file
View File

@ -0,0 +1,28 @@
/* Compatibility routine for IEEE double as long double for fadd.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#define faddl __hide_faddl
#include "nldbl-compat.h"
#undef faddl
float
attribute_hidden
faddl (double x, double y)
{
return fadd (x, y);
}

59
sysdeps/ieee754/soft-fp/s_daddl.c Normal file
View File

@ -0,0 +1,59 @@
/* Add long double (ldbl-128) values, narrowing the result to double,
using soft-fp.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#define f32xaddf64x __hide_f32xaddf64x
#define f32xaddf128 __hide_f32xaddf128
#define f64addf64x __hide_f64addf64x
#define f64addf128 __hide_f64addf128
#include <math.h>
#undef f32xaddf64x
#undef f32xaddf128
#undef f64addf64x
#undef f64addf128
#include <math-narrow.h>
#include <soft-fp.h>
#include <double.h>
#include <quad.h>
double
__daddl (_Float128 x, _Float128 y)
{
FP_DECL_EX;
FP_DECL_Q (X);
FP_DECL_Q (Y);
FP_DECL_Q (R);
FP_DECL_D (RN);
double ret;
FP_INIT_ROUNDMODE;
FP_UNPACK_SEMIRAW_Q (X, x);
FP_UNPACK_SEMIRAW_Q (Y, y);
FP_ADD_Q (R, X, Y);
#if (2 * _FP_W_TYPE_SIZE) < _FP_FRACBITS_Q
FP_TRUNC (D, Q, 2, 4, RN, R);
#else
FP_TRUNC (D, Q, 1, 2, RN, R);
#endif
FP_PACK_SEMIRAW_D (ret, RN);
FP_HANDLE_EXCEPTIONS;
CHECK_NARROW_ADD (ret, x, y);
return ret;
}
libm_alias_double_ldouble (add)

56
sysdeps/ieee754/soft-fp/s_fadd.c Normal file
View File

@ -0,0 +1,56 @@
/* Add double values, narrowing the result to float, using soft-fp.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#define f32addf64 __hide_f32addf64
#define f32addf32x __hide_f32addf32x
#define faddl __hide_faddl
#include <math.h>
#undef f32addf64
#undef f32addf32x
#undef faddl
#include <math-narrow.h>
#include <soft-fp.h>
#include <single.h>
#include <double.h>
float
__fadd (double x, double y)
{
FP_DECL_EX;
FP_DECL_D (X);
FP_DECL_D (Y);
FP_DECL_D (R);
FP_DECL_S (RN);
float ret;
FP_INIT_ROUNDMODE;
FP_UNPACK_SEMIRAW_D (X, x);
FP_UNPACK_SEMIRAW_D (Y, y);
FP_ADD_D (R, X, Y);
#if _FP_W_TYPE_SIZE < _FP_FRACBITS_D
FP_TRUNC (S, D, 1, 2, RN, R);
#else
FP_TRUNC (S, D, 1, 1, RN, R);
#endif
FP_PACK_SEMIRAW_S (ret, RN);
FP_HANDLE_EXCEPTIONS;
CHECK_NARROW_ADD (ret, x, y);
return ret;
}
libm_alias_float_double (add)

55
sysdeps/ieee754/soft-fp/s_faddl.c Normal file
View File

@ -0,0 +1,55 @@
/* Add long double (ldbl-128) values, narrowing the result to float,
using soft-fp.
Copyright (C) 2018 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
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
<http://www.gnu.org/licenses/>. */
#define f32addf64x __hide_f32addf64x
#define f32addf128 __hide_f32addf128
#include <math.h>
#undef f32addf64x
#undef f32addf128
#include <math-narrow.h>
#include <soft-fp.h>
#include <single.h>
#include <quad.h>
float
__faddl (_Float128 x, _Float128 y)
{
FP_DECL_EX;
FP_DECL_Q (X);
FP_DECL_Q (Y);
FP_DECL_Q (R);
FP_DECL_S (RN);
float ret;
FP_INIT_ROUNDMODE;
FP_UNPACK_SEMIRAW_Q (X, x);
FP_UNPACK_SEMIRAW_Q (Y, y);
FP_ADD_Q (R, X, Y);
#if (2 * _FP_W_TYPE_SIZE) < _FP_FRACBITS_Q
FP_TRUNC (S, Q, 1, 4, RN, R);
#else
FP_TRUNC (S, Q, 1, 2, RN, R);
#endif
FP_PACK_SEMIRAW_S (ret, RN);
FP_HANDLE_EXCEPTIONS;
CHECK_NARROW_ADD (ret, x, y);
return ret;
}
libm_alias_float_ldouble (add)

14
sysdeps/mach/hurd/i386/libm.abilist
View File

@ -1024,4 +1024,18 @@ GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A

4
sysdeps/powerpc/fpu/libm-test-ulps
View File

@ -79,6 +79,10 @@ ifloat128: 2
ildouble: 3
ldouble: 4
Function: "add_ldouble":
double: 1
float: 1
Function: "asin":
float: 1
float128: 1

14
sysdeps/unix/sysv/linux/aarch64/libm.abilist
View File

@ -986,3 +986,17 @@ GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F

15
sysdeps/unix/sysv/linux/alpha/libm.abilist
View File

@ -996,6 +996,21 @@ GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 __nldbl_daddl F
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.3.4 GLIBC_2.3.4 A
GLIBC_2.3.4 __c1_cabsf F
GLIBC_2.3.4 __c1_cacosf F

7
sysdeps/unix/sysv/linux/arm/libm.abilist
View File

@ -435,6 +435,13 @@ GLIBC_2.27 y1f64 F
GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A
GLIBC_2.4 _LIB_VERSION D 0x4
GLIBC_2.4 __clog10 F

7
sysdeps/unix/sysv/linux/hppa/libm.abilist
View File

@ -747,5 +747,12 @@ GLIBC_2.27 y1f64 F
GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A
GLIBC_2.4 exp2l F

14
sysdeps/unix/sysv/linux/i386/libm.abilist
View File

@ -1033,4 +1033,18 @@ GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A

14
sysdeps/unix/sysv/linux/ia64/libm.abilist
View File

@ -962,4 +962,18 @@ GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A

7
sysdeps/unix/sysv/linux/m68k/coldfire/libm.abilist
View File

@ -435,6 +435,13 @@ GLIBC_2.27 y1f64 F
GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A
GLIBC_2.4 _LIB_VERSION D 0x4
GLIBC_2.4 __clog10 F

7
sysdeps/unix/sysv/linux/m68k/m680x0/libm.abilist
View File

@ -789,4 +789,11 @@ GLIBC_2.27 y1f64 F
GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A

7
sysdeps/unix/sysv/linux/microblaze/libm.abilist
View File

@ -746,3 +746,10 @@ GLIBC_2.27 y1f64 F
GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F

7
sysdeps/unix/sysv/linux/mips/mips32/libm.abilist
View File

@ -748,6 +748,13 @@ GLIBC_2.27 y1f64 F
GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A
GLIBC_2.4 exp2l F
_gp_disp _gp_disp A

14
sysdeps/unix/sysv/linux/mips/mips64/libm.abilist
View File

@ -988,4 +988,18 @@ GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A

7
sysdeps/unix/sysv/linux/nios2/libm.abilist
View File

@ -746,3 +746,10 @@ GLIBC_2.27 y1f64 F
GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F

8
sysdeps/unix/sysv/linux/powerpc/powerpc32/fpu/libm.abilist
View File

@ -791,6 +791,14 @@ GLIBC_2.27 y1f64 F
GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 __nldbl_daddl F
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A
GLIBC_2.4 __clog10l F
GLIBC_2.4 __finitel F

8
sysdeps/unix/sysv/linux/powerpc/powerpc32/nofpu/libm.abilist
View File

@ -790,6 +790,14 @@ GLIBC_2.27 y1f64 F
GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 __nldbl_daddl F
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A
GLIBC_2.4 __clog10l F
GLIBC_2.4 __finitel F

15
sysdeps/unix/sysv/linux/powerpc/powerpc64/libm-le.abilist
View File

@ -1027,3 +1027,18 @@ GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 __nldbl_daddl F
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F

8
sysdeps/unix/sysv/linux/powerpc/powerpc64/libm.abilist
View File

@ -466,6 +466,14 @@ GLIBC_2.27 y1f64 F
GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 __nldbl_daddl F
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.3 GLIBC_2.3 A
GLIBC_2.3 _LIB_VERSION D 0x4
GLIBC_2.3 __clog10 F

14
sysdeps/unix/sysv/linux/riscv/rv64/libm.abilist
View File

@ -968,3 +968,17 @@ GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.27 ynl F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F

15
sysdeps/unix/sysv/linux/s390/s390-32/libm.abilist
View File

@ -986,6 +986,21 @@ GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 __nldbl_daddl F
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A
GLIBC_2.4 __clog10l F
GLIBC_2.4 __finitel F

15
sysdeps/unix/sysv/linux/s390/s390-64/libm.abilist
View File

@ -984,6 +984,21 @@ GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 __nldbl_daddl F
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A
GLIBC_2.4 __clog10l F
GLIBC_2.4 __finitel F

7
sysdeps/unix/sysv/linux/sh/libm.abilist
View File

@ -747,5 +747,12 @@ GLIBC_2.27 y1f64 F
GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A
GLIBC_2.4 exp2l F

15
sysdeps/unix/sysv/linux/sparc/sparc32/libm.abilist
View File

@ -993,6 +993,21 @@ GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 __nldbl_daddl F
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A
GLIBC_2.4 __clog10l F
GLIBC_2.4 __finitel F

14
sysdeps/unix/sysv/linux/sparc/sparc64/libm.abilist
View File

@ -987,4 +987,18 @@ GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A

7
sysdeps/unix/sysv/linux/tile/tilegx32/libm.abilist
View File

@ -747,3 +747,10 @@ GLIBC_2.27 y1f64 F
GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F

7
sysdeps/unix/sysv/linux/tile/tilegx64/libm.abilist
View File

@ -747,3 +747,10 @@ GLIBC_2.27 y1f64 F
GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F

14
sysdeps/unix/sysv/linux/x86_64/64/libm.abilist
View File

@ -1022,4 +1022,18 @@ GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F
GLIBC_2.4 GLIBC_2.4 A

14
sysdeps/unix/sysv/linux/x86_64/x32/libm.abilist
View File

@ -1021,3 +1021,17 @@ GLIBC_2.27 ynf32 F
GLIBC_2.27 ynf32x F
GLIBC_2.27 ynf64 F
GLIBC_2.27 ynf64x F
GLIBC_2.28 GLIBC_2.28 A
GLIBC_2.28 daddl F
GLIBC_2.28 f32addf128 F
GLIBC_2.28 f32addf32x F
GLIBC_2.28 f32addf64 F
GLIBC_2.28 f32addf64x F
GLIBC_2.28 f32xaddf128 F
GLIBC_2.28 f32xaddf64 F
GLIBC_2.28 f32xaddf64x F
GLIBC_2.28 f64addf128 F
GLIBC_2.28 f64addf64x F
GLIBC_2.28 f64xaddf128 F
GLIBC_2.28 fadd F
GLIBC_2.28 faddl F