def run_simulation_CL1(tauE_, wE_, tauI_, wI_, Group, Synapses, Target, time, dt_, delta, rate_correction): restore('initialised') # set the parameters Group.set_states({'tauE': tauE_ * ms, 'tauI': tauI_ * ms}) print(f'taueE: {tauE_} - tauI {tauI_}') Synapses.set_states({'wE': wE_ * nS, 'wI': wI_ * nS}) print(f'wE: {wE_} - wI {wI_}') _, model_spike_monitor = nc.add_monitors(Group) target_spike_monitor = SpikeMonitor(Target, name='CL1_target_spike_monitor') run(time) gf = metric.compute_gamma_factor(model_spike_monitor, target_spike_monitor, time, dt_=dt_, delta=delta, rate_correction=rate_correction) print(f'Gamma factor: {gf}') return gf
# params_TL2['tauI'] = 1 * ms # synapses_TL2['wI'] = 300 * nS print(params_TL2) print(synapses_TL2) # Neuron group G_TL2 = nc.generate_neuron_groups(N_TL2, eqs, threshold_eqs, reset_eqs, params_TL2, name='TL2_test') # Add monitors STM_TL2, SPM_TL2 = nc.add_monitors(G_TL2, name='TL2_test') # Connect heading to TL2 S_P_HEADING_TL2 = nc.connect_synapses(P_HEADING, G_TL2, W_HEADING_TL2, params=synapses_TL2, model=synapses_model, on_pre=synapses_eqs_ex) # Run simulation run(T_outbound * time_step * ms) cx_spiking.plotting.plot_rate_cx_log_spikes( cx_log.tl2, TL2_spike_rates,
# params_TN2['tauI'] = 1 * ms # synapses_TN2['wI'] = 300 * nS print(params_TN2) print(synapses_TN2) # Neuron group G_TN2 = nc.generate_neuron_groups(N_TN2, eqs, threshold_eqs, reset_eqs, params_TN2, name='TN2_test') # Add monitors STM_TN2, SPM_TN2 = nc.add_monitors(G_TN2, name='TN2_test') # Connect heading to TN2 S_P_FLOW_TN2 = nc.connect_synapses(P_FLOW, G_TN2, W_FLOW_TN2, params=synapses_TN2, model=synapses_model, on_pre=synapses_eqs_ex) # Run simulation run(T_outbound * time_step * ms) cx_spiking.plotting.plot_rate_cx_log_spikes( cx_log.tn2, TN2_spike_rates_max,
# Neuron group G_TL2 = nc.generate_neuron_groups(N_TL2, eqs, threshold_eqs, reset_eqs, TL2_neuron_params, name='TL2_test') G_CL1 = nc.generate_neuron_groups(N_CL1, eqs, threshold_eqs, reset_eqs, neuron_params, name='CL1_test') # Add monitors STM_TL2, SPM_TL2 = nc.add_monitors(G_TL2, name='TL2_test') STM_CL1, SPM_CL1 = nc.add_monitors(G_CL1, name='CL1_test') # Connect heading to TL2 S_P_HEADING_TL2 = nc.connect_synapses(P_HEADING, G_TL2, W_HEADING_TL2, model=synapses_model, params=TL2_synapses_params, on_pre=synapses_eqs_ex) S_TL2_CL1 = nc.connect_synapses(G_TL2, G_CL1, W_TL2_CL1, model=synapses_model, params=synapses_params, on_pre=synapses_eqs_ex)