[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] min/brent.c: Implementation of Brent's algorithm
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] min/brent.c: Implementation of Brent's algorithm |
Date: |
Thu, 28 Apr 2005 13:02:20 +0100 |
Peter Brommer writes:
> In your algorithm, the checks that lead to C and D are only executed if
> (f_u == f_z), whereas they are executed if (f_u > f_z) in the original
> implementation.
Thanks for the bug report. It was due to an error in transcribing the
code from the original fortran version of Brent's routine. I will fix
it in the next release.
--
best regards
Brian Gough