1 /* Fused multiply-add of _Float128 value, converting the result to _Float64x.
2 Copyright (C) 2021-2022 Free Software Foundation, Inc.
3 This file is part of the GNU C Library.
4
5 The GNU C Library is free software; you can redistribute it and/or
6 modify it under the terms of the GNU Lesser General Public
7 License as published by the Free Software Foundation; either
8 version 2.1 of the License, or (at your option) any later version.
9
10 The GNU C Library is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 Lesser General Public License for more details.
14
15 You should have received a copy of the GNU Lesser General Public
16 License along with the GNU C Library; if not, see
17 <https://www.gnu.org/licenses/>. */
18
19 #include <math.h>
20 #include <math-narrow.h>
21 #include <tininess.h>
22
23 /* math_ldbl.h defines _Float128 to long double for this directory,
24 but when they are different, this function must be defined with
25 _Float128 arguments to avoid defining an alias with an incompatible
26 type. */
27 #undef _Float128
28
29 #if __HAVE_FLOAT64X_LONG_DOUBLE && __HAVE_DISTINCT_FLOAT128
30 _Float64x
__f64xfmaf128(_Float128 x,_Float128 y,_Float128 z)31 __f64xfmaf128 (_Float128 x, _Float128 y, _Float128 z)
32 {
33 NARROW_FMA_ROUND_TO_ODD (x, y, z, _Float64x, union ieee854_long_double, l,
34 mantissa3, TININESS_AFTER_ROUNDING);
35 }
36 libm_alias_float64x_float128 (fma)
37 #else
38 /* Defined as an alias of fmal. */
39 #endif
40