def setUp(self):
"""We need a Prior object and an IVP object (with derivatives) to run the
tests."""
y0 = pnrv.Constant(np.array([20.0, 15.0]))
self.ivp = pnd.lotkavolterra([0.4124, 1.15124], y0)
self.prior = pnfs.statespace.IBM(ordint=2, spatialdim=2)
self.evlvar = 0.0005123121
def load_lotkavolterra():
"""Load LV system as a basic IVP."""
initrv = Dirac(np.array([20, 20]))
return lotkavolterra(timespan=[0, 0.55],
initrv=initrv,
params=(0.5, 0.05, 0.5, 0.05))
def setUp(self):
"""We need a Prior object and an IVP object (with derivatives) to run the tests."""
y0 = Dirac(np.array([20.0, 15.0]))
self.ivp = lotkavolterra([0.4124, 1.15124], y0)
self.prior = IBM(ordint=2, spatialdim=2, diffconst=1.7685)
self.evlvar = 0.0005123121
def ivp():
y0 = random_variables.Constant(np.array([20.0, 15.0]))
return diffeq.lotkavolterra([0.4124, 1.15124], y0)
def lv():
y0 = pnrv.Constant(np.array([20.0, 20.0]))
# tmax is ignored anyway
return pnde.lotkavolterra([0.0, np.inf], y0)
def case_lotkavolterra():
y0 = np.array([0.1, 0.1])
ode = diffeq.lotkavolterra([0.0, 1.0], y0)
return setup_solver(y0, ode)
