[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Bug-gsl] [bug #52336] odeiv2/msbdf crashes with method order 3 and maxi

From: Tuomo Keskitalo
Subject: [Bug-gsl] [bug #52336] odeiv2/msbdf crashes with method order 3 and maximum step length increase
Date: Sat, 4 Nov 2017 03:50:15 -0400 (EDT)
User-agent: Mozilla/5.0 (X11; Ubuntu; Linux x86_64; rv:56.0) Gecko/20100101 Firefox/56.0

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

(file #42326)

Additional Item Attachment:

File name: benson_debug.txt               Size:171 KB


Reply to this item at:


  Message sent via/by Savannah

reply via email to

[Prev in Thread] Current Thread [Next in Thread]