�
�-cc < @ s$ d d l Z d d l m Z m Z m Z m Z m Z m Z m Z d d l m
Z
m Z m Z d d l m
Z
m Z m Z m Z m Z m Z m Z d d l j Z d d l m Z d d l Z d d l m Z m Z m Z m Z m Z d d l m Z m Z d d l! Z! e" Z# e# r_e! j$ j% d j& e! j' � � � Z( e( e j) j* j+ j, _( e- e j) j* j+ j, _. n d d
d d d
d d d d d d d d d d d d d d d d d d d d d d d d e- e- d e- e- e- e" e- e" e- e" e" e- e" e" e" e" e" e" e" e" e- e- d d e j d d! d" d d# d$ �; Z/ d S(% i����N( t NeuronGroupt Synapsest PoissonInputt PoissonGroupt SpikeGeneratorGroupt Equationst
TimedArray( t StateMonitort SpikeMonitort PopulationRateMonitor( t Hzt nSt nFt Mohmt pAt mVt ms( t time( t randomt expt sqrtt logt pi( t shufflet samples ~/.cython/brian-pid-{}g��&� .>i g p�@id i
g >�g
ףp=
�?g g @g @gH�z�G�?g �?g �?g333333�?gffffff�?g @g333333�?g �R�g Y@t lognormali! t ERg �?t flatg{�G�z�?g�������?t eulerc; z C s�
t � }; |7 t t j _ d
t }< d t }= d }> d }? d }@ d }A d }B d t }C d t }D d
t }E d t }F d t }G d
t }H d
t }I d t }J t j t j d � � }K |+ t k r d t
t j d | � � }L t j d |L d | � | }M t j d |L d | � | }N n< d }L t j d |L d � | }M t j d |L d � | }N |9 t j d � | }O d |M }P |O |P }Q |K d d � d | � f \ }R }S }T }U }V }W }X | t k r>|M | t j | � |N | t j | � |O | t j | � |P | t j | � |Q | t j | � f \ }M }N }O }P }Q n | t k sV| t k r |R | t j | � |S | t j | � |T | t j | � |U | t j | � |V | t j | � |W | t j | � |X | t j | � f \ }R }S }T }U }V }W }X n |' t k r�|, t k r>t j d | d t
d � �| }Y n |/ t k rnt j d | � t |Y d � }Z n |. t k r�t j d | d! | � t |Y d � } q�n |* t k r#|0 t k r�t j d" |L d# | � j t
� }[ n |1 t k r#t j d$ |L d% | � t |[ d � } q#n | t } | t } |
t }
| t } | t } | t } |R t }R |S t }S |T t }T |U t }U |V t }V |W t }W |X t }X | t } |O t }O |M t }M |Q t }Q |P t }P |N t }N | t } |' t k r;| t } | t } | t } | t } |/ t k r;|Z t }Z q;n | t } | t } | t } d& | }\ d& | }] d& t |\ t |] |\ � |] |\ � t |] t |] |\ � |] |\ � }^ |" t k r�d' d& g }_ t |_ d( | d) t j �}` n d* }a d+ }b |* t k r|b d, 7}b n |' t k r'd- }c n d. }c d/ }d |! t k rL|d d0 7}d n | t k re|d d1 7}d n |" t k r~|d d2 7}d n |# t k r�|d d3 7}d n |$ t k r�|d d4 7}d n |d d5 7}d d6 }e |' t k r�|e d7 7}e n |* t k r�|e d8 7}e n d9 |e d: }f |a |d |c |b |f }g d; |g GHd< }h t | |g d= d> d? |h d@ |: �}i | |i _ |N |M |O |Q f \ |i _ |i _ |i _ |i _ |R |S |T |U |V |W |X f \ |i _ |i _ |i _ |i _! |i _" |i _# |i _$ |% t k r�|N dA t t j% j& t |N � � |i _' dB d t j% j( t |N � � |i _) dC d t j% j( t |N � � |i _* dB d t j% j( t |N � � |i _+ dB d t j% j( t |N � � |i _, n1 dD t |i _' dE |i _) dF |i _* dG |i _+ dH |i _, |& t k r d }j dI |j GHdJ t |i j' |j <dK |i j) |j <dL |i j* |j <dM |i j+ |j <dN |i j, |j <n |' t k r� dO GH|( t k r@ dP }k dQ }l n dR }k dS }l dT |k GHdU |l GHt- |i |i dV |k dW |l d@ |: �}m |, t k r� |m j. dX |Y d dY |Y d � n� |4 dZ k r� |m j. d[ | d\ d] � n� |4 d^ k r_
t
t j | | � � }n xo t/ | � D]U }o g t/ | � D] }p |p |o k r
|p ^ q
}q t0 |q � |m j. dX |q |n dY |o � q
Wn d_ GHd_ S|, t k r�
|/ t k r�
|Z |m _1 n+ | | | t j% j( t |m � � t |m _1 |( r�
d& |m _2 | |m _3 | |m _4 n |- t k r� |, t5 k s�
|. t5 k r' t j% j6 | | d` t |m � �t } n |( t k ri |) t5 k rR | |m j3 |m _7 q� | | |m j3 |m _7 q� |) t5 k r� | |m _7 q� | | |m _7 q� |( r� | |m j3 |m _7 q� | |m _7 n |* t k r8 da }r db |r GHt- |i |i |r d@ |: �}s |0 t k r |s j. dX |[ d dY |[ d � n |s j. d[ |
d\ d] � | |s _8 n dc GHt9 |i � }t t: |i dd t �}u |2 t k r� de df g }v |* t k r� |v j; dg � n |' t k r� |v j; dh � n t< |i |v dd t/ | � d( |8 t j �}w n di GHt j= | t j � |t j> t j |t j? dj |6 dk |5 � t j g }x t j@ |x d dl � t |x d dl � }y dm |y GHdn t � |; do GH|2 t k r}
|x |u |w f S|3 t k r�
|x |u |m f S|x |u f Sd S(p Ng I@g �V�g�������?g�~j�t��?g�������?g����MbP?g{�G�z�?g @g @g (@g *�g 4�g !�g T�s+ params/intrinsic/active/params_actRoman.dats mod_gL/sigma04_higg12/gLmin%d/g $@s params/intrinsic/passive/gLs/s gLs%d.dats params/intrinsic/passive/ELs/s ELs%d.dats original/s gLsRoman.dats ELsRoman.dats4 params/intrinsic/passive/taums/taumsRomanUniform.datg @�@s) params/chem_syns/conn_mat/p036/syns%d.datt dtypei s) params/chem_syns/delays/p036/delays%d.dati s params/chem_syns/gms/p036/s
/gms%d.dats params/gjs/conn_mat/s synsgj%d.dats params/gjs/ggaps/s ggaps%d.datg �?g t dtg @s�
alpham = -((v-thm1-eps)/sigm1)/(exp(-(v-thm1-eps)/sigm1)-1.)/ms : Hz
betam = km2*exp(v/sigm2)/ms : Hz
dm/dt = alpham*(1.-m) - betam*m : 1
alphah = kh1*exp(v/sigh1)/ms : Hz
betah = -kh2*(v-thh2)/(exp(-(v-thh2)/sigh2)-1.)/(ms*mV) : Hz
dh/dt = alphah*(1.-h) - betah*h : 1
alphan = -(v-thn1)/(exp(-(v-thn1)/sign1)-1.)/(ms*mV) : Hz
betan = kn2*exp(v/sign2)/ms : Hz
dn/dt = alphan*(1.-n) - betan*n : 1
alphaa = -(v-tha1)/(exp(-(v-tha1)/siga1)-1.)/(ms*mV) : Hz
betaa = ka2*exp(v/siga2)/ms : Hz
da/dt = alphaa*(1.-a) - betaa*a : 1
INa = - gNa*m*m*m*h*(v-ENa) : ampere
IKv3 = - gKv3*n*n*n*n*(v-EK) : ampere
IKv1 = - gKv1*a*a*a*a*(v-EK) : ampere
IL = - gL*(v-EL) : ampere
s�
Ip : ampere
EL : volt
taum : second
gL : siemens
Cap : farad
gNa : siemens
gKv3 : siemens
gKv1 : siemens
thm1 : volt
thh2 : volt
thn1 : volt
tha1 : volt
s3
Igap : ampere
s�
Isyn = - (v-Es)*norm_syn*(g_fall-g_rise) : ampere
dg_fall/dt = -g_fall/tau_fall : siemens
dg_rise/dt = -g_rise/tau_rise : siemens
t s Iext = s - gext*(v-Ese) s + I s + Ip*stimulus(t) s- - gsin*(1.+sin(2.*pi*fsin*t-pi/2.))*(v-Ese) s% + Isin*(1.+sin(2.*pi*fsin*t-pi/2.)) s : ampere
s INa + IKv1 + IKv3 + IL + Iext s + Isyn s + Igap s
dv/dt = ( s )/Cap : volt
s Neuron equations are:
%s
s v>=v_thresht thresholds
v>v_thresht
refractoryt methodg @g{�G�z�?g�G�z��?g >�g:�6U���?gxak���?g�{*�=��?g�drjg��?s" The perturbed neuron is kneu=%d.
g�"j��KM�g�"�dT�?g5@i�QH�?g���ۂ%�?gy! ���?s Setting synapses...s�
# Fraction of synaptic neurotransmitter resources available:
dxs/dt = (1 - xs)/tau_d : 1 (event-driven)
gmax : siemens
U_SE : 1
tau_d : second
s
rs = U_SE * xs
xs -= rs
g_fall_post+=gmax*rs
g_rise_post+=gmax*rs
s gmax : siemenss# g_fall_post+=gmax;g_rise_post+=gmaxs Synapse equations are:
%s
s Synapse actions are:
%s
t modelt on_pret it jR t pt conditions i!=jt FIDs ERROR: Unknown connectivity typet sizes�
g_gap : siemens # gap junction conductance
Igap_post = g_gap * (v_pre - v_post) : ampere (summed)
s! Gap-junction equations are:
%s
s Setting monitors...t recordt vt Iextt Igapt Isyns Running the simulation...t windowt widthi�c��s
nu = %.2lf Hzs Done (took %.2lf mins)g N@(A R R t b2t defaultclockR R t npt transposet loadtxtt Truet intt roundt onest typet lent astypeR R R R
R
R R R R t Ipt ELt gLt taumt Capt gNat gKv3t gKv1t thm1t thh2t thn1t tha1R t randnR, t randt mt ht nt aR t connectt rangeR t delayt xst U_SEt tau_dt FalseR t gmaxt g_gapR R t appendR t runt tt smooth_ratet sum(z t epst connectivity_matrixt sim_timet Nt Nrect v_thresht epsilont It gextt Ipst Isint gsint fsint pgjt ggapst gLmint dmint dmaxt gmst gmmeant gmstdt gm_scalet tau_fallt tau_riset Est Eset taudmt USEmt gm_distributiont homo_neuronst homo_actt kneut
const_currentt
const_condt step_currentt sin_condt sin_currentt randicst
one_perturbedt act_synst stdt scale_gmt gjst mod_gLt read_synst dist_synst read_gmst read_delayst read_gjst
read_ggapst return_statet return_synst connectivityt window_widtht
windowTypeR t dt_rect factorTauR" t stt ENat EKt km2t kh1t kh2t kn2t ka2t sigm1t sigh2t sign1t sigm2t sigh1t sign2t siga1t siga2t
params_actt foldert gLst ELst taumst Rinst Capst gNast gKv3st gKv1st thm1st thh2st thn1st tha1st synst delayst synsgjt c_fallt c_riset norm_synt intervalst stimulust gat_varst param_eqst syn_eqst ext_eqst currst volt_eqst equationsR! t nett kpertt synapses_eqst synapses_actiont St Ct kt kpt prest gap_eqst Sgjt rMont spMont
state_varst stMont ratest nu( ( sE /mnt/beegfs/home/rbarav/GuillemCodeNetwork/FinalCodePublished/FSIN.pyt gewnet/ s� >
"
1 | � & $ / ' /
H
! *E '###&
! ! +
! + ( ! + /)
(0 t brian2R2 R R R R R R R t brian2.monitorsR R R R
R R R
R R R t matplotlib.pyplott pyplott pltR t numpyR4 R R R R R R R t osR7 t TigerFisht patht
expandusert formatt getpidt cache_dirt prefst codegent runtimet cythonRV t multiprocess_safeR� ( ( ( sE /mnt/beegfs/home/rbarav/GuillemCodeNetwork/FinalCodePublished/FSIN.pyt <module> s� 4 4 ( !
| 
