def buildsim(self, ):
"""
Setup the assimulo IDA simulator.
"""
# Create an Assimulo implicit solver (IDA)
imp_sim = IDA(self.imp_mod) # Create a IDA solver
# Sets the paramters
# 1e-4 #Default 1e-6
imp_sim.atol = self.p.RunInput['TIMESTEPPING']['SOLVER_TOL']
# 1e-4 #Default 1e-6
imp_sim.rtol = self.p.RunInput['TIMESTEPPING']['SOLVER_TOL']
# Suppres the algebraic variables on the error test
imp_sim.suppress_alg = True
imp_sim.display_progress = False
imp_sim.verbosity = 50
imp_sim.report_continuously = True
imp_sim.time_limit = 10.
self.imp_sim = imp_sim
#imp_sim.make_consistent('IDA_YA_YDP_INIT')
#ta, ya, yda = imp_sim.simulate(0.1,5)
##
#imp_mod.set_iapp( I_app/2. )
#imp_sim.make_consistent('IDA_YA_YDP_INIT')
#tb, yb, ydb = imp_sim.simulate(0.2,5)
#imp_mod.set_iapp( I_app )
#imp_sim.make_consistent('IDA_YA_YDP_INIT')
## Sim step 1
#t1, y1, yd1 = imp_sim.simulate(1./Crate*3600.*0.2,100)
imp_sim.display_progress = False
imp_sim.verbosity = 50
imp_sim.report_continuously = True
imp_sim.time_limit = 10.
### Simulate
t01, t02 = 0.1, 0.2
imp_mod.set_iapp(I_app / 10.)
imp_sim.make_consistent('IDA_YA_YDP_INIT')
ta, ya, yda = imp_sim.simulate(t01, 2)
imp_mod.set_iapp(I_app / 2.)
imp_sim.make_consistent('IDA_YA_YDP_INIT')
tb, yb, ydb = imp_sim.simulate(t02, 2)
print 'yb shape', yb.shape
# Sim step 1
