[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Nelder Mead bug
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Nelder Mead bug |
Date: |
Tue, 06 May 2008 21:53:56 +0100 |
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/22.1 Mule/5.0 (SAKAKI) |
At Sun, 4 May 2008 21:23:27 -0400,
Michael Kuklik wrote:
> Nelder Mead , i.e. simplex search has a bug. In the function
> static int
> nmsimplex_iterate (void *vstate, gsl_multimin_function * f,
> gsl_vector * x, double *size, double *fval)
>
> in multimin/simplex.c
> a piece of code that determines the highest , second highest and lowest
> vertex fails when first value in y1 is the highest.
> It returns second highest equal to the highest. Later that influence
> decision in the algorithm.
>
> FIX:
> add lines:
> ds_hi = gsl_vector_get (y1, 1);
> s_hi = 1;
> before the loop
>
Thanks for the bug report, you're right. I will commit a fix for
that.
--
Brian Gough
GNU Scientific Library -
http://www.gnu.org/software/gsl/