The initialization for the init_diam variable in nldrmd.c is currently
if (init_diam == 0) /* initialize diam. for psi convergence test */
for (i = 0; i < n; ++i) init_diam = fabs(xl[i] - xh[i]);
The loop has the same effect as
init_diam = fabs(xl[n] - xh[n]);
I imagine that the loop was intended to be
for (i = 0; i < n; ++i) init_diam += fabs(xl[i] - xh[i]);
_______________________________________________
NLopt-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/nlopt-discuss
