On line 199 in siman/siman.c is:
} else if (gsl_rng_uniform(r) < exp (-(new_E - E)/(params.k * T)) ) {
which is the line that decides to take a step in the direction of higher
objective function energy sometimes. The problem is that my objective
function sometimes returns very large new_E for a small change in parameters,
leading to exp(-750) for example, and errno is set to ERANGE. Should this
condition be checked for?
_______________________________________________
Bug-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/bug-gsl
