def get(): l = [] solution, s_mul, s_equal = get_solution_slow_GP_striatum_2() for y in numpy.arange(1.0, -0.1, -0.25): for z in numpy.arange(1.0, -0.1, -0.25): x = 2.5 d = {} for keys in s_mul: update(solution, d, keys) misc.dict_update(d, {'nest': {'GI_GA_gaba': {'weight': y}}}) misc.dict_update(d, {'nest': {'ST_GA_ampa': {'weight': z}}}) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) l[-1] += pl(d, '=', **{'name': 'GIGA_' + str(y) + '_STGA_' + str(z)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) for nodes in [['M1', 'M2'], ['M1', 'M2', 'FS'], ['M1', 'M2', 'FS','ST', 'GA', 'GI','SN']]: x=2.5 d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate']*=0.7 for node in nodes: misc.dict_update(d,{'node':{node:{'edge_wrap':False}}}) l[-1]+=pl(d, '=', **{'name':'mod_ST_beta_'+str(x)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution() # Decrease/increase C1,C2,CF,CS, EI and EA. Test idea that at # slow wave cortical and thalamic input is decreased. mod=numpy.arange(0.5, 1.5, 0.05) for x in mod: d={} for keys in s_mul: update(solution, d, keys) misc.dict_update(d, {'node':{'CF':{'rate':x}}} ) misc.dict_update(d, {'node':{'CS':{'rate':x}}} ) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) d['node']['C1']['rate']*=x d['node']['C2']['rate']*=x d['node']['EI']['rate']*=x d['node']['EA']['rate']*=x l[-1]+=pl(d, '=', **{'name':'mod_C1_C2_CF_CS_EI_EA_'+str(x)}) return l
def get(): l = [] solution, s_mul, s_equal = get_solution() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) v0 = numpy.arange(1.25, 3.75, 0.5) v1 = numpy.arange(2.5, 5., 0.5) v = [[x, y] for x in v0 for y in v1] for x, y in v: d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST': {'beta_I_GABAA_1': f_beta_rm(x)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(y)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(y)}}}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': ('mod_ST_beta_' + str(x) + '_' + str(y))}) return l
def get(): l = [] solution, s_mul, s_equal = get_solution_slow_GP_striatum() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) v = STN_ampa_gaba_input_magnitude() for _, x, y in v: # x=2.5-1 d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) # misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) # misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) # d['node']['EA']['rate']*=0.7 misc.dict_update(d, {'node': {'CS': {'rate': x}}}) misc.dict_update(d, {'nest': {'GI_ST_gaba': {'weight': y}}}) l[-1] += pl(d, '=', **{'name': 'mod_ST_inp_' + str(x) + '_' + str(y)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) x=2.5-1.25 d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate']*=0.7 misc.dict_update(d,{'node':{'CS':{'rate': 830.}}}) misc.dict_update(d,{'nest':{'GI_ST_gaba':{'weight':0.35 }}}) l[-1]+=pl(d, '=', **{'name':'mod_ST_beta_'+str(x)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution_slow_GP_striatum() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) x=2.5 d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'GI':{'beta_I_GABAA_1': f_beta_rm(2.)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':'mod_ST_beta_'+str(x)}) return l
def get(): l = [] solution, s_mul, s_equal = get_solution() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) x = 2.5 - 1.25 d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate'] *= 0.7 misc.dict_update(d, {'node': {'CS': {'rate': 830.}}}) misc.dict_update(d, {'nest': {'GI_ST_gaba': {'weight': 0.35}}}) l[-1] += pl(d, '=', **{'name': 'mod_ST_beta_' + str(x)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution() # d0=0.8 # f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) # x=2.5 for x in numpy.arange(0.7, 1.5 ,0.2): d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) # misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) # misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate']*=x l[-1]+=pl(d, '=', **{'name':'mod_EA_'+str(x)}) return l
def get(): l = [] solution, s_mul, s_equal = get_solution_slow_GP_striatum() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) for y in numpy.arange(0.9, 0.0, -0.2): x = 2.5 d = {} for keys in s_mul: update(solution, d, keys) misc.dict_update(d, {'conn': {'nest': {'GI_GA_gaba': {'weight': y}}}}) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': 'mod_ST_beta_' + str(x)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) v0=numpy.arange(1.25,3.75,0.5) v1=numpy.arange(2.5, 5.,0.5) v=[[x,y] for x in v0 for y in v1] for x, y in v: d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_GABAA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':('mod_ST_beta_' +str(x)+'_'+str(y)) }) return l
def get(): l = [] solution, s_mul, s_equal = get_solution() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) d = {} for x in numpy.arange(3.4, 6, 0.4): # x=3.4 y = 2.25 for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) print f_beta_rm(x) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_GABAA_1': f_beta_rm(y)}}}) # d['node']['EA']['rate']*=0.7 l[-1] += pl( d, '=', **{'name': 'mod_CS_STbeta_' + str(x) + '_GA_STbeta_' + str(y)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution_slow_GP_striatum() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) x=2.5 d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':'mod_ST_beta_'+str(x)}) return l
def get(): l = [] solution, s_mul, s_equal = get_solution() # Decrease/increase C1,C2,CF,CS, EI and EA. Test idea that at # slow wave cortical and thalamic input is decreased. mod = numpy.arange(0.5, 1.5, 0.05) for x in mod: d = {} for keys in s_mul: update(solution, d, keys) misc.dict_update(d, {'node': {'CF': {'rate': x}}}) misc.dict_update(d, {'node': {'CS': {'rate': x}}}) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) d['node']['C1']['rate'] *= x d['node']['C2']['rate'] *= x d['node']['EI']['rate'] *= x d['node']['EA']['rate'] *= x l[-1] += pl(d, '=', **{'name': 'mod_C1_C2_CF_CS_EI_EA_' + str(x)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution_slow_GP_striatum() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) v=STN_ampa_gaba_input_magnitude() for _, x, y in v: # x=2.5-1 d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) # misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) # misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) # d['node']['EA']['rate']*=0.7 misc.dict_update(d,{'node':{'CS':{'rate': x}}}) misc.dict_update(d,{'nest':{'GI_ST_gaba':{'weight':y }}}) l[-1]+=pl(d, '=', **{'name':'mod_ST_inp_'+str(x)+'_'+str(y)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution_slow_GP_striatum() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) d={} for x in numpy.arange(3.4, 6, 0.4): # x=3.4 y=2.25 for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_GABAA_1': f_beta_rm(y)}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':'mod_CS_STbeta_'+str(x)+'_GA_STbeta_'+str(y)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) y=2.5 v=numpy.arange(0.1,1.2,0.1) for x in v: d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}}) misc.dict_update(d,{'nest':{'GA_M1_gaba':{'weight': x}}}) misc.dict_update(d,{'nest':{'GA_M2_gaba':{'weight': x}}}) misc.dict_update(d,{'nest':{'GA_FS_gaba':{'weight': x}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':('mod_GA_striatum_1' + str(x)) }) return l
def get(): l = [] solution, s_mul, s_equal = get_solution() # d0=0.8 # f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) # x=2.5 for x in numpy.arange(0.7, 1.5, 0.2): d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) # misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) # misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate'] *= x l[-1] += pl(d, '=', **{'name': 'mod_EA_' + str(x)}) return l
def get(): l = [] solution, s_mul, s_equal = get_solution() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) y = 2.5 v = numpy.arange(0.1, 1.2, 0.1) for x in v: d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(y)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(y)}}}) misc.dict_update(d, {'nest': {'GA_M1_gaba': {'weight': x}}}) misc.dict_update(d, {'nest': {'GA_M2_gaba': {'weight': x}}}) misc.dict_update(d, {'nest': {'GA_FS_gaba': {'weight': x}}}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': ('mod_GA_striatum_1' + str(x))}) return l
def get(): l = [] solution, s_mul, s_equal = get_solution() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) # v=STN_ampa_gaba_input_magnitude()[0:5] for _, u, w in v: v0 = numpy.arange(1.25, 3.75, 1.) v1 = numpy.arange(1.5, 4, 1.) v = [[x, y] for x in v0 for y in v1] for x, y in v: d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update( d, {'nest': { 'ST': { 'beta_I_GABAA_1': f_beta_rm(x) } }}) misc.dict_update(d, {'nest': { 'ST': { 'beta_I_AMPA_1': f_beta_rm(y) } }}) misc.dict_update(d, {'nest': { 'ST': { 'beta_I_NMDA_1': f_beta_rm(y) } }}) misc.dict_update(d, {'node': {'CS': {'rate': u}}}) misc.dict_update(d, {'nest': {'GI_ST_gaba': {'weight': w}}}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl( d, '=', **{ 'name': ('mod_ST_beta_' + str(u) + '_' + str(w) + '_' + str(x) + '_' + str(y)) }) return l
def get(): l = [] solution, s_mul, s_equal = get_solution() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) # betas MS-MS v = [0.1, 0.25, 0.5, 0.75] for x, y in [[v1, v2] for v1 in v for v2 in v]: d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': { 'ST': { 'beta_I_AMPA_1': f_beta_rm(2.5) } }}) misc.dict_update(d, {'nest': { 'ST': { 'beta_I_NMDA_1': f_beta_rm(2.5) } }}) d['node']['EA']['rate'] *= 0.7 misc.dict_update(d, {'nest': {'MS': {'beta_I_GABAA_2': f_beta_rm(x)}}}) for conn in ['M1_M1_gaba', 'M1_M2_gaba', 'M2_M1_gaba', 'M2_M2_gaba']: misc.dict_update(d, {'conn': { conn: { 'beta_fan_in': f_beta_rm(y) } }}) l[-1] += pl(d, '=', **{'name': 'mod_MSg_' + str(x) + '_MSc_' + str(y)}) #beta CTX-STN return l
def get(): l=[] solution, s_mul, s_equal=get_solution_slow_GP_striatum() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) x=2.5 d={} for y in [12]+ list(numpy.arange(5,55,5)): # x=3.4 # y=2.25 for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) ratio=12./y dd={'nest':{'GA_M1_gaba':{'weight':0.8*ratio}, 'GA_M2_gaba':{'weight':0.8*ratio}}} misc.dict_update(d,dd) # Decreasing from 2 leads to ... # Increasing from 2 leads to ... dd={'nest':{'GA_FS_gaba':{'weight':2.*ratio}}} misc.dict_update(d,dd) # Just assumed to be 12 ms dd={'nest':{'M1_low':{'GABAA_3_Tau_decay':12./ratio}, 'M2_low':{'GABAA_3_Tau_decay':12./ratio}, 'FS_low':{'GABAA_2_Tau_decay':12./ratio}, }} misc.dict_update(d,dd) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':'mod_GAtau_'+str(y)}) return l
def get(): l = [] solution, s_mul, s_equal = get_solution() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) y = 2.5 w = numpy.arange(0.5, 1.2, 0.1) v = [0.1, 0.25, 0.5] ll = [[v1, v2] for v1 in v for v2 in w] pp(l) for x, z in ll: print x, z d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(y)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(y)}}}) misc.dict_update(d, {'nest': {'MS': {'beta_I_GABAA_2': f_beta_rm(x)}}}) for conn in ['M1_M1_gaba', 'M1_M2_gaba', 'M2_M1_gaba', 'M2_M2_gaba']: misc.dict_update(d, {'conn': { conn: { 'beta_fan_in': f_beta_rm(x) } }}) misc.dict_update(d, {'nest': {'GA_M1_gaba': {'weight': z}}}) misc.dict_update(d, {'nest': {'GA_M2_gaba': {'weight': z}}}) misc.dict_update(d, {'nest': {'GA_FS_gaba': {'weight': z}}}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': ('mod_GA_str_' + str(z) + '_MS_' + str(x))}) return l
def get(): l = [] solution, s_mul, s_equal = get_solution_slow_GP_striatum_2() d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) l[-1] += pl(d, '=', **{'name': ''}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution_slow_GP_striatum_2() d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) l[-1]+=pl(d, '=', **{'name':''}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) y=2.5 w=numpy.arange(0.5, 1.2, 0.1) v=[0.1,0.25, 0.5] ll=[[v1,v2] for v1 in v for v2 in w] pp(l) for x, z in ll: print x,z d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}}) misc.dict_update(d,{'nest':{'MS':{'beta_I_GABAA_2': f_beta_rm(x)}}}) for conn in ['M1_M1_gaba', 'M1_M2_gaba','M2_M1_gaba','M2_M2_gaba']: misc.dict_update(d,{'conn':{conn:{'beta_fan_in': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'GA_M1_gaba':{'weight': z}}}) misc.dict_update(d,{'nest':{'GA_M2_gaba':{'weight': z}}}) misc.dict_update(d,{'nest':{'GA_FS_gaba':{'weight': z}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':('mod_GA_str_' + str(z)+'_MS_'+str(x)) }) return l
def get(): l=[] solution, s_mul, s_equal=get_solution() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) # v=STN_ampa_gaba_input_magnitude()[0:5] for _, u, w in v: v0=numpy.arange(1.25, 3.75, 1.) v1=numpy.arange(1.5, 4, 1.) v=[[x,y] for x in v0 for y in v1] for x, y in v: d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_GABAA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}}) misc.dict_update(d,{'node':{'CS':{'rate': u}}}) misc.dict_update(d,{'nest':{'GI_ST_gaba':{'weight':w }}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':('mod_ST_beta_' +str(u)+'_'+str(w)+'_' +str(x)+'_'+str(y)) }) return l
def get(): l = [] solution, s_mul, s_equal = get_solution() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) y = 2.5 v = [1. / 10., 3. / 10, 5. / 10, 7. / 10, 9. / 10] for x in v: d = {} for keys in s_mul: update(solution, d, keys) misc.dict_update(d, {'conn': {'GA_M1_gaba': {'fan_in0': x}}}) misc.dict_update(d, {'conn': {'GA_M2_gaba': {'fan_in0': x}}}) misc.dict_update(d, {'nest': {'GA_M1_gaba': {'weight': 1. / x}}}) misc.dict_update(d, {'nest': {'GA_M2_gaba': {'weight': 1. / x}}}) l += [pl(d, '*', **{'name': ('mod_GA_MS_' + str(x))})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(y)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(y)}}}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': ''}) for x in v: d = {} for keys in s_mul: update(solution, d, keys) misc.dict_update(d, {'conn': {'GA_FS_gaba': {'fan_in0': x}}}) misc.dict_update(d, {'nest': {'GA_FS_gaba': {'weight': 1. / x}}}) l += [pl(d, '*', **{'name': ('mod_GA_FS_' + str(x))})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(y)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(y)}}}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': ''}) for x in v: d = {} for keys in s_mul: update(solution, d, keys) misc.dict_update(d, {'conn': {'GA_FS_gaba': {'fan_in0': x}}}) misc.dict_update(d, {'conn': {'GA_M1_gaba': {'fan_in0': x}}}) misc.dict_update(d, {'conn': {'GA_M2_gaba': {'fan_in0': x}}}) misc.dict_update(d, {'nest': {'GA_M1_gaba': {'weight': 1. / x}}}) misc.dict_update(d, {'nest': {'GA_M2_gaba': {'weight': 1. / x}}}) misc.dict_update(d, {'nest': {'GA_FS_gaba': {'weight': 1. / x}}}) l += [pl(d, '*', **{'name': ('mod_GA_str_' + str(x))})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(y)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(y)}}}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': ''}) return l
def get(): l = [] solution, s_mul, s_equal = get_solution() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) # betas MS-MS v = numpy.arange(1.5, 3, 0.25) for x in v: d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}}) l[-1] += pl(d, '=', **{'name': 'mod_ST_beta_' + str(x)}) v1 = numpy.arange(0.1, 1, 0.2) for x in v1: d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) d['node']['EA']['rate'] *= x l[-1] += pl(d, '=', **{'name': 'mod_EA_' + str(x)}) v2 = numpy.arange(1.1, 2, 0.2) for x in v2: d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) d['node']['EI']['rate'] *= x l[-1] += pl(d, '=', **{'name': 'mod_EI_' + str(x)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) d={} x=2.5 for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':'mod_ST_beta_'+str(x)}) for y in numpy.arange(1,11,1): d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate']*=0.7 misc.dict_update(d,{'conn':{'GI_FS_gaba':{'lesion': False}}}) misc.dict_update(d,{'conn':{'GI_FS_gaba':{'fan_in0': y}}}) misc.dict_update(d,{'conn':{'GI_M1_gaba':{'lesion': False}}}) misc.dict_update(d,{'conn':{'GI_M1_gaba':{'fan_in0': y}}}) misc.dict_update(d,{'conn':{'GI_M2_gaba':{'lesion': False}}}) misc.dict_update(d,{'conn':{'GI_M2_gaba':{'fan_in0': y}}}) l[-1]+=pl(d, '=', **{'name':'mod_GI_M2_'+str(y)}) for y, EA_rate in zip(numpy.arange(25,125,25), numpy.arange(0.9,1.3,0.1)): d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate']*=EA_rate misc.dict_update(d,{'conn':{'M2_GA_gaba':{'lesion': False}}}) misc.dict_update(d,{'conn':{'M2_GA_gaba':{'fan_in0': y}}}) l[-1]+=pl(d, '=', **{'name':'M2GA_'+str(y)+'_EArate_'+str(EA_rate)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) # betas MS-MS v=numpy.arange(1.5,3,0.25) for x in v: d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) l[-1]+=pl(d, '=', **{'name':'mod_ST_beta_'+str(x)}) v1=numpy.arange(0.1, 1, 0.2) for x in v1: d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) d['node']['EA']['rate']*=x l[-1]+=pl(d, '=', **{'name':'mod_EA_'+str(x)}) v2=numpy.arange(1.1, 2, 0.2) for x in v2: d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) d['node']['EI']['rate']*=x l[-1]+=pl(d, '=', **{'name':'mod_EI_'+str(x)}) return l
def get(): l = [] solution, s_mul, s_equal = get_solution_slow_GP_striatum() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) for y in numpy.arange(1.0, -0.1, -0.5): for z in numpy.arange(1.0, -0.1, -0.5): x = 2.5 d = {} for keys in s_mul: update(solution, d, keys) misc.dict_update(d, {'nest': {'GI_GA_gaba': {'weight': y}}}) misc.dict_update(d, {'nest': {'ST_GA_ampa': {'weight': z}}}) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': { 'ST': { 'beta_I_AMPA_1': f_beta_rm(x) } }}) misc.dict_update(d, {'nest': { 'ST': { 'beta_I_NMDA_1': f_beta_rm(x) } }}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': 'GIGA_' + str(y) + '_STGA_' + str(z)}) for y in numpy.arange(1., 6.): for z in numpy.arange(1., 6.): x = 2.5 d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST_GA_ampa': {'delay': y}}}) misc.dict_update(d, {'nest': {'ST_GI_ampa': {'delay': y}}}) misc.dict_update(d, {'nest': {'GI_ST_gaba': {'delay': z}}}) misc.dict_update(d, {'nest': { 'ST': { 'beta_I_AMPA_1': f_beta_rm(x) } }}) misc.dict_update(d, {'nest': { 'ST': { 'beta_I_NMDA_1': f_beta_rm(x) } }}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl( d, '=', **{'name': 'STGPdelay_' + str(y) + '_GPSTdelay_' + str(z)}) for y in numpy.arange(2.5, 21., 2.5): x = 2.5 d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'CS_ST_ampa': {'delay': y}}}) misc.dict_update(d, {'nest': {'CS_ST_nmda': {'delay': y}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': 'CSSTdelay_' + str(y)}) for y in numpy.arange(2.5, 21., 2.5): x = 2.5 d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'C1_M1_ampa': {'delay': y}}}) misc.dict_update(d, {'nest': {'C1_M1_nmda': {'delay': y}}}) misc.dict_update(d, {'nest': {'C2_M2_ampa': {'delay': y}}}) misc.dict_update(d, {'nest': {'C2_M2_nmda': {'delay': y}}}) misc.dict_update(d, {'nest': {'CF_FS_ampa': {'delay': y}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': 'CSTRdelay_' + str(y)}) for y in numpy.arange(1, 8., 1.): x = 2.5 d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'GA_M1_gaba': {'delay': y}}}) misc.dict_update(d, {'nest': {'GA_M2_gaba': {'delay': y}}}) misc.dict_update(d, {'nest': {'GA_FS_gaba': {'delay': y}}}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': 'GASTRdelay_' + str(y)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution_slow_GP_striatum() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) for y in numpy.arange(1.0,-0.1,-0.5): for z in numpy.arange(1.0,-0.1,-0.5): x=2.5 d={} for keys in s_mul: update(solution, d, keys) misc.dict_update(d,{'nest':{'GI_GA_gaba':{'weight':y}}}) misc.dict_update(d,{'nest':{'ST_GA_ampa':{'weight':z}}}) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':'GIGA_'+str(y)+'_STGA_'+str(z)}) for y in numpy.arange(1.,6.): for z in numpy.arange(1.,6.): x=2.5 d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST_GA_ampa':{'delay':y}}}) misc.dict_update(d,{'nest':{'ST_GI_ampa':{'delay':y}}}) misc.dict_update(d,{'nest':{'GI_ST_gaba':{'delay':z}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':'STGPdelay_'+str(y)+'_GPSTdelay_'+str(z)}) for y in numpy.arange(2.5, 21., 2.5): x=2.5 d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'CS_ST_ampa':{'delay':y}}}) misc.dict_update(d,{'nest':{'CS_ST_nmda':{'delay':y}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':'CSSTdelay_'+str(y)}) for y in numpy.arange(2.5, 21., 2.5): x=2.5 d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'C1_M1_ampa':{'delay':y}}}) misc.dict_update(d,{'nest':{'C1_M1_nmda':{'delay':y}}}) misc.dict_update(d,{'nest':{'C2_M2_ampa':{'delay':y}}}) misc.dict_update(d,{'nest':{'C2_M2_nmda':{'delay':y}}}) misc.dict_update(d,{'nest':{'CF_FS_ampa':{'delay':y}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':'CSTRdelay_'+str(y)}) for y in numpy.arange(1, 8., 1.): x=2.5 d={} for keys in s_mul: update(solution, d, keys) l+=[pl(d, '*', **{'name':''})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'GA_M1_gaba':{'delay':y}}}) misc.dict_update(d,{'nest':{'GA_M2_gaba':{'delay':y}}}) misc.dict_update(d,{'nest':{'GA_FS_gaba':{'delay':y}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':'GASTRdelay_'+str(y)}) return l
def get(): l = [] solution, s_mul, s_equal = get_solution() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) d = {} x = 2.5 for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': 'mod_ST_beta_' + str(x)}) for y in numpy.arange(1, 11, 1): d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate'] *= 0.7 misc.dict_update(d, {'conn': {'GI_FS_gaba': {'lesion': False}}}) misc.dict_update(d, {'conn': {'GI_FS_gaba': {'fan_in0': y}}}) misc.dict_update(d, {'conn': {'GI_M1_gaba': {'lesion': False}}}) misc.dict_update(d, {'conn': {'GI_M1_gaba': {'fan_in0': y}}}) misc.dict_update(d, {'conn': {'GI_M2_gaba': {'lesion': False}}}) misc.dict_update(d, {'conn': {'GI_M2_gaba': {'fan_in0': y}}}) l[-1] += pl(d, '=', **{'name': 'mod_GI_M2_' + str(y)}) for y, EA_rate in zip(numpy.arange(25, 125, 25), numpy.arange(0.9, 1.3, 0.1)): d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}}) misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}}) d['node']['EA']['rate'] *= EA_rate misc.dict_update(d, {'conn': {'M2_GA_gaba': {'lesion': False}}}) misc.dict_update(d, {'conn': {'M2_GA_gaba': {'fan_in0': y}}}) l[-1] += pl(d, '=', **{'name': 'M2GA_' + str(y) + '_EArate_' + str(EA_rate)}) return l
def get(): l = [] solution, s_mul, s_equal = get_solution() d0 = 0.8 f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0)) for z0 in numpy.arange(1, 6., 1): for z1 in numpy.arange(1, 6., 1): y = 2.5 x = 2.5 d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'ST_GA_ampa': {'delay': z0}}}) misc.dict_update(d, {'nest': {'ST_GI_ampa': {'delay': z0}}}) misc.dict_update(d, {'nest': {'GI_ST_gaba': {'delay': z1}}}) misc.dict_update(d, {'nest': {'C1_M1_ampa': {'delay': y}}}) misc.dict_update(d, {'nest': {'C1_M1_nmda': {'delay': y}}}) misc.dict_update(d, {'nest': {'C2_M2_ampa': {'delay': y}}}) misc.dict_update(d, {'nest': {'C2_M2_nmda': {'delay': y}}}) misc.dict_update(d, {'nest': {'CF_FS_ampa': {'delay': y}}}) misc.dict_update(d, {'nest': { 'ST': { 'beta_I_AMPA_1': f_beta_rm(x) } }}) misc.dict_update(d, {'nest': { 'ST': { 'beta_I_NMDA_1': f_beta_rm(x) } }}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': 'CSTRdelay_' + str(y)}) for y in numpy.arange(1, 8., 1.): z = 2.5 # ms x = 2.5 d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'C1_M1_ampa': {'delay': z}}}) misc.dict_update(d, {'nest': {'C1_M1_nmda': {'delay': z}}}) misc.dict_update(d, {'nest': {'C2_M2_ampa': {'delay': z}}}) misc.dict_update(d, {'nest': {'C2_M2_nmda': {'delay': z}}}) misc.dict_update(d, {'nest': {'CF_FS_ampa': {'delay': z}}}) misc.dict_update(d, {'nest': {'GA_M1_gaba': {'delay': y}}}) misc.dict_update(d, {'nest': {'GA_M2_gaba': {'delay': y}}}) misc.dict_update(d, {'nest': {'GA_FS_gaba': {'delay': y}}}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl(d, '=', **{'name': 'GASTRdelay_' + str(y)}) for y in numpy.arange(1, 11., 2): for z in numpy.arange(1, 11., 2): x = 2.5 d = {} for keys in s_mul: update(solution, d, keys) l += [pl(d, '*', **{'name': ''})] d = {} for keys in s_equal: update(solution, d, keys) misc.dict_update(d, {'nest': {'C1_M1_ampa': {'delay': y}}}) misc.dict_update(d, {'nest': {'C1_M1_nmda': {'delay': y}}}) misc.dict_update(d, {'nest': {'C2_M2_ampa': {'delay': y}}}) misc.dict_update(d, {'nest': {'C2_M2_nmda': {'delay': y}}}) misc.dict_update(d, {'nest': {'CF_FS_ampa': {'delay': y}}}) misc.dict_update(d, {'nest': {'CS_ST_ampa': {'delay': z}}}) misc.dict_update(d, {'nest': {'CS_ST_nmda': {'delay': z}}}) misc.dict_update(d, {'nest': { 'ST': { 'beta_I_AMPA_1': f_beta_rm(x) } }}) misc.dict_update(d, {'nest': { 'ST': { 'beta_I_NMDA_1': f_beta_rm(x) } }}) d['node']['EA']['rate'] *= 0.7 l[-1] += pl( d, '=', **{'name': 'CSTRdelay_' + str(y) + '_CSTdelay_' + str(z)}) return l
def get(): l=[] solution, s_mul, s_equal=get_solution() d0=0.8 f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0)) y=2.5 v=[1./10., 3./10, 5./10,7./10, 9./10] for x in v: d={} for keys in s_mul: update(solution, d, keys) misc.dict_update(d,{'conn':{'GA_M1_gaba':{'fan_in0': x}}}) misc.dict_update(d,{'conn':{'GA_M2_gaba':{'fan_in0': x}}}) misc.dict_update(d,{'nest':{'GA_M1_gaba':{'weight': 1./x}}}) misc.dict_update(d,{'nest':{'GA_M2_gaba':{'weight': 1./x}}}) l+=[pl(d, '*', **{'name':('mod_GA_MS_' + str(x))})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':'' }) for x in v: d={} for keys in s_mul: update(solution, d, keys) misc.dict_update(d,{'conn':{'GA_FS_gaba':{'fan_in0': x}}}) misc.dict_update(d,{'nest':{'GA_FS_gaba':{'weight': 1./x}}}) l+=[pl(d, '*', **{'name':('mod_GA_FS_' + str(x))})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':'' }) for x in v: d={} for keys in s_mul: update(solution, d, keys) misc.dict_update(d,{'conn':{'GA_FS_gaba':{'fan_in0': x}}}) misc.dict_update(d,{'conn':{'GA_M1_gaba':{'fan_in0': x}}}) misc.dict_update(d,{'conn':{'GA_M2_gaba':{'fan_in0': x}}}) misc.dict_update(d,{'nest':{'GA_M1_gaba':{'weight': 1./x}}}) misc.dict_update(d,{'nest':{'GA_M2_gaba':{'weight': 1./x}}}) misc.dict_update(d,{'nest':{'GA_FS_gaba':{'weight': 1./x}}}) l+=[pl(d, '*', **{'name':('mod_GA_str_' + str(x))})] d={} for keys in s_equal: update(solution, d, keys) misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}}) misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}}) d['node']['EA']['rate']*=0.7 l[-1]+=pl(d, '=', **{'name':'' }) return l