From 8d97178881f103a81c0b4f6361cd9769883aa379 Mon Sep 17 00:00:00 2001 From: Ian Jonsen Date: Thu, 11 Jan 2024 12:13:21 -0400 Subject: [PATCH] track_simulation speedup --- vignettes/Track_simulation.Rmd | 5 ++--- vignettes/data/track_simulation/simdat.rda | Bin 0 -> 33100 bytes 2 files changed, 2 insertions(+), 3 deletions(-) create mode 100644 vignettes/data/track_simulation/simdat.rda diff --git a/vignettes/Track_simulation.Rmd b/vignettes/Track_simulation.Rmd index 2ca40616..029f25c4 100644 --- a/vignettes/Track_simulation.Rmd +++ b/vignettes/Track_simulation.Rmd @@ -25,7 +25,7 @@ The `sim` function is useful for situations where simulated random tracks are re #### Simulate tracks from each of the 3 process models -```{r sim} +```{r sim, eval=FALSE} set.seed(pi) sim.rw <- sim(N = 200, model = "rw") sim.crw <- sim(N = 200, model = "crw", D = 0.5) @@ -44,8 +44,7 @@ sim.mp <- sim(N = 200, model = "mp", sigma_g = 0.3) ![](images/track_simulation/plot_sim1.jpg){width="95%"} #### Fit SSM to simulated lon,lat locations (with Argos error) - -```{r SSM, warning=FALSE, message=FALSE, fig.width=7, fig.height=6} +```{r SSM, eval=FALSE, warning=FALSE, message=FALSE} # coerce simulated RW data to format expected by fit_ssm d <- with(sim.rw, data.frame(id = 1, date, lc, lon, lat)) diff --git a/vignettes/data/track_simulation/simdat.rda b/vignettes/data/track_simulation/simdat.rda new file mode 100644 index 0000000000000000000000000000000000000000..d634b70036e5ff4d36e4ed109d4c313d598483dc GIT binary patch literal 33100 zcmV(pK=8l)H+ooF0004LBHlIv03iV!0000G&sfahufu^4T>vQ&2UKVgRpfklJ zX`5N??lWQ-aN$^%=S23Cid~p$T44^fsRJS2Y<*>SmEy7;O00LA?;3$tFVcIw6&U7M zA}hMt%n?z&a1i1=#cRp>DQD3lNC8k(S{Ol%B*dkpn;9%}oc{Fj=n|9rJzX|o+0`A1 zCD_~7!B&iJ6J}LMEmqmWsf#q}w0GPi43n-4i(HoQ zZ&H!x_@yZSm<)CN>UQ#}Dq^IQw&KM;S=1L^Zy}5Mw@|+Llpej5^Lc=)8HSvhJHg0U zwv*7MQ(?^3EqwksMMOsaqRGsU33|jU=4k4`i*y%OD*fZ-ftHi+wJpzEdPzK09mUg{ zP1ly>)isr>qBfIb#E|OiFRMqws=eb-H9XEWA{d>4O+rRmC{bN~jCo!m@20)YKvqtV z8xAa)sxvM6?L+_J7mTuCw6$#X4~ zw10Gh%mmuiAcsXU%vY@@mtH(v#{l^`Xng7w{A+@ZfC>0hS>w<@E0jNL9x6-r7&1N{ z#5ZT9!O!1&I=7B7*!H~L7I=_MI)q$6>JKqDfN&)(Zo>6P)&k3(1L0XLCM-+*a#}Bf z;HWG^iSNxpf|+q3du=0nE{jY2(9?<%sT`-ueeNB`AmBHYNw36*NY%lGedzo zv8EF>o!YxlaQxgm$hCQdQc@TjTkL153=3H*NN$rm6mBj#FZ!tJ*OLkk3vF zfFbEJmb0V$69+8VKUPz$^dS5zs3m+Mz{R#yUZnTnM5`^|_OB`jc#qd-aQ76_t&ZND zY97BmPkMQdg3?w*(n;2&N{>=pF0QrM0K8dJTQ8s$>34pqUuh(VzDN87L_;`MrG!rG z4V~W?*3|s`Y{zgz5Vj&ovq~C6>gK!kgls~XYZX3MC`p+wo5gYa>X0*bBi_2c?01k) zJf*$18&CQ2_vh^{V?yY=BoDtdn_^f;(hPP#!%YyEEfxy6=G#2;C_1qwn+UaEH`U*U z6Yb7TTI!d_*%(f8X3_VQ{C_}(3hifIcHdoy^gDpb%!uV7&kwDWd(C4n3DvecKl)*t z;y&JF3RUo#-?%BuO!USf)H&N0LmBiG%Gj8#;`)UwpSKzrtdZzh^w*2 z`jU~;V)*49{y~lcN_{p;Q7Ck_K{{nHDg#3X>BA#y*;+vE+$0+jmoLKn+d!r$3at%6 zA9NfoBbQ2TN)`Wqp>GmP5s&y{HsG5QG9|SyWYj^9-gM&+CVF>3CE zqv|vsXZye^_+z_;(%oLAZJs__8ZT8;pltO?6MeQD9;KDiH5IqMvmlwi4D~kxTY1{j zdH;F_MZ(ep8@8oP*ur*1>v(TNrMA5KF$E*DG>pUGMUwpL3>7*0a%?Sne(N1|lC z&t&5EF!Ujz>k})ukwz3IlC@F=(HlFM%T_lGRTUc}zcHWxL`((B1jD;cS0#GYox>Nk z&+2s0`=lWrBk;ITLmsIB(tCD*#K>c_wSdKnRwxeOS+-}AYp=DnJQ(=7lcHL+aMm4E zQZMw1f?fY42?#Suf$zh%p^dDCSQH8=aE1&s{>w3u5y`>ezscYv4|MFbRu|`i_ipM( z%+P7p{n)O%9y>nYTh(gWq$txH?uyhqjDG3XI2&x8cXexPPGH=d+T*GAAnNn;@I_L3yDie zea`bFJMAuUF_%cO%Zq8>;}YK&uwXN<5x7KetBcNvf0zj=kuc_~Tgp>#u~}t$R(h@{ zw}x+z0SDA>DOmKdHF1(wK(U0MO2`|IdPlP5JgpQ}nlsuFMVv7(_n!#PW1em+?lumP zEzrZcI^j&o;M@RpPi>B7s|cN4kwai`a6C1uX+>1;)jCwsN)YZtWS7)BRSncsdm>0i zNQrAb8eZCC9Jm2Rw_$Rn54O18KaDR|PfzEg!{H8KuRcFNlK-n6LvZ7lp*nOYqX1%Z zNcW6SS!?Q}@|$hp4z<9ZXcFEAiGH z-gmhbCf|byndZJzJnipI@>vgyj)U#;%o{4MrtQx&3&-c%a05zZ=c>jaczZ^u#z4DX z^HGbj>TaA-Q%|T7s<^ zv@Pqkcf9g^5iP<1r>5{n| zi2CM!(}B!jyOne^idsarTYHqz9sP-Y-dML$m5ipYFX1KB0_u?1;@40Ksp0OFTE6Ks zsAi`-9<(u@0kSjOZh3u}c6gqSg8oG@VX(fWAUm`Pmpa_#uGxGXb zS&mOb)8a<6UuKDnHWgZfX3Gez@GYJ90{TI=0~>{9CxvJOqr}=U)fC`U@3LA*3i!bM z&Y*S9PEASoZZVEE*`(AR-h+ zPA{n5!3JgCRO5A@At%yb%11?mf_xTy;jAji&W3aL;N20^#0QiPwIg@liq4)_C3L+e zXS-`fPnvADsQ^emVk^olj{~eE8q7%2hhcQgUa{tN0B-~g?CQU>`vb^){o7JaO0zZNi0#H)De6>-9M5a~ zS$Qp#-V-D1FW_VS^}zs};QN*5@hRdXe+5xfgl%DHI!nLWE}A+B;=6)Afp6X?uOx2c zaVrwNEXgmniDiq~&t|f0m?vo0xmKj)A6_co4GsG^T7tXh3ct?v=*JSzg4HYq*OkHe z54;YVKa$7)bLonl859~xk$$|{-xrD&Iz60ey2H?H-UOaCNmGvEkbjCi6p}^@kHRL% zz3RbEwr?Z=-1TSVYu`|T+R!|dGsoMVzbEk;lMlh%iQ2H#df;$qa*QNrg;+ZRId{^? zuQB_=vt;yC$r2UjhqR!eR;$>+#Dpe)-J1JxSFy*%p)M2H*tAjO>HL)jNdS%YVgF%CJh9YnTXo)U!g_n zZR0$}M)Y*9xW==E)gl>jhLM%zt-j?@BzYIPA0>(zctW(7qX!;;!Zz`mb#G3sp_q}V zd4w>Jow~U45|E<8N@>`sqR8jHm83)(T_*&&ZrL{v+T50iX)eM}wj7o?yn%%cka_5K z_&=}e1}+|K2DV?fRnqm&OB5gLraXa*<0Lp6922E$u_lJ>gyR76~$@8 zH_Kh(4po}tD|rdd|F$oD9n9#d8+eI%!&w5EJJhS#FlJNIsMhLdvnZ%T!@BqA_t=Zk zJEiPrYzJY|T74D%C&lU`fTdy^7k0-TrpY!o&~JD-!yxw{*kbqOt3NBKVnGZkP9g$0 zwv~T9>1Y92xXU3I(BR}jlDzfHLndatx&c}C-Xx*s8%`9ZuwZb|PmSe5Ki zqj-!^z10Yk*wEu(+wX?f=0D`AeI5%#F@j%~*C{e)PlxSnL3$-Fb8b<1YdlTx9WMSk z*tdKnQOgObY}*+VR$5u`)-nPj#&Y`>$v~II33lTJ`z@-!sQRZk<(83T;zPW}u7!OU zn-qa`viYVW_e%zl(5kbzpFw|?5*fCRo6IwjnhW@IAc0!$$GjK{p>jk#3 zAxw@u2*j{~Ek^DA8FblAyB+1r+(MCkVrX!5hUVKPU$v$`r*3q1%GNzncIiZ@aUA?^lTKdPBi5_0h;!xz@IcG};IPfN@qPC$+{tCZQTuhIxM zqr)pQ@0WC@2&v(cSjN@olK?Oyg1OxZqD&Jj8|CrAfoGNHsPp4}0mKjF4U1u|dQ%H0 zEx^1mjFr~2Z?nJIVuTeBO5G*?IwZ9USZ3Ayz#8eLz|v&%UHd)&KDh^!4L~DQ;ZJnb zvUCPfw`FYl?_H<=(<^l4j&OTK;Xy~O8YUYX%Quh#ZM(gUrT{HhcD2$T_&0Df=_q}F zJN71<=1jMY(#7xJU)rYxUwiq_2Ft`ryw#ONcE`0EqNnH!ShjN7lJhX{1MHIcm#3qw z?7_${q1!=YffrzM4{I+mU!;%R-o-FWb<~XyN5oXLc6mZ0{4zCjSRL-2tO~+A3%B=$ z^u(#o4mr;!Vbvaocd{_Co3?C_!aor5g~O$41;*9mt89z}tRY3XaA2r%gy`RDYI7bs z6J!UY2|%>WS;84JL^&lgF7RpHVNnAS6kcg(S5FLFUaK`(DL|kd=CArwg~@>CB2aV( zgdrKy`5hA{^3XO*1i9h+XTU39plR^NG&J|Zbi|W!zUAbPEE_#(>$?Iu&`>AGg{E2(D6cg5@cNGx`y7xzw zteB+}3H9%37`t^N2jd#eSmBiLnPDV$gfY)0Nu2vx5a|37)vwPr-q*;15_yld?41XW z^czqM2yJDQeh%*`GDs%yy1Te%CHd}fhPv!#X33u;$3(qREKHAn_@KyImFa) z@g7>uIMxqTa)l%L%3m8QF&RkyFx zTSrCC)YNkbc7y2uBecG~XwrAn+R~h6+wWRFxUEf1*XE#vKj}m2G$ER3B+y~llr-}U za{jhJ%)SC0Roc|95vE^>) zdqIRSK8O!q&Y(b@ISBRWmND~Aj8<#WASO6kTa47&LE^JXamfz_K7=bD3Ky1sVxXA!4v3lRJt2IhfVo+mz%vQF=M8`p=nCmV!F9g~V3Le`@C z{1XBW^{MZy2Q9~mSS<{Cnqff9}-i zwvRVdDL_5u^?CIS_+bL3G%-#cofbd$ov=HaoTc=E?1>`JOO;p9+^KC-85a#JeK~AC$?SNq@O{re@btr zUtwjIUJTwHxptzus{h~P;2RsEg;LUERK$-tL(b}9Buq%F>i!6xU$!ba{@*@cW&G6> zemmO&P$#ECLntA9!(ysTRQKKip=lT7&=636zJwQm;5twBP}Y;H7c`N_03o{-{sJNr zsR^H!Dr`u0%s-c(^Oc_d$6qQ$t?_zJ4sgX4cLOMvRevUXwdSomJKpXi@lkV}^f7N} zkGw|*B`r3fQh~1X5!hb>8iD9fIg@b`8tr7mk$uy~^f|#jf+}&zP$+H*2yzZsne1Y) z1{8JeOk1fV3)M0Cm^@+Ty$oR@JOMKA!ur~uwjPDDqsT=m1+-^ESN@AH^TR)ACaMy{Vy>vGqqgp{U z3(O5nW1F}+nqVu`DLDPhVYhql%I&yD_>CnP{YY~e_x=wcE3oG{B2+U4{f9^%J9}4A zU_D&xR?Ze|l)8RZFeJ22XrpOmvSqqVO8N@=9_Z9a1E=^ZT;_X=Nss=rFvQ6%SH@w8tjTZ;-H5v)e0Q zPHi(dS+3fn6SI8bvWE8rOO@A|rQYYv)=d#XEzs%c(>D7t=Iq7ud27bAAo^MPB{B_= z0RZMllLXQb?hNSRkF0_ZW`>t$wEPa&u?dAeoFH-APund!fdk~+7-f(j_!0=ozlT`Mw#*(-nPs`q&EDTP@v5 zZ5kcgPDO~#!d^;BLy=n>dNPjlQ5#LnNiR4;Zp;ygKJ^|SBHNCwbIHJc%-RWWKhN5fEHp7?JT2J%`!^< zfz#KoT=^-T+Wi6e%x?A}AUgyXe|`m4?>NP~E!UAt1?5>h2E4vb2+ET-+|rbQm= z=1Bgbk|gkVTwqQx;s;3(Hs?(t4!`R3k#bpY>?cR{2<_W{zNPf0dQVDoKK^Y`Mbvg>Q}sZ4u=Z9TVC0s9MgUA;}CTW9UiKD9aaw>tUtcBMZ41e0u! zp=jmuH37I(05LB2?s|xcrV5N;{-u@6gZ7Z>S5l`vPNfW0$9=BIY16ah({ZS1Hm-8k$azk!Dyk5`% zpL8042RR9kqRpIF_g>E90ZJHF;Rt)g8~Wkwj+9JnMj)Bf4jzVb_q@!lEp zK~P*cJaJ~pf+RT;CjoNtp>v(E_~fx%AcomN)Tbf#OUV)@lb&rPHmp#_l!vHL&zA7szxj#`5s)k*;NbyU zy`y|bmejg#`ud$uW37}oe+E&yHMY-xqfzps0da$45I;El{((Yg%GAe`T!Z91lx37V zPntv4=8D)Y@5zF2@#$oJJqZEhl+vZF5cb79Z*@}pP+5jn#6_?C1Pd1J**Ba$@*d(D zr5R~3M19SE>*nx~uF)3F^UP_#lM{;t;cR+Xr>Z#URTqzm)1H&MF=gbKHGR3p(^Jq) zv7iVT`U+crqACGx$!V?XI5@0F(Id?iBGl#RGt9F{E(w0vFkjU)cW*rb#(^mJ1%a{y zFw`oUkpq$2mJwPU(JFf#r7qTwPrc$?)y3*gcO2XQje0(u_Eu@>Jl4rbHJzUU*WCJ4 zCg1zb4}V)sko7NYR1CBcpKw{2SGeB>_Z!v@yp`-@a{b!QMmMBA#q%fJe__ESK$`l1 zhsnC+y+;3+t|MmMWuQ;0INX0S0dPTjfWQ+9-}o?OfUj^9EUe$+Aw!^#bG=@ZdzRzI zQ0L!;pJ;J-G0n@ay84$(u-{{Cs}3~HNtN_(S&(Pc91mlzwLeg0$dZ`lA_1UOV;y^y zY1^>_;&N0yox;2s)B^hU-33-4?3Zy7x15!I~^UQJ+xrgW^ zhwNaj(*OrAW6=L?j_o&cG5L-fwjr89Q|Xe$s$tUhE1|3t=0f5$&=}Fd=w9~q^#`K} zb?uS>c-xfY6s^_l|=4176PX?ni_ z2W+hrO(kOIy>9pxVq*2%=O^eHWydb!S0Y3z5M{P{#NP3)w-0=+yqIMY7MltCXF!H0 zs|7C=4@N$$KR(LhIa4sw4Ee-|b3L>N9xzz0A#SPgn1CDyQ+;RuSeS)BaApQC2|NdJ zWnZpD34CR=r_6)2nR9ve#cldCI#XI#kaBY35z+F--Ys}-AGpTSTf2wxz!tsJv zF^vSZq&SP*_uEwGwX0C*!Ig{2d~jJ8Jpm{&#k|al4F29T^edBJT|!psdikiH7Aq(@ zLrVO(N?P^w0&JCK+I3Isz$d{XcOHJZdpr(4GbbHy`K(Eq7 zV6mOrqc@+HeISiCk#Z~oDlL8&^Fg)==39nkTdp0 zIr0b1k(50()8}yNr^A{tIxQu@<77U3%D7ViP0|Hhg{xDAuLmJ4O)XU%Bg?I~V4H){ z&CsSf4+{DO(#iZ6q)7%od1RudV%QJHPV=Zt7I2*@Dw^}aBGJ-!ZQ=>NPfg1Mu23%{ z`{lb~DqUqo1wPc=9|X{6QN?=&&vmm63a4Rr+m(;29{Ux&nnEkdUq#ze)<~94S&4&} zDuvRjLbMK=zqoOJ`%V-*1mG%3rGPTiZ5wp>C1X;3AvcAe806vkJb$% zVy3ShFa|%YIkl@UcODJj0NtMoL_lralmf&074!A!=)7N>*9t!~Nnm!B@@PU|8azBS zwjf`wMlHp|$Ap^^3(h2bflc5FTW9^#;LM#0v%>-dbE zhC7h(c&WH@uao?T@fJt~q}hCT;!GNwxxd3xh02`Be$Zm1J_??oS+j!fzQ!Q&xjW_} zVELo!^OyLa{O_Sp9jD>=fXp1oEWmA?%Vu>lFi%(>`*_*QeE3@(aB~Fxz{5Ik5I`qI zIs~+sySIHf?ZWByFST1cd(7~(ffDw^bahXh`!^xYe z%SFYp!V7s{%4^{6%q|KG+?(6Mv$_Yq)bzLiY^@pPTE@%ps92(6yfb7bpelu}Yz>mI zARr?7{R~MNwkZQmu8}9)AVe^c6=Hr<7fp%x5Z7yA{Y9#TR2clKFTx*Z)WJVY_+76` zf_|kydp~y{wy>F3=wpfy19|{fVhY+qjGNhWEpH&;iX>VVZ@{cC10O*UcCoN0?h~rq zrHx#jG8=l350^7k;$XXL*(68<8&5Xi1h6W}JwxwOz;dfktv=*7LAcfi4%AYuLO!6X zTpIA?^D$QXVMcZGFxBk*S-D6?X1f8hU%r3T)TXgi86LUpg_VPtq=j;MceR(T^_Lta zVo#mcVn$+()UIn4wzcq*Z{oQwwM<>6SI{15cZkDB8;jNt1n)`qj}~%|sFG%!rXX15 zv|^QOFkd$|#J&Y>Xn7d?*i$#C(azl%{#s@JLX~ANe5uCAmpf zM!@TzTqf66;$6^;@R?+|{x}R$vmrF|@EZ}^HCkWT|Td3AFAe~uK$qT(&j*t-D` z@;zQ=*FyV$h8*STSlz{)F=i!;#LW_XWYt+9*OEWJW+H$*~2JX7W~*>@1dj zO`YwOpI83-MtjWL2(Suf(1W=Xnc>Aty7-uM3va-^?W-qX!G%Mkk*02b&I)8-RI;vV z@43!+hgzNNDfoAQb=XMKNK>o;G!t`39X4jMa^k)_pSw-n1XdnfQdC&j4;l(wfMr%7 zGxQxy{hZ$UfZGDR#e)<>7JQc_a~~1j2O%qWz#X>RsQ;js&-If^fr3u8#y;LkIx$mN zPj$rJnslhHk8Dt_dluCM(xvyq)rA|oC>J@W)c{&?C8uD2jXHkTa%k#4y(sx#@VZP59CMrTJzEBf2WDtv z-ON($PH=-@GSAA(q3nncLlHxM`P=c4;Flb7TcZvBLw7-5HNy~=pb#z806#G*ZuPWF z%s6I2C-y%BYM#xRxRWB`c65npCM*`j@q^ zMCV{o#+yR8#T>ga#sZ<=;?lpH`2XC>(Bbo)Bg!n~yljP$`%r67-LvFPzm?n}s-bB@{> zG&&NO)SglzWM=_b=$GgVK0>IO!V~Do zx2on;Zdw9|bNf8&otPf_|E_0=x}(C1I^x+6x0tPhn!0~>Y$%K|uOe`n#o;&qV~oF{ zt6#H*dzswjAf7ELQ6kkH2xh5tpmxRjaOEEX!GEWOk4D@Ju zEhm&)!XpXqWRn-i+GQin#2HvP=>ksh;qHXfR~>`X{Mgp6A@;WTk*gqlfQ;JHDqC03 z;uHJgVkJhB_YXWJwL9;T*(PS`+01PXG^A?OoIf~qQI?MthZf!f1Hk>N5F6U~TY4>5 z2ltUY+mX-`DeCaq;aVjn23>XOI7DP9HsD*3;i$Ht5Jk;Po&8^)sWO?1!!yb9sPmC* zo1aB1{sM~$6L)8Y5HF*T%Rambl?KPU)aG2wxL@I-M`-H8R2&eFCv6tRUQ0`koK^j;_Xy`V z1T1Y7h%b_RBtY?)-5wg04^b!Az~Mis2thz9{V9$OqGctxT|0fqXM@Aryej(u4wW3A z4Auj6%W>!aE_eg5Zg-Y?d>p|+L+wiwRbUkYF<{o4^p?qL>92GEOVyxA2dkwJaqo#l zvGg!ysz&#RO(&_!6!ipQ>%qn^XrbeO*d##vW>TGtRgZ3P!ipfsniuzjDTR*sqU1SJ zu?h+!bC%ncR{=4T*#?cKfkM!- z_~w0BRvWJUH}M6DN#rkgCi@%(-y^3G3W&|G@ODEZl>dDE*)ijh8VdQok3?1rZV+cr z?j!N2e%Vp#gnxlkVYy;Q_)b=jk%#~8rkx(pNfF# z4HLqR9`bRHfy~P5S4Z)?Z5KWF%835{GT09X7rP-(^9~D_4agXA!4<}2YHzIhyo+2a zSk1S0`xj=K?rOXXN6i)tg{2|Bcb{Ew`||W#cDpgL({eFa8}{hY)#g2Up^em(6V~+8 zMNQC~7~@5jU**1R;*3bm0V5jwt2SxpgWf__L`Oo3SZ{2J22xXUcdru}zOa>O;@LL6 zi}DrZE&(eJ-C~KgyaAvSa8|WWrKV!JUqwj)8-R_NwYaAm>?ptk>IZ-6xg%C>9`a}A z(})klOj1svNh=1xf9AE&48>AP-aN9?~yWsgO$&NNwHZ5u}cU~DrH0r$VKX4If%gq*$g#L#FO2W+P{w7I@RGV^B zg|6(RGS0>7H~b%PsojLMs^Cg~*2GUJi6F#!nCf;~{(z}68cS>x1DYgW0{EuAsU0fp zZfu#{9jbqn;)^&S&*`Hd>M-AY8kDo2@@wo3EvsDI)RgH#t!@%YSk+r@ws6-xf&5SC zR{eQM=AI8Xzc@i@ip{(<2MeFdbFnacppwKvZOc|a$6*8fj!kZZ<2%&SL}4tyLMA4y!>LfeG|ruw+6jk%-o&rub?2x^J}06(<83$oQRbDWlKky8KQ5w zX&5c&8bON=MmiaGPfd`_F>XCoU(h*ok(0@s7K-H2d89h*$8OMdeKksyNzgbE*+nIf zjIkiSVDT;5&a-K1r~ihgN05SksD0kiqGPgW)PG1{O>&C0^~mJQ71{QN;*8TDSme~r zVgdxk*8l~v7@RYbaCG^Sw{#ls3ix`-qtIeD^WT|ljmZTtAG+zDqbU$yX}J{!q)7^D zLp>Y7stsq%n0cA$!@q5pj*)H74%Ixm8vUK4eT-pqLJ{D;S5Y)K8j>!XH_ZR5jdINj zi$%dP8hT%A=~dD8@M5|}*Yf6O3|5Svmp2pCacTMdDxz#Ta00Bk)5&(rpCC9i1e#O2 zUp%`ZM@er!dXsFqhNE^Ias?sokrY+CHej^M(21e)S^V>@T337@4r{052S@OgTbFY_ z=ZGm4>;mLQY`%lJL_%m~pFfd^T7@AC?>5Embnyvs!~#may!|kRYGDw&V`KFI)I9~- zBoa>k*zbW1cin-r?a{*qaCwdS09)&dk1+jGbI%Nz)~4oJY#FN=vCuL%@N88VJc$_6 zDe#Gt++Ix@XrYnfxoH>bw8D(CBF){_-3|$mjI&pxmOf&dnBIPmNE*3C*iFLVC$MKF z@H)47LiOGZp1}jWDhnI`Ro!>%gbSf8Bv`Unn-^VojW1|yTCq zw+xT>A$)Z5u&QZ}8f_w126Mh(bAc;RybssOx}}i~1$;kfw-zHB_HJ=PK9s>!^1^{6 zl9_1FpQndnF~*A?l~dbp+1!>Nh);ijVEQA2RHE=O{u22}$9{P7T0Q0IPrfL=ZpaOR z7{_jP`Ns(NUl7v&T!`X7HXu;w0*Y#vrKlb!z?}ehdbGkz5=lqhPIAmHp`CcZJb=+q zAxQ;qX(O4r1^Sc92Lz#IgCc^(*O9z84KUs4uq~wGgA3c;UQ2%hL`MoY(mcPG9=gAW z)qlpRsh(;Up7K+4)thg9+ijC&$Am&}&5Lg-i)$jj!>fK^$ja0DqM^hBe_n_oZ2V$p zO+ul;lcn2*Oa~BJV}t_uil9H*WyFPyIcW@MUJ&5aenL$1Xu-9y5;4_2yhrs9B^zK{ z2w(At*1m~^1=es$V1noFRC|3Ffh0LTB8+$%C`8lN%XBkA^^t-J(qggn;THF%%M!-} zt_k2jgk`8EAq6f=#B67xrr6FFD=kwgY*v{T$k(ebr8UzQ3o%%dS`m)IFu|)D-VzHl zcf}7{APC9?84--iCw_HV)`1n*6<Uj>X+FhB2DAe4sB!AQ{ zkve;sqz(}Gp-7Y^seM>N#>`JD&Ufx=u)G9HPc}VVD`lrccnKbt9%jZ+^HH`74Fb6s ze+OkS+0aZ+k0uC;ZL4JY$wF`v&ahR9=!dy5ep->>#SiR5LVOkHtL?Vj<#TvR++Z80 zdf_6$VDl728yKSfXxpl&oTGNZ&Cn#UR0&u66`8fqaL$@Dg&#!cp)FlGG{iYN>WR9N zwt8dWm1{}HfQ`rmd72{vrrZX@){(nwnHDnKGC&z#*bEa%8>z$e+dI9V`inwojKu&! zvX)J0`nxsWVq*M%>YX(uS#GOGvko(=ip zB%JNjZN!9XYQ?q>AWmnCg0U9kbnMHpwS!j;&27O7I7#q-vPhL^5$U41@&8vOtL&ecWjg6rj?tAMlE~M(@co%7*f)2(g7j^*3t_{i0S! zRHsf!$6pvctH_A>jz_1pYZZu0@?>1O;kk(P0@8S5uYEUNXDRY#uz@0G_4HFH+_rbw zy+87;Ciw&mOn)MEM#tG13(RF>!c=xAv5vH9h+lG>;9>*d#+7g6bPHt}Qwerir?gMt zn7AH%TFy$Cq#wn|J3QzukFEH2(udQUpd*scfDHttgJULt z6~d6Dgr6TI=&HKQj;LU_XocO%RXy|?e&2Ngm+#&`DVq(h;_$L2OgucH{Ho+3Dk%$V zR&hLAMp_~5`-p(}y|#xpoTn?k0)s#hEVuiUFdEPLzcPhlfY;=VY2nP?)k#z6UJ^MFwhXi13kjsBS2 zJP+xn1_ruJJL;TfN2ugmZX&N-K2FXOxmp2{bR*&=tJw-fi**Pq+;L60b}j@L$k6S7 ze4EGtDSL`yRcvf1jgaMCD0!qsIe*Wx2wPLtU~Kzk!#kI1s-sQ72Vu~@#dXGYPBj?c zUvpegg#)n{9Ae2JEUWK0cijY1sgA~xly>w;?zH5vw-NfG4W(XN@A(R}2bqUO;pDJ2 z362abq7~tiwLy@=Duzudzwu8XT>2zA?GWSqNBfOcS^sddjs}{;a3b=U*E~zQtPQi7 z+f7-EHRYDKc3!GVL>Zy+1a@DQ^9m0{4Da_sK{F5zIecu1Z5h3=@rjG0sEnQNn?n~B zA2WrY)qwMN>$0l084O5IB(f4-w~76;XG8)35$B@;1H%j+MLU4`2){}bi>+^G*gfNB zU^7Oa!X;tT`AG+OgGaRPntqXQsaoI4L@^c9?`vV!cW0%DK0M5He^(UxRGLS2x@q0Y zNc=4;B(xTpseC$X9bzU{5z>-%1BJS*2{fHD8|^Pr19xfdrcNJlSG64zg!4{igtHHr z^ZTT7XszcI;CT@~6#0L}g-DqUvpu|Vk6WsKk(9C8#OtepD z@mm$m1WytMS(T-0`*fvY6vtEVqi~$z42yN#bJJUOfbE8)*F5z^B}2;pT-4-)ncGo~ z?mhp&lPc;~5By)x*Y~-VFov_^>6hR)(vn~2&z33#Q#r2YA2Rf2Yi9^ zw_-BSUyl7h{xpi5zh!<}U2z>blOR&&YDvMU=iH`!pA}JM`-u%{L;}c} z+ejG*DJy*3EZW4=qE#ZNr-&d4~OUV&{L??hT+=X27RvMjOr1<#MVo8-GeC->id+#l3k|~eCj#QF#{j!(`A1A zHjuW1_cS&&@r_`zV9_bC=mdO2U39a;8#qx#Bo4ntEiKG5W^jw8YDcaaEi-3S9&EB^ zV#aYA?vR&b&$fV5>rg*-v0h zx(am`Ap%YY*b8wN$NuQ2{4w2h^$GM9q}Lbs?5(&6mXXtOvPfjgTX7qIdKd|YmssfA zM{W-(!!g8)mo+-EnsyA6Ou!y z$#9(VoA~k4%sUeWNS&8roFTDK@O4Gf6?^5T=ysu-K9P&%eX@_a>m!#VZ|iw ziz;Z8pgiALw+(qrr2GGOH@A}vXj?3r=3!4FtKu8}rYWlfD5#0T)(qXth-3NVjMdske2}(OY`q_(ycy z9TrbXzJP+)cz$GS@A;rOBvLfbFu^dHmQ01dg7Sf{&-TI6&8BeX%Q1#S!cNaOmeG&r zeLaMmqud}QvC067;!p%_pWvYY+2(e@8L$~C)?E4gbCTT!`V-O>pi;v#vdeKSF5xtC7=V?l0N`NH`8R&FwsbT}e}c)o_?3>68E z4dI>H@q)({*#f!vTw?~*kiV46W`{>(V>MmB2Y(3-{~Xzda$^)pV8;TKI9toYOHeEF z~U*d6I2KGlOU2M`Q$>ZQ2M9kTGfwNn6lX6 zCCHdR3Gq2QX3M#@IQk@j&}Y;|KS@jFJ;d(y`)4mySnAcME|7_4DbL->)k&335xncs zbI5{~tmzkke=t5@%CHE~{r`Gci-ai_nNay9vR9snF3$OqhZPlB5Q$g{U5TMkNY=5XS!qZYM3KG;Kbim zR9vsi*%h1q{%_H^iBD8tbX9~;{X`ICeiJIFa5YwqpWkh|&P^lDf<`WoAdcP|Q5-7g zz!=gaa57U{Zv>3~lRm1SuYNXBl&2-`Qo}#5yu#@+BRX!hqql*~HK45&&q*8)=IjYQ zQ{f^GvrbWpKScetdm8b}EDAN81c5P@6oV0h3ZzewOmaurENX6JG2TPrHX=%xj^P%3 zE%x4nn;!6>15*o_3*FUfQNUuSMvy?xLaa21Vi^k^X^g zJk-I}INpGvuf;*fbefm0FHOgSs2n>}7TI7MqXs&Dr;esp_Xc48;yuc3NyRtV^v6?n zhp*jBW@-k>Fa|G*Korvg?Eq~VjynE6K+->{+ihB1`b!SNJKe?PIC9mHjazJELi}M; zrtaN_17JAVs5cckETfEf$ARnplDthiB@p6fNje3sY21W(BX%E{)y5dmV=hav*YF(tkiv0mRfo3-%FcG?VLnMGfmn$c zG|KO-pnyCzgOO3-s@(U$TPQ#rA49219t0yV$lB|>=?9Cmz?NO+X40G`LQZ>YjA6|I zx|7lSs9!jWJ5euV>cM`FF`Y^=Rxb;i6RavfVvYLH;IKI`U#rDM1by}snJ9$t3g+`k zkd|lS$XuJxmiEY-dA+8N=N4ztin=6nB7TBFp)=#n9Q?x#OJC>c`zusuL-FN1%4>B_ zN-XJXw#9)3wX`?)bjzT-Sj$mrcKdCjj-InVs2N~?zg)Wje)i3LbglG<_#e!NqRtQ~ z!dd(ly}DR*uSvG$h$`f#2q_}jUY*7_cSuRahD02is zc7Jj=qKUY{Z14Af^1=kdKQcX`2fw%-e75RP2_bDzD|o~>LK^bl%hSR{*piwRYbcl! zP;y-8Mi!KOdB~wfg42w196=cZR007O;$3b$z-$Yas5HQ4-`k)s&Hf!!$>{aMrw-H8 zR(V;sHgQw3->bv%DtlA4l)Tn@&^?>PP3&xz2=6s^T2^tF)@>;YLmnn*CEJtyD+r+n zM^qFMRFki-piaQ9i22}HDb2T1&*nka-hP@IQSyPOig66(d1g$&?0}h@D-bV`b<42x`F>c)tpMC3Ey1 z;lgU!+n2+Txe4|<1V9J)tgt4ffW4v)b)ByK!7qAnwr}vTMfJ5|0ZVDk(OFrhha_DT z3wftaa4o0Zm6jB!A~YFh1j6svM%Fm^fKvcFK*YcC$nX=FDm@V7<+A){ZZp?0JE$Hc zPM$OXeiKv92iW3Y@Ie()>EwujI}7$Dh0^}{WnC)-u(ApotYb{3f1|qA82djaQKXtS zqxwtfvdnGNi;O7PyH{3{?%yXm8j+SkhI}%0QcJHJ+)xrhs~TzuRl%FGXBQrO3|hv> zH=~cvP0CC)a(lqZ1Mjz|d1+g0Q5FP{=<(fxFg960_&9q5AFO_uA+cDC0fofZ5J3U~uo$qGMNDjg@g%T88Z-`EM`D|SH~41yyBbBayrlPQmMV8X(b3}DX<7#7Ql1s8kx55lrsbV77{Y2h2I<{xux%EWe-GFvNLc`Z=#i#h zY10$Cg|=R$m+_hoZ?-=B#Nsz{!J9H=T|0x&5{XaIFD>-~zT^#J`A0gWNgqI~%=GE4 z^g9lPIQ#eiT8eQRHffnH1na+Y&PY*aBnhh$4^qFIZ$0CO+kUihK-O{5>Vuf+m3FYa z0gEeQz@3(4VJqHc=E5&lN|8q}$H{$nIw5hAJyB2$sBHbchjj`}t85jU(jOZ@DUsK^ zFnz&IS7%O?``|Slt_K9Ib`XGpaBpp!*m1~t^=*Azk90V|u`vJcvd`9K>Id*RNHgJe z=RM%@@YXg!CxwQsIa4B5E}ku*VSdfsU8p_Pw5qT95g^RtfGi1<;*2ecgl(hzj59Qm zrj;)Q@-}ePrpXbP1>zVW<*$c1MIR*R)q+Dd+A*-ISscW4#>l3?jHCL7!|oGfLyfWi znX=(k+6~eCf$pW>e_L&u_>zhbGKo>;@G!E5E5;~U8-xm$c^)8g)&=L1$*kNZ%Tnwa z*uni;Dtj3N#)1AjbCO7XarW}y56Qlv`jt9Vqxx*U%DVK;#`Fk{^X5xKiZsF6=$J_^ z2r?@~YnVdT-#9}sT`LdHJC;9ok9$~n8He#a-eG$FxnP|wY(W~*fNcp4y9IYyCZ^PVH<{hGR@mIg&5AIyLARINlo{+=vM%Vk#X?-mHDX`PPw zP!pfba5HV4r*by98Zf^x*FH* zh9=%y&8O!WefoKDj|2f1fUEITyX$X*eH~vq;)H8+LGF7&C$r-5M5L^M;Xlc{1h<

&JUYYhYkb2BN|ajK}3 z6kezPkK&I}kd3}9N??$$SWs0{I}A3Kh#7_%DF|CJ00+Qt=2B_mF#lt)!d&{R0gGlu z9yMv-RaU7eG1-_h8sV>Q-S5k5mJ4{wk6TWE%!}=Jy_Ktih7NTSg#bEjQ=;^{L@c_K zKb4Th%d@tQ1f6skME42oUEKY&whdMTLvLEGA03XbutD5O$>lF+qT7>{<$K!0KKn6J z9(jRMZ3P|)sRVpviD7WR{m2BOYlEnBtEFs=YFB5F`yelTI-8UP`H*=*ic6G5dOh#f+x0#74s>Cbzh-iqZC7`;t$xdbWnol;`geH2Afxti zTrFKJPatm&ZZ~^*QIy68W9$}E9vFBQ$C6?3cvo}MbV9t_--rJ=&rK~_&qhR7)a*?& zeA5o|FjYR%uy-z@#(GwR&6nFM*a)rSTpBa95_i*;bO>P3&+o77{5<3S>WXKX6P8M6 ztQerW9Nf1*5G`lN6=b9Cg186r0BBUKi z9Wv%3H9ClE;axxr>ZwCl#FuD>rdJeXPA9$R@FruV(K)XuC~%3OTQ$J|DL$|W9CFPLQl=`%U?2nr5iEh95Nhx2oV%z(TEZxegzhq@W`Ey z6Qf%_Y6uX9sKY$MWcuN1gAwEwsK46`^co8Uv3Ms%}C7R6$Vfx${|ksI8hl;zfu%p~j9l5OJYMCme20Bo9pB2+CJh5YCmJp;6~M zC3VTEp6YT&vu@=gZP|||ew8NN=1lQX zTlyEOZ@)-`3-B%VvGs9>o`&OA!bHD2J{SQ3#Su@zwJ*BIWBtR@B|EG1t;RH+A;2<71$We)t2?a~S#KXL{3K3T#kY*DQ4^l19&qpRC)ro|##L*YG>BMy0F`)@NF=_o%^Lhi`=06*vC zDUgbOD0+jb4VkLj4%td3D9QC-{`SXC743=`v?7V0t3m7o_eWdaZ@>m;OH#2={OMrd zz_2p9dS#ep`&6cj+9jkRpUr6ie|5K0K9x#L`PZFEF~!`&mMq$}F}`kgf3|7%ZT<;G zXbc^`A6br%xAi=_8MvvM6Ri04JfFMie74Buj3G= z)z&gR0Qipp@TgsOxhAM%+jo|&($K`w57r@~^|X*Ej{vWDnrC$|W$aqjp(337Hae^r+7 zs2CRfp@f|$_ObIA!D{a#jaLxc(G<>?P1nVHlkK#FG61(SQRqr?tZ)0)v7~(qfX4)W zkC|d#&w1g~+qu*tlkDBOb9xol-xk>9J3tOF^h;!tWe@EYQq*^phy{rWyUO{V-YOxh zPH2@I`5{JWG!ivB4;CKJf#Sa$&AK#3W*u5dtCrJI4p+lacm5Bg#^M1u_AvumwQDmL7#2|0k4Y_$o)}B?s&H2rUH@(5kQTpecHm2VZu+<( z^R5;SL!!+omu%J)Dt(P(MQR^uUkl9i7IKLsgKp*hbU;2W6vky?fM8twjT#!{d4$6M z^;wXz6g5tt=6XrjNr18}{1$W_-h3OC`M9}VTw~L+kb`trWP*3W0dCgPXWsHhY<_UR zs|K&zYrT*)YHb;MGB@%J6#*9vxQnb$_mm~9%7HD;hBs5!x?htcb|CXy8{;YRj&dq! zjXIbV)XF&j)RI|;hO$^;pRS*VIIN*0V_gPAq0aqy!N950 zSx<%upj8q&Fu7)?krc-ZGm`G^5sIERrWe^pbh2$>!=BGP$-53yUlK`S;R2D3jXq*= z1&6l{6|7amf(|J87N0KJDg&u5JBxFA4dfs|EIb!@6YK?wiZZATQg-c?bVm_DLb>O( zB~1xOlL7_FE4f&s(~d~LP__qdkvouDlOiw8W9_s{pUhsJ|Hah(#X)!S3WMNBZ5LAhS>MMPp9fIvRMy=StclnBi4&V- zgs3+dK+4HXq2KN6)Skyp!&AW~3b9-RF_+{{SMs)ioHBMBQo zcUSEfD3vh`;3`lTG!kfI@~tZPxDhD48 zh?4!wfx+4)fqj63@~!>yj`QH5V^E(Hmyon^mD}pS7zdfr47ZT2@h>-vYH^ z0S)xa_=E8`duB#)Z|~M~t?V2Jke8SN>E?iV*0I`8B5YXMXsNRQqg4hIC{H*kv*p5U zuDNU2)aCo$R7X`M5Idsz`8bZg*^Jm=|2Ocm9vsagyh1^PG)}fE9=vSF0V|8mPFk)G z!E|@j{voT0W<|-zZw$Ih?1sWF1?!;FuA8js-F1CnqSkUZx#<-?W1GCCC>_!9s0FcC8S9&RDKhV3or1M*nz>oTEi4M8C1^W;^DOvXz==`Ew7{%W? zTTe%0{u)!~Ba?+AKLZr6LF%?FZuLN|I5SO^JD&yD2RWpE^6j^oJdOP}qtzY$4u|;O zZPz?Q?ZEPzu7)qDaicAJP5@fIq(jSf^LT$7K|!6e`s3sshLj4!DqKHNw8qeUn-k8m zy5|?awy{?1n|@LimEcY?JAn`o=PO6V7*)W@; zBgJEIZXC#?T)+G!MCuPF_kU*&wxo;)PFigzK+O)xEGRhmCcQG^ciITe4)JI>1cJy5uNf&qXJqzaV%3_esfiJ$@>%BL z0bQ>a!V!sx6cL+y4hN5H^wtlMaxm9Rl-Gyfj?u{kVB*jboch(O!Z_J!zu%ZkSsTEUS;ebjfZQay7&z>#)KJxbK|nI=|6*Q+8dr z<)h8Ociu77ly24{KfA|I||ONU-j5qv3*o0yY{j(sly6ohvFm5O?7>E`?;fb ztLkZ>3C#Mr=ib;Uy<`3#7`&>zx77EMCrP6l`N;`X3oC{$m~48d9ZUjRR%o%I=&!>f zG-PvGn!`u}`XWpTc={>fWK(S>JTRqY?AglL0PfiJNoGsaftmM2SJL}vZwBr7#FgEC z&wAASnt68C!nCM~Zz0&0j&TV7(Z!sS;EV_+CecDdClfyH;(-<&!&B!&9f`g_eRVwb@PYLG27q=PRZ_I~r6 zuSbOv5zg*p^Y}!GVASNA4X=)unTq{#Rg!W)IXOmXm*ZtD_UyG1pVyFnRs67buEh62y>Twl3eW-l$Sx3 zFbZI*+froQNV3CZ)JT1&%qciXi?=?%P9<#71x@MAYzD$88H&N&&O#_2lh4~i3ZOlP zD2n*VRaUe62fW2ioBDEu*ZJpU-L|bQA7_x2yBXx0VpafT$JNr&MuK4n<9~1khm5`S zp~ZCCKO4eJ07!9o0Poo2VUz?I5T8j`u)6`zZwq(9--Z<%WZaWRcfocErQx=ak3!gh zi@ewWw`v#UDf*1kIL;oeoA7L3pnu_=cZPSk)-e05Pwa!Nn zp|306ND_1v{qDBTlC7^v=OLvoFgbvcGqV=Oz35EaEP-FXFHsvcfRV<$ncTm4tWPWy zqruM>C&W3xWkHEOfaXw71u1JxSr9Q)_!E)D*+0VeYcsdQ+`mRIYt9OoeSzXZxQ~Le z>PBpPmr$v)@MYs|-r8?g6|zg&<+A&W(1=@J^fXD)&D5s&+Y(3r+Gx>&vJ!ZBmi29 z_Tpr)+BLE@FQQo>WG|M)cfqD+72NB#AwYL5&a5nVa}#?Hh(lALoklsb2cbYQjIE%BT$<+kivNebNEn!q31Gv`dy-c%+*{ z1Ig{fZBbCt2E=o6#;zzpTQF`QZ(n{!#RHilX-*lNN2qLJA1ej+tmDQzBtP1fAB;4~6G0|XaIg-OgNqzX1;Z4aA6DjT?buF)$S%tq~ zv4p+t)BG2K+OPRwM2w&!Vr%~Hu^5p*fRUxE zXb+~Twg^#I9cj^WBw=LW8^abdudl^f#IkRZi!Ct^!JV29$5V80yMV$ayb<;US-z$wHJHWSb}d2y<>J4A(YcydXjz8WFCWEBdDwWrNQBj_ z!-@*O+v8izr5}O$-EuY)nR*;RKc}@WK?V-lS4&Lx>?I1dGqkO1qO=l?d6%Y$`Bet_ z&b{BqSQ>H}mXYI#Uy3liUo*~58Rc?--l$?whC;w3)w)gh#H@LjsF@F#m}~}r6dKG0 zPgjEFJ5*iIwxF>WNF=RwyO3fD!vld!*B9|xRd%t-1uw@;Whl)wAC=v(_0MXx!TlRY`{P zcMd;Qe!8^LUS}DwQ7s=7w%cyecma(}$oIzTU-TRG$w2?PYebuTJI0H43?!zTCt^Tq z6qin#5OM-9Zwr*t8B*L;KYQi1p`728BX#X{mC(NZtF!uQg$e_9*F<4SXl7!ZS2u8Q zE~YI0B!3Fe>u#RT+8=WzrlrdW3Zb|MPpa=?`lVh9z(~zr*OvlC@52|Pt{OMt-_n3H z3K6~?yF86$qh-PT{onI@6%bILeV#~_ML(%Bo+MWPXZ~74c(#C8r{&L)3fFyj#9QJjypRNRi-YVh z3b)SN(tff!DkD-F^6>YFbiP1mz*bxcrX}h-0cZw1kN5y;$^ z_4czLZp}bq&ahp!?KLcMHam3|K$ZP5MWNQOxE!Aux_Ox1?G6SwVy5SdA(U z-=9!^|0f>W<-fF=DS}g8>;$#7kP6@l%hgK8*Rx!XzSDg*WrtOy2{?dC3=PWkqg8(e z6N@$Om(8`HX+98SDK&1I6$Z_V{p$&$tOqlMo1Ewix33m-X|L#T*X(98)WJ_&ekaG55oi#)wY7IHg#xwS=MDh{iXi$5JF8y8 z*fMjmd7A%6gB`O(+DBkHj51)BsJ7NT<2H>Qb=$}u1A$Mzq1%t*pS|sMr-%ekY1XY9 z48CC1JBB+_i5O6CsUj*7=9&C%sV?dyRLcWTgT{ zscq3 z2VIJ6_&KnqLcuXNYxl9=0KiQ;Yg<SKZ+V|+fX7u3AxwYFI6Z$m#f#rqF=j9$eUJGJi-?$#| z&SnQ*7L~$mql?gVSgds9)rU^a1r)^shAX&3IUh?_^r%jfwuo6Q-YLG>8CS|W%>ZWZ zBe>m%Aa>TYYg$BBGAaFKBN)!nY3DjjVxxuJEWNl0`)XoxkNcQGj58j7L2zKx?O`B8 zG8aun6rH>t8=_fso|m%$7cY=gTEp613Y>AzvIH}LRX){(OP*isnh0gW_HghH~&^* z=!fW+J}hS^KX|u?q`*F-=%H}l3|={oVtQ^>teLP=Y02>TCF&*(An_G7@;C4S?(i?# zSPwRmxD7!^`XoJ)Lr0K6t*0j5&o`>FFZ4&`0CFjz^3ysx%P(n2E-0MfTaxlxhq_UR~Kz5CGUOur^GH~GF{Z{g1Oi3 zFUwzSr&Vih(xUI~FVAP$^zgpnxw5Y*b%M1T_cqYS{2ODQ5}vEjs{_`pQ67- zo8^cg=oE$`RThD%=r@}alARL<5Mb!r2Eu$K<=ch=GS?N5ej zij&N4k=ft>@e2;GCjMG@fr-WJRUAR76{^93dXdWpi-(`O7cF-jgdhiwS8IZlWL+xl zD7nM2Q1Z6D&Knf5z!tK6{_{2oK#A|PgaZkdk_P}^N^3sKts)zY^!5c0`xJQs5E&5i zWC5c=J}8rfPH|-BB;@6Ye>f1m@M2Kw3^3p`tL@sv0hYiM-hxWg+Kg*QdSsTB2dEr8 zc#mVDtwQJ zFTgzbOi?E7cQACLQlA{p=M2abv|+>wgMdOSLy*{~CGp`h4fpl6LS;5FK!1L>eOfzm zR@de~(cDkt>qnOSI6s}({2}`^`M&Qvj|Y+?8+#Pr$spV8ui2B!Vw!bD^z|l1zAjM} zvEW4`5lsiSgo-CBB_`tC2_=ob)|hZE0Hnofx|pi}`y$rJUIGbOiBYN9@zKBtaXkldN6_|pL;PlEigRMfc;kglkd1qjY~ zXuDdd*sTFHShxM%29X$oHcU{O)es#rzZK_r(y`%i=b|wM%x3m-vTKF_7CSMI$=B1h zq@j1mGgSyB%maK_fKulvzkVMQN{xAh8Zn~eEy(l2S_HfjWX%}q+1;r2P8d|HChkJ2 zVpnop$YMt75^mo=AOJsMFoJQ6g|M8&Cb?UssVc}iWA_@%U!BA1v+CE~swh|@N1J=JygH?|CBJ-!6$rZt8pwbFO zOMG)0ZWFAP*De<-OYj^2$M`xbb4&sV=LikdByh9n2{T|#l-giWd8>1YTh0AKpV#4{ zHdw=cr+PSom$qDfPtBYrLQjf4Kp55DQoWADQ@j+sqptCFY1fQ_A|(0jh5oIDA6f}p zqle+rV(l9m@`$tSjLeIUfqepIT}QxrFAuyBlfb*z;wf*BRQ4S=1RQZ|5}rv7mNH=2JWskP1( z_ECeYlcKdT&{^%42VQIcBPmR@P4mqe<~3Y~LUalt*JPfqQ{+oh3 zOL}{tRjjSueCV5HSe4wUyKS?1Rr7GoASQ+(t9#}W0^q0jmLMV)dKqxsM)8RY0%Dom zf%fka!IS%%>VWd7HcZ^yr)NjW3Mu1$Pza2dvD7F=kgQ{Ko0@^HvU=T)zB~xjS&R~M zI-DS04VH-e&f^TyX{B9BJm5vr)>-lj9$D!a-Q=hcG>kvdPZRYJ9WREi=QFmGmljcz zydu!fbz=x2ksi52@u4D+DG?>*GXOE{{3@o1HJ``Rh3`Yu!=^0K=9YJW;QomgzOtlG z+)ED{{B@Kf-ex9Et1hLvgkxDCGmB5_mho~Vqd)qv(dMWu7{!oTnLwS{XZhd`{D@Lj zJEyf~8(<*9echUDFCjoU_ZYJlXZM0`5c$Ms3w#>yL!lHt>SHR24jrgd4xdX^lFo@& zrXoBio`-wlREMEzD-z#T68`mf@OY~cgFfkDO@T$b(y+5%I=-c6c1h-932M75$j~Y& z;bnJZIJ>8?7l6dXZlyP%X%-i;cK!@H4Kf{Q=5v4Fy>f{4(xa~l3ZxNJVI*_4+>_+L zl@|m_kt9uALao@Fj#UC7rt~PP&K$~4$Fqm&_$R5ul*wPhEo2qzlhIg4wH41Wpv$wJ zjHJKce?vJr2Fp73(u>uPpK!p9F&7`SZ(EJ-KhBNA0x$#xT6#)mIQ0>PA8!R�#2Tix%W)ncpYiOw^xe$z!193<9Wu)6x4S=KIt(UvVq0 z5Ik&h_pE-mnVUQ&7Sfl#af}!AIStEKuS2kDtW5FqcL$y<#95=W25h#Q3ZidE)`iH1 zQCg*&i;e!?W61!0?bietWzccM-nW^bfoUbn3;OMHwXA6*Q_E+XIQZg#5ff=;aUOwE z4j2mNN-%o!z4C^R!mIAZR~clxFZK#39T@DR5PMS0#g!td-L#}5{{!j{!JJ^Z4)d!1 zDM#61iG}x%#w(D3j}HXm96i7AFfRX9Hsu_aJ2wu!jFu!^tY#XFv87MRhczDZEW*s> zBx4aV1-YYtKwrk-u}kK?CykAup4i1I%E;Yut^fYV2-u0un*3kb%8&2va#**Gb+RlQ z5U0-HLc(!sQbHl0i$t@ASDJR2@Y7ep$~Sw-w3BRsDRBgGaEBW3M)SkyB3F{7rHX3F zb7(~A0giC$%nZm=O+zWv14Len)O%>HEyl;pSW0M#b<`J*{Ami&<~*#085`5!5tF;; zg6&K#a?MY_ua9VB?EUGG{f*Q-01Q|PkR|)qo-Xvl6}+#r#!;(Ru&Ce(y#~8ZmLi{i ztCG=Jxm!vH^q6RByrvtH9DQK&3M%A@hZe*C>O8c9*B-}+VG!nq=MI-&5$d{mLq+7b%ngn8Uet-b-)m%&G*yLPB1Xg}eJd3<;1F zXs6?0H;)}x+Wz{Y!1q)_md#QL#3IfG&0Oq1i){1v5((j# zr3o6PUOM5+lRm8y%9*geLo+*B8e8JG1ANeJMkxn5~u5$$O0+QJQN_qD7kY#+P z7Oq_;r#!3PVhA`MInBJxJcymHS8~<3i8bjcUDD&km zw^F!XacUfnV|zICRDEcmlO`;*8C(WHeEC;LAXd(x@7EwW!BpU{rl!!C)*vmbN&2($ zZ6zyFOce#0SLQ%3uC(#nx?h&66HP%yG5GlZ?@G*^hwzcf(uh3f@nvExl-vHgen_MYWVjZ z)!K5@{0(%8@sZ0G*R;tZd|cgPC81?&HsvS+lO_9};&)-PxBYHU@YEQZIoUkD^49a~ z)0ZR{&?_HJFeQIsa3U><+z~;mgCfoTC{lJcXp?M4X#c4Azzor00iL~!PfwtrzHoeZ zaOe>ntm12`uJa0ta7VW~)*Bnw*swZtrnuJ!1h+#Ex=NU#JAt8Qv8Jvxo)j=A?L5fV znodA#{e)aBr4d4cWG53wNSkH}h$Bs0-_JZSB!ETG0`z;&SYXMlt;)9u@Gbk{v!ahl zXj-+CLat70pSdpz);dSu!n@HB2vgGXi7YLm)fP2b-_qPd{>8y>Y?ssYCkdAI4mG0+ zPW7`yTWwHGKGw!*@*A;YIQBrnIMnC&J8h-aKEyC4iUz8xM9nZ0kYFFu^GEp!kG~?M zhHCb*3nX2Go7Z=@_IK0{0gppf$)x2i1-4S+b`*%!{r)26OZ7*jfj&}$pFzUkmP-n1 zctBT>yD%XBppKe&MZzC{nwHIkH4l$K?^=NeUe7WLw{dPZ=!;hZ-L?TOk5xR~oVE#n z2B)7o;x^aP%Xl@?ziH%)Y-eygx>uhsqNVTeLW}ShV~tK~g=jH;n9kx2SN4{WC+YW( z=h`^l2y`wX+gYiG*t?b!g8!&(BuF|I*S8ecyX9&Gkw?T6EL2m5W9It8TGhbtjKH{y zIfNU!AG0W!8?9B6c-#G%`EYNKee_!F7*Vgy>Ij`3`+L(uXxl!dlWI6Jm<(fZU#?jj z#`LY>0%!e#TZ`hQqppRqB7(Dit@%B4n$>B6g@|{MYtzt%CFudn0A2L+V9?)naRMz9 z0kk^l0O$?JK8*IyI#o@xQ8zldgp~vzz9q*%N#%NbiZI3IIm$lD_%{+!)5)HBi44t< zv?J_+Nb;SRC(CU6YP5{LIqq@?tDP-%LdpmPUnSf-zNO;~%}D*mvmg7U5_VFKpzvkc za_s4`ZTK(CWR1FJ54!QdLsq1gaN$|?Qg_@~x23{bD%aNaURpwW#@6NUkJ45D@rK6v zBV>Gxpuv`cEfFtAZS62CY#m*lv7ViT+!YMnVo{$p4zM4q*igLK9R}>BSMiX1nntfT z41I3mFA_?jzGdfw4eiWWHZA_4a%%33v3#OXO;P9h77mM~Bo)w{fO-0HaEn9#n@L>l z;XLl}3ki7$Y_b{UsDT|U!@1BWn}h7u4Na$ zlnNJ``JSs$^euiJ1|5xf3>S&3KgnPZNg4g%^Bu%x6)b!V*U_LNDyk*EuZ-v{@Yps$ zMIkgi5d8-}*&Bb#T^d_(ZyGo=fcC-poUE|nRMS!`d8&~9JSp3jzH?P(RN$$kcH)zm znUynfO*1KTSq4!@L{%!>G*0$Y_GtTc9#nuEn6zdJrp+vcKp99d-FB&hK8_t6nTT8CnsZ*9FLv~n0{ZLte67AeET9^ z^1$#SZ#N?()-N;h7Clb_m`{<}$D zf@YR=zr~6+j-_sgRxgjY-Sf0`>@O6LZ)OgZisPg}emPcA56ekNAr=w!%ZOgfJ+wQi z-xY&IDHc>t^kFSLuh5{16`)E+61fn71@uud{mSMdY2DPIP~Bjz&>p&k`MoQ0UWqkp zMRTBo-)XwYN5x(vaS?FG0gu@?GpC)y7-_er(jyOuR!k+!cs5ZWCL$^0dDF=zrrUcs zDvBlK?S7Z6lg)W~X|Q{)*Y_$Sfr~?U`!qwJ%)I$_&LB2QtciUGExP& z^V>cizty+@RRSbPaOC1%?y@7DPa_myq`u}YFy1aWb7vo4Hrp<;7p;wzk|uvsi4#*n zfV0%O?9kMze)23Zgm??h=|6TE1v`PRZAM)`6ucL-dA4Bpx^+?5tD$uYna^NX9an`k zhl~J;r5Y7EnLKC2`f#t(!{NAfiD7KOjz9O(y0(ZcbfYC zB)jA=C(Df`awX&prZ-l09dM5AneME^;CEM_szoujGwc)Z8aUMqbLY@BK%WMu$xsXo{C^9D;NEb`j-0DbcSz3Esp?@?IObaIX`DG>N3X|ct}C~2An>#8C3$}B zh~QGj_*nxQ;wF~-3}KLV<4rUjQ>=Km(&{v(6mAIR9>?;VmySb5Y4?nnJ6)$NBdug$ z_|BvPA}g0U`bQ+s%Yi}TH+1;!W!-iRXBh>iD#i!*XW`THBhH{tLZjYjlYgh^-`JO! z0UY7Gyj>+0wYS}bWsNK#Ix^=fX2<}We^|D_{!OqNXH}%M-wzTm8Ex`~GC}PCDjMhO z>I6G<3!97ad1Y53R+Kw~k(<i8p6Iuz?I;+{+f}yaHuVe#UoO@W zfv+d0F>P^Cc62ojCXZ=LdXr7!JQTOzFi&G~A zmKn+C4y5?;fpA7r4#w&{74q)twharJS5+BZbcSq}-5n0qkk3&)_9D5;Rg=;DlAAtI z*0MLRG$9&0T9J@-!htlb=p++6I}0Fo6%wSon#dll+-USyVg~oE;fTz5c^m#7Yr$ha zDx=P`&#r$(0I}i(g!raUGdLLUJ+|sXU(?E&!9CYXGq z==9*!{Ox&sg@ly zktU6Z1d@a(HP5El8furpv$?1-+xzwS6Lg{(C??@hIe;-B==GAMF~E2;f}T8>)fKHf zAe@yAG;^ND;zec8k|$xKo`WEmM)wo`|5mct0gNM#(i+I`0Bj;zx7fntR z6VqxsHWFNbD7yG`U*|FLq&E95nxDPVr+rc)`uSZ zQFt?WZ_vAiK-cd%jRC zmh~9{(WRoh8ROGyYIRe8%!NVD!FTKT<0M@MJXTFxu6%wtYAN!#n{+uKiIO8S{rD#x zAGk2__I33Z_53Tt7w`QdqS^G!KJ;VAh?tOk-0$Mq)IZE)VIMU=$FZnf(;3q$j}Po$ zf5L9xH3H~%JBXEgqjQo9R{54q2cMg+Sdy|AiR4c2O&h7YbUfLhW7LnnDXM`$d;5)@ zr?a`jRr(UzTNQ50liOBxWP#C(|#^H+Ka;++llf{aftXydMh|+8POB7k+?w-sy zw|YNDHx1T3aG-RQBG<6ID^hRKw-v8cE|m^t$YH-;4uC5(zZ{k4y(qMr7X(Z1l*n&+ z+0^j+9BqVjW1O4>_ty}k=T9WaZ^Nj|o9iJ4f>YTa{IOr~Oc{mtYQ~*{#glxr)1ULf z9*$FK@Zr!djg@YzV$my-XJJN|!EnK-Cn(X;oG{PjpHE9`%jM%uEEaw1VID48ETKWw zWfExhGD{uhFPr8`ii;vv^#Le#iG${2`re|&+*0KDK>NN!@-6V3Ak8GkB6UYAzG!=? z^-uV5*O#KDC+}W&ZANhR>d+4kQ*CS-p{1wzqf(#XnUn>+3=Z*SD4Qrq@K*jinxHCPr?McB>Xyi?(TI6`v;LD1A4F*uMHq!x$+VZnK?7ucZ83&bntvj z!=xl`)9^HXLpFk7^wG73T2hem6D$&?JWg{Mw!%wt$A@&7q~*@KA&PfN9@_6WHUzcz zo$#qJ@81*E$fprD=6<&1B~r);bmOA|k8H9v{-SzhmnC)mBk{AK>$F)i#9-|`0!{~K8@JSf*!2Pn%5J?bw%D6z&i_Y^5wQq>Oowp}X!y5!2M zRe|-;p%HdPcE=F0>pb|*l|EO!M3*!t-|?!2(bJ;Evc4jBBgagz>EIqS?M!@pIhTV2#U+ifdqNn4Z~+u|RNNom##ZgYy95Bs z0KbPnnF<7;?4I+tXc_DwoG?_=_Q*T|i-+lkXfD-meH6@C3p4b5b0sa1x(%rbW zBV6g7=;kzkLC^chUqNAWNW?>pxxRi6n5hL&* zo7;Y