Created
May 20, 2025 20:25
-
-
Save asoplata/c1635c3f6debbe00a4f1b548372f65d2 to your computer and use it in GitHub Desktop.
Test of HNN-Core network complexity limit (2025-05-20)
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| # 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