def make_sim(sim_len=1000.0):
sim = simulator.Simulator(
connectivity=connectivity.Connectivity.from_file(),
model=models.MontbrioPazoRoxin(),
integrator=integrators.HeunStochastic(
dt=0.1, noise=noise.Additive(nsig=np.r_[0.001])),
monitors=[monitors.Raw()],
simulation_length=sim_len)
sim.configure()
return sim
def test_infinite_theta(self):
model = models.MontbrioPazoRoxin()
self._validate_initialization(model, 2)
model = models.CoombesByrne()
self._validate_initialization(model, 4)
model = models.CoombesByrne2D()
self._validate_initialization(model, 2)
model = models.GastSchmidtKnosche_SD()
self._validate_initialization(model, 4)
model = models.GastSchmidtKnosche_SF()
self._validate_initialization(model, 4)
model = models.DumontGutkin()
self._validate_initialization(model, 8)
def _tvb_sim(self, dt):
import numpy as np
from tvb.simulator.simulator import Simulator, models, connectivity, monitors
sim = Simulator(
connectivity=connectivity.Connectivity.from_file(),
model=models.MontbrioPazoRoxin(I=np.r_[1.0],
Delta=np.r_[1.0],
eta=np.r_[-5.0],
tau=np.r_[100.0],
J=np.r_[15.],
cr=np.r_[0.01],
cv=np.r_[0.0]),
integrator=self._integrator(dt=dt),
monitors=[monitors.Raw()],
).configure()
assert sim.history.buffer.shape == (513, 2, 76, 1)
sim.history.buffer[:, 0] = 0.0
sim.history.buffer[:, 1] = -2.0
sim.current_state[:] = sim.history.buffer[0]
return sim