Re: [Bug-gsl] Simulated annealing range problem

[Brian Gough]

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/