It seems to just shrink the simplex and then exit. I think it doesn't
manage to find the narrow trench in configuration space in which
possible improvements lie.
Also tried gsl_multimin_fminimizer_nmsimplex but this is even slower
than siman. So we will go back to siman now.
Signed-off-by: Ian Jackson <ijackson@chiark.greenend.org.uk>
oc._dbg(repr(fc_input))
- findcurve_epsilon = 0.0001
+ findcurve_epsilon = 1E-10
cl = ['./findcurve', '%d' % (nt+1), '%.18g' % findcurve_epsilon]
oc._dbg('STARTING FINDCURVE %s' % cl)
gsl_vector_set_all(step_size, 0.01);
gsl_multimin_fminimizer *minimiser =
- gsl_multimin_fminimizer_alloc(gsl_multimin_fminimizer_nmsimplex2, NX);
+ gsl_multimin_fminimizer_alloc(gsl_multimin_fminimizer_nmsimplex2rand, NX);
for (;;) {
/* NINPUT + 1 doubles: startpoint, epsilon for residual */