Skip to content

Instantly share code, notes, and snippets.

@asoplata

asoplata/test_network_cx_limit.py

Created May 20, 2025 20:25
Show Gist options

  • Select an option

    Learn more about clone URLs
  • Save asoplata/c1635c3f6debbe00a4f1b548372f65d2 to your computer and use it in GitHub Desktop.

Select an option

Learn more about clone URLs
Save asoplata/c1635c3f6debbe00a4f1b548372f65d2 to your computer and use it in GitHub Desktop.
Download ZIP
Test of HNN-Core network complexity limit (2025-05-20)
Raw
test_network_cx_limit.py
# ATTENTION: See which simulations fail and succeed at the bottom of the script.
import os.path as op
import matplotlib.pyplot as plt
from hnn_core import jones_2009_model, simulate_dipole
from hnn_core.viz import plot_dipole
import matplotlib.pyplot as plt
from copy import deepcopy
import pandas as pd
import pickle
import numpy as np
from hnn_core import MPIBackend
net = jones_2009_model()
weights_ampa_d1 = {'L2_basket': 0.006562, 'L2_pyramidal': .000007,
'L5_pyramidal': 0.142300}
weights_nmda_d1 = {'L2_basket': 0.019482, 'L2_pyramidal': 0.004317,
'L5_pyramidal': 0.080074}
synaptic_delays_d1 = {'L2_basket': 0.1, 'L2_pyramidal': 0.1,
'L5_pyramidal': 0.1}
net.add_evoked_drive(
'evdist1', mu=63.53, sigma=3.85, numspikes=1, weights_ampa=weights_ampa_d1,
weights_nmda=weights_nmda_d1, location='distal',
synaptic_delays=synaptic_delays_d1, event_seed=274)
net.add_evoked_drive(
'evdist2', mu=63.53, sigma=3.85, numspikes=1, weights_ampa=weights_ampa_d1,
weights_nmda=weights_nmda_d1, location='distal',
synaptic_delays=synaptic_delays_d1, event_seed=274)
net.add_evoked_drive(
'evdist3', mu=63.53, sigma=3.85, numspikes=1, weights_ampa=weights_ampa_d1,
weights_nmda=weights_nmda_d1, location='distal',
synaptic_delays=synaptic_delays_d1, event_seed=274)
weights_ampa_p1 = {'L2_basket': 0.08831, 'L2_pyramidal': 0.01525,
'L5_basket': 0.19934, 'L5_pyramidal': 0.00865}
synaptic_delays_prox = {'L2_basket': 0.1, 'L2_pyramidal': 0.1,
'L5_basket': 1., 'L5_pyramidal': 1.}
# all NMDA weights are zero; pass None explicitly
net.add_evoked_drive(
'evprox1', mu=26.61, sigma=2.47, numspikes=1, weights_ampa=weights_ampa_p1,
weights_nmda=None, location='proximal',
synaptic_delays=synaptic_delays_prox, event_seed=544)
# Second proximal evoked drive. NB: only AMPA weights differ from first
weights_ampa_p2 = {'L2_basket': 0.000003, 'L2_pyramidal': 1.438840,
'L5_basket': 0.008958, 'L5_pyramidal': 0.684013}
# all NMDA weights are zero; omit weights_nmda (defaults to None)
net.add_evoked_drive(
'evprox2', mu=137.12, sigma=8.33, numspikes=1,
weights_ampa=weights_ampa_p2, location='proximal',
synaptic_delays=synaptic_delays_prox, event_seed=814)
net.add_evoked_drive(
'evprox3', mu=137.12, sigma=8.33, numspikes=1,
weights_ampa=weights_ampa_p2, location='proximal',
synaptic_delays=synaptic_delays_prox, event_seed=814)
nets=[]
for isi in [90]:
# AES
for t in [15]:
r_net = deepcopy(net)
# copy evprox 1
r_net.add_evoked_drive('evprox4',
mu=r_net.external_drives['evprox1']['dynamics']['mu'] + t,
sigma=r_net.external_drives['evprox1']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox1']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox1']['synaptic_delays'])
# copy evprox 2
r_net.add_evoked_drive('evprox5',
mu=r_net.external_drives['evprox2']['dynamics']['mu'] + t,
sigma=r_net.external_drives['evprox2']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox2']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox2']['synaptic_delays'])
# copy evprox 3
r_net.add_evoked_drive('evprox6',
mu=r_net.external_drives['evprox3']['dynamics']['mu'] + t,
sigma=r_net.external_drives['evprox3']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox3']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox3']['synaptic_delays'])
# copy evdist 1
r_net.add_evoked_drive('evdist4',
mu=r_net.external_drives['evdist1']['dynamics']['mu'] + t,
sigma=r_net.external_drives['evdist1']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist1']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist1']['synaptic_delays'])
# copy evdist 2
r_net.add_evoked_drive('evdist5',
mu=r_net.external_drives['evdist2']['dynamics']['mu'] + t,
sigma=r_net.external_drives['evdist2']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist2']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist2']['synaptic_delays'])
# copy evdist 3
r_net.add_evoked_drive('evdist6',
mu=r_net.external_drives['evdist3']['dynamics']['mu'] + t,
sigma=r_net.external_drives['evdist3']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist3']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist3']['synaptic_delays'])
#---------------
# copy evprox 1
r_net.add_evoked_drive('evprox7',
mu=r_net.external_drives['evprox1']['dynamics']['mu'] + 2*t,
sigma=r_net.external_drives['evprox1']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox1']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox1']['synaptic_delays'])
# copy evprox 2
r_net.add_evoked_drive('evprox8',
mu=r_net.external_drives['evprox2']['dynamics']['mu'] + 2*t,
sigma=r_net.external_drives['evprox2']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox2']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox2']['synaptic_delays'])
# copy evprox 3
r_net.add_evoked_drive('evprox9',
mu=r_net.external_drives['evprox3']['dynamics']['mu'] + 2*t,
sigma=r_net.external_drives['evprox3']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox3']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox3']['synaptic_delays'])
# copy evdist 1
r_net.add_evoked_drive('evdist7',
mu=r_net.external_drives['evdist1']['dynamics']['mu'] + 2*t,
sigma=r_net.external_drives['evdist1']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist1']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist1']['synaptic_delays'])
# copy evdist 2
r_net.add_evoked_drive('evdist8',
mu=r_net.external_drives['evdist2']['dynamics']['mu'] + 2*t,
sigma=r_net.external_drives['evdist2']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist2']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist2']['synaptic_delays'])
# copy evdist 3
r_net.add_evoked_drive('evdist9',
mu=r_net.external_drives['evdist3']['dynamics']['mu'] + 2*t,
sigma=r_net.external_drives['evdist3']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist3']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist3']['synaptic_delays'])
# ------- SECOND BURST -------
# copy evprox 1
r_net.add_evoked_drive('evprox10',
mu=r_net.external_drives['evprox1']['dynamics']['mu'] + 2 * t + isi,
sigma=r_net.external_drives['evprox1']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox1']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox1']['synaptic_delays'])
# copy evprox 2
r_net.add_evoked_drive('evprox11',
mu=r_net.external_drives['evprox2']['dynamics']['mu'] + 2 * t + isi,
sigma=r_net.external_drives['evprox2']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox2']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox2']['synaptic_delays'])
# copy evprox 3
r_net.add_evoked_drive('evprox12',
mu=r_net.external_drives['evprox3']['dynamics']['mu'] + 2* t + isi,
sigma=r_net.external_drives['evprox3']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox3']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox3']['synaptic_delays'])
# copy evdist 1
r_net.add_evoked_drive('evdist10',
mu=r_net.external_drives['evdist1']['dynamics']['mu'] + 2 * t + isi,
sigma=r_net.external_drives['evdist1']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist1']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist1']['synaptic_delays'])
# copy evdist 2
r_net.add_evoked_drive('evdist11',
mu=r_net.external_drives['evdist2']['dynamics']['mu'] + 2 * t + isi,
sigma=r_net.external_drives['evdist2']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist2']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist2']['synaptic_delays'])
# copy evdist 3
r_net.add_evoked_drive('evdist12',
mu=r_net.external_drives['evdist3']['dynamics']['mu'] + 2 * t + isi,
sigma=r_net.external_drives['evdist3']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist3']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist3']['synaptic_delays'])
#---------------
# copy evprox 1
r_net.add_evoked_drive('evprox13',
mu=r_net.external_drives['evprox1']['dynamics']['mu'] + 2 * t + t + isi,
sigma=r_net.external_drives['evprox1']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox1']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox1']['synaptic_delays'])
# copy evprox 2
r_net.add_evoked_drive('evprox14',
mu=r_net.external_drives['evprox2']['dynamics']['mu'] + 2 * t + t + isi,
sigma=r_net.external_drives['evprox2']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox2']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox2']['synaptic_delays'])
# copy evprox 3
r_net.add_evoked_drive('evprox15',
mu=r_net.external_drives['evprox3']['dynamics']['mu'] + 2 * t + t + isi,
sigma=r_net.external_drives['evprox3']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox3']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox3']['synaptic_delays'])
# copy evdist 1
r_net.add_evoked_drive('evdist13',
mu=r_net.external_drives['evdist1']['dynamics']['mu'] + 2 * t + t + isi,
sigma=r_net.external_drives['evdist1']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist1']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist1']['synaptic_delays'])
# copy evdist 2
r_net.add_evoked_drive('evdist14',
mu=r_net.external_drives['evdist2']['dynamics']['mu'] + 2 * t + t + isi,
sigma=r_net.external_drives['evdist2']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist2']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist2']['synaptic_delays'])
# copy evdist 3
r_net.add_evoked_drive('evdist15',
mu=r_net.external_drives['evdist3']['dynamics']['mu'] + 2 * t + t + isi,
sigma=r_net.external_drives['evdist3']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist3']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist3']['synaptic_delays'])
#---------------
# copy evprox 1
r_net.add_evoked_drive('evprox16',
mu=r_net.external_drives['evprox1']['dynamics']['mu'] + 2 * t + 2*t + isi,
sigma=r_net.external_drives['evprox1']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox1']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox1']['synaptic_delays'])
# copy evprox 2
r_net.add_evoked_drive('evprox17',
mu=r_net.external_drives['evprox2']['dynamics']['mu'] + 2 * t + 2*t + isi,
sigma=r_net.external_drives['evprox2']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox2']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox2']['synaptic_delays'])
# copy evprox 3
r_net.add_evoked_drive('evprox18',
mu=r_net.external_drives['evprox3']['dynamics']['mu'] + 2 * t + 2*t + isi,
sigma=r_net.external_drives['evprox3']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox3']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox3']['synaptic_delays'])
# copy evdist 1
r_net.add_evoked_drive('evdist16',
mu=r_net.external_drives['evdist1']['dynamics']['mu'] + 2 * t + 2*t + isi,
sigma=r_net.external_drives['evdist1']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist1']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist1']['synaptic_delays'])
# copy evdist 2
r_net.add_evoked_drive('evdist17',
mu=r_net.external_drives['evdist2']['dynamics']['mu'] + 2 * t + 2*t + isi,
sigma=r_net.external_drives['evdist2']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist2']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist2']['synaptic_delays'])
# copy evdist 3
r_net.add_evoked_drive('evdist18',
mu=r_net.external_drives['evdist3']['dynamics']['mu'] + 2 * t + 2*t + isi,
sigma=r_net.external_drives['evdist3']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist3']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist3']['synaptic_delays'])
# ------- THIRD BURST -------
# copy evprox 1
r_net.add_evoked_drive('evprox19',
mu=r_net.external_drives['evprox1']['dynamics']['mu'] + 4*t + 2*isi,
sigma=r_net.external_drives['evprox1']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox1']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox1']['synaptic_delays'])
# copy evprox 2
r_net.add_evoked_drive('evprox20',
mu=r_net.external_drives['evprox2']['dynamics']['mu'] + 4*t + 2*isi,
sigma=r_net.external_drives['evprox2']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox2']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox2']['synaptic_delays'])
# copy evprox 3
r_net.add_evoked_drive('evprox21',
mu=r_net.external_drives['evprox3']['dynamics']['mu'] + 4*t + 2*isi,
sigma=r_net.external_drives['evprox3']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox3']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox3']['synaptic_delays'])
# copy evdist 1
r_net.add_evoked_drive('evdist19',
mu=r_net.external_drives['evdist1']['dynamics']['mu'] + 4*t + 2*isi,
sigma=r_net.external_drives['evdist1']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist1']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist1']['synaptic_delays'])
# copy evdist 2
r_net.add_evoked_drive('evdist20',
mu=r_net.external_drives['evdist2']['dynamics']['mu'] + 4*t + 2*isi,
sigma=r_net.external_drives['evdist2']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist2']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist2']['synaptic_delays'])
# copy evdist 3
r_net.add_evoked_drive('evdist21',
mu=r_net.external_drives['evdist3']['dynamics']['mu'] + 4*t + 2*isi,
sigma=r_net.external_drives['evdist3']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist3']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist3']['synaptic_delays'])
#---------------
# copy evprox 1
r_net.add_evoked_drive('evprox22',
mu=r_net.external_drives['evprox1']['dynamics']['mu'] + 4*t + 2*isi + t,
sigma=r_net.external_drives['evprox1']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox1']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox1']['synaptic_delays'])
# copy evprox 2
r_net.add_evoked_drive('evprox23',
mu=r_net.external_drives['evprox2']['dynamics']['mu'] + 4*t + 2*isi + t,
sigma=r_net.external_drives['evprox2']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox2']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox2']['synaptic_delays'])
# copy evprox 3
r_net.add_evoked_drive('evprox24',
mu=r_net.external_drives['evprox3']['dynamics']['mu'] + 4*t + 2*isi + t,
sigma=r_net.external_drives['evprox3']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox3']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox3']['synaptic_delays'])
# copy evdist 1
r_net.add_evoked_drive('evdist22',
mu=r_net.external_drives['evdist1']['dynamics']['mu'] + 4*t + 2*isi + t,
sigma=r_net.external_drives['evdist1']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist1']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist1']['synaptic_delays'])
# copy evdist 2
r_net.add_evoked_drive('evdist23',
mu=r_net.external_drives['evdist2']['dynamics']['mu'] + 4*t + 2*isi + t,
sigma=r_net.external_drives['evdist2']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist2']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist2']['synaptic_delays'])
# copy evdist 3
r_net.add_evoked_drive('evdist24',
mu=r_net.external_drives['evdist3']['dynamics']['mu'] + 4*t + 2*isi + t,
sigma=r_net.external_drives['evdist3']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist3']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist3']['synaptic_delays'])
#---------------
# copy evprox 1
r_net.add_evoked_drive('evprox25',
mu=r_net.external_drives['evprox1']['dynamics']['mu'] + 4*t + 2*isi + 2*t,
sigma=r_net.external_drives['evprox1']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox1']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox1']['synaptic_delays'])
# copy evprox 2
r_net.add_evoked_drive('evprox26',
mu=r_net.external_drives['evprox2']['dynamics']['mu'] + 4*t + 2*isi + 2*t,
sigma=r_net.external_drives['evprox2']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox2']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox2']['synaptic_delays'])
# copy evprox 3
r_net.add_evoked_drive('evprox27',
mu=r_net.external_drives['evprox3']['dynamics']['mu'] + 4*t + 2*isi + 2*t,
sigma=r_net.external_drives['evprox3']['dynamics']['sigma'],
numspikes=1,
location='proximal',
weights_ampa=r_net.external_drives['evprox3']['weights_ampa'],
weights_nmda=r_net.external_drives['evprox3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evprox3']['synaptic_delays'])
# copy evdist 1
r_net.add_evoked_drive('evdist25',
mu=r_net.external_drives['evdist1']['dynamics']['mu'] + 4*t + 2*isi + 2*t,
sigma=r_net.external_drives['evdist1']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist1']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist1']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist1']['synaptic_delays'])
# copy evdist 2
r_net.add_evoked_drive('evdist26',
mu=r_net.external_drives['evdist2']['dynamics']['mu'] + 4*t + 2*isi + 2*t,
sigma=r_net.external_drives['evdist2']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist2']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist2']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist2']['synaptic_delays'])
# copy evdist 3
r_net.add_evoked_drive('evdist27',
mu=r_net.external_drives['evdist3']['dynamics']['mu'] + 4*t + 2*isi + 2*t,
sigma=r_net.external_drives['evdist3']['dynamics']['sigma'],
numspikes=1,
location='distal',
weights_ampa=r_net.external_drives['evdist3']['weights_ampa'],
weights_nmda=r_net.external_drives['evdist3']['weights_nmda'],
synaptic_delays=r_net.external_drives['evdist3']['synaptic_delays'])
nets.append(r_net)
# ATTENTION: The below MPI run FAILS
# ATTENTION: Comment out this MPI run to use the subsequent serial run.
with MPIBackend(n_procs=5, mpi_cmd='mpiexec'):
dpl = simulate_dipole(nets[0], tstop=450, n_trials=1)
# ATTENTION: However, the below serial run SUCCEEDS
dpl = simulate_dipole(nets[0], tstop=450, n_trials=1)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment