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)