1/* long double square root in software floating-point emulation.
2 Copyright (C) 1997-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 <stdlib.h>
20#include <soft-fp.h>
21#include <quad.h>
22#include <shlib-compat.h>
23
24long double
25__ieee754_sqrtl (const long double a)
26{
27 FP_DECL_EX;
28 FP_DECL_Q(A); FP_DECL_Q(C);
29 long double c;
30 long _round = 4; /* dynamic rounding */
31
32 FP_INIT_ROUNDMODE;
33 FP_UNPACK_Q(A, a);
34 FP_SQRT_Q(C, A);
35 FP_PACK_Q(c, C);
36 FP_HANDLE_EXCEPTIONS;
37 return c;
38}
39
40/* ??? We forgot to add this symbol in 2.15. Getting this into 2.18 isn't as
41 straight-forward as just adding the alias, since a generic Versions file
42 includes the 2.15 version and the linker uses the first one it sees. */
43#if SHLIB_COMPAT (libm, GLIBC_2_15, GLIBC_2_18)
44compat_symbol (libm, __ieee754_sqrtl, __sqrtl_finite, GLIBC_2_18);
45#endif
46

source code of glibc/sysdeps/alpha/e_sqrtl.c