[thirdparty/glibc.git] / sysdeps / x86_64 / fpu / e_sqrt.c

/* Square root of floating point number.
   Copyright (C) 2002-2018 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>

#undef __ieee754_sqrt
double
__ieee754_sqrt (double x)
{
  double res;

  asm ("sqrtsd %1, %0" : "=x" (res) : "xm" (x));

  return res;
}
strong_alias (__ieee754_sqrt, __sqrt_finite)
Commit	Line	Data
2b904739	1	/* Square root of floating point number.
688903eb	2	Copyright (C) 2002-2018 Free Software Foundation, Inc.
2b904739 AJ	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/>. */
2b904739 AJ	18
	19	#include <math_private.h>
	20
0ac5ae23	21	#undef __ieee754_sqrt
2b904739 AJ	22	double
	23	__ieee754_sqrt (double x)
	24	{
	25	double res;
	26
0ac5ae23	27	asm ("sqrtsd %1, %0" : "=x" (res) : "xm" (x));
2b904739 AJ	28
	29	return res;
	30	}
0ac5ae23	31	strong_alias (__ieee754_sqrt, __sqrt_finite)