diff options
author | Richard Maw <richard.maw@codethink.co.uk> | 2012-01-17 14:43:55 +0000 |
---|---|---|
committer | Richard Maw <richard.maw@codethink.co.uk> | 2012-01-17 14:43:55 +0000 |
commit | 2de9abc5c9d40b3c716307d67d16146f823fd554 (patch) | |
tree | 6979db67934ddc8b564150b465846a383b428ff8 /gnulib/lib/sqrtl.c | |
parent | 33cc1c6fda6e72a7bae1401e9b2cec495a4d3ff1 (diff) | |
download | patch-baserock/bootstrap.tar.gz |
add the output of bootstrapbaserock/bootstrap-pass2baserock/bootstrap
Diffstat (limited to 'gnulib/lib/sqrtl.c')
m--------- | gnulib | 0 | ||||
-rw-r--r-- | gnulib/lib/sqrtl.c | 71 |
2 files changed, 71 insertions, 0 deletions
diff --git a/gnulib b/gnulib deleted file mode 160000 -Subproject 443bc5ffcf7429e557f4a371b0661abe98ddbc1 diff --git a/gnulib/lib/sqrtl.c b/gnulib/lib/sqrtl.c new file mode 100644 index 0000000..6795d66 --- /dev/null +++ b/gnulib/lib/sqrtl.c @@ -0,0 +1,71 @@ +/* Emulation for sqrtl. + Contributed by Paolo Bonzini + + Copyright 2002-2003, 2007, 2009-2011 Free Software Foundation, Inc. + + This file is part of gnulib. + + This program is free software: you can redistribute it and/or modify + it under the terms of the GNU General Public License as published by + the Free Software Foundation; either version 3 of the License, or + (at your option) any later version. + + This program 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 General Public License for more details. + + You should have received a copy of the GNU General Public License + along with this program. If not, see <http://www.gnu.org/licenses/>. */ + +#include <config.h> + +/* Specification. */ +#include <math.h> + +#if HAVE_SAME_LONG_DOUBLE_AS_DOUBLE + +long double +sqrtl (long double x) +{ + return sqrt (x); +} + +#else + +# include <float.h> + +/* A simple Newton-Raphson method. */ +long double +sqrtl (long double x) +{ + long double delta, y; + int exponent; + + /* Check for NaN */ + if (isnanl (x)) + return x; + + /* Check for negative numbers */ + if (x < 0.0L) + return (long double) sqrt (-1); + + /* Check for zero and infinites */ + if (x + x == x) + return x; + + frexpl (x, &exponent); + y = ldexpl (x, -exponent / 2); + + do + { + delta = y; + y = (y + x / y) * 0.5L; + delta -= y; + } + while (delta != 0.0L); + + return y; +} + +#endif |