[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Bug-gsl] bug in gsl_root_test_delta()

From: Brian Gough
Subject: Re: [Bug-gsl] bug in gsl_root_test_delta()
Date: Thu, 04 Jan 2007 11:46:13 +0000
User-agent: Mozilla/5.0 (X11; U; Linux i686; en-US; rv: Gecko/20060815 SeaMonkey/1.0.4

Brian Gough wrote:
Dave Spiegel writes:
 > I think that I have found a bug in the function gsl_root_test_delta.
> When x1 and x0 are both 0.0, gsl_root_test_delta returns -2, indicating > failure to converge.

Thanks for the bug report.  The result is a bit unintuitive for x1==0
but it is following the inequality in the documentation.

For roots close to zero, add a small epsabs term like GSL_DBL_MIN to
work around it.  I'll consider changing the behaviour in future to
make it more obvious.

This will be changed in the next release.  Thanks for the bug report.

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]