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'])
