Rosenbrock method of order (3)4 with step-size control and continuous output.
Based on the FORTRAN code RODAS by E.Hairer and G.Wanner, which can be found here: http://www.unige.ch/~hairer/software.html
Details about the implementation (FORTRAN) can be found in the book,:
Solving Ordinary Differential Equations II,
Stiff and Differential-Algebraic Problems
Authors: E. Hairer and G. Wanner
Springer-Verlag, ISBN: 3-540-60452-9
Import the solver together with the correct problem:
from assimulo.solvers import RodasODE
from assimulo.problem import Explicit_Problem
Define the problem, such as:
def rhs(t, y): #Note that y are a 1-D numpy array.
yd = -1.0
return N.array([yd]) #Note that the return must be numpy array, NOT a scalar.
y0 = [1.0]
t0 = 1.0
Create a problem instance:
mod = Explicit_Problem(rhs, y0, t0)
Note
For complex problems, it is recommended to check the available examples and the documentation in the problem class, Explicit_Problem
. It is also recommended to define your problem as a subclass of Explicit_Problem
.
Warning
When subclassing from a problem class, the function for calculating the right-hand-side (for ODEs) must be named rhs and in the case with a residual function (for DAEs) it must be named res.
Create a solver instance:
sim = RodasODE(mod)
Modify (optionally) the solver parameters.
Parameters:
atol
Defines the absolute tolerance(s) that is to be used by the solver.backward
Specifies if the simulation is done in reverse time.clock_step
Specifies if the elapsed time of an integrator step should be timed or not.display_progress
This option actives output during the integration in terms of that the current integration is periodically printed to the stdout.fac1
Parameters for step-size selection.fac2
Parameters for step-size selection.inith
This determines the initial step-size to be used in the integration.maxh
Defines the maximal step-size that is to be used by the solver.maxsteps
The maximum number of steps allowed to be taken to reach the final time.num_threads
This options specifies the number of threads to be used for those solvers that supports it.report_continuously
This options specifies if the solver should report the solution continuously after steps.rtol
Defines the relative tolerance that is to be used by the solver.safe
The safety factor in the step-size prediction.store_event_points
This options specifies if the solver should save additional points at the events, \(t_e^-, t_e^+\).time_limit
This option can be used to limit the time of an integration.usejac
This sets the option to use the user defined Jacobian.verbosity
This determines the level of the output.
Methods:
RodasODE.interpolate
Simulate the problem:
Information:
RodasODE.get_options()
Returns the current solver options.RodasODE.get_supports()
Returns the functionality which the solver supports.RodasODE.get_statistics()
Returns the run-time statistics (if any).RodasODE.get_event_data()
Returns the event information (if any).RodasODE.print_event_data()
Prints the event information (if any).RodasODE.print_statistics()
Prints the run-time statistics for the problem.