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