1 /* Square root 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
22 /* math_ldbl.h defines _Float128 to long double for this directory,
23 but when they are different, this function must be defined with
24 _Float128 arguments to avoid defining an alias with an incompatible
25 type. */
26 #undef _Float128
27
28 #if __HAVE_FLOAT64X_LONG_DOUBLE && __HAVE_DISTINCT_FLOAT128
29 _Float64x
__f64xsqrtf128(_Float128 x)30 __f64xsqrtf128 (_Float128 x)
31 {
32 NARROW_SQRT_ROUND_TO_ODD (x, _Float64x, union ieee854_long_double, l,
33 mantissa3);
34 }
35 libm_alias_float64x_float128 (sqrt)
36 #else
37 /* Defined as an alias of sqrtl. */
38 #endif
39