From 9075fe1be283a6ecf6c94c2fc3657e872c788692 Mon Sep 17 00:00:00 2001 From: Jan Benda Date: Mon, 26 Nov 2018 19:47:06 +0100 Subject: [PATCH] [pointprocesses] update exercise --- pointprocesses/exercises/pointprocesses01.tex | 2 +- pointprocesses/lecture/dattomat.py | 21 ++++++++++++++++++ .../lecture/p-unit_baselineactivity.mat | Bin 0 -> 42016 bytes 3 files changed, 22 insertions(+), 1 deletion(-) create mode 100644 pointprocesses/lecture/dattomat.py create mode 100644 pointprocesses/lecture/p-unit_baselineactivity.mat diff --git a/pointprocesses/exercises/pointprocesses01.tex b/pointprocesses/exercises/pointprocesses01.tex index b05f216..2672158 100644 --- a/pointprocesses/exercises/pointprocesses01.tex +++ b/pointprocesses/exercises/pointprocesses01.tex @@ -15,7 +15,7 @@ \else \newcommand{\stitle}{} \fi -\header{{\bfseries\large Exercise 8\stitle}}{{\bfseries\large Point processes}}{{\bfseries\large November 27th, 2017}} +\header{{\bfseries\large Exercise 8\stitle}}{{\bfseries\large Point processes}}{{\bfseries\large December 3th, 2018}} \firstpagefooter{Prof. Dr. Jan Benda}{Phone: 29 74573}{Email: jan.benda@uni-tuebingen.de} \runningfooter{}{\thepage}{} diff --git a/pointprocesses/lecture/dattomat.py b/pointprocesses/lecture/dattomat.py new file mode 100644 index 0000000..f2f31d4 --- /dev/null +++ b/pointprocesses/lecture/dattomat.py @@ -0,0 +1,21 @@ +import sys +import numpy as np +import scipy.io as io + +file = sys.argv[1] +with open(file, 'r') as sf: + x = [] + for line in sf: + line = line.strip() + if len(line) > 0 and line[0] != '#': + x.append(float(line)) + spikes = np.array(x)*0.001 + # break the single spiketrain in trials of 5 seconds: + trials = [] + w = 5.0 + for t0 in np.arange(0.0, spikes[-1], w): + trials.append(spikes[(spikes>=t0)&(spikesnA>pWq_O4xZPFL-{)^E+V>a?TIIrXTc zD-Ybed*x|IowL`K{jR;DvgfsZ23&c=P1hY-+5d){`wi??dCW~$^cir)<&`&Fe`w`t zH(XzN`VF^M?srh-zK0yT|9*$=dr0Me-S_R0`2X=Q>!j1laoCVI6N!9m3-LcI4wLtP z?f>_4_PM!#ziY0zIg#LTy|&nN9NYNs|32>f>`mhpp{*nt*Gx#qwJ<_%4+r+ zJlgJlj6Fu%8)FV`nRxlfcHu|gta|nEF?Q6w3x6E)<)=36z4vEbdh2KQ(u@4}x&3tO z*RO58&sR3#=hu(>wDBuz_}Ao%Pnq-RSHt=oUx{m>umu&mSWFFWH$JLtV5&%5xWA7lGZR@`@| zYX{!(lf8P$?!P2|{Mk-fk$XVupvktsE_AZ(bN!oH&*~-4-YIZ13`iy^5=O;O#%`d-a=YJO9;s_toR4+NX^j(*`b|Y9}tZ zW3vxV_|x9zWj+0;6_f`QoG*w)e&zJDP8q z5yzPk+h^LsUR!?lM`D)UdLw`6S+RYVO-=lrUG&Io8#wfhoj-4%ZTh*p&9Srd-hBGd ztiSBwt*bArJnAny`}rN8xaG#bYp>?Uyw_eS7qLyMtfnm-&W+JCE5gKlZD!ZTaUP zR%1=qo_p5^pVru-gLvWvw%0M+u3UJ~f;i3sJM#8kBMMo+UQp+*SL0 z`H!X6F!i*{lD(GNO5E*jzbv%_sy_337cRAT>)Ei(ZpK|O$d=h_zrOqX^mCV4`CcEt zvi9v|cG%Ou{QJajmRa_A-ng~4^T~JKx6i4yb}b*5Pil>$t@)`oUdM9V{h}u)uIsuy zJ~zv4-xv4j^4zrLHach7V_jD-H{Evot%&#k3ftU%m^=3V6|sGV;U=jWnp$a<_$wWK zQLtKxZ9#phy`&AR92R@;Gyv7EgcmAbo?9spaCT_TAjopdaVcV5!ESHyCwl=n}wZ_psuX*XV zwT8*zKNtMDHuhUAl^~R_2=JV?9UVg0y z>TM_eK=rmA?!K+7>+O6#K7ZH8`?lUTtX^4t$j%KmWjg-whc(!q`1OnWHrO?%H9m3n zfQIsC-t)mN8ttINzdGyr9UD#4=DCfw zeD)|jQH^09yuC4AN23k7K7C*H&y9B0KZZR#dO>V&Z4C3sKbow*>E!bt=+R`=cWsAn z+7!oWvP({xRB*)0P1eL|<*lZ0o*A)yb(5`r_3&=fwq6&nW1UUrW4hP6IG?Pue`&&8 zXT8tinb+A?x#!}RU1#^d${To{y)*oD&n#GHlZJfx{NRe_aK3Gt?P~lrembW)wD)PY zFXlYbam8)Tb}S#~mzwQ9JY7|9HQUj>mrU8|uV%aIO8i^eY`m*4s#qVccenL%ob_RT z8noUxE1h=#dc&lBOZH3a?OEO?4x>jzgBipOK3m5#kdM>abk;o zsq1L5`86faymN1hWpVm@t;NppOrLN6Y%xt^^IJ5VmS(rw4wy9>%Uk0(t)bu9t#Mv$ zwOgL)l6&qGtzljt(i-NOUt8_d`v*<=WOl1HVphDgV1rE@$scM%Sci_@U`se%Ub-Q) zKe!>zTN_N%+m{pR)5GXfrHMSaI8I zP1)DoUwpUCcI0C)ug#7>;ilG#tRp-pZ+3+9OzsHJNo|Mi ztA)A4u6^k?tUR4ovP;*yA3LNo&WoL)eL$z(I0`dHRcD+ZJMBlzVuPl1#^Utc9 zZ`Nrt54PB7i#OZurzN)0v|OIN(O%?jdB;Xe>p|OSm%TpYfVJZ`8aiQu+i2_f*v!~y z_dIvZqt`Cm7=NFHf5d5M+k`J%{@ll<`zBoLo+bRjEjZyM{02UbeG(q*nS^V+yA!Tu z@7{R)(1f3fVTQhu@Gy@jd?n3>@d@AWpWXL(=g)-i@_GA%O@Adk%nJz*_D{mcoWJh# zImKBX?1ijwz8$l|IQwV$KRErJoaM5cddBuXS^grdgp;q%axKqyXL$!kpLtJ~>v?@H z%Vm$ep5=ev_{bBBKhAQkH{ZnZe$Mi4(_UJ--Lfp#Jg_>;r^7nCwk0ckzs@Wl`p?(i zd%sJz%RJdJJGA%6_6OfS@8|uF&5qZT9nRAy+c8-!f9|Gi-*4#$yDq*j+ZW%5iQ~y^ zzk>$W%h|pHKktpB;yB}C`|sIaj@x&SS=nAf!)$rB_v5m;F56|V6zBMPTy86JT=qs~ zPB_ngIesE<^MiBzI4%b#<+$dz3v&E5OmiZ@_Bocs9pd_r(Pc z&vDsb@8tNMu)Y?2k>j)mKb#iFnUmvpb9!B#<0J4i_sPliy0?Eg<=~=Rzq9}r$kw?| zt7i3{xt@nV^WjJ2dQX}S$Hevva{XOSq`h;!`EyQ0xgPHCT!*pyZtrJvJ;gg;3GV}BD{TdqlnS9@kH^3hGp6YeTm`Mla`y4KJ$K|`$ z!*lcfxho#}*Tom*JFSnQH|M*a@7weJG~CuJ?$3AG5wGRD?5a2Ny(i7OPxAfK%lWVGE(QL;Pv=Y+xl4h|emuOuWuF{X z;KeiF?0>|`1^(3lPVfbO)%0ao?|Ey1%MZD;z-Jz?w#RV~7R2YNz#ESG{mSWY7r4yW zaRn~BenLFYtOA!Eyr96J-~zL@z%QqpQcxJLr_gzsd+k{0DcY_7DD-Uy^}!TT=(77x zD)e@mt!EWF?viuPx~edoXJDbrA9=9Q%FSTWl!Bv*8YX_49t; zeRO?scncG<_z7W<$E8koiYQsFkQUh+<{%U&B(?6P-%D-O@$v|{gsxwKn-ak%cL zc%DQ_=vQ76o|~;oTy{rgi3k6&#OWg^yu@YaTu|cI(Ve-hBs_P6N}S&2E%%oAmvW0s z{M3UVKC1VK68{C4yY~krK8fzpHzmF+f979GTz=Nv66b20y|TpL$@uhssBkgr!|gKSmv^yGG(p@ciS?T9lKXq=yyn& z(`sG*&oYS?%-iT*eq83qb7CG>=GQH2 zd#8SCS@_;F%Hnes+gr-~sS2zton=0o?%WpT{sP_BZOZ+^qmQhe-L2f;dilG`(fgJA z8k+AXm;0$SkN;Wj!UC7ae%F?}?CF8!{u14)hs%9E9-mX6D-Y*;z1$C{`}=jdbJO{4 zd3a8LDtDQqi^`pUzDHM;JO23pp50OI8!Gr&fB zuT?l!<%S*3s)*NB;j-JWjQy^y@KTyt4_Acs@5u_6{r!4u9~H;>IF2*5!Z~d<&y4M> zD!c{W_>88Cu)b}maQSbANvHYsc}X(ZA6q5;W8#Z#lTPcZV)vvUMe}&iWLSUpO@?{o zu%ye4J38s7(2P1R>GF?zCVhYYI~)D=#gZRTTwn+uMZl_dO zzjjT9BPYUQ~or}o&)br`C7PT?>`*Jc`WvOE*0ndl#5%umGW9X9`C1I{_ck<&v0QH8~c43 z`~8se--sV3r(Aq#YAV>fb5kDtpE%Cyls^vd{Oh`u^RItJYbwm=9jSOfr=8Q+l4Lr3 zuT0u8d;R;w9n!%L-YFgK(|yt|zihvBxE~KrhvSb)yZFZmah#s%VBemdcA8E1^-g;~ z(1P46(&4`Emv&O3mK)Q-P8gK-t?2IElXhXa`(yi47?(fb3+Zs5zZAzAkq-9TJLpfl zZ%o?ZgwOu%vvjbpC#1tX_+#4XozMO)9qjeFY1i^RKkexeUrpF@X*#Y`>0oEHr2UUi zV6bmUdyg64Uw1$v)!>LxZY*LbzhYUcHNDcVDH_Q3G2`unXpb(WnBK>qZz0B zkr|wE(#$W0#Orz`<6qI>9Fg(wD*-G&%Ea|C6W67T%Z-|l@z;T2KKMSK?~iz%Iho)G zEY1YGU|Ggx#;nNr>?7dWu8Zeu$vF9P&HuueXZ#PoJR1EIqX^RaTLR`_?g95$VmPpC z`dLe8;d|Z^0o+*PbgWN@47VyxC~sRrj>EfFwaW_|2jBgkRdXIXeY91>J^8rwBda3+ znfb9*i`R^`D!NTKe`*OWbf25^KA-Z1B_3kKmsU+Xb?(=eAP#xu8>^u^QZm6RfthZo zn_$D_cTTh#nt4ZjYgNRcGro)0{kjryX|rFR}ess}i98)e^^&MEz=E!u{QPYayFrqh*IpvF`G-r&y(c?jKeugX|Bh z<~qL1RO>IkIMup~ga2uD;wOJv&qX7!Qcknc#O-HHGnZX3%?8W=oo zV0hUFy*IN zt;z3dthIrPCoH#G+G`gsw_Y@N|GnIL%79#M!{~P}T^{?bupV%mPPl1BbEQ?ufLdvj=znyqv|2u{H?FeoaxB58iPwHxWs}6$R$Cpp*B+~5 zJGQznR$CWhn?#^967k_EtdhX>#;w83(orRyuH5Oh5z~{Z$H6}5ohcjaS%!62ZLn)`dvwWd zi{rG}Kt!EJ9MNXO`8c1`7S4B7n?0-NugxAKFZ{U86wdhxpKF{p>%ss{W1H1*nyG9z z{mebu?HqEhOWUm*14CD~$M4sU2pR~|6YVya^Vhf7^17$DTctQvdmN|3#J%?Humfm7 z9^Dbg>4?X7Sl|U6R;Rdfhdspg@6QgCpRuIFdebi6s?%B+PTRH9F49EXX_}s|?zC<+ zBX8@pcHAxlp6RsF^sB~q+O-0(opB!MG?_7*ZL~__k?l9yXwcGEPum#djT__oxzT#l zo_TMh&E+~fcB5^}d2GT)lUdWY(aw?mkno=sk4g9w1k{y@VDEQJxWbG_BwS&F;}gDy z`0j#)_Yi{{X8Ayd>$l4a{BzeVpUOD!p;_LaoZn-=^Rj#w#@{Z=@>>bOuFvxB zy5F(__jok656yC!BX4GTr2@yX-}hOO=VbZb3K(X&rpcA@d>gU?*DTKtymjmBcs=0M zw3iOZc37I=IN7fF@afr(+2@6MXJ5l%DC*$7@y|!lpNQ3 zF+cWOnd5hG9&XF=YOSlek(1{7-2!yEq5Uv$Yt1j<)to=h%?*5{Pj29@*TsH!=K4;= zFpuZ@-js(tpX-ZdujKk1-N)eGAc`)znbSx z9xy5|URUh*Q(oYiv+{fw`thss`~m!V&R>`3xdRco>dbRm9bGoh_ahZR$q(n-FF){t z!}4A6vQzWJ`OeDs%V^j3%lD=53J$m--{p=}#rCK3UGe2N^Idl6JNZt0HR_AluR7oR zGQK)3-$Pt4-(@$f2lrMyq9C+yR^WFm#|_u5AaI`p3VbK7|3?(~RQjtI6u6d;z6Ft& z#`e1k!hQclLFA|hk3JYE9VhS8v zaq~i#b$>{qH*sD)vCxCvROsVqUiU2wJnEW4-;58|fWmOThho1c3SHs$Hws<5KD zfQH34;N}G4QwjrDSrXgV6neGvL4{r;JEh zD_&RPHHtfygzxuxiOW8pUJ|(8iV`1B|7abKXL!1;#APRJR_d1{ymNb(QV;RsQeVtJ z*pa3Fn0RTaOHOrOsmo69SL)kR?{Pz^w~)KuR~nv&XVG5)+ftX^KC0AtS$lmP`+bWo z7rg1E9ymg&2RpAca-FjHe3iNU_8rT7GxDI_%DflCwT~Q8=A0h79aHAxs6RL#T6?)aMA=t%g(n#xDy(qf z^&Ki){`Ic-Tyeh&5B^d`xc|?taQS)XSH$b82=T>1ah!)LTz=jY6~2oC$Q4dG(3_)i zJQuc!6)wHpkJy(0DfU}Y;RzZj8!Mb&w>~Es?u&|Kn4i~XGwHH#A5I3l?6Ek`bIA~w9+vc;3}n3to-KPT=?o8_^KmkA%Q)UR@Os9h ze@ce9&96xxPr1X4q#w^4Y<3)HQPPF8mL^?~MP2OI49-n`NJs3KlM3S%rCeNQi&R)A zwoWrL<8h;&lyaJNi8E4uuk>1}@OS@L zDy}=JU{~Ih3hUQxDZhZv(Os!vM?RDa@$^UII8UZRd}(OP59c4^)l`V1zMTr|(3q4T zO8#1%a>{yk_%>cAIQK&8Q~rqkrlrEVHV6Ej`md$n;nY2WcPA-KSf2{Z_iAT9mJV^P zr_(-2x{7p&GrgMjIpPj!|ABt?=yaHG$Hw!3cgtS-4!oXz*VNcPJsqBlS!o{%s&rFr z+QkW0ro(#Nln!K!H-qIy(3yv)vLTm@;{*HL5 zdnP;&hh=;z^-u<0%z)A@dEhx+#>e)rtx5sPYO59)t+p!8XQxiEDmgk6EHRk6g>S8z z>q6movEO&G{d=pHzXZN5f$5L2eUeqv?p-^{s>BC=wkpPZs(-d>>0~EcHSy-vlVkgD zhF=+=^EXRuNf7Y6RVz?D#S$_Yr&zUk+!U)KU-$}~o3hyze?)FE)oLg^&HXcS9&FYA zr`5>8on|$}{l`zUYVlWaZT();WBUwq36^JA4Rs#NW>}>b`k8SYaBk`i_MR2v1mN59 zcV}BQ{kZdITW`*59kZ=Q{Bn-!5=93FaWoqti}`irbaymV2FPb{+j z;%+hpK{=Z*r)zBe~I;%ZeocgHE}G7aeyT@SOP2Pu$S-# z0_X0azw3hSxA3%1lT z<*qcvG51|*NeK^E+A#94X)7%u;rgoBzAEtIlULbju0Kz#vMF3we_9pCS#3S3r`;Ny zTl{9V$-cT39Grf_h}G7g@}=@MRxKUVns_~Hf^PZCHP&BxiZwP^{?A$)t>@EZ0X5=uD2xD+uC{y@zw?# zjj-uvdpB4e!!U)R8>}zo@3jrq!s(*1!Q|IK zm#yjIxJFBS!U?R=y2$Tq4DHW1T7S+~ug8AhH=5S-RgLj_noREEPE9sM4#p;3w&tPx z!MzojYO<~p)-^?LitX71#yE~)sG@atB7g3#>ueVNy+hW;_{h4@4n4M>|H135C(}%z z%dV9pi+$j^9?RBQ4|2SmW~*b^q+7Ffm2R}zTBJX0j`5ggt0f17E?fOZ;d9hp>Du+H8G{v#k&DvU9+-dH9vsSG;Atjh0}1eT@ID59eF5-UiE21^?#! zwnIz2u9g^Q0rzIS@B;8`hV2Hmgt*y@EioR{Vv|n61Xj~x@?Yv(Ov}#>t?~YCjq$wJ z(0*fUh_4I+@7D8!e)LZ!VynE5*ss1d#)Vqre6hiL$ne?_=F2NK#CQ_+876*uLx_tE z-(UloSMuuy8$#LCtPL>^)n={)j z?XX^4;7;wZDb#iL?XbRB1s;DPwvWV?`A}bU#JE|9%~HNjhYe=fJq6y)aP|S6As*(P zG5?}7v_I8p{l#NCLte$q*q*o1dh+&8ZVdD0{@~jJ$Qx}4=k@D0T6f5z-nxBbxXzL2 zFGF~v&C-01{^6f1ELj4+O#3gB2>P~N6W+r7p@R~xc=O?j$dwWd{$$qM@LH)6lh zS-zM7hOe`N|NJMm*_J5&W>t*bBXO8!x z4*W^4y%@bEY*6eIw5w;=(y5FW%Sje80r@g?Z8M%kvZI0Jr9O3)3~B zr{NV5HogX;T!TCO!f!tH_eVFpw@&o^bo?HH2|NMAe&~KA_J_BAYJx9Km zN`ROjc+EHYQCFTHxX|o;SKM;6qT?%>B4 zy7URd!M7RU`?}B-PyMdYd+EV03_8-Kg)V<&RiSsIo~s4zP?g-;S?FD<2ZxS(C;{m| zz^lm}x)phM+6B;Y(_H%Eq@vL8!lLlJUI{KOoko#Qlg_=!l~4CrkynXF#c@UzdC(u> zc={7RfuqyUsEPfS7P;cO9ng0xUjlrZ=`LFq2miUUIK&(G1;?hla#C?P&#AG!cd-W@ zMX^fJ zx%`omcs(Wl81u)V<5nDbw-PU;9Ocv!C;sSrW{GpUe*8S}Yv~V5TxNGw?DtsId%si? z0QkXxfB*ZnL=e}Gz_Yz;FJXY{_x)0gl+&s^grJ>)pr5^ma zQXen79Nb!d4EQtw*QKQ~K2jRiuIcZh>&jhu+4q#YbP9vfUx_BzXB_dv z^6=b1*R6HxyYiqz{|g*ldN6#>c;!0qZu!OKK7jULS%tUL4c)OK^xGXATJ|HjHvQS- zEBsmM!z$uBRuQiM>I$!+A5<0lJx~$WpVuq=A_e>^0(btR!n-NK0eWzz5zne{nO|!w zyeAV1)`5o$v?jwkl@I<)KdUI|BN5j5E|m=T&ftY565}wkAibEEH`&-(zlo13Vd46buw^R=(#1M`y=UcN9HB{E%5_z zc-phI@p>AQ{wg`h`ee`@fOD(85L}pf=q0J3M=wuBT|&ww`vT|I_v)H*F3$QNUVje=g-rwRl;N8J) z#&HNgt*lA;W%75SD;Gye`IXeAHm8EU+`dWIolXC~BJF31Uo&K@k zt>Es=p8)TkM!W7|aBRk92ZK}7UVJVc;<2x#!+rgF+IOPk_EJcoB?g8t%uaA@kb9t3x%;r}$wNBiWtjE}^!`ta+q-*D*0 zX@|TU+s9^n4<*FK>zM#f&iKT4;PKQ?{T9cYnsJ%ee`kC@I;xAo-bB_s@a&Jy~kc!gjKdhU?;W7P`$K-Vq9;Cri*Uit^C=6u=&-YsDVI5*|J&-`RAeK@wn`e*%Y z2?dfSn@b-)*%G0^g$2O}c(=~?yVZy*f_E$N6MAk19>KY#Px!+Us-rN~dei>BZmJDq znta)x=F)xtX;q9{)?h1rIP~0_U#7)=)2&MW?{rHjF=~b-$k+OUcT1m+E#qVFL)RVk z+{08a1-x79>ns~a`S+2tqwX8rn>_9V@NNm==2+q!3d-l0OOFojtqK0Gcs=0V8fUHz z(*iKpYE;i-UOewS8%AB~<$qg9hy(Z5_XqdZIP-0=^wHqjvUkC`6_Ba1{*YztH5eQ` z>bk2{AEPF6)diN201AAYWf-1VU{y+J!spW&&{$|amFT?CMk@dT-W_$_!=xhz@7DP? z>AJ%>i{dzotf%yNi=#dqdTtqbi%ovOpNlPw1FkLm0KA*?#bZk>7%EF_7;m>N!LwEG z1w32Vvoz|%v9Ix$#kw-fEEqb=q7EJVGI*B7ezkEt@NMc>2h>`f_)%@_2VJ+o1$egZ z_vLXsaBb$JKZ4KoJb;6%-wLay-*pvux8gx7Y?wIj3QH()6FggjYGwbvBxM?@@K( zeBaiYF!CJezIET?b7Uzzxo^FNiMBpm$LaNPoO-KNf^dDzbEuEwfPX74gU>Y%K35`O zL+H0>gH_7GYOre6uV^rZ2cLzmoA%9TQRh9e!MZbHcMd-1dRowE-K8(WR&{h5!}vEg zT9Uf+L9zWA=)9>L_%!OgCqmb)1-a3Nk#lX;WC`Yn@7WY}*G(~w)nr;<2jO!i1jqKV zO+ha{5nFw)CaY5&3~*`c&bMD@RTN|$ge@2N^P%gO4t!m>?gyj3dn~q^kHNVaURb-% zMr$5#4)eEzPm5z>%kWRH<`8Fk8tuxLS z;M?M3E#Z7K@i}!Ki(0I^;)Sh2-~A8huxX#3f;uDOD5x*Oyp2ckx$-^mIrHN`X|*8e zgZNS&)6!C z*&fHizT(>LG5@69M$4dqUYmaRz8%(|VXY%N!g|&lTg9n5Y#6!a%N;Qe0$nz7{Cw=o z@WpWgIGvVY-I%RAEy+CElR87cOHoIJ`1hJltC4QJGtTp!7U~a!_E{$>2RdxU z#W#jLj1$1KrHkJf=JR_uS{=E=FzC6d%N>n=GGw9SmJSfd^FGK2uV!94|Fey{rzuJ|hS+5M#>!{@A*F*Xsn$rs?- zcx%9~nK+S`<6X&3O2CsDkJu8tnf_+ioN&GapueX3c_O&9 z9K+b&7oRJ>f#Vogz8U?Q|1mhn!}E-~B#iTojQ!pRpH`v z--2^f-}QNJtOtO7=J8L@4LYp3vAs1n^2XRN7xhF`P$SP(k8H0z??FDWf1Ya{g+5#3 zK%YI1_Vb0PCqiBM)u<;T9tz$}UHSds&$Kh12bZRO_;#KrIT4J*xT=SS&jo()x$3>; zc^&I#tw!Av@z4M2#Phv3(=fLIS61GAzDrjDUG^-2!+h@|y?B1m_w|nbu8i#u#`fp1 z&-nLn^pnm6T$*;mx7f1q&m4>^a})ikD{X*&oB1EmWfyY3DJpR3Xf`hh`oZ17nb~hY zaBA5b;K-_*Rp7EKpwHIvw-vbT3FxvpZ5)3ew$d>bxMVk@FfQW)AAv{H9{Q=kmH+Y^ z#-ZD>JnFmGV4qxaeSud}*4zP(&1rQv@MGE+dqQW;`SjqZw}viTcE&k{F5r1Vp)01t_H;l{teCV>J16f(w51Vmb>1v8y0hz6$?s`}36F?je&aC`AY-ylE zpFMs79zZQo$AN1MEUtEFDnJ3g< z9O8uqCEh}FJyQ~|2X#hh$969X@wWqG``IO~yuM3Hd^G*@8%pB!fIri3co4i>eoTo= zm+*dxKO}p*#AToTi2l-vqYjDk_u}}_V@vk~J+|sArc1q${`dB!-k)4+pHf#|>w)0X zD$oX9HSMOepu3h`7Ta$DZ>FyMPVBSp9Ym=QR{j_2kBE1oo!n@0smt#B z8{=_ZTZhjX-`jxkXpWUaf319F@MGTZyO)LbL!qOluJx{WJBX`QX%S-ws@w{@2dn z#|l8j_T$m6I9YkTuJSNWzj6=yuGsH^sM~(DJm|hh#D4E!T;@Z4AKQO{ep~wRc>EgZ zyU9%x6)~O&&P;oBXVm{-ysc}67pguLbkR)TJGvsoD=)4H&*e4fM?K8V6)~<^5#o8z zS9q=R5Gz7F@;&fr+QHz{vR8irho;{$zalrVn z#X)z7y6LV-pG13hZ}4LIEznJKoj5ER;*h8_Qi}3DyPuMD>3+bcHBK*ZV#YTw0|#dw z))mPRAMBs>c731Nes3J-KXDvzX(gn*ob<``PhUv}9s3)}5N8|-9X9=5@aV(&{Er7u z*ZLfF;=hA8)6e=d8FT>9U59xk8RM4V&*BO=A316(__*?Zp~t4ZyAgFp7*_$0*7pIA z4tXKiQs0g`BSTrQ2)gSK-;C}1r(F6Mi~Wv0X@*`9R>%T4(M^g?xvq*!~E3HT9&>TdQu}Gbw+FIy-P^Vfo?Mr+(s-c%DyV`-FI& z?@}R;;rG}+Bjru1$A#->fk)IC(RHjyMg2rp>rmWJ0~6DAm_#7`#{&MdUvZYKpY2eFmqHU)H^yBTv zUC&ifhdm&Ub4wiO9&lvki9&x(zvanzT~9;5ObOTa$V5)y>H4*fL^Cg`k{AoC6?LNb6jD(ZH@ zpH)x^Tb9B2^%F~QJve?Wjw6uARy<^!C6p-gW&Hdr=%8qrfd% zRf_L{J4+u8eYGA8aOaSKi|b|p`3I}weONZh5-KQ-t;PX=R=rN}X7z*KTKO!{S4*dY zt@O>mS(WO3LT|0Y-Oy8Wy*+S~zB=3x;LSI2qQ_Qt{#2`$zHX|y?BS_V?+xB8zk6C7 z2i#fpmZn*v|D_ZYgx*>@7i`6Q7DnEj0+F92;AKaM!>n_ko^E?NFE7KA02X0I~!<}eXK6HZxJs-Dd?hgd{flt zH`y?*(*wYZsWO@9`37Fhc?LW=#7(eIUBvU?#*CAF4X#YyFdd&u7YW`>efdUw zuJOQ+8F%Z`V*R;3UyFJlia&r8%TEG7rY`SuaA(F}mVqm?4r(K|Tpx>CLqBk2@sXb3 z#nLN+AB!Kf#<*B($P1VOeysTwbwa4SZpU%TQveqhpV$zt7j-|xUj}Y4EtfB0tGeRY zGC^_LhH$;8`%$SlYg>rNRHDv@bY{>s>-&HwD=(tW21~z%E&ccL=r3J1__5YY=$|FD zYPTumLi>Ru%ijh^W}M_k@M7w!9%_%*6WjlU9-0D*`Orx#E(#8;xC1yb$3MLzeExWc z)hIs>yqNy{M6}EQ?TEhvbknReQ~>U){6}zL&M)VK6U(0NjCFOerT_mT+EuU^`x*yz zK(rohs_Rk1x(D6CfffJV81hg0Z;WvUaAeg*07usG&_hdSi2ho?pHS#`lpG=M-qVteFol4odxvG%72`R&zT-E z2mD#jT_SLw4GFJQJ@l-Qr{6WphtSd5J1cOYGqYT2MC8Iu6Sx@t6tBwikng7p zTveYLnjNn1Q0R{J z{NZyw$IvPBfx9LUc2jG6VW$>Y$f_ zAIpEscIAO=jN=r8ze>-Z;~{?rI%KXJhk)10Kg|ie=KP%Cmt6(!OMmWO@Lutl9G}GU zF>j#$26@hC)c4>z^LcFl5#v*ju^8N!e%lK4r!EBCS9v74;I*tLRSd36zO)PKKP%q} z+?RE(P@hBX;J+ciG&ks;Z^{kVc{}vV)WbZI8@W(!;4H6#CLw|p3! zm}T&u%!@n-JXdxKK3Dt#AG&{1y&eoOuLw&1*skL+FOg)~1{ zcZ2yHXBPTM@e+KlI^cyd{~n*~ItnBA!Fciy3WF~7GwhQef$xr^V>cTdm*zz+bj;F2 zgU3?03a%UK{(!eqP`wlCY%mV8e^KO5;IfJX7kL-y!BCe&*HPrkheEv#t_yEs`K;v1;5F_;`<|6;FXx+%DBsIvne z%QAzP7Y9B8o~t4I+Nmp8 z4K6Hu4#(+wz+Vv#D*@LXCVh5^hkAo0t~!E8K~KzjP~f?tP93{DtGFw^#t4t2hb9m3|eRm~oUM)VJU`;JT9Q>;bOJ zG|&B@KUO^*@L9%jF9e^JP6eM64}#|op%C%0(h!FkTpHse*rzUaEZXU3eo^Ym_p6TQ znHt9d&kc3`O5?hQdKS{VKo3lNb7$y{S^p3mSNFYR%e<+xpeH8B=!2~awwDE7`~bKv z*QY1STzQGFpk4=MOYcO^`vpFyo%l;U4>)efM}!WUfkD*MP@ZTBxT>zBJgg7kx5`I6 z0QEJL2U70JmpZRJ#9vS^I;^wk&wo|mv8;#iN_mWbLRYN3!*UOKRp`&S$sgr$y@ifh z@q3(C@eq8@^uj#oiK*w4tX=+wA7iSPV@nczd}zO>inQ?1`X(=l0i=nUaRw+3C^klUP*t9?(wC` zz>C0d#dog(r=@`M7I0wlBGk{QR(vNJ=9!0ZJ=CW@4Sg})ZSY&Qk3fA3&O6|?i>asm zC>ipsK95{?BKE0|L49c5Pg5h81-A`(kKnePA8V6Ar@lI#rvaRp{%Sk;G|PsA+bZ6G zdKsgr`>O!Ir5&|h$|ZNJOoh5Xd&GLmd!@p2bQm}%*M0C>;e}(tO9|4?NQL|Tf>>wy z;@A#uEB*X+;H-3$$!+EL#c}RU1>Nl6RIEz`eye#Yw!eWbb+e;W@qG!XM?=OD$n&ya9%$5=f-g^#OL&{E<=3_&i6OO^W2ni z3BgeZL-_~bv*S3wKb(p0U5Pvxd{*=LtFhnf;J=)IKFGN05RL&C<^wVw^)P5pqaKFh z63_!H-wJgwLfxjwdFSK$8JJ!g+v~t{l@AFH&HR$~SO?>O;j^Fp4?a7Nd{#oq`=GBC zn6SiVbezFukLAqrGU}X?d%bB1y|@cpR)TBr*dPRlzLN$Mc&uJX2p+2!@Qgv8Jsq{r zV*Lp4Sm`-Y-`Qv=jkmqTg{`;11k^@F01!`fy?SJ=zwJ> zq3(rVwDhAT^r9Q^SQR1#e^uZCJXUq9v88|dDmbj_?4!Q3#)tk_?budO82)?YWm7Dn z7f4`B8S(RfL|rcGI_rfyQ)4|1a99=K2Zt3m`_sa_4qdSF;Zg5`I;G0#vETHlg9d+9 z92Y#6e&$-#zffK-^uHPh99DX|nNim}%c{h`z+Z{6-vW;nKY-p>=K+V6t{Hk?`TgLr z;vd-3?ixSWg5f&Xs%RIT1`aEMJ?dSEcYwp{ePe%H;y7kVg2U=O^Q|}iqg`sEUK#wA z^?$0cFTWq0RrRSC#5&yA3e+r!`O(+k!m97LFR$K;KuE)ng*Q>l? z@K+^Lqt1mUB-C+M9c0w65Ko4_m-goy*q2Tfd{#PDa8}lFIt%<&6Y3JHl%It<&gus~ zs~1v%x2oUL`1f5J`3|@&W!-0hztSDP0KAoDS4M%u(j590999M;_^R?#p!XFJs*SuD z{8fdpP|reoS?GG{FAQ29+m~B!#V5dJmA|wi)*T0bRX!5-Rc{%5R&`2HpIPI8yQ(f1 zwi+K>)iH&>m%7q}P^X#e=|kA3-8dFo)!q1S`|7AC1%Fju{OZVKz+ZKqHIef|->ddD zQ9rpR)~!KZ3+W7?>y^L@+*Nf+P{)Gi-Wt@epue?QomJC+J`#LY90^=ioTM)5a=~Zm zcf5ywdan^UtqcWx&g(C%54>pCdJA!L?5my&>M^t5p!!fpeFXHrCy^&@()FreJ^IJ` z%xXUxb(fi6dP1zfd_Fj<<^k|nhK)uv#CQUDE$iW|*re-~uCP0Fx;h^Gm2uMkjj=up zxT}tDj64~i^Y@vBeT@U1FZKVZ#~kc=@Koj%+|p!0KM4Lx-taCst$s)7c=_jBgnrc3 zZU>Gkz_BjsUa_xr2s&NrE$#t_)xsXz$6#OeJJ!Ygh1f3(yj2CnnnOGXb(p)W?hbe> z^}-LL?u2-9)caO%((wk448BSbn~M5ea98O@jt6Je^`I_=jz>FXc^{y^{C99x<~P)$ zJ_UJk>w2r>W3+8c)ZK!&N_UHT%gTcQXBF^`I^QwS^(t=?^(QzVfV+ly0^F9oi27Xu zv_qiVCFZ@VHPrXNy*21+-`=F#tx{n%)T7Wi_?+Pv)LjmBOg6;#XM(GWXF;#4xD@og zJkQ(cC%rB9$%~=il`bR?{FGtIZLy8|-TvfceZf~%r?D-}i*G>B%f|zCmQ{ylRhvx_ z0E6C^+@%}%sp{#2r!s$}H|k0-9(5PC%Dczsir+z}t8vg@@4bQUSM`OU&n0i%8oX6{ zThyDNE(^Nd?$pKI2JWhHHtBX%r)I_`{u(%R)Z^}pEyJ8ggR8zxlyD~AbnDae_3gAP{-q2Q_Vo4{o`U)=}}D@O|ZT8~g?LNDk7k5#@o_L)}% zk1P1=;GojSqCN!YvjY>s4?Z>OdBI=hcV7tIF8%U;(DhP}i1!1^&VMrMd51*3?}yO! z>b+R_oYTe+_*@03pwDH3-U9Gf`mxJV&w_Sn7TycUJkTy#!H(NLE9iLl0cWMTgF4C@ z=XmVPUq(9%I$xd@IMt1)hs;EP67cg+5pQGw+~$|(~D+8|4Tb#Vbtw5;e5Id=xCMCnG?A4cF^52p3x0DT(&!S zD&qvuHq43QLk8W9Pm+I?@qbFpWYw!9i&eIkL7#53_4u;E$4yXlCxiieV*@D?2DI! zyDF{#J+5>Vxn8LFG3rx@Q%8RLGx~9vT^Y~Q7|#P;ZnzHUY6aq;i)G?r_q7q|^A7a366ipOt9s_(y3D`+4qJw+X2gEW z;(1U%BJgU|g`j~y^1o&*w1?NYdJ9Me1AyxVt(M>f5qb)(2wj5fSW2_4-Tt%ebmvO z1br>_7nc+Sf9+~;Sj9=fX&H{XA39sT&l2^O6$b-{Wg6L7^k-Vcm-w9i$aL^uEi@RP zyrKd6UtI_IF4xN}a8m+>&7r5|IMCyYBOFl}uct8ZF}(jzyrCaH=QyYvAzksqg+Vv` zEVwJz`61X+w=ov{R(2uAV_AK?|4-lNZyZnFjJgq_t`l^)tiRs@j;jK+MS+X$4jrxP z;enr0hk0sI;A$5ag}m|`p_i3Cjk*!~La28f0+3O!h`JEWKl}ommvP?eqQJ$c#eOwK zF8{n9oS1Q<_1H2!vk3JbXpV1O2R(0!9&zK9xqtNX#-|&qR5AQbt zuO;`w`~GH8&w%<6A?}3jcQ80=QFobk%bM7KJ;qZ$7kH}d71VQJ9?<^asT{|^Wmz8+ z@BLHzMbOJK-vRIaQ$M`-Pv7fa@LlOCz;(%QUc|oY>0)2&AL=Qq?g_>fPsi~LPcHx; z7QcqRmpX;Zi&`yCFBOF!_~vY3Am`+?hre2Ca@5b7z@jeiDw z7AlhaUPb)~^7r9oVZECGj!S#+M~owV4f-)oinE^?iA)AF&o(Sb8VuWZ7N{T`hSz-uD;k=AjNU+mFS*bY{@a^0g4XV!unkYk8hq z@Hyk<(V!y92F8%Bul0Jxe0h4fjeEz3lpZ4Y)@MqeI3!&3xe18det^Cbc z2LtsWLcJB}RoR{g-pRNe>OzF~8>f5>?SierVF~cSbM;=WJ;6z-A3!~b-ozq3z-8$_ z9}aHHf1V2S%qcjI_Q2_}&MhQyYl!{8b@dvm524dl-3+ud&I29paM}|SQ%)7(WYj+vub&1zF4yhp z;M2;F1HaXRjO(M`Xbm_v^LOiU9QFU5DG%=zg3eZX)8Jbx>@E$TnU{lFC5(WL@vCygR?3gkO^@O)Pd0Z1|G@8cmS?TdI!{9mfisRUSgyV zWBVuIwOj{2&xGG^BEFB}n0PNCamcUO*E#_1OMy0YxZ?ESy{gBw7WFIiBF0$%Vq@UF zz5bUv5a0g~-rJRRAOw!?MjZ$yY&>WQ1sa~QgkDheyd}CY(GNP@8C)<&T0)>1I$YH= z0QWtVg%>}xg!t*lmXIzI+*j`lgZ@@8_WaZmdciO1KZu`VE1fU6uiBx*)eDTFzhxY> z0bA8S0QXfL*l#VN!ZP5#Ixlp%T4?cJKfQ1rTgHF;fdA^o2mk#yCwlN-)!Y5as-$}W z2NqZT#S+vzcftGo^g?fJH4Zqi+EMSAx~&}Oam5*-$5p*)?DO?Uw?c<|J_|mA_li?u zs{$V2z-kBYRozu=bzRWoDz6V)>Ak00wF-EE|Ell|c(3v#!GVts^^U8hO9k(}iwhIB zI?pT%dH`&7Jh-oT=IqEb!Gom(onuvcj~#TmjGMg%?yG*NcdUGT@LzF9@L$<4;J+$x z1^z3?1l(7K9Jc)P{{bE>kc~RWvKLYJff#(pn)rR7zZD0^R^JmkT|Va}3#?l8IiSat zeii%DFQD!N?aevR>nh)5VSMi$xG&eGyRgsvw%OR{`#Cp5-3Q{L3$T@84E`&DAoRGZ z@3AQIUT|PNm(bxdAL?o7aP>lHaAD~=QP)^`qTs&z9Z>gyvdmw>hjl)*Q&+zhoLJ|B z9#`Xm1FPTCsGG&-43DJ2g*6Viuyo_7Yb+f%xG&e`uGmU94*sil?8}dXK3B(s`^sO& zd;1t4eGoiY-wS;BGA?-doN=#RHH&#l0BrSC(1V)?P)!g`)i_d)XyKIe^}A>%r=NiS(hFkC0e1%f zRsC7$ZuR#9_vLziNn_-*_*@*lG1TYy0sYBc=QUd3Oz6+)cT3cBVBCKfa9`E&fDV^- zF6F%4VfD7|{JA)UCn}83~PwERUEL}izm>=JVK39Gx__BBr zwi*Zhbv^4teOmD1phH-1J(ULtKCE%D)qRc6_1?c&A7bXFdJny&`vVVF9USmr8U{VV zjdeZX##&dfPrvAOaAEOhY#HuqLfvCJCY{jp>URMrrhU=_9GKkVpQr;t`|*0{a<%Sk zs`n6hFZ$E%T@Agi{4a1~zCU6c)Nj!IiG5uU>J;-ax(@Xmc)hImzH=ip}r6}GQam0(C2c!zj&~xC&0u$}=`QKFcLH#=ZTR&n`-G{+i&!N-RICw9g>MTIVOa0wP;LEBv z8{0c#dnLFqmz_PqgOz}Y`ozp{JEJq!k3jtgjnf(G&Af-tRWB*lg;*K;t?!I+BJXvh)%v>)_%0nMaO4`j z*9iOMN@svGYu*C)y%(Y%b?Y}_%ir@p)OlcikB3pu zfjY65QO7~^9NLwK6Z?G&4t*FM_{gIdqrdc~`@ zwHD`-y#?;e`=lar<4W*f9)B?Q^>+eirmp%V)D@OKE+^{Ip!ZdMZnTU0;B(El;Kqud zLoZB+?IY-ZnU6CO@7bgIH$BId$1^7<>Pv7Q0byJR^IP(u-&MRGI$iml(A_GZ2RvEV zfqKNMQ;PZw)Dc{h8}uu;K^H8Yd^~;_biw3UBd}%NllRbHoEy4d#zm&#bJ-uU9b7ro z4Z!&Hv)VB(U!%Gi^u6TumC*C*eG|~@l8>K=dd0G9P@h40U8uvrdRjN4KlN>cpbw@F zb1;rm+y?!b|N0)b9Dh9ebM^lf{8@1}j3fInwzq>jORoVwEbb3oFWYwo7bYOy6ZIZA z&f(z0{N0X29b#VZIcTSyavkbB5EI@CU9ajOVI0-fK|7bXq2SbfUPogbhK(lTc;z|8 zaq!+f9p8+t^7O&2#luj?LArwiPq2RJ9;mZGJydsm&i13gff>I$8Fhz=?JmQX%iUGr z#=0KVf1q9XAo?@U0QHDvm%Ii}O&uVoUxowJ=;abilGJ*Wq|xcPq9ESQL4vyHw<|N45tCRz4B~0?;Ka=P zI13z^b`a_i>;3G#OT0fh|1I%-d3S;v^F1QZLf1=h&w9kP_dksN#^Q6u0nlI9f%ot! zp9))rkMTLbSU$L~@=UO${f&CWnio;$LH!N|M^@ezxG~Lr)OiT+|A_tY-o5aC9>6iPkVXQ=Z)eb67J9u`ELj}yY0*slru@M`W$K;O%LDR5%dF#r!%9Ruin zV||C9yD4+!g`HRCvhz^qA=H}yN9H(Hs87uJKXk&X&j5XJ_FU9UbzRu1fNXgf=Tz`!%5TryRM$auiw59x z<+p(w%kGWGzlqP4XNKdc@0$o6Fm)$*-=5Al6a95Ps0X2V3AXZo!G#qMsBl-lBkC^b zeJHX1!r|bQR?R)#C_vZQyg3xij&&DyML*VofKIrM@~wkWA3=Hv@L<+&I6AhU zf_<$cvHg7LgZcOF1MW<{LBCih;(GL_-ww{KL6n9mM$b8^_{&wCBdhIu76C`)EA^_vQlxy>M8U@cul;h30Rn@1TCr z3u`@TMt{n*vQck=e{XQ+|JHZV@4N+cwbW&Am5T3eOoiulALwhjfTO;H>fIj%ZcCZQ zG2q530ED{3T%UTPz5{gwy;32554~_$S1=y!%qvnMUUW11F+UM}Svq{wdFaJDp^v9r z?U^Y1XlTo>KU2 z#c?v>_x=ZTxzug!4Q|Z3LHpxSB zyy|)5c-pDIfKQW~Kvzuj?Y+P8{it6=o#AoHZ;kEH6)Uf>Bk<+w{Qq~ohpE8HS@_R? zMfg7${S#Fbgq?@?$gD>K#*KJ?9bZ^K7Rlo2fa3~sb;LLi_@)%1i3K3D6Zwr}Dd1$sYwzR{kw?#mW!DzVsE)7b^||&MX}>`0^VxG{Bd+ zp6vanC8SHiR&`&e#eUen!vGMtGC5a&=#145Ta5$0tn0_}AQ6_X2!b^$F%hJvFv^{!srx zbuhq}bz{zt*E`=5dT}T8#tM|yM4kUW0TSlk#qS^dD16<^1`>K=K{lJ-39~=9skASVlgU(p*X9QQ49gD5FGjzzhZqyOheFwg*I0yE{ zpYi@YeLw7L9Mlokd9ao48#-jw+rU=!u%I{A^FqsS>xbyO!v+CsPCXaB)GEv-r&t@$5!WA73)WUBkMe<^Ps;k_SFyVDu9oDz0d|c zc@F)I)mE*1hfVb#Lct>FjHO#g-3Q{vr?4--0$bIIz*cdgwefevRzg4MkhK6}pY~if z>Ith(88|cb4rgOuzbp1tf3z;%56~k^PmOl@Rj5B4Idc``GSD9n(*gp%Ouz1O=#4c$ zV=Mm=`egb~L$KvM^IrY`R}Vt|&GOj48C&V+!I7DVy)U@3#%qXi1n_3YrG_+we#228 zLiy>aLo8h+K9_#}|Fm?*HVnc*6s5~h6jqvuDWFb9Ad94D0D2aPlufAPh=L~XKe