def test_jansen_rit(self): """ """ model = models.JansenRit() history_shape = (1, model._nvar, 1, model.number_of_modes) model_ic = model.initial(dt, history_shape) self.assertEqual(model._nvar, 6)
def test_jansen_rit(self): """ """ model = models.JansenRit() self._validate_initialization(model, 6) self._test_steady_state( model, numpy.r_[1.105518e-02, 4.790847e+00, 8.826628e-01, 1.308187e-03, 2.367387e-01, 6.615113e-02])
def test_jansen_rit(self): """ """ model = models.JansenRit() history_shape = (1, model._nvar, 1, model.number_of_modes) model_ic = model.initial(dt, history_shape) self.assertEqual(model._nvar, 6) assert_array_almost_equal( model_ic, numpy.array([[[[-0.01381552]], [[-0.30892439]], [[-0.09769047]], [[-0.03384097]], [[-0.06178488]], [[-0.30892439]]]]))
def test_jansen_rit(self): """ """ model = models.JansenRit() history_shape = (1, model._nvar, 1, model.number_of_modes) model_ic = model.initial(dt, history_shape) self.assertEqual(model._nvar, 6) self.assertTrue(ArrayAlmostEqual(model_ic, numpy.array([[[[ 0.40556541]], [[ 2.52434922]], [[ 3.73943152]], [[ 3.04605971]], [[ 0.85500724]], [[ 4.27473643]]]])) is None)
def test_jansen_rit(self): """ """ model = models.JansenRit() self._validate_initialization(model, 6)
import tvb.simulator.simulator as simulator import tvb.simulator.models as models import tvb.simulator.noise as noise import tvb.simulator.coupling as coupling import tvb.simulator.integrators as integrators import tvb.simulator.monitors as monitors import tvb.datatypes.connectivity as connectivity ##----------------------------------------------------------------------------## ##- Perform the simulation -## ##----------------------------------------------------------------------------## LOG.info("Configuring...") #Initialise a Model, Coupling, and Connectivity. jrm = models.JansenRit() nsigma = 0.022 white_matter = connectivity.Connectivity() white_matter.speed = numpy.array([4.0]) white_matter_coupling = coupling.Linear(a=0.0) #Initialise an Integrator adding noise to only one state variable hiss = noise.Additive(nsig=numpy.array([0., 0., 0., 0., nsigma, 0.])) heunint = integrators.HeunStochastic(dt=2**-4, noise=hiss) #Initialise some Monitors with period in physical time momo = monitors.Raw() mama = monitors.TemporalAverage(period=2**-2)