{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "0b85a598",
"metadata": {},
"outputs": [],
"source": [
"from brian2 import *\n",
"\n",
"from adex_sine import *\n",
"\n",
"defaultclock.dt = 10 * us\n",
"\n",
"\n",
"class Model:\n",
" # Model type flags\n",
" SINE = 0\n",
" EXP2SYN = 1\n",
"\n",
" # Noise flags\n",
" HIGH = 0\n",
" LOW = 1\n",
" OFF = -1\n",
"\n",
" # Noise parameters\n",
" EXCITATORY_NOISE_VARIANCE = {HIGH: 0.5 * nS, LOW: 0.25 * nS, OFF: 0 * nS}\n",
" INHIBITORY_NOISE_VARIANCE = {HIGH: 1.25 * nS, LOW: 0.625 * nS, OFF: 0 * nS}\n",
"\n",
" # Noise mean conductance\n",
" EXCITATORY_CONDUCTANCE = 1 * nS\n",
" INHIBITORY_CONDUCTANCE = 4 * nS\n",
"\n",
" DEFAULT_PARAMETERS = {\n",
" \"sigma_flux\" : 6.75*pA, \n",
" \"c\": 85 * pF,\n",
" \"tau_w\": 18 * ms,\n",
" \"b\": 0.25 * nA,\n",
" \"a\": 1.3 * nS,\n",
" \"v_T\": -45 * mV,\n",
" \"v_thresh\": 0 * mV,\n",
" \"DeltaT\": 0.2 * mV,\n",
" # EQUILIBRIUM POTENTIAL\n",
" \"e_l\": -65 * mV,\n",
" \"e_ex\": 0 * mV,\n",
" \"e_in\": -70 * mV,\n",
" # CONDUCTANCES\n",
" \"g_l\": 3 * nS,\n",
" \"mu_ex\": 0 * nS,\n",
" \"mu_in\": 0 * nS,\n",
" # EXCITATORY NOISE\n",
" \"sigma_ex\": 0 * nS,\n",
" \"tau_noise_ex\": 3 * ms,\n",
" # INHIBITORY NOISE\n",
" \"sigma_in\": 0 * nS,\n",
" \"tau_noise_in\": 10 * ms,\n",
" # SINE INPUT\n",
" \"f\": 100 * Hz,\n",
" \"A\": 0 * pA,\n",
" \"i_injected\": 0 * pA,\n",
" \"v_reset\": -70 * mV,\n",
" # m current\n",
" \"g_adapt\": 10 * nS,\n",
" \"e_k\": -90*mV,\n",
" \"beta_z\": -35*mV,\n",
" \"gamma_z\": 4*mV, #5\n",
" \"tau_z\": 100*ms,\n",
" }\n",
"\n",
" def __init__(\n",
" self, n, *, stim=None, noise=None, resistance=None, additional_vars=()\n",
" ):\n",
" if resistance is None:\n",
" raise ValueError(\"Resistance must be specified\")\n",
"\n",
" if noise is None:\n",
" raise ValueError(\"Noise must be specified\")\n",
"\n",
" self.stim_type = stim\n",
" self._input_resistance = None\n",
" self._noise_level = None\n",
" self._duration = 0\n",
" self.recorded_vars = (\"v\",) + additional_vars\n",
"\n",
" self.neurons = self.set_default(n_neuron=n)\n",
" self.set_resistance(resistance)\n",
" self.set_noise(noise)\n",
"\n",
" self.spikes = None\n",
" self.spiker = None\n",
" self.synapses = None\n",
" self.inhib_synapses = None\n",
" self.smon = None\n",
" self.network = None\n",
" self.build_network()\n",
"\n",
" def create_model(self):\n",
" return ADEX_MODEL, self.DEFAULT_PARAMETERS\n",
"\n",
" def set_default(self, n_neuron):\n",
" model, parameters = self.create_model()\n",
"\n",
" neurons = NeuronGroup(\n",
" n_neuron,\n",
" model=model,\n",
" method=\"Euler\",\n",
" name=\"neurons\",\n",
" threshold=\"v > v_thresh\",\n",
" reset=\"v = v_reset; w += b\",\n",
" )\n",
"\n",
" for parameter, value in parameters.items():\n",
" neurons.__setattr__(parameter, value)\n",
"\n",
" neurons.v = neurons.e_l # remove most of transient\n",
"\n",
" return neurons\n",
"\n",
" def set_resistance(self, level):\n",
" if level == self.LOW:\n",
" exc_conductance = self.EXCITATORY_CONDUCTANCE\n",
" inhib_conductance = self.INHIBITORY_CONDUCTANCE\n",
"\n",
" else:\n",
" exc_conductance = inhib_conductance = 0\n",
"\n",
" self._input_resistance = level\n",
" self._set_variable(\"mu_ex\", exc_conductance)\n",
" self._set_variable(\"mu_in\", inhib_conductance)\n",
"\n",
" def set_noise(self, level):\n",
" if level == self.HIGH or level == self.LOW:\n",
" exc_noise = self.EXCITATORY_NOISE_VARIANCE[level]\n",
" inhib_noise = self.INHIBITORY_NOISE_VARIANCE[level]\n",
"\n",
" else:\n",
" exc_noise = inhib_noise = 0\n",
"\n",
" self._noise_level = level\n",
" self._set_variable(\"sigma_ex\", exc_noise)\n",
" self._set_variable(\"sigma_in\", inhib_noise)\n",
"\n",
" def set_injected_current(self, amplitude):\n",
" self._set_variable(\"i_injected\", amplitude)\n",
" self._set_variable(\"A\", 0 * pA)\n",
"\n",
" def set_stimulus_current(self, amplitude):\n",
" self._set_variable(\"A\", amplitude)\n",
" self._set_variable(\"i_injected\", 0 * pA)\n",
"\n",
" @property\n",
" def f(self):\n",
" return self.neurons.f\n",
"\n",
" @f.setter\n",
" def f(self, new_f):\n",
" self._set_variable(\"f\", new_f) # this will reset smon\n",
" if self.stim_type == self.EXP2SYN:\n",
" self.spiker.T = 1 / new_f\n",
"\n",
" def run(self, duration, report=\"stdout\"):\n",
" self._duration = duration\n",
" self.network.run(duration, report=report)\n",
"\n",
" def build_network(self):\n",
" self.smon = StateMonitor(\n",
" self.neurons, self.recorded_vars, record=True, name=\"smon\"\n",
" )\n",
" self.spikes = SpikeMonitor(self.neurons, name=\"spikes\")\n",
"\n",
" self.network = Network(self.neurons, self.smon, self.spikes)\n",
"\n",
" def _set_variable(self, name, value):\n",
" self.neurons.__setattr__(name, value)\n",
" self.reset_recording()\n",
"\n",
" def reset_recording(self):\n",
" try:\n",
" self.network\n",
" except AttributeError:\n",
" return # network not yet initialized\n",
"\n",
" self.network.remove(self.smon, self.spikes)\n",
"\n",
" self.smon = StateMonitor(\n",
" self.neurons, self.recorded_vars, record=True, name=\"smon\"\n",
" )\n",
" self.spikes = SpikeMonitor(self.neurons, name=\"spikes\")\n",
"\n",
" self.network.add(self.smon, self.spikes)\n",
"\n",
" @property\n",
" def spike_train(self):\n",
" return self.spikes.spike_trains()\n",
"\n",
" @property\n",
" def firing_rate(self):\n",
" return self.spikes.count / self.duration\n",
"\n",
" @property\n",
" def duration(self):\n",
" return self._duration\n",
"\n",
" @property\n",
" def input_resistance(self):\n",
" if self._input_resistance == self.HIGH:\n",
" return \"HIGH\"\n",
" else:\n",
" return \"LOW\"\n",
"\n",
" @property\n",
" def noise_level(self):\n",
" if self._noise_level == self.HIGH:\n",
" return \"HIGH\"\n",
" elif self._noise_level == self.LOW:\n",
" return \"LOW\"\n",
" else:\n",
" return \"NO\"\n",
"\n",
" def __repr__(self):\n",
" return f\"{self.neurons.N} Neurons with {self.input_resistance} input resistance and {self.noise_level} noise\"\n",
"\n",
" def __str__(self):\n",
" return self.__repr__()\n",
"\n",
" def store(self, name):\n",
" self.network.store(name)\n",
"\n",
" def restore(self, name):\n",
" self.network.restore(name)\n",
"\n",
" @property\n",
" def v(self):\n",
" return self.smon.v\n",
"\n",
" @property\n",
" def t(self):\n",
" return self.smon.t\n",
"\n",
" @property\n",
" def injected_current(self):\n",
" return self.neurons.i_injected\n",
"\n",
" @property\n",
" def stimulus_amplitude(self):\n",
" return self.neurons.A\n",
"\n",
"\n",
"class CurrentModel(Model):\n",
" def __init__(self, **kwargs):\n",
" super().__init__(stim=self.SINE, **kwargs)\n",
"\n",
" def create_model(self):\n",
" model, parameters = super().create_model()\n",
" model += CURRENT_INPUT\n",
"\n",
" return model, parameters\n",
"\n",
"\n",
"class SineModel(CurrentModel):\n",
" def create_model(self):\n",
" model, parameters = super().create_model()\n",
" model += SINE_INPUT\n",
"\n",
" return model, parameters\n",
"\n",
"\n",
"class SawModel(CurrentModel):\n",
" def create_model(self):\n",
" model, parameters = super().create_model()\n",
" model += SAW_INPUT\n",
"\n",
" return model, parameters\n",
"\n",
"\n",
"class SynapticModel(Model):\n",
" def __init__(self, **kwargs):\n",
" super().__init__(stim=self.EXP2SYN, **kwargs)\n",
"\n",
" SYNAPTIC_PARAMETERS = {\n",
" \"tau_input_1\": 0.4 * ms,\n",
" \"tau_input_2\": 4 * ms,\n",
" \"offset_A\": 1.48793507e-11,\n",
" \"offset_B\": -2.66359562e-08,\n",
" \"offset_C\": 1.77538800e-05,\n",
" \"offset_D\": -8.05925810e-04,\n",
" \"offset_E\": -3.51463644e-02,\n",
" \"offset_switch\": 0,\n",
" }\n",
"\n",
" def create_model(self):\n",
" model, parameters = super().create_model()\n",
" model += EXP2SYN_WAVEFORM + SUMMATION_OFFSET\n",
" parameters = {**parameters, **self.SYNAPTIC_PARAMETERS}\n",
"\n",
" return model, parameters\n",
"\n",
" def build_network(self):\n",
" super().build_network()\n",
" self.spiker = NeuronGroup(\n",
" self.neurons.N,\n",
" \"\"\"T : second (constant)\n",
" lastspike : second\"\"\",\n",
" threshold=\"timestep(t-lastspike, dt)>=timestep(T, dt)\",\n",
" reset=\"lastspike=t\",\n",
" )\n",
" self.spiker.T = 1 / self.neurons.f\n",
" self.synapses = Synapses(\n",
" self.spiker, self.neurons, on_pre=\"input_aux += 1\"\n",
" ) # connect input to neurons\n",
" self.synapses.connect(\"i==j\") # one synapse goes to one neuron\n",
"\n",
" self.network.add(self.spiker, self.synapses)\n",
"\n",
"\n",
"class SynapticCurrentModel(SynapticModel):\n",
" def __init__(self, offset=True, **kwargs):\n",
" self.offset = 1 if offset else 0\n",
" super().__init__(**kwargs)\n",
"\n",
" def create_model(self):\n",
" model, parameters = super().create_model()\n",
" model += CURRENT_INPUT + SYNAPTIC_INPUT_CURRENT\n",
" parameters = {**parameters, **{\"offset_switch\": self.offset}}\n",
"\n",
" return model, parameters\n",
"\n",
"\n",
"class SynapticConductanceModel(SynapticModel):\n",
" FLAT = 0\n",
" ACTIVE = 1\n",
"\n",
" CONDUCTANCE_PARAMETERS = {\n",
" \"A\": 0 * nS, # overwrite A to be conductance\n",
" \"g_i\": 1 * nS,\n",
" }\n",
"\n",
" INHIBITION_PARAMETERS = {\n",
" \"tau_inhibition_1\": 1 * ms,\n",
" \"tau_inhibition_2\": 10 * ms,\n",
" }\n",
"\n",
" def __init__(self, offset=ACTIVE, **kwargs):\n",
" self.offset = offset\n",
" super().__init__(**kwargs)\n",
"\n",
" def create_model(self):\n",
" model, parameters = super().create_model()\n",
" if self.offset == self.FLAT:\n",
" model += CONDUCTANCE_INPUT + SYNAPTIC_CONDUCTANCE_FLAT\n",
" parameters = {\n",
" **parameters,\n",
" **self.SYNAPTIC_PARAMETERS,\n",
" **self.CONDUCTANCE_PARAMETERS,\n",
" **{\"offset_switch\": 1},\n",
" }\n",
"\n",
" elif self.offset == self.ACTIVE:\n",
" model += CONDUCTANCE_INPUT + SYNAPTIC_CONDUCTANCE_STIM\n",
" parameters = {\n",
" **parameters,\n",
" **self.SYNAPTIC_PARAMETERS,\n",
" **self.CONDUCTANCE_PARAMETERS,\n",
" **self.INHIBITION_PARAMETERS,\n",
" }\n",
"\n",
" return model, parameters\n",
"\n",
" def build_network(self):\n",
" super().build_network()\n",
" if self.offset != self.ACTIVE:\n",
" return\n",
"\n",
" self.inhib_synapses = Synapses(\n",
" self.spiker, self.neurons, on_pre=\"input_inhib_aux += 1\", delay=2 * ms\n",
" ) # connect input to neurons\n",
" self.inhib_synapses.connect(\"i==j\") # one synapse goes to one neuron\n",
"\n",
" self.network.add(self.inhib_synapses)"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "c82769f2",
"metadata": {
"scrolled": false
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"WARNING Cannot use Cython, a test compilation failed: Microsoft Visual C++ 14.0 or greater is required. Get it with \"Microsoft C++ Build Tools\": https://visualstudio.microsoft.com/visual-cpp-build-tools/ (DistutilsPlatformError) [brian2.codegen.runtime.cython_rt.cython_rt.failed_compile_test]\n",
"INFO Cannot use compiled code, falling back to the numpy code generation target. Note that this will likely be slower than using compiled code. Set the code generation to numpy manually to avoid this message:\n",
"prefs.codegen.target = \"numpy\" [brian2.devices.device.codegen_fallback]\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"Starting simulation at t=0. s for a duration of 12. s\n",
"0.33678 s (2%) simulated in 10s, estimated 5m 46s remaining.\n",
"0.6669 s (5%) simulated in 20s, estimated 5m 40s remaining.\n",
"1.00401 s (8%) simulated in 30s, estimated 5m 29s remaining.\n",
"1.30555 s (10%) simulated in 40s, estimated 5m 28s remaining.\n",
"1.66319 s (13%) simulated in 50s, estimated 5m 11s remaining.\n",
"2.01629 s (16%) simulated in 1m 0s, estimated 4m 57s remaining.\n",
"2.37158 s (19%) simulated in 1m 10s, estimated 4m 44s remaining.\n",
"2.72623 s (22%) simulated in 1m 20s, estimated 4m 32s remaining.\n",
"3.08984 s (25%) simulated in 1m 30s, estimated 4m 20s remaining.\n",
"3.44787 s (28%) simulated in 1m 40s, estimated 4m 8s remaining.\n",
"3.81075 s (31%) simulated in 1m 50s, estimated 3m 56s remaining.\n",
"4.1707 s (34%) simulated in 2m 0s, estimated 3m 45s remaining.\n",
"4.53086 s (37%) simulated in 2m 10s, estimated 3m 34s remaining.\n",
"4.89271 s (40%) simulated in 2m 20s, estimated 3m 23s remaining.\n",
"5.24906 s (43%) simulated in 2m 30s, estimated 3m 13s remaining.\n",
"5.59608 s (46%) simulated in 2m 40s, estimated 3m 3s remaining.\n",
"5.95208 s (49%) simulated in 2m 50s, estimated 2m 53s remaining.\n",
"6.30408 s (52%) simulated in 3m 0s, estimated 2m 43s remaining.\n",
"6.66648 s (55%) simulated in 3m 10s, estimated 2m 32s remaining.\n",
"7.02333 s (58%) simulated in 3m 20s, estimated 2m 22s remaining.\n",
"7.38178 s (61%) simulated in 3m 30s, estimated 2m 11s remaining.\n",
"7.73919 s (64%) simulated in 3m 40s, estimated 2m 1s remaining.\n",
"8.0986 s (67%) simulated in 3m 50s, estimated 1m 51s remaining.\n",
"8.45184 s (70%) simulated in 4m 0s, estimated 1m 41s remaining.\n",
"8.81481 s (73%) simulated in 4m 10s, estimated 1m 30s remaining.\n",
"9.17027 s (76%) simulated in 4m 20s, estimated 1m 20s remaining.\n",
"9.52854 s (79%) simulated in 4m 30s, estimated 1m 10s remaining.\n",
"9.88859 s (82%) simulated in 4m 40s, estimated 1m 0s remaining.\n",
"10.25651 s (85%) simulated in 4m 50s, estimated 49s remaining.\n",
"10.60198 s (88%) simulated in 5m 0s, estimated 40s remaining.\n",
"11.0361 s (91%) simulated in 5m 10s, estimated 27s remaining.\n",
"11.46854 s (95%) simulated in 5m 20s, estimated 15s remaining.\n",
"11.91097 s (99%) simulated in 5m 30s, estimated 2s remaining.\n",
"12. s (100%) simulated in 5m 32s\n"
]
},
{
"data": {
"text/plain": [
"[<matplotlib.collections.EventCollection at 0x29ba054c7f0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba055fac0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba056ec40>,\n",
" <matplotlib.collections.EventCollection at 0x29ba057ae50>,\n",
" <matplotlib.collections.EventCollection at 0x29ba058f0d0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba059b1f0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba05a73d0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba05b56a0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba05c3940>,\n",
" <matplotlib.collections.EventCollection at 0x29ba05d0af0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba05dcbb0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba05ead60>,\n",
" <matplotlib.collections.EventCollection at 0x29ba05f7f10>,\n",
" <matplotlib.collections.EventCollection at 0x29ba0605f10>,\n",
" <matplotlib.collections.EventCollection at 0x29ba0617370>,\n",
" <matplotlib.collections.EventCollection at 0x29ba0624610>,\n",
" <matplotlib.collections.EventCollection at 0x29ba06317f0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba063f8e0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba064c940>,\n",
" <matplotlib.collections.EventCollection at 0x29ba0659a90>,\n",
" <matplotlib.collections.EventCollection at 0x29ba0666b80>,\n",
" <matplotlib.collections.EventCollection at 0x29ba0674df0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba0680fa0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba06911c0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba06a1130>,\n",
" <matplotlib.collections.EventCollection at 0x29ba06ae0a0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba06b5fd0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba06c7310>,\n",
" <matplotlib.collections.EventCollection at 0x29ba06d52e0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba06e2310>,\n",
" <matplotlib.collections.EventCollection at 0x29ba06f03a0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba06fc520>,\n",
" <matplotlib.collections.EventCollection at 0x29ba0709670>,\n",
" <matplotlib.collections.EventCollection at 0x29ba07157f0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba0725970>,\n",
" <matplotlib.collections.EventCollection at 0x29ba0731bb0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba073fd60>,\n",
" <matplotlib.collections.EventCollection at 0x29ba074cd30>,\n",
" <matplotlib.collections.EventCollection at 0x29ba0759ee0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba076a1c0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba077a2e0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba0786430>,\n",
" <matplotlib.collections.EventCollection at 0x29ba0793580>,\n",
" <matplotlib.collections.EventCollection at 0x29ba079f670>,\n",
" <matplotlib.collections.EventCollection at 0x29ba07ae5b0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba07bb5e0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba07c7880>,\n",
" <matplotlib.collections.EventCollection at 0x29ba07d5a00>,\n",
" <matplotlib.collections.EventCollection at 0x29ba07e3ac0>,\n",
" <matplotlib.collections.EventCollection at 0x29ba07efc10>]"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXEAAAD4CAYAAAAaT9YAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAnCUlEQVR4nO2db6xm11Xen1VPQmhgGrsZW7d21EkkK+qIURN6G2WaqkoxximNsFs1dSICU9WtR2orhVKJjJsPvXzC0AqhqlIZC9KOSgjjQlJbFu1gGaIKaeRkDCHjjG1sICRubjOXpMKBfiiG3Q/37DvrrrvW/nPe971z9/Xzk0b3vPvsP2vtc95tv89Zex1JKYEQQsiY/LkbbQAhhJD5cBEnhJCB4SJOCCEDw0WcEEIGhos4IYQMzJH9HOzNb35zOn78+H4OSQghw/PMM8/8QUrpmHduXxfx48eP4/Lly/s5JCGEDI+I/H50jnIKIYQMDBdxQggZGC7ihBAyMFzECSFkYLiIE0LIwDRFp4jIlwB8E8CfAng1pbQuIrcAuADgOIAvAfiHKaX/sxozCSGEePT8n/jfTim9I6W0Pn0+C+CplNKdAJ6aPhNCCNlHFpFT7gVwfjo+D+C+ha0hhBDSResingD8iog8IyIPTmW3pZQ2AWD6e6vXUEQeFJHLInJ5a2trcYsJIYTs0Lpj8z0ppa+KyK0AnhSR51sHSCk9AuARAFhfX+cbKAghZIk0LeIppa9Of6+JyKcBvAvA10RkLaW0KSJrAK6tysj7z13aOb66+QpOrB3ddd6WXThzKmx34cypXeWa6NzVzVcAYM+4Xr2abRF5bG1nJrKp1K/1Rde3/rTaOIeWvqM6eg6A8jz0XvPWMs+mljF6r3trf3Pt9mzS9xyAPf3b+Qfi72LPPVSywyuP2umy0nfX2ljzO9eNznm21Xz15nJZVOUUEXmjiHx7PgbwPQCeBfA4gNNTtdMAHluVkYQQQnxa/k/8NgCfFpFc/+dTSv9DRD4H4FEReQDAlwF8YHVmEkII8agu4iml3wXwV53yrwO4axVGEUIIaUP282336+vraU4qWk+Psue0JlrT6lp0VD1WpIfN0fAiPA2uRQdt1Rdr2LmI5tP6GR1rn6xfrZR88K5hj0Zbe2ZQen6gP7fYHLWZoyN717fF797nPbXnTPY72ONjyRevXcv5PPbc+9ij5/vd8sxpEV1cRJ5Re3R2wW33hBAyMFzECSFkYLiIE0LIwAylibdqUhEtumIUQ5rtaIk3t7TEsLYyVwOv6f+98dyt8fo9Mf6tOnZJY2yNjc/95Tat8xfVb9VXa887ajHzLbp2q0Ye9VeyPbepafK939XW+PherbxkV8s6EM2FN35tXwM1cUIIIXvgIk4IIQPDRZwQQgZmKE0caIvj1vV6NDBvLGBefG+UA6UlD0uPXlmyqZb/oUZvHG9LToqSvSXfcn9z6NXabZsWHTbbV9K859x/PX2U7GuJZ57bv9duTr6ROc91bNuW3EjRd7PXVmDec4g5UBMnhJBDChdxQggZGC7ihBAyMMNp4kBdM+7Jxb1oXGt0PrKjNTa5lruiZkOLfhjRk1e7JQ9G1GeUu7pnblrjgltyW9TGaInJ742Vn5MDJ9sf+VDT7W0fujyy2daJfPbm2Wvbuj+gZT/CIsx99lHr09ZbNJ84NXFCCDmkcBEnhJCB4SJOCCEDM+Qifv+5SzuandX29HlP171w5hSubNyz8zkf2/50nRNrR3eNozXSC2dO7dK+cr2sK1o7Tm5cbI4/tX2d3Li4x4bIxsj/XDf3m/9p7j93aadM26DJtuRxcxs9b5mS9m3tLeWbsDbk8XK5jv3Vx/YanVg7umNn9CzBm9tsm+7buy+8edJ96fp53vQ10XMTadiejfZYX7voecCJtaN77h3vPrP3Y80OPWakaet+rH/edfCurbXH3tP2+lhqc2htzvMQ2Ru17d2j0cOQizghhJBtuIgTQsjADBliCJS32bb8fLf1cp+Av33ajmn7awkRaw2XmrMlPdttaU1V2zJGyxbwnr6izy0hii1hgV6/3uu7avdP6bglxDETjePVtW1a05zWtpz3nptTd06K10xLaGfL50xr2uDWFBMt9nrjMMSQEEKICxdxQggZGC7ihBAyMMNp4iXNqictrcfcrfZztn7XbGtNdVna4tu6pb0lDapnZ+/2bN2mRae0fc7RxXu191Ifkc+Rhp/pSb/aq9fWbI7uj5bnAC3fm7nPXBbR4VvGatG2W22Y49uyoSZOCCGHFC7ihBAyMFzECSFkYA6FJg7UU7y2pvSMWDR1ph6jxQcPrT96/fWmh/VsK/lRsnGuvtiaYrfnuUJvPH5PzHdrauOI3nmK5kfT+zyo9FzDO5/Hn/MMQtvn+dFz37Xq9D2U5rmFnvj8RbRyauKEEHJI4SJOCCED07yIi8hNIvKbIvLE9PkWEXlSRF6c/t68OjMJIYR4NGviIvLDANYBHE0pvV9EfgLAN1JKD4vIWQA3p5Q+WupjWZp4phZDDfgx0y05VGp5LzRz48c9ey1RzHdv3DDQ9nqtEtF8LGpHz9yU4up7Yqh7cuq04PXXkz8HiF/nFWnZ1vZa39F3Yc790evHnPj+uc84LLUY+Ix3rTKL7BlZRtz4wpq4iNwB4O8C+BlVfC+A89PxeQD3LWAjIYSQGbTKKT8F4EcA/Jkquy2ltAkA099bvYYi8qCIXBaRy1tbW4vYSgghxFBdxEXk/QCupZSemTNASumRlNJ6Smn92LFjc7oghBASUNXEReTHAPwAgFcBvAHAUQCfAvDXAbw3pbQpImsAPpNSenupr7maODBPgy2da8lN0ho/7LWp0ZJnpNa+ltuk1C6ilNsEWCxfdW8Mfc910H14PntjlJ4v9MZd27mZo+O2+qvH1NRi2CNqMenRvPTGzfc+H/HoyWeejy097wEo3fdz5nYuC2niKaWHUkp3pJSOA/gggF9NKX0YwOMATk/VTgN4bCErCSGEdLNInPjDAO4WkRcB3D19JoQQso8c6amcUvoMgM9Mx18HcNfyTSKEENLKEDs2S/rW1c1X9uSSiLDaXE2fzeNe2bhn1/n898rGPTs2RXq0tjkf5/Fzv3Z8i25nsTae3Ljo1vNs9ubO9pc5uXER95+75NpydfOVnXa2XPeZ59XOw4m1o7vmXtuax7bn9LWzc3nhzCl3DMv95y7t+OT5a23V4+r4avv5wplTe3zV/tjyOWTb7z93aefaaDu9MTVeDLR3ja2Gnf0sXXN7baJrps9ZW3Jfek6zjZ4PGmuj/af71v7l+vl+sM9MLLatZhU5xSOGWMQJIYT4cBEnhJCB4SJOCCEDM1w+caCcW7gUYxvVy8zJbaKpxfX25DvxKMUmL9p/b96Nmr/W5trYeZxa+1ounFLceYnWeOCemOrI15Z49MjG6JnL3Dznte9OT18t9ra2r43XGn9eq6frLuprpiVPSy/MJ04IIYcULuKEEDIwXMQJIWRghtHE5+jNpfcK6r41NQ23VX+PjltyXsx5D2PrHOm2QPn9jHacVvt69ExNSf/ueTbRak+v9hvNM9D3PGVuLvToXMtzh6hNiWW+d7J1LlrGj/oH4vfO1qg9Q4vsjZ63afsORD5xQgghBxMu4oQQMjDDyCmZ6KeO99O31kfpOPfRGrbXk9a09NPN+wy0/1yv/Yys2TtHConmPrpudtyW+cnHeS5aZC9twxzJoebjIil/bT+tIW81ezzZrjWEsDZ2j4yZxy1JT3aMufJN1FevDFtKs1sKg205n/ucC+UUQgg5pHARJ4SQgeEiTgghAzOsJm4paYBe21KI3Zwt5HO2Eufxe/uJ+tuvMLVF7ADi7fVA/EwhnyvVnTu2HmeOht7al7ZlWde4J7y1p9+aPdEzm0V87AnvLfXTOv4iIaA94by5bBGoiRNCyCGFizghhAwMF3FCCBmY4RdxrY/a1y+VXsV0Yu3666Hsa6b069h0/K73mrCS1qXb6uPcR+5H95VtrGloVoezml30Srea3Vazb9HyrB22jaeb2tfW6Wtnsef0Z288O4a+3tYmfWzHsDZ7/UfzrMe3rxTz5iG6V0rXOJdZ//V9FNXx9imU7LDXMJof+x3SPto6uTz3ob9f1h89pmdfbueN71HbV6LLW+c948WXr5LhF3FCCHktw0WcEEIGhos4IYQMzFCLuKdP2bwTUb0oplPraJ7eF2mHkTYYaWC6L62lnlg7uqPfeZqs7tNqvB56/DwPkd0XzpzaNbbXv7Yj0jbtsdaA85jRHNt449pzAT13ur0dR9thy+fEE3vX9eTGxTAeWJ/T/VuttqVONL7V8r1rqW312ls9Wo9d0uOj62WfM2n7dDtb7o1h63n3Q3RdT25cdL831s6Mnp/o2F6D6L7TY3rfo2Uz1CJOCCFkN1zECSFkYLiIE0LIwAyXOyVTymfQkmMB2B3rWcvdXRoXqOfK0OMtI/9KS77kln5q52p5LaK87V4um9wf0Je32isr5YKu+VezpTam51/Nvppfc3Kq12zuuTdqr0lr+X7MyTcTzYkex7OvxR9g8VcfemOVcvnoMexztEVg7hRCCDmkcBEnhJCBqS7iIvIGEfmsiPyWiHxRRH50Kr9FRJ4UkRenvzev3lxCCCGaqiYuIgLgjSmlPxKR1wH4dQAfAfD3AXwjpfSwiJwFcHNK6aOlvlaVT7ymE2p6dTtPB62927NXc/PazdW+a/pplHc9f27RXPP56DlALbd7yzVozfXsUdND5+SCrs1Nj6a9TD279vyg9HzDo/U69Wj9QPk9uT3PYGpz0KPHe37PvT6Z2rzPZSFNPG3zR9PH103/EoB7AZyfys8DuG+2hYQQQmbRpImLyE0i8nkA1wA8mVJ6GsBtKaVNAJj+3hq0fVBELovI5a2trSWZTQghBGhcxFNKf5pSegeAOwC8S0S+o3WAlNIjKaX1lNL6sWPHZppJCCHEoztOXET+DYA/BvBPAbw3pbQpImsAPpNSenup7TI18VqcrqYl18icmG3PpihWXJd5up+1oVdjLcWmt9jcMn8tPvXi6YR2TubGxLfOWWvMfst1y/Ws7VG/tecHluhcLaa9FC/d+hynhqeDz9X252riXtmccb3j7FOJ0v6FG6aJi8gxEXnTdPytAL4bwPMAHgdweqp2GsBjsy0khBAyiyMNddYAnBeRm7C96D+aUnpCRC4BeFREHgDwZQAfWKGdhBBCHKqLeErpCwDe6ZR/HcBdqzCKEEJIG8Ps2NQ649XNV3bes+i9l+/q5u73Vnrv28t5qa0e1qNblXTnlhwb0TsFr26+sif/c5Tn2fqk62h9Useylvoo1V2WHu7lW8n5t71c7l7O7XwfRDZbdNs8l3b+bT712jsUtY3aDq88Ivsc5Ve3vtm879o/ff/qerlv208us/nHvZztLbTo4brv0jWzOdIjO+x3ODofjWWvvc0lntvpaxTZou223/Fl5E+JGGYRJ4QQshcu4oQQMjDDLOK1LaxWStB4P0/zzyP9k8/+nG+1SY9vfyqWXkFVOrY/AW0/+ZyWBrzwp1L6WP1T2rPFzlntJ6l3PpIF9M9q+5PX/uTW5d64+tjeB57d3k91Xd7yUziaiysb9+zYX0ql4JGvR0QpzE3f07muPtbX2MpGus9SP9YOPXbG+17lebG2R99pfb1L4X6l+bLhffm6WHu13aV+rL+1UEP9/VyllAIMtIgTQgjZCxdxQggZGC7ihBAyMEO9nq1nm3RpW3FrGs3ofI1FU4qWQhejrcCt6TAXGas1Lam2o7aFPNsfjefVK9GTYqEn1UDp/uoda05KBz0uMC/VRCm9Q5TOwuuvtkW+JdVAz71Ua1tKf1Czp+f7WRqjVLYMTZyvZyOEkEMKF3FCCBkYLuKEEDIww2riQD3Nqj5exhbxFo22Z7xS3Vbdbq5GWtMtc1nG0x9b+qjZUatT0zVb0vtG16slrWvJX10/1/H6iXTn6H6qfS6N4VGK9e65RzVRWomovT22dnk+znmOE41d8id6xlTyO7KvNK83LBUtIYSQgwsXcUIIGRgu4oQQMjDDauKRtjhX967FApfis3s1yR6dvLW85EtPLHyvvhvZkGm5Hr3aaM0+297izWWL7lzytRbD3hJX7tlUymHSM0btu9Laf+6j9Z7q2dtRord9KcY88q3l+UDus3WtWVbeFGrihBBySOEiTgghA8NFnBBCBmYoTRxozz+iy4B6Hg+vnx6NvaYZZxtKdkY6raZkT03b79GNvX6trbXxazHE0Vi2TY82XLMtqhftQejJmaHtXzQHTC9z9hhom1v1cKAtfr1mp51voC2Hf9Rnbl/L/dL6rAkoP19qrW/Pz4WaOCGEHFK4iBNCyMBwESeEkIEZfhHXOpPWoTxtM78zUaM1Od3m6qb/Hj1dX5flcfIYuY98bDVWi20PXH9fYbajlo8h19Nja+3O8yuyR58r5SvX6HGs/d48R5q/9tnTlC+cObWrrf3s9W3fl2k1Sz0n+tyVjXvc+8q+n1OPc/+5SzvjnNy4uMv3PBfefFg8v6J2Xp4PWx7VKT2Dsbq3tt9eW3tNPN/yOydzf3Z+Ix+9/rznFfefu7TrfbO2TnTve/eo9/2339OW7593HZfJ8Is4IYS8luEiTgghA8NFnBBCBmaYOPGWeNc5tPY7N89DKR9LK725yTWl/BiL5rOwNtrxWvuLtN1S/u9SnK7Xf+161PKp1GKFW3KEzJ3/KFbdo5bbpmfPQ63MG7e3r9b9FC15WVritHvqt9jg+ViKqZ8L48QJIeSQUl3EReQtIvJrIvKciHxRRD4yld8iIk+KyIvT35tXby4hhBBNy/+JvwrgX6WU/gqAdwP45yJyAsBZAE+llO4E8NT0mRBCyD5SXcRTSpsppd+Yjr8J4DkAtwO4F8D5qdp5APetyMY9lGKubVkU+xrF1No+ckxsy5i5X62PRbG5Nla4ZJ+OqY3ihm2srI61tu28Ni26axSjrGNqL5w5tTOenYtSjHP+l3VKL6bfxlxbvFhf+6zA+pf79DRrPYc27t+OZ2PktY9e/YhavLcXi63HyzHqLfp56XrkslY9PtoL4N2z2n5vXj0bS/Mf7dGoxX3r74c3pqW210P36+2PWBVdmriIHAfwTgBPA7gtpbQJbC/0AG5dunWEEEKKNC/iIvJtAH4JwA+llJr/8yIiD4rIZRG5vLW1NcdGQgghAU2LuIi8DtsL+CdSSp+air8mImvT+TUA17y2KaVHUkrrKaX1Y8eOzTbU+6mVfwZpOUP/7ALg/izXP/1zv9H2ff05klXsNnO7Vdj+RM/9RD9RrZygf7Lrvjy5ItoOfHLj4q6fmxpv3mwfuV3+Sa/nXv+ctVKO7qv001rbqm3UEkI0Z/baaRnFSzVgt3B7W7I10VZ7O1/6J7ydUy8tgPezP/Ip26fvC3tvl9JFeP1r+7w0AlHYXMlWDxtmZ+UPPbaVIkuhf7lOSXLJ40b1tH12LrywXU+q0t8Le2/lvlaRfjjTEp0iAH4WwHMppZ9Upx4HcHo6Pg3gseWbRwghpMSRhjrvAfADAK6IyOensn8N4GEAj4rIAwC+DOADK7GQEEJISHURTyn9OgAJTt+1XHMIIYT0MMy2+0xtW7Nt07rNvRZK1bN9uNb3nC383pbt3OeyXgFW2o4c2eG11bS+Sq22Hd6r4/UBtL0GL7JHt++l9brWtvpH1zlq3zJ2T4oAO15Ln7YP776xRPdGSwqH0veiNlc1Stv2rd2912Au3HZPCCGHFC7ihBAyMFzECSFkYIZbxEtbgqPtrtFWdR3TGcWktsbRetui9TlbN4oftjHDUZxpFAOtY7SjMaJxo3hpHWerXz1W2wJv+7H1dSx4FEPeoq3qui2x6HnM0li5jvU1KtPtM9H8RNvF7bHts/S8o7Yl3PtOePe7vq8i8nfP3mf52LvXa/eotdXTwz3bvRh0fY3y/Rq9Tk/bXTtvXxuox9OU7plVMNwiTggh5DpcxAkhZGC4iBNCyMAMFSdei1/26mpKmmJPPHKtz1LdltjYErW499Z47ta2es6B9jj66Hq0xOu2XotSjHBvfHnPnGgfo1e5AXt1au9zzz6G2jz02Nr6CrfW70IUn22Po/a5HrB3DpdxnWvrQbSOeP602hB9J+bAOHFCCDmkcBEnhJCB4SJOCCEDM5QmnpmbK6SUZ8Prt5Z/pVXb9HK62DFbddsebdOrU9IGtX0t+T16cplE+uOc3Co9unqP/tvrW8vzgV6dttZPZF9vPpU5OVI8bb+k80fPAqJxF81zVHtWEn0PS37ZPlqoaf5zoCZOCCGHFC7ihBAyMFzECSFkYIbTxJehk0ZlNZ1c15+rxXt9trTXds2Nce7RWFvqzI27L8Xi5vJWW3u1SmC5+cJrMeEtfXj03iMtbVrvw2X4WbJvzvWrxWlnW2ox8q1jLvuZyKK5U6iJE0LIIYWLOCGEDAwXcUIIGZihF/FII87/PGz5ibXruZOvbNyz00/WsE5uXNzVn87vnfU926fOJ+xpd1E7z75sY7az9EwgIuuBOqeyp+d5c6Pb2zZ5DqJ8yXbOgOv6pM6brecll3v2zdFTI59K9fKx9Ulrm1qHtfboe8r2Z++v6B7I81Crp6lp6NYum7tc17NjaD+9z9Yee39rn/QcWJsitG+2nbbFmwObS92jtC7oMm987z7JfZbsWhZDL+KEEPJah4s4IYQMDBdxQggZmOEXcfuuQE+rq53TupUXw3xi7eiOrmt1RK2jA9s6mNaQrXZXssXTTD37swZt+7YaXSlfhB1D26PL81ilOFercVtfM54uqu303vsYvQtSvzexpBe3PAewsby5XqTbR1q81Uwj7D0S9ZPHy2Wezp7Hs9cnmhM777X63mfv/vB8sL558dL2/ZdeX/a5gH7GY22JfOl9jlSidO1uBMMv4oQQ8lqGizghhAzMMNvuoy3Di2y9b+nTK2/Z2uuN51F6vVtLn7W6nm3LSKtZq1vbtqztWvYYQP2nbrRtO2JOmoWeNqU+aimRl9GHt62+9R5rmcvSlv/W1NLeFnYvpaz1o+aXR8896bWJ5n4u3HZPCCGHFC7ihBAyMNVFXEQ+LiLXRORZVXaLiDwpIi9Of29erZmEEEI8Wv5P/D8DeJ8pOwvgqZTSnQCemj6vFG/rei63ulekPdW2YEc6XA49zGPZerrMjqHDp/J5XSdvXffsatEFvfBIYO8r1+zW71K6gJr9nk22bi10DsCucMNSqoSoDy/ETIdhluz1tsrnPizRcwtb34Z25tDLKDTQ25qu/bBj1vqyx14fpTQJGk83tmGkuR993XPIrRfuqil91+w2+Uhnt1pzvo7etbX3aE3zrqUDaA1PbNnyvyjVRTyl9D8BfMMU3wvg/HR8HsB9yzWLEEJIC3M18dtSSpsAMP29NaooIg+KyGURuby1tTVzOEIIIR4rf7CZUnokpbSeUlo/duzYqocjhJDXFHMX8a+JyBoATH+vLc+kMi1aa6umBsTb2+1x1gO9LcIt23hrscglrVnru9HWZRvDqzVyr523LdpLkVtKAWq3/3vPCjyfIn2yJy5X66eez9HW6CiWV6cezr7pfnV5pFnb8T2NPtr2rol0VK3velvoS/eY7u/kxkX3/rDPGEro+ddauad1l7bp6/s2+n5692lpT0TG6y9KY2DbttLyfGHVzF3EHwdwejo+DeCx5ZhDCCGkh5YQw08CuATg7SLysog8AOBhAHeLyIsA7p4+E0II2WeO1CqklD4UnLprybYQQgjpZIgdmyXNKUrB2aqvWh04ii+2GqiXRyLSmXX7UipMq+PbMfVnGz/svc7La2fx/Izmzksj6+Wv0X1lrFYZ5b6JsNeopCdHvnlpi3Mba4uN+9e6v6dZWy0c2B1XbTV1a4O12cat23n1NF3dlzdfLd+NaP5q6XW1rbUxdF0v3a/+bK9DpKnb/SN63qJnSZZov4Znm713vHh5r59VMMQiTgghxIeLOCGEDAwXcUIIGZhh8olnSnmqrT7Wk9u6tcwbo5aHvHTs1S2VWZsim7X/LUQ5mhfBy2PeYlfvdfByXXtae+u4NtY8uq8iH1vyVWta8r3Puc6l8Vu/Iy11vPnPfgHl50glormMnrvUbOzBW1O8/QW1fPjefTsH5hMnhJBDChdxQggZGC7ihBAyMMNp4kCs30ZaldbmojKgrP3lcUuaXqtG2aqrt+q6LXp6qd8eTbVFj+/VhGu+ztE059rjja0/e8zVtBetH93XkQa9yLxG81DSxEtjtHyXSteg1L6kO89516gd29pn23r1a3bVoCZOCCGHFC7ihBAyMFzECSFkYIZbxKM8xsDeHCE234WtozWxlnhZ4HruhiinhX0npddHHje/jzCXlfKI6HwZOi+G7s/OUT4u6cI6p7fF5puIcnDYnBklHVnnwPDG0jk6bN6S1r40NodFrZ88z7lcf7Z4797UvrT4EJ3TZfYaRO8WtT7ZvCT2/rRtrM16DH2vRr7afjzf8v3ovZM2mp8SXpucL937l4nGt9c7uv4t39NF/OphuEWcEELIdbiIE0LIwHARJ4SQgRl+Ebd6U6S16vM97xHMfQDX37VoNXZPn7Z2ZU2xVNfqlVp307boc54vJd3T6rX6XZCejqfzrXsxx1EssldP51y273jMZdE7GiONWNvYi9Wrs03e+zA930pxzic3LrrPH2r56q0OncfJdnlz5Gm99h7JfXl55/X8azttbu4I+w5S3c7m2faeOVh/cz17z+c2eUz7rMd7FmOfH1nfLS1atneNtM1RnpVVMfwiTgghr2W4iBNCyMAMue0e6E+Ladu0bnVu2VoOxK9fm7tdvEZp67VnF1BPV1CyuXdOa9vVo23i2raW7ey1c6U0srpOHnfuPRD9fG5JvVq7b1v8aJnvfN7bmm9tKG3j7/3utaagbbm+LelfW/vL8zGHnnsCWGzLPcBt94QQcmjhIk4IIQPDRZwQQgZmmEW8tC1ZlwF79V5v63NJC6ttG9aap91mq8O6vLC4lq3i1l+vD60P5vAuL5xQo9MVRKGA0Vb5SNO0W5R1uGRpq3kpjC0f67A1u7W8dE2sz7Xt815ahlLIYin9gL4HozmzYW/efattsPdjdE9GW9yB3dfDhiR6era2yQsxtTbaEEY9vg4njeYgH9e+H1Fob/Qdifq8snHPrvDAWuhxLc2BLfO+e6timEWcEELIXriIE0LIwHARJ4SQgRkmTrwnhjQ6Hx2XYmOj8QA/5lbbrMta4n9bY39L5V5saksccMu82HOWmu7XE5/s+dl6fbz5aL0Han1oe3rG8eqWxvLw5ixTiyHviam297em9ZVmLXslansFFtlzUdpjUItzL/k553ouY8s948QJIeSQwkWcEEIGZqFFXETeJyIviMhLInJ2WUYRQghpY7YmLiI3AfhtAHcDeBnA5wB8KKV0NWqzaO6USBNcVL9q1VhbiLS4ljwgczT56HOkmbY8A2jJC5Fp0UE9LTzqt1X3LD1D0HbZ+svIx1Iaf5H5KOnykU+2vOdeLj3XaMlRMkf/j+aipKv39Ndj3zL67RlzEW18VZr4uwC8lFL63ZTS/wPwCwDuXaA/QgghnSyyiN8O4Cvq88tT2S5E5EERuSwil7e2thYYjhBCiGWRRVycsj3aTErpkZTSekpp/dixYwsMRwghxLKIJn4KwEZK6Z7p80MAkFL6sajNXE2cEEJey6xKE/8cgDtF5K0i8noAHwTw+AL9EUII6eTI3IYppVdF5F8AuAjgJgAfTyl9cWmWEUIIqTJ7EQeAlNIvA/jlJdlCCCGkE+7YJISQgeEiTgghA8NFnBBCBoaLOCGEDAwXcUIIGRgu4oQQMjBcxAkhZGC4iBNCyMDs6zs2RWQLwO/PbP5mAH+wRHNuJIfFl8PiB0BfDir0ZZu/nFJyMwju6yK+CCJyOUoAMxqHxZfD4gdAXw4q9KUO5RRCCBkYLuKEEDIwIy3ij9xoA5bIYfHlsPgB0JeDCn2pMIwmTgghZC8j/Z84IYQQAxdxQggZmAO/iIvI+0TkBRF5SUTO3mh7PETkLSLyayLynIh8UUQ+MpXfIiJPisiL09+bVZuHJp9eEJF7VPlfE5Er07l/LyLeC6lX7c9NIvKbIvLE4H68SUR+UUSen67NqYF9+ZfTvfWsiHxSRN4wii8i8nERuSYiz6qypdkuIt8iIhem8qdF5Pg++/Jvp3vsCyLyaRF50776klI6sP+w/dq33wHwNgCvB/BbAE7caLscO9cAfOd0/O0AfhvACQA/AeDsVH4WwI9PxycmX74FwFsnH2+azn0WwCkAAuC/A/g7N8CfHwbw8wCemD6P6sd5AP9kOn49gDeN6AuA2wH8HoBvnT4/CuAfjeILgL8F4DsBPKvKlmY7gH8G4Ken4w8CuLDPvnwPgCPT8Y/vty/7+qWaMWGnAFxUnx8C8NCNtqvB7scA3A3gBQBrU9kagBc8P7D9ntJTU53nVfmHAJzbZ9vvAPAUgO/C9UV8RD+OYnvhE1M+oi+3A/gKgFuw/UrFJ6aFYxhfABw3C9/SbM91puMj2N4VKfvlizn39wB8Yj99OehySr55My9PZQeW6efPOwE8DeC2lNImAEx/b52qRX7dPh3b8v3kpwD8CIA/U2Uj+vE2AFsA/tMkDf2MiLwRA/qSUvpfAP4dgC8D2ATwhymlX8GAviiWaftOm5TSqwD+EMBfXJnlZf4xtv/PepddEyvx5aAv4p5ed2BjIkXk2wD8EoAfSim9UqrqlKVC+b4gIu8HcC2l9ExrE6fshvsxcQTbP3v/Y0rpnQD+GNs/2yMOrC+TXnwvtn+S/yUAbxSRD5eaOGUHwpcG5th+IPwSkY8BeBXAJ3KRU23pvhz0RfxlAG9Rn+8A8NUbZEsREXkdthfwT6SUPjUVf01E1qbzawCuTeWRXy9Px7Z8v3gPgO8TkS8B+AUA3yUiP4fx/MBkw8sppaenz7+I7UV9RF++G8DvpZS2Ukp/AuBTAP4GxvQls0zbd9qIyBEAfwHAN1ZmuYOInAbwfgDfnyYtBPvky0FfxD8H4E4ReauIvB7bQv/jN9imPUxPln8WwHMppZ9Upx4HcHo6Po1trTyXf3B6Ev1WAHcC+Oz0s/KbIvLuqc8fVG1WTkrpoZTSHSml49ie619NKX14ND8mX/43gK+IyNunorsAXMWAvmBbRnm3iPz5yYa7ADyHMX3JLNN23dc/wPZ9u5+/YN8H4KMAvi+l9H/Vqf3xZT8eaiz4EOF7sR3t8TsAPnaj7Qls/JvY/snzBQCfn/59L7a1rKcAvDj9vUW1+djk0wtQEQIA1gE8O537D1jhA5qKT+/F9QebQ/oB4B0ALk/X5b8BuHlgX34UwPOTHf8F2xEPQ/gC4JPY1vL/BNv/p/nAMm0H8AYA/xXAS9iO+njbPvvyErZ17Pzd/+n99IXb7gkhZGAOupxCCCGkABdxQggZGC7ihBAyMFzECSFkYLiIE0LIwHARJ4SQgeEiTgghA/P/AY7SpJsQ8V6NAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"#from models import Model, SynapticConductanceModel\n",
"from brian2 import *\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"\n",
"\n",
"model = SynapticConductanceModel(resistance=Model.LOW, # Model.LOW, Model.HIGH\n",
" noise=Model.LOW, n=50, # Model.OFF, Model.LOW, Model,HIGH\n",
" offset=SynapticConductanceModel.ACTIVE)\n",
"\n",
"\n",
"# model = SineModel(resistance=Model.LOW, # Model.LOW, Model.HIGH\n",
"# noise=Model.OFF, n=1, # Model.OFF, Model.LOW, Model,HIGH\n",
"# )\n",
"\n",
"\n",
"model.f = 200 * Hz\n",
"model.set_stimulus_current(400 * nS) # current should be scaled by 100x for Active Offset so (500nS is actually 5nS)\n",
"model._set_variable(\"i_injected\", 90 * pA)\n",
"\n",
"\n",
"model.run(12*second)\n",
"\n",
"spike_times = [s/ms for s in model.spike_train.values()]\n",
"\n",
"plt.eventplot(spike_times)\n"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "7f22e21e",
"metadata": {},
"outputs": [],
"source": [
"#print(spike_times)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "20e3191c",
"metadata": {},
"outputs": [],
"source": [
"import scipy.io\n",
"import numpy as np\n",
"\n",
"\n",
"file_path = 'data0.mat'\n",
"scipy.io.savemat(file_path, {'data0': spike_times[0]})\n",
"\n",
"file_path = 'data1.mat'\n",
"scipy.io.savemat(file_path, {'data1': spike_times[1]})\n",
"\n",
"file_path = 'data2.mat'\n",
"scipy.io.savemat(file_path, {'data2': spike_times[2]})\n",
"\n",
"file_path = 'data3.mat'\n",
"scipy.io.savemat(file_path, {'data3': spike_times[3]})\n",
"\n",
"file_path = 'data4.mat'\n",
"scipy.io.savemat(file_path, {'data4': spike_times[4]})\n",
"\n",
"file_path = 'data5.mat'\n",
"scipy.io.savemat(file_path, {'data5': spike_times[5]})\n",
"\n",
"file_path = 'data6.mat'\n",
"scipy.io.savemat(file_path, {'data6': spike_times[6]})\n",
"\n",
"file_path = 'data7.mat'\n",
"scipy.io.savemat(file_path, {'data7': spike_times[7]})\n",
"\n",
"file_path = 'data8.mat'\n",
"scipy.io.savemat(file_path, {'data8': spike_times[8]})\n",
"\n",
"file_path = 'data9.mat'\n",
"scipy.io.savemat(file_path, {'data9': spike_times[9]})\n",
"\n",
"file_path = 'data10.mat'\n",
"scipy.io.savemat(file_path, {'data10': spike_times[10]})\n",
"\n",
"file_path = 'data11.mat'\n",
"scipy.io.savemat(file_path, {'data11': spike_times[11]})\n",
"\n",
"file_path = 'data12.mat'\n",
"scipy.io.savemat(file_path, {'data12': spike_times[12]})\n",
"\n",
"file_path = 'data13.mat'\n",
"scipy.io.savemat(file_path, {'data13': spike_times[13]})\n",
"\n",
"file_path = 'data14.mat'\n",
"scipy.io.savemat(file_path, {'data14': spike_times[14]})\n",
"\n",
"file_path = 'data15.mat'\n",
"scipy.io.savemat(file_path, {'data15': spike_times[15]})\n",
"\n",
"file_path = 'data16.mat'\n",
"scipy.io.savemat(file_path, {'data16': spike_times[16]})\n",
"\n",
"file_path = 'data17.mat'\n",
"scipy.io.savemat(file_path, {'data17': spike_times[17]})\n",
"\n",
"file_path = 'data18.mat'\n",
"scipy.io.savemat(file_path, {'data18': spike_times[18]})\n",
"\n",
"file_path = 'data19.mat'\n",
"scipy.io.savemat(file_path, {'data19': spike_times[19]})\n",
"\n",
"file_path = 'data20.mat'\n",
"scipy.io.savemat(file_path, {'data20': spike_times[20]})\n",
"\n",
"file_path = 'data21.mat'\n",
"scipy.io.savemat(file_path, {'data21': spike_times[21]})\n",
"\n",
"file_path = 'data22.mat'\n",
"scipy.io.savemat(file_path, {'data22': spike_times[22]})\n",
"\n",
"file_path = 'data23.mat'\n",
"scipy.io.savemat(file_path, {'data23': spike_times[23]})\n",
"\n",
"file_path = 'data24.mat'\n",
"scipy.io.savemat(file_path, {'data24': spike_times[24]})\n",
"\n",
"file_path = 'data25.mat'\n",
"scipy.io.savemat(file_path, {'data25': spike_times[25]})\n",
"\n",
"file_path = 'data26.mat'\n",
"scipy.io.savemat(file_path, {'data26': spike_times[26]})\n",
"\n",
"file_path = 'data27.mat'\n",
"scipy.io.savemat(file_path, {'data27': spike_times[27]})\n",
"\n",
"file_path = 'data28.mat'\n",
"scipy.io.savemat(file_path, {'data28': spike_times[28]})\n",
"\n",
"file_path = 'data29.mat'\n",
"scipy.io.savemat(file_path, {'data29': spike_times[29]})\n",
"\n",
"file_path = 'data30.mat'\n",
"scipy.io.savemat(file_path, {'data30': spike_times[30]})\n",
"\n",
"file_path = 'data31.mat'\n",
"scipy.io.savemat(file_path, {'data31': spike_times[31]})\n",
"\n",
"file_path = 'data32.mat'\n",
"scipy.io.savemat(file_path, {'data32': spike_times[32]})\n",
"\n",
"file_path = 'data33.mat'\n",
"scipy.io.savemat(file_path, {'data33': spike_times[33]})\n",
"\n",
"file_path = 'data34.mat'\n",
"scipy.io.savemat(file_path, {'data34': spike_times[34]})\n",
"\n",
"file_path = 'data35.mat'\n",
"scipy.io.savemat(file_path, {'data35': spike_times[35]})\n",
"\n",
"file_path = 'data36.mat'\n",
"scipy.io.savemat(file_path, {'data36': spike_times[36]})\n",
"\n",
"file_path = 'data37.mat'\n",
"scipy.io.savemat(file_path, {'data37': spike_times[37]})\n",
"\n",
"file_path = 'data38.mat'\n",
"scipy.io.savemat(file_path, {'data38': spike_times[38]})\n",
"\n",
"file_path = 'data39.mat'\n",
"scipy.io.savemat(file_path, {'data39': spike_times[39]})\n",
"\n",
"file_path = 'data40.mat'\n",
"scipy.io.savemat(file_path, {'data40': spike_times[40]})\n",
"\n",
"file_path = 'data41.mat'\n",
"scipy.io.savemat(file_path, {'data41': spike_times[41]})\n",
"\n",
"file_path = 'data42.mat'\n",
"scipy.io.savemat(file_path, {'data42': spike_times[42]})\n",
"\n",
"file_path = 'data43.mat'\n",
"scipy.io.savemat(file_path, {'data43': spike_times[43]})\n",
"\n",
"file_path = 'data44.mat'\n",
"scipy.io.savemat(file_path, {'data44': spike_times[44]})\n",
"\n",
"file_path = 'data45.mat'\n",
"scipy.io.savemat(file_path, {'data45': spike_times[45]})\n",
"\n",
"file_path = 'data46.mat'\n",
"scipy.io.savemat(file_path, {'data46': spike_times[46]})\n",
"\n",
"file_path = 'data47.mat'\n",
"scipy.io.savemat(file_path, {'data47': spike_times[47]})\n",
"\n",
"file_path = 'data48.mat'\n",
"scipy.io.savemat(file_path, {'data48': spike_times[48]})\n",
"\n",
"file_path = 'data49.mat'\n",
"scipy.io.savemat(file_path, {'data49': spike_times[49]})\n",
"\n",
"\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "d353eff4",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.8.8"
}
},
"nbformat": 4,
"nbformat_minor": 5
}