search

KineticPreProcessor / KPP

The KPP kinetic preprocessor is a software tool that assists the computer simulation of chemical kinetic systems
GNU General Public License v3.0
19 stars 11 forks source link

SDIRK solver convergence #79

Open maqsoodrajput opened 1 year ago

maqsoodrajput commented 1 year ago

I am trying to check if the SDIRK solvers in kpp/int/ (with fixed step-size) give correct order of convergence. In all cases (second/third/fourth order solver) I am getting an order ~1 which is strange.

RolfSander commented 1 year ago

I'm not familiar with SDIRK. I almost always use the Rosenbrock integrator. Anyway, how do you determine the order of convergence? Is there a variable in sdirk.f90, or are you calculating it yourself?