I want to port a model from use with the GUI to a batch mode controlled by Python. Everything seems to work except for the variable timestepping settings which were done using the AtolTool (see viewtopic.php?f=28&t=3232). I have the following code in my .ses file:
Code: Select all
objectvar save_window_, rvp_
objectvar scene_vector_[6]
objectvar ocbox_, ocbox_list_, scene_, scene_list_
{ocbox_list_ = new List() scene_list_ = new List()}
//Begin VariableTimeStep
{
ocbox_ = NumericalMethodPanel[0]
}
{object_push(ocbox_)}
{
atol_ = 0.001 CVode[0].atol(atol_)
restore(412, 1)
atoltool_ = new AtolTool()
ats("v", 10)
ats("vext", 10)
ats("mp_axnode", 0.1)
ats("m_axnode", 0.1)
ats("h_axnode", 0.1)
ats("s_axnode", 0.01)
ats("Vector", -1)
atoltool_.scales()
}
{object_pop()}
{
ocbox_.map("VariableTimeStep", 12, 990, 287.1, 342.9)
}
objref ocbox_
//End VariableTimeStep
objectvar scene_vector_[1]
{doNotify()}
Code: Select all
cvode.atol(0.001)
cvode.atolscale("v", 10)
cvode.atolscale("vext", 10)
cvode.atolscale("mp_axnode", 0.1)
cvode.atolscale("m_axnode", 0.1)
cvode.atolscale("h_axnode", 0.1)
cvode.atolscale("s_axnode", 0.01)
//cvode.atolscale("Vector", -1)
So what is the correct way to port it? A negative tolerance does not make sense, but what does "-1" mean in this case?
Thank you for hints and explanations...
Matthias