[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] gsl_odeiv_evolve_apply can only decrease step size, not increa
From: |
Chris Rogers |
Subject: |
[Bug-gsl] gsl_odeiv_evolve_apply can only decrease step size, not increase |
Date: |
Thu, 13 Oct 2016 11:16:00 +0100 |
User-agent: |
Mozilla/5.0 (X11; Linux x86_64; rv:45.0) Gecko/20100101 Thunderbird/45.2.0 |
I am using gsl v1.15
I note that ODE evolver, gsl_odeiv_evolve_apply, defined in
gsl-1.15/ode-initval/evolve.c, can only decrease the step size in the
integration. If I define a control function that results in an
increasing step size, the evolver does not increase the step size.
For example, if the integration goes from a region where the derivative
changes quickly to one where the derivative changes slowly, the step
size should get larger. This is not possible, even if a control function
is defined that yields a large step size.
Reviewing the source code, there is a code block like
if (hadjust_status == GSL_ODEIV_HADJ_DEC) {
<some stuff>
}
but no code block like
if (hadjust_status == GSL_ODEIV_HADJ_INC) {
<some other stuff>
}
Let me know if I can help in any way resolve the issue.
Regards,
Chris Rogers
--
ASTeC Intense Beams Group
tel: +44 (0)1235 44 6983
Building R2,
Rutherford Appleton Laboratory,
Harwell Science and Innovation Campus,
Didcot.
OX11 0QX
- [Bug-gsl] gsl_odeiv_evolve_apply can only decrease step size, not increase,
Chris Rogers <=