[MATH-563] resetting relative and absolute tolerances in ODE adaptive step size integrators

• Type: Improvement
• Status: Closed
• Priority: Trivial
• Resolution: Fixed
• Affects Version/s: 2.2
• Fix Version/s: 3.0
• Assignee: Luc Maisonobe
• Reporter: Luc Maisonobe
• Created: 25/Apr/11 09:11
• Updated: 24/Mar/12 16:16
• Resolved: 25/Apr/11 09:46

Description

There is a limitation of the current implementation of ODE integrators with
adaptive step size. For now, the tolerances that are used to adjust the step
size are specified only at construction time and cannot be changed afterwards.
However, these tolerances are highly problem-dependent and in fact the dimension
of the problem (which is related to the dimension of the vectorial version of
the tolerances) is specified only at integration time, not at construction time.

So it would be interesting to add at the top-level hierarchy (abstract class
AdaptiveStepsizeIntegrator) a few setters to allow users to change these
tolerances after the integrator has been built. It seems the integrators by
themselves were not documented as immutable (I first thought they were), so this
change is probably harmless.

There was a discussion on the developers mailing list about this feature here: http://markmail.org/thread/djrjpm5pwcgg3ev4