]>
Commit | Line | Data |
---|---|---|
f819dbea PM |
1 | /* soft-fp sqrt for _Float128 |
2 | Return sqrt(a) | |
688903eb | 3 | Copyright (C) 2017-2018 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 | |
27 | <http://www.gnu.org/licenses/>. */ | |
28 | ||
1c09524e TMQMF |
29 | #include <math.h> |
30 | ||
31 | /* Unavoidable hacks since TFmode is assumed to be binary128 when | |
32 | -mabi=ibmlongdouble is used. */ | |
33 | #if __HAVE_FLOAT128_UNLIKE_LDBL | |
34 | # define TFtype KFtype | |
35 | # define TF KF | |
36 | #endif | |
f819dbea PM |
37 | |
38 | #include <soft-fp.h> | |
39 | #include <quad.h> | |
40 | ||
41 | __float128 | |
42 | __ieee754_sqrtf128 (__float128 a) | |
43 | { | |
44 | FP_DECL_EX; | |
45 | FP_DECL_Q (A); | |
46 | FP_DECL_Q (R); | |
47 | __float128 r; | |
48 | ||
49 | FP_INIT_ROUNDMODE; | |
50 | FP_UNPACK_Q (A, a); | |
51 | FP_SQRT_Q (R, A); | |
52 | FP_PACK_Q (r, R); | |
53 | FP_HANDLE_EXCEPTIONS; | |
54 | return r; | |
55 | } | |
56 | strong_alias (__ieee754_sqrtf128, __sqrtf128_finite) |