From 8400b18bfb630399e4b028cabe72ba55b580ef66 Mon Sep 17 00:00:00 2001 From: Roland McGrath Date: Wed, 1 Feb 2006 03:14:01 +0000 Subject: * sysdeps/unix/sysv/linux/alpha/fpu/Implies: New file. * sysdeps/alpha/fpu/s_ceil.c: Include math_ldbl_opt.h, add compat_symbol if LONG_DOUBLE_COMPAT. * sysdeps/alpha/fpu/s_copysign.c: Likewise. * sysdeps/alpha/fpu/s_fabs.c: Likewise. * sysdeps/alpha/fpu/s_floor.c: Likewise. * sysdeps/alpha/fpu/s_rint.c: Likewise. * sysdeps/alpha/soft-fp/e_sqrtl.c: New file. * sysdeps/alpha/Implies: Add ieee754/ldbl-128. * sysdeps/unix/sysv/linux/alpha/bits/wordsize.h: New file. * sysdeps/unix/sysv/linux/alpha/Implies: Add ieee754/ldbl-64-128, ieee754/ldbl-opt. * sysdeps/alpha/bits/mathdef.h (__NO_LONG_DOUBLE_MATH): Remove. * sysdeps/unix/sysv/linux/alpha/nldbl-abi.h: New file. * sysdeps/unix/sysv/linux/alpha/Makefile [$(subdir) = math] (libm-routines): Add multc3, divtc3. --- sysdeps/alpha/soft-fp/e_sqrtl.c | 37 +++++++++++++++++++++++++++++++++++++ 1 file changed, 37 insertions(+) create mode 100644 sysdeps/alpha/soft-fp/e_sqrtl.c (limited to 'sysdeps/alpha/soft-fp/e_sqrtl.c') diff --git a/sysdeps/alpha/soft-fp/e_sqrtl.c b/sysdeps/alpha/soft-fp/e_sqrtl.c new file mode 100644 index 0000000000..a1d09725b7 --- /dev/null +++ b/sysdeps/alpha/soft-fp/e_sqrtl.c @@ -0,0 +1,37 @@ +/* long double square root in software floating-point emulation. + Copyright (C) 1997, 1999, 2006 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, write to the Free + Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA + 02111-1307 USA. */ + +#include "local-soft-fp.h" + +long double +__ieee754_sqrtl (const long double a) +{ + FP_DECL_EX; + FP_DECL_Q(A); FP_DECL_Q(C); + long double c; + + FP_INIT_ROUNDMODE; + FP_UNPACK_Q(A, a); + FP_SQRT_Q(C, A); + FP_PACK_Q(c, C); + FP_HANDLE_EXCEPTIONS; + return c; +} -- cgit v1.2.3