[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] msadams.c:973 aborts with GSL_ESANITY
|
From: |
Tuomo Keskitalo |
|
Subject: |
Re: [Bug-gsl] msadams.c:973 aborts with GSL_ESANITY |
|
Date: |
Thu, 28 Sep 2017 19:08:57 +0300 |
|
User-agent: |
Mozilla/5.0 (X11; Linux x86_64; rv:52.0) Gecko/20100101 Thunderbird/52.3.0 |
Hello,
On 26.09.2017 21:36, Michael Kaufman wrote:
gsl: msadams.c:973: ERROR: msadams_apply too large order change
Default GSL error handler invoked.
I'm running on OSX 10.11.6 using a gsl-2.4_1 from macports.
This might be compiler optimization issue. Can you please try with
non-optimized (-O0) version of GSL? I have no idea how to get that with
macports... Alternatively, can you try to run your program on Linux?
Strangely, if I solve the ODE with and abstol of 1e-8 and then
immediately afterward solve again with an abstol of 1e-10, (for
convergence testing) then the abort happens (call this case 1).
How exactly do you "immediately afterward solve again"? Do you call
gsl_odeiv2_driver_reset in between calls to gsl_odeiv2_driver_apply?
BR,
Tuomo
--
address@hidden
http://iki.fi/tuomo.keskitalo