def test_simulate(self): # Make system 1 x = SX.sym('x') y = SX.sym('y') t = SX.sym('t') tau = SX.sym('tau') ode = SX(1) alg = y - (5 - t) sys = DAESystem(x=x, y=y, ode=ode, alg=alg, tau=tau, t=t) # Test res = sys.simulate(x_0=0, t_f=5, t_0=0, y_0=1) self.assertAlmostEqual(res['xf'], 5) self.assertAlmostEqual(res['zf'], 0)
def test_simulate_with_options(self): # Make system 1 x = SX.sym('x') y = SX.sym('y') t = SX.sym('t') tau = SX.sym('tau') ode = SX(1) alg = y - (5 - t) sys = DAESystem(x=x, y=y, ode=ode, alg=alg, tau=tau, t=t) # Test res = sys.simulate(x_0=0, t_f=5, t_0=0, y_0=1, integrator_options={'abstol': 1e-10}) self.assertAlmostEqual(res['xf'], 5) self.assertAlmostEqual(res['zf'], 0)