X-Git-Url: http://erislabs.net/gitweb/?a=blobdiff_plain;f=lib%2Flogl.c;h=52adaecc17d08469521f4180efeb227c23601ac6;hb=cc7f5298d74867d113d32db1ddfb5f3667cbaeff;hp=c4738ebc6a8cb93e9a4f96ef9adab75f66652cf7;hpb=a03ce664037ae2d0565aedc481fdd8150276b233;p=gnulib.git diff --git a/lib/logl.c b/lib/logl.c index c4738ebc6..52adaecc1 100644 --- a/lib/logl.c +++ b/lib/logl.c @@ -1,3 +1,23 @@ +/* Copyright 2001 by Stephen L. Moshier + + 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 . */ + +#include + +/* Specification. */ +#include + /* logll.c * * Natural logarithm for 128-bit long double precision. @@ -42,27 +62,6 @@ * */ -#include - -#include "mathl.h" - -/* Copyright 2001 by Stephen L. Moshier - - 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. - - 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; see the file COPYING. - If not, write to the Free Software Foundation, - 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ - /* log(1+x) = x - .5 x^2 + x^3 l(x) -.0078125 <= x <= +.0078125 peak relative error 1.2e-37 */ @@ -189,12 +188,16 @@ long double logl(long double x) { long double z, y, w; - long double u, t; - unsigned int m; + long double t; int k, e; /* Check for IEEE special cases. */ + /* log(NaN) = NaN. */ + if (isnanl (x)) + { + return x; + } /* log(0) = -infinity. */ if (x == 0.0L) { @@ -205,8 +208,8 @@ logl(long double x) { return (x - x) / ZERO; } - /* log (infinity or NaN) */ - if (x + x == x || x != x) + /* log (infinity) */ + if (x + x == x) { return x + x; }