Daniel Webb writes:
> 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?
Thanks for the bug report. I will add a check to avoid that.
--
Brian Gough
Network Theory Ltd,
Publishing Free Software Manuals --- http://www.network-theory.co.uk/
_______________________________________________
Bug-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/bug-gsl
