Example #1
0
 def setUp(self):
     initrv = Normal(20 * np.ones(2),
                     0.1 * np.eye(2),
                     cov_cholesky=np.sqrt(0.1) * np.eye(2))
     self.ivp = lotkavolterra([0.0, 0.5], initrv)
     step = 0.1
     self.solution = probsolve_ivp(self.ivp, step=step)
Example #2
0
 def setUp(self):
     initrv = Constant(20 * np.ones(2))
     self.ivp = lotkavolterra([0.0, 0.5], initrv)
     step = 0.1
     f = self.ivp.rhs
     t0, tmax = self.ivp.timespan
     y0 = self.ivp.initrv.mean
     self.solution = probsolve_ivp(f,
                                   t0,
                                   tmax,
                                   y0,
                                   step=step,
                                   adaptive=False)
Example #3
0
    def setUp(self):
        initrv = Normal(20 * np.ones(2),
                        0.1 * np.eye(2),
                        cov_cholesky=np.sqrt(0.1) * np.eye(2))
        self.ivp = lotkavolterra([0.0, 0.5], initrv)
        step = 0.1
        f = self.ivp.rhs
        t0, tmax = self.ivp.timespan
        y0 = self.ivp.initrv.mean

        self.solution = probsolve_ivp(f,
                                      t0,
                                      tmax,
                                      y0,
                                      step=step,
                                      adaptive=False)
Example #4
0
 def setUp(self):
     initrv = Constant(20 * np.ones(2))
     self.ivp = lotkavolterra([0.0, 0.5], initrv)
     step = 0.1
     self.solution = probsolve_ivp(self.ivp, which_prior="ibm3", step=step)
Example #5
0
def ivp():
    initrv = Constant(20.0 * np.ones(2))
    return ode.lotkavolterra([0.0, 0.25], initrv)
Example #6
0
 def setUp(self):
     initrv = Dirac(20 * np.ones(2))
     self.ivp = lotkavolterra([0.0, 0.5], initrv)
     step = 0.1
     self.solution = probsolve_ivp(self.ivp, step=step)
Example #7
0
 def setUp(self):
     initdist = Constant(20 * np.ones(2))
     self.ivp = ode.lotkavolterra([0.0, 1e-4], initdist)
     self.step = 1e-5
     self.prior = "ibm3"
Example #8
0
 def setUp(self):
     """Setup odesolver and Lotka-Volterra IVP."""
     initdist = Constant(20 * np.ones(2))
     self.ivp = ode.lotkavolterra([0.0, 0.5], initdist)
     self.tol = 1e-1
     self.step = 0.1
Example #9
0
 def setUp(self):
     """Setup odesolver and solve a scalar ode."""
     initrv = Constant(20 * np.ones(2))
     self.ivp = ode.lotkavolterra([0.0, 0.5], initrv)
     self.tol = 1e-2
Example #10
0
 def setUp(self):
     """Setup odesolver and Lotka-Volterra IVP"""
     initrv = Dirac(20 * np.ones(2))
     self.ivp = ode.lotkavolterra([0.0, 0.5], initrv)
     self.tol = 1e-1
     self.step = 0.1