]>
Commit | Line | Data |
---|---|---|
f819dbea PM |
1 | /* soft-fp sqrt for _Float128 |
2 | Return sqrt(a) | |
dff8da6b | 3 | Copyright (C) 2017-2024 Free Software Foundation, Inc. |
f819dbea PM |
4 | This file is part of the GNU C Library. |
5 | ||
6 | The GNU C Library is free software; you can redistribute it and/or | |
7 | modify it under the terms of the GNU Lesser General Public | |
8 | License as published by the Free Software Foundation; either | |
9 | version 2.1 of the License, or (at your option) any later version. | |
10 | ||
11 | In addition to the permissions in the GNU Lesser General Public | |
12 | License, the Free Software Foundation gives you unlimited | |
13 | permission to link the compiled version of this file into | |
14 | combinations with other programs, and to distribute those | |
15 | combinations without any restriction coming from the use of this | |
16 | file. (The Lesser General Public License restrictions do apply in | |
17 | other respects; for example, they cover modification of the file, | |
18 | and distribution when not linked into a combine executable.) | |
19 | ||
20 | The GNU C Library is distributed in the hope that it will be useful, | |
21 | but WITHOUT ANY WARRANTY; without even the implied warranty of | |
22 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
23 | Lesser General Public License for more details. | |
24 | ||
25 | You should have received a copy of the GNU Lesser General Public | |
26 | License along with the GNU C Library; if not, see | |
5a82c748 | 27 | <https://www.gnu.org/licenses/>. */ |
f819dbea | 28 | |
b637306d | 29 | #include <float128_private.h> |
1c09524e | 30 | #include <math.h> |
220622dd | 31 | #include <libm-alias-finite.h> |
1c09524e TMQMF |
32 | |
33 | /* Unavoidable hacks since TFmode is assumed to be binary128 when | |
34 | -mabi=ibmlongdouble is used. */ | |
35 | #if __HAVE_FLOAT128_UNLIKE_LDBL | |
36 | # define TFtype KFtype | |
37 | # define TF KF | |
38 | #endif | |
f819dbea PM |
39 | |
40 | #include <soft-fp.h> | |
41 | #include <quad.h> | |
42 | ||
43 | __float128 | |
44 | __ieee754_sqrtf128 (__float128 a) | |
45 | { | |
b637306d PM |
46 | __float128 r; |
47 | /* For multiarch builds, this will be built twice. */ | |
48 | #ifndef _ARCH_PWR9 | |
f819dbea PM |
49 | FP_DECL_EX; |
50 | FP_DECL_Q (A); | |
51 | FP_DECL_Q (R); | |
f819dbea PM |
52 | |
53 | FP_INIT_ROUNDMODE; | |
54 | FP_UNPACK_Q (A, a); | |
55 | FP_SQRT_Q (R, A); | |
56 | FP_PACK_Q (r, R); | |
57 | FP_HANDLE_EXCEPTIONS; | |
b637306d PM |
58 | #else |
59 | asm ("xssqrtqp %0,%1" : "=v" (r) : "v" (a)); | |
60 | #endif | |
f819dbea PM |
61 | return r; |
62 | } | |
220622dd | 63 | libm_alias_finite (__ieee754_sqrtf128, __sqrtf128) |