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 _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
