import utils
import matplotlib.pyplot as plt
t = 200.0
parameter.parameter_simulation['path_result'] = './new_parameter_test_stim'
# tools_simulation.print_bistability(parameter.parameter_model)
simulator = tools_simulation.init(
parameter.parameter_simulation,
parameter.parameter_model,
parameter.parameter_connection_between_region,
parameter.parameter_coupling,
parameter.parameter_integrator,
parameter.parameter_monitor,
parameter_stimulation={
'onset': 100.0,
"tau": 10.0,
"T": 100.0,
"weights": np.ones(84)
})
tools_simulation.run_simulation(simulator, t, parameter.parameter_simulation,
parameter.parameter_monitor)
tools_simulation.print_all(parameter.parameter_simulation['path_result'], 0.0,
t, 0, simulator.connectivity)
#
# tools_simulation.print_all(parameter_PosterWeek_desperation_disconnected.parameter_simulation['path_result'],
# 0.0,t,
# 0,5)
import tools_simulation
import parameter_76connectome_second_order_sleep
import numpy as np
parameter_76connectome_second_order_sleep.parameter_simulation[
'path_result'] = './76connectome_secondorder_wakey/'
t = 10000.0
parameter_76connectome_second_order_sleep.parameter_model['b_e'] = 0.0
parameter_76connectome_second_order_sleep.parameter_model['E_L_e'] = -50.0
simulator = tools_simulation.init(
parameter_76connectome_second_order_sleep.parameter_simulation,
parameter_76connectome_second_order_sleep.parameter_model,
parameter_76connectome_second_order_sleep.
parameter_connection_between_region,
parameter_76connectome_second_order_sleep.parameter_coupling,
parameter_76connectome_second_order_sleep.parameter_integrator,
parameter_76connectome_second_order_sleep.parameter_monitor)
tools_simulation.run_simulation(
simulator, t,
parameter_76connectome_second_order_sleep.parameter_simulation,
parameter_76connectome_second_order_sleep.parameter_monitor)
tools_simulation.print_all(
parameter_76connectome_second_order_sleep.
parameter_simulation['path_result'], 0.0, t, 0, 0)
tools_simulation.print_all(
parameter_76connectome_second_order_sleep.
parameter_simulation['path_result'], 0.0, t, 0, 5)
parameter_Matteo.parameter_simulation[
'path_result'] = '/mnt/usb-WD_easystore_25FC_575856314136373548355350-0:0-part1/4TB_copy/Athens_Matteo_TVB_params/MTOO_b0_try/'
#'/media/gpi/4TB/Athens_Matteo_TVB_params/MTOO_b100_try'
Nregions = 84
parameter_stimulation = {}
parameter_stimulation['onset'] = 9000 # onset time in ms
parameter_stimulation['tau'] = 10 # stimulus duration in ms
parameter_stimulation['T'] = 100000 # interstimulus interval in ms
parameter_stimulation['weights'] = np.zeros(Nregions)
parameter_stimulation['weights'][59] = 1e-3
simulator = tools_simulation.init(
parameter_Matteo.parameter_simulation,
parameter_Matteo.parameter_model,
parameter_Matteo.parameter_connection_between_region,
parameter_Matteo.parameter_coupling,
parameter_Matteo.parameter_integrator,
parameter_Matteo.parameter_monitor,
parameter_stimulation=parameter_stimulation)
# tools_simulation.run_simulation(simulator,
# t,
# parameter_Matteo.parameter_simulation,
# parameter_Matteo.parameter_monitor)
tools_simulation.print_all(
parameter_Matteo.parameter_simulation['path_result'], 4000.0, t, 0,
simulator.connectivity)
# tools_simulation.print_all(parameter_Matteo.parameter_simulation['path_result'],
# 1.0,t,
# # parameter_one_second_order.parameter_model['a_e']=0.0
# # parameter_one_second_order.parameter_model['T']=40.0
# parameter_one_second_order.parameter_integrator['noise_parameter']['nsig']=[1.0e-8,1.0e-8,0.0,0.0,0.0,0.0,0.0,]
# parameter_one_second_order.parameter_model['initial_condition']={
# "E": [0.01,0.01],
# "I": [0.01,0.01],
# "C_ee": [0.0,0.0],
# "C_ei": [0.0,0.0],
# "C_ii": [0.0,0.0],
# "W_e": [0.0,0.0],
# "W_i": [0.0,0.0]}
simulator = tools_simulation.init(
parameter_oneRegion_second_order_SlowWaves.parameter_simulation,
parameter_oneRegion_second_order_SlowWaves.parameter_model,
parameter_oneRegion_second_order_SlowWaves.
parameter_connection_between_region,
parameter_oneRegion_second_order_SlowWaves.parameter_coupling,
parameter_oneRegion_second_order_SlowWaves.parameter_integrator,
parameter_oneRegion_second_order_SlowWaves.parameter_monitor)
tools_simulation.run_simulation(
simulator, 30000.0,
parameter_oneRegion_second_order_SlowWaves.parameter_simulation,
parameter_oneRegion_second_order_SlowWaves.parameter_monitor)
tools_simulation.print_EI_one(
parameter_oneRegion_second_order_SlowWaves.
parameter_simulation['path_result'], 0.0, 30000.0, 0, 0)
#
tools_simulation.print_bistability(
parameter_oneRegion_second_order_SlowWaves.parameter_model)
import parameter_one_first_order
import tools_simulation
simulator = tools_simulation.init(
parameter_one_first_order.parameter_simulation,
parameter_one_first_order.parameter_model,
parameter_one_first_order.parameter_connection_between_region,
parameter_one_first_order.parameter_coupling,
parameter_one_first_order.parameter_integrator,
parameter_one_first_order.parameter_monitor)
tools_simulation.run_simulation(simulator, 20000.0,
parameter_one_first_order.parameter_simulation,
parameter_one_first_order.parameter_monitor)
import tools_simulation
import parameter_Matteo as parameter_PosterWeek_desperation_disconnected
import numpy as np
import utils
import matplotlib.pyplot as plt
t=10.0
parameter_PosterWeek_desperation_disconnected.parameter_simulation['path_result']='./posterweek_apoint2_b0_disconnected'
simulator = tools_simulation.init(parameter_PosterWeek_desperation_disconnected.parameter_simulation,
parameter_PosterWeek_desperation_disconnected.parameter_model,
parameter_PosterWeek_desperation_disconnected.parameter_connection_between_region,
parameter_PosterWeek_desperation_disconnected.parameter_coupling,
parameter_PosterWeek_desperation_disconnected.parameter_integrator,
parameter_PosterWeek_desperation_disconnected.parameter_monitor)
tools_simulation.run_simulation(simulator,
t,
parameter_PosterWeek_desperation_disconnected.parameter_simulation,
parameter_PosterWeek_desperation_disconnected.parameter_monitor)
tools_simulation.print_all(parameter_PosterWeek_desperation_disconnected.parameter_simulation['path_result'],
0.0,t,
0,0)
#
# tools_simulation.print_all(parameter_PosterWeek_desperation_disconnected.parameter_simulation['path_result'],
# 0.0,t,
# 0,5)
# tools_simulation.print_EI_one(parameter_PosterWeek_desperation_disconnected.parameter_simulation['path_result'],
# 0.0,t,