{
"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*3 * nS, LOW: 0.25 * nS, OFF: 0 * nS}\n",
" INHIBITORY_NOISE_VARIANCE = {HIGH: 1.25*3 * 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.57934 s (4%) simulated in 10s, estimated 3m 17s remaining.\n",
"1.15726 s (9%) simulated in 20s, estimated 3m 7s remaining.\n",
"1.69823 s (14%) simulated in 30s, estimated 3m 2s remaining.\n",
"2.20673 s (18%) simulated in 40s, estimated 2m 58s remaining.\n",
"2.72607 s (22%) simulated in 50s, estimated 2m 50s remaining.\n",
"3.25204 s (27%) simulated in 1m 0s, estimated 2m 41s remaining.\n",
"3.77415 s (31%) simulated in 1m 10s, estimated 2m 33s remaining.\n",
"4.29701 s (35%) simulated in 1m 20s, estimated 2m 23s remaining.\n",
"4.82068 s (40%) simulated in 1m 30s, estimated 2m 14s remaining.\n",
"5.3365 s (44%) simulated in 1m 40s, estimated 2m 5s remaining.\n",
"5.85269 s (48%) simulated in 1m 50s, estimated 1m 56s remaining.\n",
"6.37219 s (53%) simulated in 2m 0s, estimated 1m 46s remaining.\n",
"6.89328 s (57%) simulated in 2m 10s, estimated 1m 36s remaining.\n",
"7.40856 s (61%) simulated in 2m 20s, estimated 1m 27s remaining.\n",
"7.92897 s (66%) simulated in 2m 30s, estimated 1m 17s remaining.\n",
"8.44672 s (70%) simulated in 2m 40s, estimated 1m 7s remaining.\n",
"8.9615 s (74%) simulated in 2m 50s, estimated 58s remaining.\n",
"9.4794 s (78%) simulated in 3m 0s, estimated 48s remaining.\n",
"9.99601 s (83%) simulated in 3m 10s, estimated 38s remaining.\n",
"10.50688 s (87%) simulated in 3m 20s, estimated 28s remaining.\n",
"11.02289 s (91%) simulated in 3m 30s, estimated 19s remaining.\n",
"11.53009 s (96%) simulated in 3m 40s, estimated 9s remaining.\n",
"12. s (100%) simulated in 3m 49s\n"
]
},
{
"data": {
"text/plain": [
"[<matplotlib.collections.EventCollection at 0x197289d9ac0>,\n",
" <matplotlib.collections.EventCollection at 0x197289f42e0>,\n",
" <matplotlib.collections.EventCollection at 0x19729093c70>,\n",
" <matplotlib.collections.EventCollection at 0x197290a5550>,\n",
" <matplotlib.collections.EventCollection at 0x197290b9100>,\n",
" <matplotlib.collections.EventCollection at 0x197290c6a30>,\n",
" <matplotlib.collections.EventCollection at 0x197290d9430>,\n",
" <matplotlib.collections.EventCollection at 0x197290e5bb0>,\n",
" <matplotlib.collections.EventCollection at 0x197290f9430>,\n",
" <matplotlib.collections.EventCollection at 0x19729107ee0>,\n",
" <matplotlib.collections.EventCollection at 0x1972911a9a0>,\n",
" <matplotlib.collections.EventCollection at 0x1972912a400>,\n",
" <matplotlib.collections.EventCollection at 0x1972913abe0>,\n",
" <matplotlib.collections.EventCollection at 0x1972914d550>,\n",
" <matplotlib.collections.EventCollection at 0x1972915afd0>,\n",
" <matplotlib.collections.EventCollection at 0x1972916db20>,\n",
" <matplotlib.collections.EventCollection at 0x19729181520>,\n",
" <matplotlib.collections.EventCollection at 0x1972918dd30>,\n",
" <matplotlib.collections.EventCollection at 0x1972919f640>,\n",
" <matplotlib.collections.EventCollection at 0x197291aee50>,\n",
" <matplotlib.collections.EventCollection at 0x197291c17c0>,\n",
" <matplotlib.collections.EventCollection at 0x197291d4070>,\n",
" <matplotlib.collections.EventCollection at 0x197291e1b50>,\n",
" <matplotlib.collections.EventCollection at 0x197291f35e0>,\n",
" <matplotlib.collections.EventCollection at 0x19729203f40>,\n",
" <matplotlib.collections.EventCollection at 0x197292158e0>,\n",
" <matplotlib.collections.EventCollection at 0x197292282b0>,\n",
" <matplotlib.collections.EventCollection at 0x19729234bb0>,\n",
" <matplotlib.collections.EventCollection at 0x19729249940>,\n",
" <matplotlib.collections.EventCollection at 0x1972925b400>,\n",
" <matplotlib.collections.EventCollection at 0x19729269ca0>,\n",
" <matplotlib.collections.EventCollection at 0x197292796d0>,\n",
" <matplotlib.collections.EventCollection at 0x19729289fa0>,\n",
" <matplotlib.collections.EventCollection at 0x1972929c970>,\n",
" <matplotlib.collections.EventCollection at 0x197292ae130>,\n",
" <matplotlib.collections.EventCollection at 0x197292bcca0>,\n",
" <matplotlib.collections.EventCollection at 0x197292d04c0>,\n",
" <matplotlib.collections.EventCollection at 0x197292deeb0>,\n",
" <matplotlib.collections.EventCollection at 0x197292ef7c0>,\n",
" <matplotlib.collections.EventCollection at 0x19729302280>,\n",
" <matplotlib.collections.EventCollection at 0x19729310b20>,\n",
" <matplotlib.collections.EventCollection at 0x19729323760>,\n",
" <matplotlib.collections.EventCollection at 0x197293371c0>,\n",
" <matplotlib.collections.EventCollection at 0x19729343c10>,\n",
" <matplotlib.collections.EventCollection at 0x197293564f0>,\n",
" <matplotlib.collections.EventCollection at 0x19729363f70>,\n",
" <matplotlib.collections.EventCollection at 0x19729378a90>,\n",
" <matplotlib.collections.EventCollection at 0x197293883d0>,\n",
" <matplotlib.collections.EventCollection at 0x19729399c10>,\n",
" <matplotlib.collections.EventCollection at 0x197293ab5e0>]"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXEAAAD4CAYAAAAaT9YAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAoc0lEQVR4nO2db6xm11Xen1VP4uCAG08ztqZx1Jsgy+qIURN6G2WaCqUYMymJcEBNnaihU9WtLbWVQqlE7OZDL59s0gqhCqnYglTTNoRxIaktQztYBqtCspxci5Bxxja2YUjcTPElU2GCJZeB3Q/37PG666619zrve++d2ZfnJ43u+56z/6y1zz5b8z57nXWklAJCCCFj8pcutwGEEEIWh4s4IYQMDBdxQggZGC7ihBAyMFzECSFkYA7sZWdvfetby8rKyl52SQghw/PUU0/9YSnlkHduTxfxlZUVrK+v72WXhBAyPCLy+9E5yimEEDIwXMQJIWRguIgTQsjAcBEnhJCB4SJOCCEDk4pOEZFzAP4YwJ8BuFhKWRWRgwBOAVgBcA7APyil/N/dMZMQQojHnP+J/91SyrtKKavT97sBPFZKuQnAY9N3Qgghe8gycsptAE5On08C+PDS1hBCCJlFdhEvAH5NRJ4SkTunYzeUUs4DwPT3eq+iiNwpIusisr6xsbG8xYQQQi6RfWLzfaWUb4jI9QAeFZFnsx2UUh4A8AAArK6u8g0UhBCyg6QW8VLKN6a/L4vIFwC8B8AfiMjhUsp5ETkM4OXdMvL2+5/Y8v3UXce2HQOAs+dfAQAcOXyt2049r/HK2vZ1vSOHr8XZ86/g1dcu4pqrD+DI4WsvlT97/pVL7dnPuq+onO7bHtdle/57fnnnovYimzy0vVG/dhw0mTFZxMd6/tRdx3B07fSl77rt9XMXcM3VB3Bm7TgAbPHDa6eOwam7jl061/O95aOHN5/s+Fe7a1vVHmt/aww9O+bYGZW3vnt2aL8s3tzvjYd33Xu2zLm/9HXs9e2Vi+bOTtGVU0TkzSLyHfUzgO8H8DSAhwGcmIqdAPDQjltHCCGkSeZ/4jcA+IKI1PK/UEr5nyLyJQAPisgdAL4G4CO7ZyYhhBCP7iJeSvldAH/DOf5NALfshlGEEEJyyF6+7X51dbUskoo2o4d5+nOkVbfK9/rXdWo93Uakv2kyujmwXW+2bXhaaEaHr1hba39ac9Wfezru0bXTePW1i1hdORj67Y1VpAPb9rTdek/CaqyRfu+Nq+3rmqsPbLNJ92tt9ezX427702Nhr4PV7Gvd6pO1t1Wnta8CAK++dhEALvmrbepp5paWhp7R+L32LK39htZeWMZ+zw+Ntje6VhndfRk9XESeUs/obIGP3RNCyMBwESeEkIHhIk4IIQMzjCYOtOO3ge06pNYNNZ6GFsWCRvqntsu23YrDjvRB64vVfLW93nevvB6XeuzM2vFUXHdWX9b22D4rUVx1LRddq158rh3nOg6ehp6hFeOu/Yj0UetTTzfv9dnTx6M2dFu2jHd9F9Gxo3leqdfTzjcPT1eu81nH8Ncx8Pao9PzXLPJMhVc2s1fV2+uiJk4IIWQbXMQJIWRguIgTQsjADKWJA/m8Hr046azmF8XpartaNvQ0z2pjSzfP+KPLaFo5MHo6ei+HTOtYJv7dtuvV8eJyrXa6unIwpeFH/ej2a3uZ2P/Wswi98am2e7Hunua+fu4CAGzR+7O6u+6/pedm8oJ4ZaMY+Vb8tqa19+GVzexxZfrv3YutPbeKN+8ivZ2aOCGEkG1wESeEkIHhIk4IIQMzlCbu5fbw4nN1HSCvUVp653WfvdhdmwNat9vTEHvx5pEt2bwSc8amp8O3roFXrqWNt+zM6KdAO4+IPe/52dqLsG1GOnqr3Za/+ry1I9qryeQ3r+da+nmrb2+ford/0suFk42Xb83RzLyJ/G89r5BpzxK1vyjUxAkhZJ/CRZwQQgaGizghhAzMMJp4JmY1G8es6cXV6nI2h7NXtqdnZnNXeG15PmS014xeWn3sxR3XPNRezvBeHL72d8710fVavkXldR27nxJpodp+L55bk4nh17bZcdB9RONo28nGpmvmPlfQi+3vxV/bOPfeGLb2Ibx5Z/cjPL3d8zWrr2fWld4eEOPECSGENOEiTgghAzOUnNKSACrLPiKdDS2z57w+e99bYXLZn5yaXkih14ctO+eR/zm227Y92z3bAGxLM9p6PVl0rbNpDXSbrZ/nLZ9boWfa1zlSF4Btj95HP/1b8pVnrz6mbczIBC0psBXGOPfeicIYPbkia28vzLJnoxc2GIX9LgvlFEII2adwESeEkIHhIk4IIQMzjCZeyT762gv38qhhUPqVUFl9UxOFjGXD0+p3G962aIpazz7rQ2u/oPeKut7jzLbPVmiatcv6lA2La+mYUZ+2nL0WXp9AP7VuLROFqPb2PnbCp6htj0j3jfYhIlp7UV5ZzwbATzGgz0Xhp/oe1Hbr74vsnUT3nvXbwhBDQggh2+AiTgghA8NFnBBCBuZAv8jlJ6MFVzw902LP2fjbo2unQ000Ewfb0iitnbVtr77u6+z5Vy5pklGbWru0ry2rVL1Pa+4aT/uL0r22dG+vbeu7/uw9kj5X/9bj7sWPR2U1NtbYK2fbtI/se5qtvna2jL32dr5pndba442BPrZ+7sKleXD7/U+ktX5PIz6zdnzbXM++0s2ej8bfm0u33//Epetiz7dSEdt4ejsf6j2lxyxjdyYFcq0T3Yc7Cf8nTgghA8NFnBBCBia9iIvIVSLyWyLyyPT9oIg8KiLPT3+v2z0zCSGEeKTjxEXkxwCsAri2lPIhEfk0gAullPtE5G4A15VSPtlqY9E4cWB7jKrG5tYA/LhTe8zT+Ho5IbxXw7ViX7Mxvt5nq116+mcmRWhUvhfPbe21bbX6tO1749TLY9Ea8x4Z7TwTY+2NUZTDJPLD2qT7quXtq8uWHYOWtu35a2P/gf5elI7FtvUz+1eRDXNy99h1waYNtuNR0c+E9J5DacXN22sWxbDbfZ+5LB0nLiI3AvgggJ9Th28DcHL6fBLAhxe2kBBCyEJk5ZSfBvDjAP5cHbuhlHIeAKa/13sVReROEVkXkfWNjY1lbCWEEGLoLuIi8iEAL5dSnlqkg1LKA6WU1VLK6qFDhxZpghBCSEBXExeRewH8CICLAN4E4FoAnwfwtwC8v5RyXkQOA3i8lHJzq62dyCdeNSud36RF9EosoJ1nxNMzM1pnK7bcazvKERHZ07Kj6oEtnS/KFdHTAiN7ra29+PyWThrtY/Q07JZfmtb+QU+DjvYMshps7R/IvaJuztxs1c1o0961n2OP9TGb8ztr56L22Wto74/WfNKfe69pa/mny182TbyUck8p5cZSygqAjwL49VLKxwE8DODEVOwEgIcWtpAQQshCLBMnfh+AW0XkeQC3Tt8JIYTsIbMeuy+lPA7g8enzNwHcsvMmEUIIyTJEPnFgvoY2Rwfs5YNu6WSWln7a06AtvTwNrXcK6jZ0+Z4m7NkUxZdnfZ8TKx+Vr9R8LDY219sr6Wme9njtu6fJ2rbseXuNtc1eP3pMomvUi/PWbXg2Z763YpxbfXlk7rvsfIjaWWYtiHLk2zw4nt+2/eh4FF++CMwnTggh+xQu4oQQMjBcxAkhZGCG0MQjLdbTnazmpmNEbV7flm7q0eqrl9si0seqP1rTzeYPaWns1pdsbHAvvrnVbutdij0NVpeJ7PI0xdb+wlwt3paxRPp0NF+8tnpt2Gtq49B7Mdct21rzV7Mb75n07svonG7Pi9fX88Bq2L38/t59ltHMveO9daN1L82FmjghhOxTuIgTQsjAcBEnhJCBGU4Tr997WlYt7+lXUY6SReOoW8d68ba92PFFNNw5Mee9GHlLa78guj49u2yZnkaZ0dW9nCY9PdnamI2xjmzvad+tmHHrd7XZyyPU2ifpzUGrF7fyY3s+tsjEUmfKZvYg5txbrblkx6Dld0s79+oyTpwQQsg2uIgTQsjADCGnVLz0q5WeLKHL6LZaMsCicoYmCs/zfubrx35tCNTcx66t3169llwQ/Ty14Zpz0th6dmbC86KfqtkUoRm5xPusw9ui/mo7+jr3XuFnw1wjevPW82duWgftc63b8zsbquj5EpVrjUlWpsvMZ+2P/Wzb1L55qa9b/kcS7qJQTiGEkH0KF3FCCBkYLuKEEDIwQ2ji2RCe+lmXyYSDWQ3Q1svo0bWc1bTXz10AsD11auszsD2czKbJzPgZhUt5PkY+9zRXS2bcvbYyYYNee16Zlk2RDuyFJuq258wnwL9+Xqhg7RPYOke8MYrGIKN3e7b3tGLA34eyWnGUvrUVRtnT7+d+tvV7Y2GJ6vTSaUT+tPbjFoGaOCGE7FO4iBNCyMBwESeEkIEZQhMH8ulogVgX9uJgs/HEmVjQ6Jg+V9vopQz1Uubq/iN90dP3e2MV6Z+tx4db8dJRPC+AUNv3fJu7dxC107Ktpyd7xyK/or2CuXHNlijWvJVCwdPkAf8+6unuUVz13DQYc3yfOw971yw6NvdaarLzcyegJk4IIfsULuKEEDIwXMQJIWRghtDEo7SaERmtqh7T5S0ZHRrY/oqoTNu2PZsCtOVTy6ZWmUyM7BytOBMbW7/34nuzGm3d4wCAa64+4PqV1Vm9OPo5seu9ZwiszdVe3f/cGO6MLa04ae+7LafJvo4sG/fvPUsRjVlkkz1vdeco/bTdi8nGj0fXOaP375Q+Tk2cEEL2KVzECSFkYLiIE0LIwAyniQPtnA/6fEujW1Q/zGjsXlu92GZdJ9KE63cv13hvnOZqqRmNM9L8WnsYrTj1zN5CLwZ5kRjvqN2eHnp07fQWTT3KwdO6dvY61PHz7NM6vu4ncx1sP4vO/daY1zxB0XMKui3v9Xm6z1753rz3iOa/JZs3fu5+FDVxQggh2+AiTgghA9NdxEXkTSLyRRH5bRH5qoj8xHT8oIg8KiLPT3+v231zCSGEaLqauIgIgDeXUr4lIm8A8JsAPgHghwFcKKXcJyJ3A7iulPLJVlvLaOKVlg7l6WyR1tjSKjWt+GZLRgOO2u3Vq755sbW2zd5YteKqe+8fbOU+sTbodjPx7baf7DsV9fFoLOfm5V425n7O3or93ppnXuy5nRNA/O5IbWe07xGNR1Zrr/1X5r4H1evfq5e9zrauR0sv7823aKyid4UuwlKaeNnkW9PXN0z/CoDbAJycjp8E8OGFLSSEELIQKU1cRK4SkS8DeBnAo6WUJwHcUEo5DwDT3+uDuneKyLqIrG9sbOyQ2YQQQoDkIl5K+bNSyrsA3AjgPSLyXdkOSikPlFJWSymrhw4dWtBMQgghHrPjxEXk3wL4EwD/DMD7SynnReQwgMdLKTe36i6bT9zG4eoc2gCasbP2WOt81RtrDvLaf0v/rfX18Uy8uldft+G94y+K123Z1Yr5tXG3mXzfXn0vhrmHbkvn1Kh2tGKio/h7/bmnveu6ldZ7H+0eivZDk9F3vbnVq2dt8s55bWW0eP2uz1ZekN4YRPtQ2TjyOecyey2ZfQ3v+ut5ad9JkCHaI1iEpTRxETkkIm+ZPn8bgO8D8CyAhwGcmIqdAPDQUlYSQgiZzYF+ERwGcFJErsLmov9gKeUREXkCwIMicgeArwH4yC7aSQghxKG7iJdSvgLg3c7xbwK4ZTeMIoQQkmOY3CmVnkbYym8QoeNtdb2MZg3kcy7rulaDq+1E/vb8mhvDbvVum+/aq2/t0XsHkc+9uPBsXL2tB2zdL7D7Jb19jUwMsOdDL7Za158zLzPPCrTGrTdWUZ3ePJ5zLT3/e7nhM2227oHMuOm6UY4l22ZkXzQ/M88iMHcKIYSQbXARJ4SQgRlCTgHmh2HVct5PHC+tpW1Hh/gBaEoGunzvsWmvjm5z7uPdrbasXXPKa7+zr2Cr3zMpbK0dWTklShOb6cvzG+inBm6N06LSRybUtCUV9OZl9h7R/WRDOrPpWKPjUQhi5KOV+TwpLyPHeD73rosXNtuS0QC4Y3lZQwwJIYRcuXARJ4SQgeEiTgghAzPUIn7k8LXbQuOqPnVm7birU1k98czacbx47we3tHP2/CuubnbqrmN48d4PbtEdddlqT9W8Tt117FK7tozni0b3X9uxxyNtr5atulvVEmvIZNRfSye95uoD2x4zPrN2fIufntZ39vwrl8rb87VPe72sv1k9HAC+855fueSv9bWW077X9r1+bFl9XWtd236rP53GwM4HG6pp54wd0/pIfGV15SBWVw5u88O7R7w2LfWaVf1ZX6vqx6uvXcT6uQs4unZ6S4oEb9zsvWFts4/j67Gr5auPuq7uQ9+X3v2vffDO6Xmqy1jNe3XlIM6sHd/SXh1fPTfqMT2WtZ3so/qLMNQiTgghZCtcxAkhZGC4iBNCyMAMESfuxXJmH1HOPorbixFuvbYs0kY1cx+dzsQz93TObOysHgd9Xsfr2nS4tt3IZ0//jcZV25WNR9bt67QJ3li04sgjP+bEeGdin6NrMid1bmusbZ/R2C8zh+bcU5FNvTQPUZ+6j0rmkft6LIpl1+WA/nVvtevFwTNOnBBCiAsXcUIIGRgu4oQQMjDDaOIVrVHV2NlWXoNeGtKeXpbRpr02o35b2rNnn35llu3fG5NI56/Ho/wYmqxuqn31Xovm+dsax4zm6uHpsHXsdG6NjF6aTY3bG6PIZjsmWs+3+ql+LVjUl/cKsFZa5NZejX09n7Y5oye38tj09nMy+zeeFq7HIMr349kajYdNaVyxr2mL6lt/a3v6mY1FoCZOCCH7FC7ihBAyMFzECSFkYIbUxDPxmi1tK2qv1Uclq6FndeqWDZaMZh9pjlqb03qsjdfNtjc3nrmnqdZ2MrHWQDsmW5ftxfBm47gj/zz7q6as6/fmiNdOK296S3/2YsOjfQtdzr5qrzfGmTGy/fWuRe+6RXsvWdu8GHpdNoqxt+ORuQfrOW/fYi7UxAkhZJ/CRZwQQgaGizghhAzMEJo4sDUmVMc6A9u1x6rtZd4NqWlppZY5sc2etmdjTiv23Yk9bT3SR62dFi/GPhqDyA/dRzTmVg9s6cWL7DdE9aJ8Ly27Pb8zzxT04qCjmHVPD87k7NDHM3l8ALjPG7TGEsjPheh4b+8qiuXOPqfgkYlVt75H89rT4O1eW8ZGe70XgZo4IYTsU7iIE0LIwHARJ4SQgRlCE7d5JjydMxNvWumd75XL6uHed6sLZ+pm9EdLRpvvtT+3nhejnBkX3U+UhznSTr3yc+K/dax8S5OOYn1tXhgv3roSxV23NPVlteVeu15um1beEdu23eOIcozYvr32WmUzcfy2rLcv0loLbDst+1rHWxr7olATJ4SQfUp3EReRt4vIb4jIMyLyVRH5xHT8oIg8KiLPT3+v231zCSGEaDL/E78I4F+XUv46gPcC+BcicgTA3QAeK6XcBOCx6TshhJA9ZLYmLiIPAfiZ6d/7SynnReQwgMdLKTe36i4bJx7p14vENms9r+YTyepmPd00ox338oq34oz1cc9HzzerD9a/Nre4tSeTBzqzp+D5kYkHnxun77URlWtp+rpML+YfaOer1+Ncz/dysFda896OkZ1fmjnvh43GMfLd7vG0NPvoPqm2aVsj/yp67vbsjuL3s3Pca1uPZ+/+uGLesSkiKwDeDeBJADeUUs4DwPT3+qWsJIQQMpv0Ii4i3w7glwH8aCkl3rbeXu9OEVkXkfWNjY1FbCSEEBKQklNE5A0AHgFwupTyU9Ox57BHckov/aT9eWp/Hmd/Ks551VWlJ6lkH5NvhTplQxf1WLUkAl1Hyy66fjQuPVt1+/Z8S1aKXq3l/RTupVXN/FxvSV6tNiKfbRkvfUJmTGxb2VfYWZtaMkG1TYdVepKIJ7f15KheitjWvdiT9yzZcMFoPO3806maW+Po2W3HLpqXi7KUnCIiAuDnATxTF/CJhwGcmD6fAPDQwhYSQghZiAP9IngfgB8BcEZEvjwd+zcA7gPwoIjcAeBrAD6yKxYSQggJ6S7ipZTfBCDB6Vt21hxCCCFzGOqxeyD3GLo+F4UXLaNLz9U25+jFunxPn4u0Wq15enb0dNWsX9bebFif96h3xQsz08frudY1tXZZW6IwPmtfJgQtGyraa9/z2RubKLTV+tLS1W3d7OvGonkS4WnIto9sGU979trp3ZtZ36JrH9kZzaeozFz42D0hhOxTuIgTQsjAcBEnhJCBGUITB/zUpJXso8L6nCXSmFs6q+3Ltt96xLnGxEavYevZZsu3/JujWXs2zH2NWMtu71ym/17scU/3rbSeOfBs9GzyxjqKqc6OubYzSgfhkdWAs5p81Fc0H7zv0TMK0f1Vie6V6NVoupwtW9uesxfg+ZvR0T27Ij8WhZo4IYTsU7iIE0LIwHARJ4SQgRlCE4/ieIF+jpNsTgpbphePHul+2sZql9dfpFVaf22bkQ49Jx42Sp/a2m9YNrbc2lvHxPbnpczt5aGw+m70+rsov07kqz3eOudpnq30r724ZF0vY5dnn52jvecW5sTzA35+mF67reORL717pvXMQes1f4tcZ3tukdw9i0BNnBBC9ilcxAkhZGC4iBNCyMAMoYkDr2tZr752EUA/LwiQ0yJrudqHxdNlozY8LdYyV3+cE8fa0lVte5nY6J52uaiWGl2rSpTT3Dvn2RTplHNfIxb5Vu1YJMa4ZbutD8T5Urw2Wv5k67Vs0Lb39Pnsns+cmHr7OrbM/LTx+/r5jEX2gLy9Fx3P3/J7GaiJE0LIPoWLOCGEDAwXcUIIGZghNPFe7pJIs+vlPan17TlPq/X01EyMaiaXBdDPAW597cX9ZsvV43WvYXXloBtf3MuHkRnbjIYexTZH53qavxerm80VHtGLb7e04pLnxNTb/CGZZwmsLqzLHV07HeZ7ae0hRPTmRSvG3hsDzx9bPjPHWn232ov2Ylr299ajRaEmTggh+xQu4oQQMjBcxAkhZGCG0MSBOLbU+16J9C3dZqUXt6q160ysaqTn6T40i8YD2/Za+TJa8ctWz4v8tLZmYq0trfHIvCfT2hyRjeX1iK5RJfP8gbXD6tK6fu/YnFj9Vpm5unS2vbkx67WMjtv27LLtAfG+Rq9undMAFsrlv8z9VO1eFGrihBCyT+EiTgghA8NFnBBCBmbIRbyljdd/9biOcz66dhq33/8Ejq6dxtG105fqnbrrGI4cvhavvnYRZ8+/cum7R9WJozJHDl+LM2vHAbyum2qb5miP2v4za8dxZu34tphtb1z05/r3zNrxLZpcbbv+q8fquKyuHAxjtb/znl8J29TjCmzVAfVn7zrVMdPfbd/2mnnjqNuqYxdp+Pq61PHV+ye6zzpvtC92nE/ddWyb9nl07TSOHL52m+5rx8O2qfu3tlo79Xnvux5Le909vLGp/6yNuk40R729ofq96tRe//X5Bc3t9z+x5frXa2L9rG3Ufl6894OXzmk9XNu9fu7CpWc2ahtnz7+C9XMXtt273hjq69V7j+dOMeQiTgghZBMu4oQQMjBDhBi2Qs5aaTK9n+Stx7Szr1rSEg2wNXwtakszJ0ypZb8OXdM/1bOpaFvlNK00ubUd/Yi27ctjTmhmNozRhqtFc6UVmtaaB/Zne73uES17o7BVO3bR/PCuiZ6P0SvqPLsWJZL1eufn3F9euoAoTUIrTLeVFsGzW6eh0DbrNluhhV7oJEMMCSGEbIOLOCGEDEx3EReRz4jIyyLytDp2UEQeFZHnp7/X7a6ZhBBCPLqauIh8D4BvAfjPpZTvmo59GsCFUsp9InI3gOtKKZ/sdbbMY/fA9kfvIw1s7qPb+rincwH+o+a6jV4d23fvkfyett/qdxH/NDXEyoZhRWk59fFWSl3v2nipVTMaa6YdjR1760drvHt6aK9u1I7WfFu6f8aWlr/a79ZrxLy2lrmfvBTDwOuPvEd7GJm52tvL0vV75z2/azm779VLGd3bj1uUpTTxUsr/AnDBHL4NwMnp80kAH17YOkIIIQuzqCZ+QynlPABMf6+PCorInSKyLiLrGxsbC3ZHCCHEY9c3NkspD5RSVkspq4cOHdrt7ggh5C8UqThxEVkB8IjSxJ8D8P5SynkROQzg8VLKzb12dioVLbD1MeXovCbSnluabFaTs2UXjY/22un5EbXb0s+zGn5kl0dGb/VssXV1fO6cV3Nl9emeTz19WbfR0rvnHLNj4PW9yN6GNy722QIvVWrrOtrr3BuPjI/R+Ng+rY3az8x4eT70yIxBy+6dYjfixB8GcGL6fALAQwu2QwghZAkyIYafA/AEgJtF5CURuQPAfQBuFZHnAdw6fSeEELLHbE8fZiilfCw4dcsO20IIIWQmQ+ROAbbrdBGZOHJdTh+zbXjnWv3V73N0Od1OVsdfRPu3GqbO+9EaJ6/dRcarFy+tP/d0y0Visu35yD4bu5yNP8/o9t74ZPXcVjstXXqObp559WB2Drau4aL7P72x1FjtPBPr3tvDydwnto61Z1GYO4UQQvYpXMQJIWRguIgTQsjADKGJR/mwvWNz44Ttd0+v1v1k45S9eFwAbp5n608LL1Y+m9O6EvlmbcnqynP07ozt0Xi2coEDW3VPrf1Xn1t2Ws2ylyenlRfHMmf/Ixv7ru2e03Z0/8y9R1r7GNVnIJ5rkY02l0yLyFbPXm8cenlborY9MjYwnzghhJBtcBEnhJCB4SJOCCEDM6QmvqjmDeRySbc0z1ascEQ2rrSer3pdlGe62h9pnS39Lqsvt/qMxjMTD2zbttpkj0yst7at95xAS9+3/ml673btxRVH87j1ntdWX1GbGu/aW3p7DXasPL90X5rWe0Z7MeH6mO2/d4/05ktvz6Dls2d3b79lEaiJE0LIPoWLOCGEDAwXcUIIGZhhNPFWjCcQa5S9eN5efLQ9Z/vq6ZVeG5FmvUheCs+PqO0au27bzWjGPW2x5UOvnejaZvTbSi/PSi3fixf3bNc5ziP/I725tx8yN79PJk6/Z0vmvG0/smHu3MnMbT02lTm51m17mTVgzpi3ns1orUPLQE2cEEL2KVzECSFkYLiIE0LIwAyniVeyupYu24sxz8T59urUsjpfd9Sf/r5+7gIAbMmzEtWfEzesbW5pzFldtzcOQP9dit6ehh2PXpx+dK1sey27re+R/RkN2tpWfayx/q0257y3snVtda4ea6/nR+S/PVbzjEd5aHq05qGmFSdfieZE7z62dmdyDkV29fbVoj0z5k4hhBCyDS7ihBAyMMPIKZasfKJ/0rZC0HrhabpM66e9127vc8837ydh5ENWOrL1s6GYc463flp6bfTa11KVxo6JveY9OaKGENoUqIu80ivypfU4fe969OyP+rL9ZueJbbNKfZE01PN/jiyl71k7HkBfYqzHeve9rWfJyHEtKS+6fxeFcgohhOxTuIgTQsjAcBEnhJCBGUYTz+im9VwmPC0TSueF/7X68+zx2vNC0WwKAe81bvaz7q8VsmaP2zItn7U2af2NwhJ74VuWjO7v9W1tjkJDW/sQrb4zYaa2He9zZGclm37WG4+ob++aANtfD5jZQ8nYoP2oePdNT0+249EK74uumbW5t/fUsjuTnsALh2zZuwjUxAkhZJ/CRZwQQgaGizghhAzMUJq4fQzY07qsluvFoQKvPxbfewzaY278tedDNn65Pu7c0kfn6NC9fQFNKzY38lP7qs/bvq1NnvY/R5Nt6d22bubZgEhPz+4H9PZwMnpwzyZd3tLT3Ws70T7EnLFs7bW07rmWb9k9sGwZu4dkbeuNae+66vj21p7EolATJ4SQfQoXcUIIGZilFnER+YCIPCciL4jI3TtlFCGEkBwLa+IichWA3wFwK4CXAHwJwMdKKWejOjuliXt6XxRXXctbMvGkrThrrw+tidn6c2OqW7GoUble/hNbr5Uq1bYZ2RW1p8drTjzynBjh7PMDPZ3VuwbWH6vr2vbn+FLP6Tmr93my8e2eb953/SxC73pm9keiOdK7R+fs2WRtitLWZnxsxXR79SytvTndxpWsib8HwAullN8tpfw/AL8I4LYl2iOEEDKTZRbxtwH4uvr+0nRsCyJyp4isi8j6xsbGEt0RQgixLLOIi3NsmzZTSnmglLJaSlk9dOjQEt0RQgixLKOJHwOwVko5Pn2/BwBKKfdGdRbVxAkh5C8yu6WJfwnATSLyDhF5I4CPAnh4ifYIIYTM5EC/iE8p5aKI/EsApwFcBeAzpZSv7phlhBBCuiy8iANAKeVXAfzqDtlCCCFkJnxikxBCBoaLOCGEDAwXcUIIGRgu4oQQMjBcxAkhZGC4iBNCyMBwESeEkIHhIk4IIQOzp+/YFJENAL+/YPW3AvjDHTTncrJffNkvfgD05UqFvmzy10opbgbBPV3El0FE1qMEMKOxX3zZL34A9OVKhb70oZxCCCEDw0WcEEIGZqRF/IHLbcAOsl982S9+APTlSoW+dBhGEyeEELKdkf4nTgghxMBFnBBCBuaKX8RF5AMi8pyIvCAid19uezxE5O0i8hsi8oyIfFVEPjEdPygij4rI89Pf61SdeyafnhOR4+r43xSRM9O5/yAi3gupd9ufq0Tkt0TkkcH9eIuI/JKIPDtdm2MD+/Kvprn1tIh8TkTeNIovIvIZEXlZRJ5Wx3bMdhG5WkROTcefFJGVPfbl301z7Csi8gURecue+lJKuWL/YfO1by8CeCeANwL4bQBHLrddjp2HAXz39Pk7APwOgCMAPg3g7un43QB+cvp8ZPLlagDvmHy8ajr3RQDHAAiA/wHg710Gf34MwC8AeGT6PqofJwH80+nzGwG8ZURfALwNwO8B+Lbp+4MA/vEovgD4HgDfDeBpdWzHbAfwzwH87PT5owBO7bEv3w/gwPT5J/falz29qRYYsGMATqvv9wC453LblbD7IQC3AngOwOHp2GEAz3l+YPM9pcemMs+q4x8DcP8e234jgMcAfC9eX8RH9ONabC58Yo6P6MvbAHwdwEFsvlLxkWnhGMYXACtm4dsx22uZ6fMBbD4VKXvlizn3QwA+u5e+XOlySp28lZemY1cs08+fdwN4EsANpZTzADD9vX4qFvn1tumzPb6X/DSAHwfw5+rYiH68E8AGgP80SUM/JyJvxoC+lFL+N4B/D+BrAM4D+KNSyq9hQF8UO2n7pTqllIsA/gjAX9k1y9v8E2z+z3qLXRO74suVvoh7et0VGxMpIt8O4JcB/Ggp5ZVWUedYaRzfE0TkQwBeLqU8la3iHLvsfkwcwObP3v9YSnk3gD/B5s/2iCvWl0kvvg2bP8n/KoA3i8jHW1WcY1eELwkWsf2K8EtEPgXgIoDP1kNOsR335UpfxF8C8Hb1/UYA37hMtjQRkTdgcwH/bCnl89PhPxCRw9P5wwBeno5Hfr00fbbH94r3AfhBETkH4BcBfK+I/FeM5wcmG14qpTw5ff8lbC7qI/ryfQB+r5SyUUr5UwCfB/C3MaYvlZ20/VIdETkA4C8DuLBrljuIyAkAHwLwD8ukhWCPfLnSF/EvAbhJRN4hIm/EptD/8GW2aRvTzvLPA3imlPJT6tTDAE5Mn09gUyuvxz867US/A8BNAL44/az8YxF579TmP1J1dp1Syj2llBtLKSvYHOtfL6V8fDQ/Jl/+D4Cvi8jN06FbAJzFgL5gU0Z5r4hcM9lwC4BnMKYvlZ20Xbf197E5b/fyF+wHAHwSwA+WUl5Vp/bGl73Y1FhyE+EHsBnt8SKAT11uewIb/w42f/J8BcCXp38/gE0t6zEAz09/D6o6n5p8eg4qQgDAKoCnp3M/g13coOn49H68vrE5pB8A3gVgfbou/x3AdQP78hMAnp3s+C/YjHgYwhcAn8Omlv+n2Pyf5h07aTuANwH4bwBewGbUxzv32JcXsKlj13v/Z/fSFz52TwghA3OlyymEEEIacBEnhJCB4SJOCCEDw0WcEEIGhos4IYQMDBdxQggZGC7ihBAyMP8ftC5NyUsXpFcAAAAASUVORK5CYII=\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.HIGH, 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\", 65 * 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
}