[thirdparty/glibc.git] / sysdeps / sparc / sparc64 / soft-fp / qp_sqrt.c

/* Software floating-point emulation.
   (*c) = sqrtl(*a)
   Copyright (C) 1997-2014 Free Software Foundation, Inc.
   This file is part of the GNU C Library.
   Contributed by Richard Henderson (rth@cygnus.com) and
		  Jakub Jelinek (jj@ultra.linux.cz).

   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 "soft-fp.h"
#include "quad.h"

void _Qp_sqrt(long double *c, const long double *a)
{
  FP_DECL_EX;
  FP_DECL_Q(A); FP_DECL_Q(C);

  FP_INIT_ROUNDMODE;
  FP_UNPACK_QP(A, a);
  FP_SQRT_Q(C, A);
  FP_PACK_QP(c, C);
  QP_HANDLE_EXCEPTIONS(__asm (
"	ldd [%1], %%f52\n"
"	ldd [%1+8], %%f54\n"
"	fsqrtq %%f52, %%f60\n"
"	std %%f60, [%0]\n"
"	std %%f62, [%0+8]\n"
"	" : : "r" (c), "r" (a) : QP_CLOBBER));
}
Commit	Line	Data
d876f532 UD	1	/* Software floating-point emulation.
d876f532 UD	2	(c) = sqrtl(a)
d4697bc9	3	Copyright (C) 1997-2014 Free Software Foundation, Inc.
d876f532 UD	4	This file is part of the GNU C Library.
	5	Contributed by Richard Henderson (rth@cygnus.com) and
	6	Jakub Jelinek (jj@ultra.linux.cz).
	7
	8	The GNU C Library is free software; you can redistribute it and/or
41bdb6e2 AJ	9	modify it under the terms of the GNU Lesser General Public
	10	License as published by the Free Software Foundation; either
	11	version 2.1 of the License, or (at your option) any later version.
d876f532 UD	12
	13	The GNU C Library is distributed in the hope that it will be useful,
	14	but WITHOUT ANY WARRANTY; without even the implied warranty of
	15	MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
41bdb6e2	16	Lesser General Public License for more details.
d876f532	17
41bdb6e2	18	You should have received a copy of the GNU Lesser General Public
59ba27a6 PE	19	License along with the GNU C Library; if not, see
59ba27a6 PE	20	<http://www.gnu.org/licenses/>. */
d876f532 UD	21
	22	#include "soft-fp.h"
	23	#include "quad.h"
	24
	25	void _Qp_sqrt(long double c, const long double a)
	26	{
	27	FP_DECL_EX;
	28	FP_DECL_Q(A); FP_DECL_Q(C);
	29
	30	FP_INIT_ROUNDMODE;
	31	FP_UNPACK_QP(A, a);
	32	FP_SQRT_Q(C, A);
	33	FP_PACK_QP(c, C);
	34	QP_HANDLE_EXCEPTIONS(__asm (
ba77ce24 AJ	35	" ldd [%1], %%f52\n"
	36	" ldd [%1+8], %%f54\n"
	37	" fsqrtq %%f52, %%f60\n"
	38	" std %%f60, [%0]\n"
	39	" std %%f62, [%0+8]\n"
	40	" " : : "r" (c), "r" (a) : QP_CLOBBER));
d876f532	41	}