Problems with simulations of stiff models
-
Allow to specify options for
IQMPsimulate(reltol,abstol,maxnumsteps), otherwise stiff models simulation may fail instead of just taking longer time. -
Do not stop on
IQMPsimulateerrors, but instead:- try to increase maxnumsteps (e.g. twice by ten)
- allow to specify a fixed cost for non-computable simulations, e.g. as an offset wrt viability threshold given by its multiplier (
infvtmby default 0, i.e. not specified).
Sample from Claude attached; Km* < 1 values are too stiff for current 1e-12 tolerance and 1e6 max number of steps between two output points, but works with default values.
topo_MAPK.m paramSpecs.txt MAPK_model5.txtbc expData_design.xls exp-input.exp exp-alpha.exp