projects / thirdparty / glibc.git / blame
| Commit | Line | Data |
|---|---|---|
| f819dbea PM |
1 | /* soft-fp sqrt for _Float128 |
| 2 | Return sqrt(a) | |
| 04277e02 | 3 | Copyright (C) 2017-2019 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 | |
| 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) |