Inline __ieee754_sqrt and __ieee754_sqrtf. Also add external definitions.

This commit is contained in:
Wilco Dijkstra 2015-07-06 12:43:31 +01:00
parent a1b85e1d6c
commit 6471190491
4 changed files with 81 additions and 0 deletions

View File

@ -1,3 +1,12 @@
2015-07-06 Wilco Dijkstra <wdijkstr@arm.com>
* sysdeps/aarch64/fpu/math_private.h (__ieee754_sqrt):
New function. (__ieee754_sqrtf): New function.
* sysdeps/aarch64/fpu/e_sqrt.c (__ieee754_sqrt):
New function.
* sysdeps/aarch64/fpu/e_sqrtf.c (__ieee754_sqrtf):
New function.
2015-07-06 Wilco Dijkstra <wdijkstr@arm.com>
* sysdeps/ieee754/support.c: Remove unused file.

View File

@ -0,0 +1,28 @@
/* Square root of floating point number.
Copyright (C) 2015 Free Software Foundation, Inc.
This file is part of the GNU C Library.
The GNU C Library is free software; you can redistribute it and/or
modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.
The GNU C Library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with the GNU C Library; if not, see
<http://www.gnu.org/licenses/>. */
#include <math_private.h>
double
__ieee754_sqrt (double d)
{
double res;
asm ("fsqrt %d0, %d1" : "=w" (res) : "w" (d));
return res;
}
strong_alias (__ieee754_sqrt, __sqrt_finite)

View File

@ -0,0 +1,28 @@
/* Single-precision floating point square root.
Copyright (C) 2015 Free Software Foundation, Inc.
This file is part of the GNU C Library.
The GNU C Library is free software; you can redistribute it and/or
modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.
The GNU C Library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with the GNU C Library; if not, see
<http://www.gnu.org/licenses/>. */
#include <math_private.h>
float
__ieee754_sqrtf (float s)
{
float res;
asm ("fsqrt %s0, %s1" : "=w" (res) : "w" (s));
return res;
}
strong_alias (__ieee754_sqrtf, __sqrtf_finite)

View File

@ -22,6 +22,22 @@
#include <fenv.h>
#include <fpu_control.h>
extern __always_inline double
__ieee754_sqrt (double d)
{
double res;
asm __volatile__ ("fsqrt %d0, %d1" : "=w" (res) : "w" (d));
return res;
}
extern __always_inline float
__ieee754_sqrtf (float s)
{
float res;
asm __volatile__ ("fsqrt %s0, %s1" : "=w" (res) : "w" (s));
return res;
}
static __always_inline void
libc_feholdexcept_aarch64 (fenv_t *envp)
{