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
IQMPsimulate
errors, 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 (
infvtm
by 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