Integrator status (ISTATUS and RSTATUS)

In order to obtain more information about the integration, KPP provides the arrays ISTATUS (integ mer) and RSTATUS (real). Each of them is an array of 20 elements. Array elements not listed here are currently not used. Details can be found in the comment lines of the individual integrator files in $KPP_HOME/int/.

ISTATUS

Summary of ISTATUS usage in the f90 integrators. Here, Y = used.

ISTATUS

1

2

3

4

5

6

7

8

9

beuler

Y

Y

Y

Y

Y

Y

Y

Y

dvode

exponential

feuler

gillespie

lsode

Y

Y

Y

radau5

Y

Y

Y

Y

Y

Y

Y

Y

rosenbrock

Y

Y

Y

Y

Y

Y

Y

Y

rosenbrock_adj

Y

Y

Y

Y

Y

Y

Y

Y

rosenbrock_autoreduce

Y

Y

Y

Y

Y

Y

Y

Y

rosenbrock_h211b_qssa

Y

Y

Y

Y

Y

Y

Y

Y

rosenbrock_tlm

Y

Y

Y

Y

Y

Y

Y

Y

Y

runge_kutta

Y

Y

Y

Y

Y

Y

Y

Y

runge_kutta_adj

Y

Y

Y

Y

Y

Y

Y

Y

runge_kutta_tlm

Y

Y

Y

Y

Y

Y

Y

Y

sdirk

Y

Y

Y

Y

Y

Y

Y

Y

sdirk_adj

Y

Y

Y

Y

Y

Y

Y

Y

sdirk_tlm

Y

Y

Y

Y

Y

Y

Y

Y

sdirk4

Y

Y

Y

Y

Y

Y

Y

Y

seulex

Y

Y

Y

Y

Y

Y

Y

tau_leap

ISTATUS(1)

Number of times that the function containing the chemical ODE system was called.

ISTATUS(2)

Number times the Jacobian matrix was constructed.

ISTATUS(3)

Total number of integration timesteps.

ISTATUS(4)

Number of accepted timesteps.

ISTATUS(5)

Number of rejected timesteps (except at very beginning).

ISTATUS(6)

Number of LU decompositions that were performed.

ISTATUS(7)

Number of forward/backward substitutions that were performed.

ISTATUS(8)

Number of singular matrix decompositions that were performed.

ISTATUS(9)

Number of times the Hessian matrix was evaluated.

ISTATUS(10) .. ISTATUS(20)

Currently not used.

RSTATUS

Summary of RSTATUS usage in the f90 integrators. Here, Y = used, s = solver specific usage.

RSTATUS

1

2

3

4

beuler

Y

Y

Y

dvode

exponential

feuler

Y

gillespie

lsode

Y

Y

radau5

rosenbrock

Y

Y

Y

rosenbrock_adj

Y

Y

Y

rosenbrock_autoreduce

Y

Y

Y

s

rosenbrock_h211b_qssa

Y

Y

Y

rosenbrock_tlm

Y

Y

Y

runge_kutta

Y

Y

Y

runge_kutta_adj

Y

Y

Y

runge_kutta_tlm

Y

Y

Y

sdirk

Y

Y

Y

sdirk_adj

Y

Y

Y

sdirk_tlm

Y

Y

Y

sdirk4

Y

Y

seulex

tau_leap

RSTATUS(1)

Texit, the time corresponding to the computed \(Y\) upon return.

RSTATUS(2)

Hexit: the last accepted step before exit.

RSTATUS(3)

Hnew: The last predicted step (not yet taken. For multiple restarts, use Hnew as Hstart in the subsequent run.

RSTATUS(4)

(Solver-specific for rosenbrock_autoreduce) AR_thr: used to output the calculated (used) auto-reduction threshold for the integration. Useful when ICNTRL(10) > 0 where the threshold is dynamically determined based on a given species.

RSTATUS(5) .. RSTATUS(20)

Currently not used.