X-Git-Url: http://erislabs.net/gitweb/?a=blobdiff_plain;f=lib%2Fsqrtl.c;h=61e9574d1a8c20ccc2d4773a72e5b9513c1e9f31;hb=beae0bdcf7fe30f2415c16f6b8c1368d469e519c;hp=bf1b2d5faa74594c771578fa6bd692e619b1ff07;hpb=55eed6bbc87934f2db04981e6f0f1773ff77974f;p=gnulib.git diff --git a/lib/sqrtl.c b/lib/sqrtl.c index bf1b2d5fa..61e9574d1 100644 --- a/lib/sqrtl.c +++ b/lib/sqrtl.c @@ -1,35 +1,43 @@ /* Emulation for sqrtl. Contributed by Paolo Bonzini - Copyright 2002, 2003, 2007 Free Software Foundation, Inc. + Copyright 2002-2003, 2007, 2009-2013 Free Software Foundation, Inc. This file is part of gnulib. - This program is free software; you can redistribute it and/or modify + 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 2, or (at your option) - any later version. + 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, write to the Free Software Foundation, - Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. */ + You should have received a copy of the GNU General Public License + along with this program. If not, see . */ #include /* Specification. */ #include -#include -#include "isnanl.h" +#if HAVE_SAME_LONG_DOUBLE_AS_DOUBLE + +long double +sqrtl (long double x) +{ + return sqrt (x); +} + +#else + +# include /* A simple Newton-Raphson method. */ long double -sqrtl(long double x) +sqrtl (long double x) { long double delta, y; int exponent; @@ -40,7 +48,7 @@ sqrtl(long double x) /* Check for negative numbers */ if (x < 0.0L) - return (long double) sqrt(-1); + return (long double) sqrt (-1); /* Check for zero and infinites */ if (x + x == x) @@ -59,3 +67,5 @@ sqrtl(long double x) return y; } + +#endif