X-Git-Url: http://erislabs.net/gitweb/?a=blobdiff_plain;f=lib%2Fexp2.c;h=23e5f07b40daebeee1b13b92d6fd777e8150184e;hb=fa1db0dd22768f09a507674a30beb5b8a87bb35f;hp=21074939402753876b9149d0c82de70d0449b448;hpb=98d50dd6077362e1bc9d9bdeb4193051e611e677;p=gnulib.git diff --git a/lib/exp2.c b/lib/exp2.c index 210749394..23e5f07b4 100644 --- a/lib/exp2.c +++ b/lib/exp2.c @@ -1,5 +1,5 @@ /* Exponential base 2 function. - Copyright (C) 2012 Free Software Foundation, Inc. + Copyright (C) 2012-2013 Free Software Foundation, Inc. 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 @@ -83,9 +83,9 @@ exp2 (double x) + 21844/6081075 * z^13 - 929569/638512875 * z^15 + ... - Since |z| <= log(2)/1024 < 0.0007, the relative error of the z^7 term - is < 0.0007^6 < 2^-60 <= 2^-DBL_MANT_DIG, therefore we can truncate - the series after the z^5 term. */ + Since |z| <= log(2)/1024 < 0.0007, the relative contribution of the + z^7 term is < 0.0007^6 < 2^-60 <= 2^-DBL_MANT_DIG, therefore we can + truncate the series after the z^5 term. */ { double nm = round (x * 256.0); /* = 256 * n + m */