[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Question about the Documentation
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Question about the Documentation |
Date: |
Sat, 18 Aug 2007 07:49:59 +0100 |
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/22.1 Mule/5.0 (SAKAKI) |
At Thu, 16 Aug 2007 10:44:22 -0400,
Karl Edler wrote:
> I believe that intrinsic to the algorithm in gsl_fcmp is that it is not
> a reasonable thing to write something like gsl_fcmp(1e-16,0.0,1e-13) and
> expect to receive a 0 indicating that 1e-16 is approximately 0.0.
>
> (|u-v|<=|u|*e and |u-v|<=|v|*e )cannot be true if either u is 0.0 and e
> and v are not.
>
> I was hung up for a moment on this point and think that the gsl
> documentation for gsl_fcmp should address this point so that novice gsl
> programmers don't make the same mistake.
Hello,
Thanks for the suggestion. I have added the following note to the
documentation:
Note that @math{x} and @math{y} are compared to relative accuracy, so
this function is not suitable for testing whether a value is
approximately zero.
--
Brian Gough
(GSL Maintainer)
Network Theory Ltd,
Publishing the GSL Manual - http://www.network-theory.co.uk/gsl/manual/