[thirdparty/glibc.git] / sysdeps / powerpc / powerpc64 / fpu / e_sqrt.c

/* Double-precision floating point square root.
   Copyright (C) 1997-2014 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.h>
#include <math_private.h>

#undef __ieee754_sqrt
double
__ieee754_sqrt (double x)
{
  double z;
  __asm __volatile ("fsqrt %0,%1" : "=f" (z) : "f" (x));
  return z;
}
strong_alias (__ieee754_sqrt, __sqrt_finite)
Commit	Line	Data
5593c5ec	1	/* Double-precision floating point square root.
d4697bc9	2	Copyright (C) 1997-2014 Free Software Foundation, Inc.
5593c5ec RM	3	This file is part of the GNU C Library.
	4
	5	The GNU C Library is free software; you can redistribute it and/or
	6	modify it under the terms of the GNU Lesser General Public
	7	License as published by the Free Software Foundation; either
	8	version 2.1 of the License, or (at your option) any later version.
	9
	10	The GNU C Library is distributed in the hope that it will be useful,
	11	but WITHOUT ANY WARRANTY; without even the implied warranty of
	12	MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
	13	Lesser General Public License for more details.
	14
	15	You should have received a copy of the GNU Lesser General Public
59ba27a6 PE	16	License along with the GNU C Library; if not, see
59ba27a6 PE	17	<http://www.gnu.org/licenses/>. */
5593c5ec RM	18
	19	#include <math.h>
	20	#include <math_private.h>
	21
8a6d5255	22	#undef __ieee754_sqrt
5593c5ec RM	23	double
	24	__ieee754_sqrt (double x)
	25	{
	26	double z;
	27	__asm __volatile ("fsqrt %0,%1" : "=f" (z) : "f" (x));
	28	return z;
	29	}
0ac5ae23	30	strong_alias (__ieee754_sqrt, __sqrt_finite)