Update of bug #59710 (project gsl):
Open/Closed: Open => Closed
_______________________________________________________
Follow-up Comment #1:
Thanks for your report. The problem is that on the first iteration,
the lambda value jumps to -12.8, and the value of exp(12.8 * 60) is too
large to represent in double precision, so infs and nans appear in the
residual vector and Jacobian matrix.
For this particular problem, the value of lambda should be > 0, however
there is no constraint placed to assure that, so when the first
iteration jumps to -12.8, bad numerical things happen.
I think the best solution for now is to adjust the model parameters a
little to prevent this type of behavior, and bring TMAX closer to 1
(instead of 40). A better solution would be to implement a constraint so
that lambda > 0, but this would require a lot of development effort.
Fixed in commit 7caac3346
_______________________________________________________
Reply to this item at:
<https://savannah.gnu.org/bugs/?59710>
_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/
