Ejemplo n.º 1
0
    get_path_rate_runs, get_path_logs, get_args_list_oscillation,
    get_kwargs_list_indv_nets, par_process_and_thread,
    pert_set_data_path_to_milner_on_supermicro, pert_add_oscillations)

from core.network import default_params
from core.network.manager import Builder_beta as Builder
from core.parallel_excecution import loop

import numpy
import sys
import scripts_inhibition.base_oscillation_beta as module
import oscillation_perturbations8 as op
import pprint
pp = pprint.pprint

path_rate_runs = get_path_rate_runs('simulate_inhibition_ZZZ8/')
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 = numpy.arange(0.8, 1.6, 0.1)
freqs = numpy.arange(0.5, 3., 0.5)
n = len(amp_base)
m = len(freqs)
amp_base = list(numpy.array([m * [v] for v in amp_base]).ravel())
freqs = list(freqs) * n
num_runs = len(freqs)
num_sims = NUM_NETS * num_runs
                      pert_add_oscillations) 

from core.network.manager import Builder_beta as Builder
from core.parallel_excecution import loop
from core import directories as dr
from core import my_socket

import eNeuro_fig_defaults as fd
import eNeuro_fig_01_and_02_pert as op
import eNeuro_fig_07_pert_nuclei as op_neuclei
import sys
import scripts_inhibition.base_oscillation_beta as module
import pprint
pp=pprint.pprint

path_rate_runs=get_path_rate_runs('eNeuro_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

ops=[op.get()[fd.idx_beta]]
amp_base=[fd.amp_beta] #numpy.arange(1.05, 1.2, 0.05)
freqs=[fd.freq_beta] #numpy.arange(0.5, .8, 0.2)
STN_amp_mod=[fd.STN_amp_mod_beta]

NUM_NETS=2
NUM_RUNS=len(op_neuclei.get([0,7])) #A run for each perturbation
num_sim=NUM_NETS*NUM_RUNS

dc=my_socket.determine_computer
CORES=40*4 if dc()=='milner' else 10
Ejemplo n.º 3
0
from core.network.manager import Builder_beta as Builder
from core.parallel_excecution import loop
from core import directories as dr
from core import my_socket

import eNeuro_fig_defaults as fd
import eNeuro_fig_01_and_02_pert as op
import eNeuro_fig_07_pert_conn as op_add
import sys
import scripts_inhibition.base_oscillation_beta as module
import pprint

pp = pprint.pprint

path_rate_runs = get_path_rate_runs('eNeuro_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 1
LOAD_MILNER_ON_SUPERMICRO = False

ops = [op.get()[fd.idx_beta]]
amp_base = [fd.amp_beta]  #numpy.arange(1.05, 1.2, 0.05)
freqs = [fd.freq_beta]  #numpy.arange(0.5, .8, 0.2)
STN_amp_mod = [fd.STN_amp_mod_beta]

NUM_RUNS = len(op_add.get([0, 6]))
NUM_NETS = 2
num_sims = NUM_NETS * NUM_RUNS

dc = my_socket.determine_computer
CORES = 40 * 4 if dc() == 'milner' else 10
Ejemplo n.º 4
0
                                              get_kwargs_list_indv_nets,
                                              pert_add_oscillations)

from core.network.manager import Builder_beta as Builder
from core.parallel_excecution import loop
from core import directories as dr
from core import my_socket

import numpy
import sys
import scripts_inhibition.base_oscillation_beta as module
from python.scripts_mohammad.old 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 = [1.1]  #numpy.arange(1.05, 1.2, 0.05)
freqs = [0.1, 0.2, 0.3, 0.4]  #numpy.arange(0.5, .8, 0.2)
ops = [op.get()[0]]
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 = [3.]  #range(1, 6, 2)
num_runs = len(freqs) * len(STN_amp_mod) * len(ops)
Ejemplo n.º 5
0
    get_path_rate_runs, get_path_logs, get_args_list_oscillation,
    get_kwargs_list_indv_nets, par_process_and_thread,
    pert_set_data_path_to_milner_on_supermicro, pert_add_oscillations)

from core.network import default_params
from core.network.manager import Builder_slow_wave2 as Builder
from core.parallel_excecution import loop

import numpy
import scripts_inhibition.base_oscillation_sw as module
import sys
import oscillation_perturbations34_slow_sw as op
import pprint
pp = pprint.pprint

path_rate_runs = get_path_rate_runs('simulate_inhibition_ZZZ34_slow_sw/')
FILE_NAME = __file__.split('/')[-1][0:-3]
FROM_DISK_0 = int(sys.argv[1]) if len(sys.argv) > 1 else 2
LOAD_MILNER_ON_SUPERMICRO = False

NUM_NETS = 2
amp_base = numpy.arange(.85, 1.0, 0.05)
freqs = numpy.arange(0.5, .8, 0.2)
# amp_base=[1.0] #numpy.arange(0.7, 1.1, 0.05)
# freqs=numpy.arange(0.1, 1.0, 0.1)

n = len(amp_base)
m = len(freqs)
amp_base = list(numpy.array([m * [v] for v in amp_base]).ravel())
freqs = list(freqs) * n
Ejemplo n.º 6
0
    'file_name': FILE_NAME,
    'freqs': [0.75],
    'freq_oscillation': 1.,
    'from_disk_0': FROM_DISK_0,
    'i0': FROM_DISK_0,
    'job_name': 'sw_dop_pert',
    'l_hours': ['00', '00', '00'],
    'l_minutes': ['15', '10', '5'],
    'l_seconds': ['00', '00', '00'],
    'local_threads_milner': 10,
    'local_threads_superm': 4,
    'module': module,
    'nets': ['Net_' + str(i) for i in range(2)],  #The nets for each run
    'op_pert_add': op_dop,
    'path_code': default_params.HOME_CODE,
    'path_rate_runs': get_path_rate_runs('simulate_inhibition_ZZZ4/'),
    'path_results': get_path_logs(LOAD_MILNER_ON_SUPERMICRO, FILE_NAME),
    'perturbation_list': [op.get()[7]],
    'sim_time': 40000.0,
    'size': 20000.0,
}

d_process_and_thread = par_process_and_thread(**kwargs)
pp(d_process_and_thread)
kwargs.update(d_process_and_thread)

pp(kwargs)

p_list = pert_add_oscillations(**kwargs)
p_list = pert_set_data_path_to_milner_on_supermicro(p_list,
                                                    LOAD_MILNER_ON_SUPERMICRO)
Ejemplo n.º 7
0
# from toolbox.network import default_params
from toolbox.network.manager import Builder_beta as Builder
from toolbox.parallel_excecution import loop
from toolbox import directories as dr
from toolbox import my_socket

import config
import numpy
import sys
import scripts_inhibition.base_oscillation_beta as module
import oscillation_perturbations0 as op
import pprint
pp=pprint.pprint

path_rate_runs=get_path_rate_runs('simulate_inhibition0/')
print path_rate_runs
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=[1., 1.05, 1.1] #numpy.arange(1.05, 1.2, 0.05)
amp_base=[1.1] #numpy.arange(1.05, 1.2, 0.05)

freqs=[ 0.6] #numpy.arange(0.5, .8, 0.2)
n=len(amp_base)
m=len(freqs)
amp_base=list(numpy.array([m*[v] for v in amp_base]).ravel()) 
freqs=list(freqs)*n
Ejemplo n.º 8
0
# from core.network import default_params
from core.network.manager import Builder_slow_wave2 as Builder
from core.parallel_excecution import loop
from core import directories as dr
from core import my_socket

from scripts_inhibition import config
import numpy
import sys
import scripts_inhibition.base_oscillation_sw as module
import oscillation_perturbations_final as op
import pprint
pp = pprint.pprint

path_rate_runs = get_path_rate_runs('simulate_inhibition_final/')
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 = [0.9]  #numpy.arange(1.05, 1.2, 0.05)
freqs = [0.2, 0.3, 0.4, 0.5]  #numpy.arange(0.5, .8, 0.2)
ops = [op.get()[1]]
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 = [1.]  #range(1, 6, 2)
num_runs = len(freqs) * len(STN_amp_mod) * len(ops)
# from core.network import default_params
from core.network.manager import Builder_slow_wave2 as Builder
from core.parallel_excecution import loop
from core import directories as dr
from core import my_socket

from scripts_inhibition import config
import numpy
import sys
import scripts_inhibition.base_oscillation_sw as module
import oscillation_perturbations_new_beginning_slow4 as op
import pprint
pp=pprint.pprint

path_rate_runs=get_path_rate_runs('simulate_inhibition_new_beginning_slow4/')
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=[0.7, 0.75, 0.8, 0.85, 0.9] #numpy.arange(1.05, 1.2, 0.05)
freqs=[ 0.5] #numpy.arange(0.5, .8, 0.2)
ops=op.get()[:4*7*8]
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=[1., 3.]#range(1, 6, 2)
num_runs=len(freqs)*len(STN_amp_mod)*len(ops)
Ejemplo n.º 10
0
    get_path_rate_runs, get_path_logs, get_args_list_oscillation,
    get_kwargs_list_indv_nets, par_process_and_thread,
    pert_set_data_path_to_milner_on_supermicro, pert_add_oscillations)

from core.network import default_params
from core.network.manager import Builder_beta as Builder
from core.parallel_excecution import loop

import numpy
import sys
import scripts_inhibition.base_oscillation_beta as module
import oscillation_perturbations13_slow_2 as op
import pprint
pp = pprint.pprint

path_rate_runs = get_path_rate_runs('simulate_inhibition_ZZZ13_slow_2/')
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 = numpy.arange(0.9, 1.25, 0.05)
freqs = numpy.arange(0.3, 1.5, 0.1)
n = len(amp_base)
m = len(freqs)
amp_base = list(numpy.array([m * [v] for v in amp_base]).ravel())
freqs = list(freqs) * n
num_runs = len(freqs)
num_sims = NUM_NETS * num_runs
import scripts_inhibition.base_Go_NoGo_compete as module

# from scripts inhibition import scripts_inhibition.base_Go_NoGo_compete as module
import oscillation_perturbations41_slow as op
import oscillation_perturbations_conns as op_conn
import oscillation_perturbations_nuclei as op_nuc
import oscillation_perturbations_dop as op_dop

import sys
import pprint
pp = pprint.pprint

# for p in [op_dop.get()[5]]:
#     print p

path_rate_runs = get_path_rate_runs('simulate_inhibition_ZZZ41_slow/')
oscillation_returbations_index = 5
# op_pert_add=[op_dop.get()[5]]
# l_op_conn=[36, 66, 97, 114, 121, 127, 138]
# l_op_nuc=[33, 49, 57]
# l_op_dop=[5,6]
l_op_conn = []  #12, 97, 108, 109, 127, 132 ]
l_op_nuc = []  #16, 33, 49, 57, 64]
l_op_dop = [-1]  #,6,6]

opc = op_conn.get()
op_pert_add = [opc[i] for i in l_op_conn]

opn = op_nuc.get()
op_pert_add += [opn[i] for i in l_op_nuc]
from core.network.manager import Builder_slow_wave2 as Builder
# from core.network.manager import Builder_beta_GA_GI_ST as Builder
from core.parallel_excecution import loop
from core import directories as dr
from core import my_socket

from scripts_inhibition import config
import numpy
import sys
import optimization_GA_GI_ST as module
import oscillation_perturbations_new_beginning_slow_fitting_GI_GA_ST_sw_0 as op
import pprint
pp = pprint.pprint

path_rate_runs = get_path_rate_runs(
    'simulate_inhibition_new_beginning_slow_fitting_GI_GA_ST_sw_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 = 2

amp_base = [0.9]  #numpy.arange(1.05, 1.2, 0.05)
freqs = [0.3]  #numpy.arange(0.5, .8, 0.2)
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 = [1.]  #range(1, 6, 2)
ops = [op.get()[0]]
num_runs = len(freqs) * len(STN_amp_mod) * len(ops)
                      par_process_and_thread,
                      pert_set_data_path_to_milner_on_supermicro, 
                      pert_add_oscillations) 

from core.network import default_params
from core.network.manager import Builder_beta as Builder
from core.parallel_excecution import loop

import numpy
import sys
import scripts_inhibition.base_oscillation_beta as module
import oscillation_perturbations_control_slow_beta as op
import pprint
pp=pprint.pprint

path_rate_runs=get_path_rate_runs('simulate_inhibition_control_slow_beta/')
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=[1.05] #numpy.arange(1.05, 1.2, 0.05)
freqs=[0.7] #numpy.arange(0.5, .8, 0.2)
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=[3.]#range(1, 6, 2)
num_runs=len(freqs)*len(STN_amp_mod)*len(op.get())
num_sims=NUM_NETS*num_runs
Ejemplo n.º 14
0
                      get_kwargs_list_indv_nets)

from core.network.manager import Builder_Go_NoGo_with_lesion_FS as Builder
from core.parallel_excecution import 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=5

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 core.network.manager import Builder_beta as Builder
# from core.network.manager import Builder_beta_GA_GI_ST as Builder
from core.parallel_excecution import loop
from core import directories as dr
from core import my_socket

from scripts_inhibition import config
import numpy
import sys
import optimization_GA_GI_ST as module
from scripts_inhibition import oscillation_perturbations_new_beginning_slow_fitting_GI_GA_ST_beta_0 as op
import pprint
pp=pprint.pprint

path_rate_runs=get_path_rate_runs('simulate_inhibition_new_beginning_slow_fitting_GI_GA_ST_beta_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=2

amp_base=[1.1] #numpy.arange(1.05, 1.2, 0.05)
freqs=[ 0.3] #numpy.arange(0.5, .8, 0.2)
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=[3.]#range(1, 6, 2)
ops=op.get()
num_runs=len(freqs)*len(STN_amp_mod)*len(ops)
Ejemplo n.º 16
0
    get_kwargs_list_indv_nets, par_process_and_thread,
    pert_set_data_path_to_milner_on_supermicro, pert_add_oscillations)

from core.network import default_params
from core.network.manager import Builder_beta as Builder
from core.parallel_excecution import loop

import numpy
import sys
import scripts_inhibition.base_oscillation_beta as module
import oscillation_perturbations13 as op
import pprint

pp = pprint.pprint

path_rate_runs = get_path_rate_runs('simulate_inhibition_ZZZ13_dopMSGP/')
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 = [1.]  #numpy.arange(0.9, 1.25, 0.05)
freqs = numpy.arange(.8, 3.5, 0.1)
n = len(amp_base)
m = len(freqs)
amp_base = list(numpy.array([m * [v] for v in amp_base]).ravel())
freqs = list(freqs) * n
num_runs = len(freqs)
num_sims = NUM_NETS * num_runs