bug-gsl
[Top][All Lists]
Advanced

[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/




reply via email to

[Prev in Thread] Current Thread [Next in Thread]