[Bug-gsl] Simulated annealing range problem

[Daniel Webb]

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?