Re: gsl_multiroot_test_delta

[Brian Gough <bjg@network-theory.co.uk>]

At Wed, 7 Nov 2007 11:04:51 +0100 (CET),
franky.backeljauw@........ wrote:
> When checking the one-dimensional case, the gsl_root_test_delta code 
> contains the test
> 
>       if (fabs(x1 - x0) < tolerance || x1 == x0)
> 
> The problem above can be eliminated by replacing the test with
> 
>       if (fabs(dxi) < tolerance || dxi == 0)
> 
> which is the generalization of the one-dimensional test to the 
> multidimensional case.

Thanks for the bug report.  I'll change the routine to allow dxi == 0.

-- 
best regards,

Brian Gough
(GSL Maintainer)

Network Theory Ltd,
Publishing the GSL Manual - http://www.network-theory.co.uk/gsl/manual/