|
Packit |
6c4009 |
/* soft-fp sqrt for _Float128
|
|
Packit |
6c4009 |
Return sqrt(a)
|
|
Packit |
6c4009 |
Copyright (C) 2017-2018 Free Software Foundation, Inc.
|
|
Packit |
6c4009 |
This file is part of the GNU C Library.
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
The GNU C Library is free software; you can redistribute it and/or
|
|
Packit |
6c4009 |
modify it under the terms of the GNU Lesser General Public
|
|
Packit |
6c4009 |
License as published by the Free Software Foundation; either
|
|
Packit |
6c4009 |
version 2.1 of the License, or (at your option) any later version.
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
In addition to the permissions in the GNU Lesser General Public
|
|
Packit |
6c4009 |
License, the Free Software Foundation gives you unlimited
|
|
Packit |
6c4009 |
permission to link the compiled version of this file into
|
|
Packit |
6c4009 |
combinations with other programs, and to distribute those
|
|
Packit |
6c4009 |
combinations without any restriction coming from the use of this
|
|
Packit |
6c4009 |
file. (The Lesser General Public License restrictions do apply in
|
|
Packit |
6c4009 |
other respects; for example, they cover modification of the file,
|
|
Packit |
6c4009 |
and distribution when not linked into a combine executable.)
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
The GNU C Library is distributed in the hope that it will be useful,
|
|
Packit |
6c4009 |
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
Packit |
6c4009 |
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
Packit |
6c4009 |
Lesser General Public License for more details.
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
You should have received a copy of the GNU Lesser General Public
|
|
Packit |
6c4009 |
License along with the GNU C Library; if not, see
|
|
Packit |
6c4009 |
<http://www.gnu.org/licenses/>. */
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
#include <math.h>
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
/* Unavoidable hacks since TFmode is assumed to be binary128 when
|
|
Packit |
6c4009 |
-mabi=ibmlongdouble is used. */
|
|
Packit |
6c4009 |
#if __HAVE_FLOAT128_UNLIKE_LDBL
|
|
Packit |
6c4009 |
# define TFtype KFtype
|
|
Packit |
6c4009 |
# define TF KF
|
|
Packit |
6c4009 |
#endif
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
#include <soft-fp.h>
|
|
Packit |
6c4009 |
#include <quad.h>
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
__float128
|
|
Packit |
6c4009 |
__ieee754_sqrtf128 (__float128 a)
|
|
Packit |
6c4009 |
{
|
|
Packit |
6c4009 |
FP_DECL_EX;
|
|
Packit |
6c4009 |
FP_DECL_Q (A);
|
|
Packit |
6c4009 |
FP_DECL_Q (R);
|
|
Packit |
6c4009 |
__float128 r;
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
FP_INIT_ROUNDMODE;
|
|
Packit |
6c4009 |
FP_UNPACK_Q (A, a);
|
|
Packit |
6c4009 |
FP_SQRT_Q (R, A);
|
|
Packit |
6c4009 |
FP_PACK_Q (r, R);
|
|
Packit |
6c4009 |
FP_HANDLE_EXCEPTIONS;
|
|
Packit |
6c4009 |
return r;
|
|
Packit |
6c4009 |
}
|
|
Packit |
6c4009 |
strong_alias (__ieee754_sqrtf128, __sqrtf128_finite)
|