[Bug-gsl] [bug #52336] odeiv2/msbdf crashes with method order 3 and maximum step length increase

[Tuomo Keskitalo]

Follow-up Comment #1, bug #52336 (project gsl):

Added debug output from Andrew to bug report.

It looks like you have found a rare collision. The errors of your function
just before crashing are suddely all smaller continuously, so controller is
accelerating and finally reaches maximum acceleration, which hits the
singularity of the solver at ord=3 as you described.

Also you are right about the easy cure: r could be limited in cscal.c to
something else than 5.0. I think that 4.9 is better than 5.1, because it
eliminates possibility for hitting 5.0 accidentally. The maximum value 5.0 is
anyway more or less arbitrary, so I'm not worried about changing it in any
case.


(file #42326)
    _______________________________________________________

Additional Item Attachment:

File name: benson_debug.txt               Size:171 KB


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?52336>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/