Re: [Bug-gsl] msadams.c:973 aborts with GSL_ESANITY

[Tuomo Keskitalo]

Hello,

On 05.10.2017 21:44, Michael Kaufman wrote:

> t:5115.96 count:41 failed:10 order:5
> t:5116.00 count:43 failed:11 order:5
> t:5116.02 count:46 failed:13 order:5
> t:5116.04 count:48 failed:14 order:5
> t:5116.05 count:50 failed:15 order:5
> t:5116.06 count:52 failed:16 order:5
> t:5116.07 count:53 failed:16 order:5
> t:5116.08 count:54 failed:16 order:5
> t:5116.13 count:55 failed:16 order:5
> t:5116.26 count:56 failed:16 order:5
> -- order change -2
> t:5116.26 count:63 failed:22 order:1
> t:5116.27 count:64 failed:22 order:1

Thanks! Looks like the solver hits a bump while accelerating on a high 
gear. That should not be a problem, though. I'll take a look at the 
code, but it will take time.

Meanwhile, if you have time, you can add "#define DEBUG" in the 
beginning of msadams.c and recompile GSL. This should make msadams print 
out _lots_ of debug info and data, that would help decipher this issue. 
Interesting stuff happens between t=5116.13 and t=5116.27.

--
address@hidden
http://iki.fi/tuomo.keskitalo