This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: simplex minimization
- From: Tim F <fenn at agora dot dhs dot org>
- To: ivo dot alxneit at psi dot ch
- Cc: gsl-discuss at sources dot redhat dot com, <brian dot gough at network-theory dot co dot uk>
- Date: Mon, 6 Jan 2003 00:49:41 -0500 (EST)
- Subject: Re: simplex minimization
On Fri, 3 Jan 2003 ivo.alxneit@psi.ch wrote:
> > // set up solver
> > T = gsl_multimin_fminimizer_nmsimplex;
> > s = gsl_multimin_fminimizer_alloc(T, p);
> > ss = gsl_vector_alloc(p);
> > gsl_multimin_fminimizer_set(s, &f, ¶m.vector, ss);
> >
> you did not initialize ss! ss should be initialized by a reasonable step
> size for each of the parameters. this is needed to calculate n more
> corners (n parameters) for the initial simplex. the elements of ss seem
> to be all zero thus the initial size of the simplex is zero as well and
> the minimizer stop right away
>
doh! I knew it would be something dumb I was doing wrong... Thanks -
that fixed it.
Regards,
Tim F