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)
