step {rodeo}R Documentation

Numerical Integration Over a Single Time Step

Description

Performs integration over a single time step using a built-in ODE solver. At present, a single solver is implement with limited options. The interface of this method may change when support for other solvers is added.

Arguments

t0

Numeric. Initial time.

h

Numeric. Length of time step of interest.

hmin

Minimum tolerated internal step size. The default of NULL sets this to 10 times the value of .Machine$double.eps.

maxsteps

Maximum tolerated number of sub-steps.

tol

Numeric. Relative accuracy requested. This is currently a global value, i.e. one cannot set the accuracy per state variable.

method

String. Currently, 'rk5' is the only method implemented. This is a Runge-Kutta Cash-Karp solver adapted from Press et al. (2002), Numerical recipes in Fortran 90. It is designed to handle non-stiff problems only.

check

Logical. Can be used to avoid repeated checks of arguments. This may increase performance in repeated calls.

Value

A named numeric vector holding the values of state variables and process rates in all boxes.

Note

This method can only be used after a call to initStepper has been made.

Author(s)

david.kneis@tu-dresden.de

See Also

Use deSolve for advanced solvers with more options and capabilities to handle stiff problems.


[Package rodeo version 0.8.2 Index]