def setUp(self):
self.file = tempfile.NamedTemporaryFile()
self.path = self.file.name
self.family = os.path.basename(self.path)
self.prefix = os.path.dirname(self.path)
self.file.close()
from odelab.scheme.classic import ExplicitEuler
self.params = params = {
'family': self.family,
'system': System,
'system_params': {'f': f},
'solver': Solver,
'scheme': ExplicitEuler,
'scheme_params': {
'h': .1,
},
'initialize': {
'u0' : np.array([1.]),
'name': self.name,
},
}
scheme = ExplicitEuler(**params['scheme_params'])
s = Solver(system=System(f), scheme=scheme, path=self.path)
s.catch_runtime = False
s.initialize(**params['initialize'])
s.run(1.)
def setUp(self):
self.file = tempfile.NamedTemporaryFile()
self.path = self.file.name
self.family = os.path.basename(self.path)
self.prefix = os.path.dirname(self.path)
self.file.close()
from odelab.scheme.classic import ExplicitEuler
self.params = params = {
'family': self.family,
'system': System,
'system_params': {
'f': f
},
'solver': Solver,
'scheme': ExplicitEuler,
'scheme_params': {
'h': .1,
},
'initialize': {
'u0': np.array([1.]),
'time': 1.,
'name': self.name,
},
}
scheme = ExplicitEuler(**params['scheme_params'])
s = Solver(system=System(f), scheme=scheme, path=self.path)
s.catch_runtime = False
s.initialize(**params['initialize'])
s.run()
def test_rkdae(s):
sys = GraphSystem(sq)
u0 = np.array([0.,0.,1.])
scheme = RKDAE(.1, tableau=RK.RadauIIA.tableaux[s])
sol = Solver(scheme, sys)
sol.initialize(u0=u0)
CompareExact('RadauIIA-{0}'.format(s)), sol, u0, 2
