diff options
| author | Roland McGrath <roland@gnu.org> | 2006-02-01 03:14:01 +0000 |
|---|---|---|
| committer | Roland McGrath <roland@gnu.org> | 2006-02-01 03:14:01 +0000 |
| commit | 8400b18bfb630399e4b028cabe72ba55b580ef66 (patch) | |
| tree | a8c8c9563007f8c8e54718ced8ff2776b2a9349c /sysdeps/alpha/soft-fp/e_sqrtl.c | |
| parent | fea3f995bd004c0e31898eed82d5e15113d50e70 (diff) | |
* sysdeps/unix/sysv/linux/alpha/fpu/Implies: New file.cvs/fedora-glibc-20060201T0846
* 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.
Diffstat (limited to 'sysdeps/alpha/soft-fp/e_sqrtl.c')
| -rw-r--r-- | sysdeps/alpha/soft-fp/e_sqrtl.c | 37 |
1 files changed, 37 insertions, 0 deletions
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; +} |