def simulate(nest_model, neuron_type): par=default_params.Inhibition(**{'perturbations':op.get()[0]}) kw={'model':nest_model, 'n':1, 'mm':{"withgid": True, 'to_file':False, 'to_memory':True, 'record_from':['V_m'] }, 'params':par.dic[neuron_type]} n=my_population.MyNetworkNode(**kw) model=kwargs.get('model', 'iaf_neuron') n=kwargs.get('n', 1) params=kwargs.get('params',{})
from core.network.manager import Builder_Go_NoGo_with_lesion_FS_ST_pulse as Builder from core.parallel_excecution import get_loop_index, loop from core import directories as dr from core import my_socket import scripts_inhibition.base_Go_NoGo_compete as module import sys import fig_01_and_02_pert as op import pprint pp = pprint.pprint from copy import deepcopy path_rate_runs = get_path_rate_runs('fig_01_and_02_sim_inh/') ops = [op.get()[0]] FILE_NAME = __file__.split('/')[-1][0:-3] FROM_DISK_0 = int(sys.argv[1]) if len(sys.argv) > 1 else 0 LOAD_MILNER_ON_SUPERMICRO = False NUM_NETS = 1 NUM_RUNS = 1 num_sims = NUM_NETS * NUM_RUNS dc = my_socket.determine_computer CORES = 40 * 4 if dc() == 'milner' else 10 JOB_ADMIN = config.Ja_milner if dc() == 'milner' else config.Ja_else LOCAL_NUM_THREADS = 40 if dc() == 'milner' else 10 WRAPPER_PROCESS = config.Wp_milner if dc() == 'milner' else config.Wp_else kwargs = { 'Builder': Builder,
from toolbox.network.manager import Builder_Go_NoGo_with_lesion_FS as Builder from toolbox.parallel_excecution import loop from toolbox import directories as dr from toolbox import my_socket import scripts_inhibition.Go_NoGo_compete as module import sys # from scripts inhibition import scripts_inhibition.Go_NoGo_compete as module import fig_01_and_02_pert as op import pprint pp=pprint.pprint path_rate_runs=get_path_rate_runs('fig_01_and_02_sim_inh/') ops=[op.get()[0]] FILE_NAME=__file__.split('/')[-1][0:-3] FROM_DISK_0=int(sys.argv[1]) if len(sys.argv)>1 else 0 proportion_connected=[0.1, 0.5, 1.] NUM_RUNS=len(proportion_connected) NUM_NETS=5 num_sims=NUM_NETS*NUM_RUNS dc=my_socket.determine_computer CORES=40 if dc()=='milner' else 10 JOB_ADMIN=config.Ja_milner if dc()=='milner' else config.Ja_else LOCAL_NUM_THREADS= 20 if dc()=='milner' else 10 WRAPPER_PROCESS=config.Wp_milner if dc()=='milner' else config.Wp_else kwargs={ 'Builder':Builder,
import scripts_inhibition.base_oscillation_beta as module import fig_01_and_02_pert as op import pprint import fig_defaults as fd pp = pprint.pprint path_rate_runs = get_path_rate_runs('fig_01_and_02_sim_inh/') FILE_NAME = __file__.split('/')[-1][0:-3] FROM_DISK_0 = int(sys.argv[1]) if len(sys.argv) > 1 else 0 LOAD_MILNER_ON_SUPERMICRO = False NUM_NETS = 2 amp_base = [fd.amp_beta] #numpy.arange(1.05, 1.2, 0.05) freqs = [0.2, 0.3, 0.4] #numpy.arange(0.5, .8, 0.2) ops = [op.get()[fd.idx_beta]] n = len(amp_base) m = len(freqs) amp_base = list(numpy.array([m * [v] for v in amp_base]).ravel()) freqs = list(freqs) * n STN_amp_mod = [fd.STN_amp_mod_beta] #range(1, 6, 2) num_runs = len(freqs) * len(STN_amp_mod) * len(ops) num_sims = NUM_NETS * num_runs dc = my_socket.determine_computer CORES = 40 if dc() == 'milner' else 10 JOB_ADMIN = config.Ja_milner if dc() == 'milner' else config.Ja_else LOCAL_NUM_THREADS = 20 if dc() == 'milner' else 10 WRAPPER_PROCESS = config.Wp_milner if dc() == 'milner' else config.Wp_else kwargs = {
pert_add_inhibition) from core import my_socket import config import fig_01_and_02_pert as op import scripts_inhibition.base_inhibition_striatum as module import pprint pp=pprint.pprint FILE_NAME=__file__.split('/')[-1][0:-3] FROM_DISK_0=0 NUM_NETS=1 ops=op.get() NUM_RUNS=len(ops) #A run for each perturbation num_sim=NUM_NETS*NUM_RUNS dc=my_socket.determine_computer CORES=40 if dc()=='milner' else 4 JOB_ADMIN=config.Ja_milner if dc()=='milner' else config.Ja_else LOCAL_NUM_THREADS= 20 if dc()=='milner' else 4 WRAPPER_PROCESS=config.Wp_milner if dc()=='milner' else config.Wp_else kwargs={ 'Builder':Builder, 'cores':CORES, 'file_name':FILE_NAME,
import sys import scripts_inhibition.base_oscillation_sw as module import fig_01_and_02_pert as op import pprint pp=pprint.pprint path_rate_runs=get_path_rate_runs('fig_01_and_02_sim_inh/') FILE_NAME=__file__.split('/')[-1][0:-3] FROM_DISK_0=int(sys.argv[1]) if len(sys.argv)>1 else 0 LOAD_MILNER_ON_SUPERMICRO=False NUM_NETS=2 amp_base=[fd.amp_sw] #numpy.arange(1.05, 1.2, 0.05) freqs=[ 0.4, 0.5, 0.6, 0.7] #numpy.arange(0.5, .8, 0.2) ops=[op.get()[fd.idx_sw]] n=len(amp_base) m=len(freqs) amp_base=list(numpy.array([m*[v] for v in amp_base]).ravel()) freqs=list(freqs)*n STN_amp_mod=[fd.STN_amp_mod_sw]#range(1, 6, 2) num_runs=len(freqs)*len(STN_amp_mod)*len(ops) num_sims=NUM_NETS*num_runs dc=my_socket.determine_computer CORES=40 if dc()=='milner' else 10 JOB_ADMIN=config.Ja_milner if dc()=='milner' else config.Ja_else LOCAL_NUM_THREADS= 20 if dc()=='milner' else 10 WRAPPER_PROCESS=config.Wp_milner if dc()=='milner' else config.Wp_else kwargs={
''' Created on Aug 26, 2015 @author: mikael ''' import fig_01_and_02_pert as op from core.network import default_params import pprint pp = pprint.pprint par = default_params.Inhibition(**{'perturbations': op.get()[0]}) pp(par.dic['nest'].keys()) #Look at MSN D1 (high or low gabaa reversal potential E_rev) print('MSN D1') pp(par.dic['nest']['M1_low']) #I am using low # pp(par.dic['nest']['M1_high']) print('MSN D2') pp(par.dic['nest']['M2_low']) #I am using low print('FSN') pp(par.dic['nest']['FS'])