From 3ab0a57b4c3b9992bcdaee35c2f6fbbf11c71cae Mon Sep 17 00:00:00 2001 From: jaimerzp Date: Wed, 17 Apr 2024 09:51:53 +0100 Subject: [PATCH 01/10] gg --- src/data_utils.jl | 31 +++++++++++++++++++++++++++++-- 1 file changed, 29 insertions(+), 2 deletions(-) diff --git a/src/data_utils.jl b/src/data_utils.jl index 9dce28d..8b9b69f 100644 --- a/src/data_utils.jl +++ b/src/data_utils.jl @@ -49,6 +49,33 @@ function _apply_scale_cuts(s, yaml_file) return s end + +""" + make_data(sacc_file, yaml_file) + +Process `sacc` and `yaml` files into a `Meta` structure \ +containing the instructions of how to compose the theory \ +vector and a `npz` file with the neccesary redshift distributions \ +of the tracers involved. + +Arguments: +- `sacc_file` : sacc file +- `yaml_file` : yaml_file + +Returns: +``` +struct Instructions + names : names of tracers. + pairs : pairs of tracers to compute angular spectra for. + types : types of the tracers. + idx : positions of cls in theory vector. + data : data vector. + cov : covariance of the data. + inv_cov : inverse covariance of the data. +end +``` +-files: npz file +""" function make_data(sacc_file, yaml_file; kwargs...) kwargs=Dict(kwargs) @@ -88,7 +115,7 @@ function make_data(sacc_file, yaml_file; kwargs...) # build struct instructions = Instructions(names, pairs, types, idx, - cls, cov, inv_cov) + cls, cov, inv_cov) # Initialize files = Dict{String}{Vector}() @@ -119,4 +146,4 @@ function make_data(sacc_file, yaml_file; kwargs...) end return instructions, files -end \ No newline at end of file +end \ No newline at end of file From 4e97b0b57acfbd1e7f26734cc6630bb5442b1b48 Mon Sep 17 00:00:00 2001 From: Jaime RZ Date: Wed, 17 Apr 2024 12:21:52 +0100 Subject: [PATCH 02/10] chi_max --- src/core.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/core.jl b/src/core.jl index 2df0f60..324ccbf 100755 --- a/src/core.jl +++ b/src/core.jl @@ -277,7 +277,7 @@ Cosmology(cpar::CosmoPar, settings::Settings; kwargs...) = begin else @error("Pk mode not implemented") end - Cosmology(settings, cpar, chii, zi, ti, chis[end], + Cosmology(settings, cpar, chii, zi, ti, chii(3.0), chi_LSS, Dzi, fs8zi, pki, Pki) end From 10e98d773445a79e26d0b4b473c38fa5c9468130 Mon Sep 17 00:00:00 2001 From: jaimerzp Date: Wed, 17 Apr 2024 15:05:39 +0100 Subject: [PATCH 03/10] pk_mopde --- src/core.jl | 20 ++++++++++---------- test/runtests.jl | 4 ++-- test/test_output.npz | Bin 69358 -> 69358 bytes 3 files changed, 12 insertions(+), 12 deletions(-) diff --git a/src/core.jl b/src/core.jl index 324ccbf..30fb9c2 100755 --- a/src/core.jl +++ b/src/core.jl @@ -14,7 +14,7 @@ Kwargs: - `cosmo_type::Type=Float64` : type of cosmological parameters. - `tk_mode::String=:EisHu` : choice of transfer function. - `Dz_mode::String=:RK2` : choice of method to compute the linear growth factor. -- `Pk_mode::String=:linear` : choice of method to apply non-linear corrections to the matter power spectrum. +- `pk_mode::String=:linear` : choice of method to apply non-linear corrections to the matter power spectrum. Returns: ``` @@ -38,7 +38,7 @@ mutable struct Settings cosmo_type::DataType tk_mode::Symbol Dz_mode::Symbol - Pk_mode::Symbol + pk_mode::Symbol end ``` """ @@ -63,7 +63,7 @@ mutable struct Settings cosmo_type::DataType tk_mode::Symbol Dz_mode::Symbol - Pk_mode::Symbol + pk_mode::Symbol end Settings(;kwargs...) = begin @@ -87,11 +87,11 @@ Settings(;kwargs...) = begin cosmo_type = get(kwargs, :cosmo_type, Float64) tk_mode = get(kwargs, :tk_mode, :EisHu) Dz_mode = get(kwargs, :Dz_mode, :RK2) - Pk_mode = get(kwargs, :Pk_mode, :linear) + pk_mode = get(kwargs, :pk_mode, :linear) Settings(nz, nz_chi, nz_t, nk, nℓ, xs, zs, zs_chi, zs_t, ks, ℓs, logk, dlogk, using_As, - cosmo_type, tk_mode, Dz_mode, Pk_mode) + cosmo_type, tk_mode, Dz_mode, pk_mode) end """ @@ -215,8 +215,8 @@ the primordial power spectrum is calculated using: Depending on the choice of power spectrum mode in the settings, \ the matter power spectrum is either: -- `Pk_mode = :linear` : the linear matter power spectrum. -- `Pk_mode = :halofit` : the Halofit non-linear matter power spectrum (arXiv:astro-ph/0207664). +- `pk_mode = :linear` : the linear matter power spectrum. +- `pk_mode = :halofit` : the Halofit non-linear matter power spectrum (arXiv:astro-ph/0207664). Arguments: - `Settings::MutableStructure` : cosmology constructure settings. @@ -267,17 +267,17 @@ Cosmology(cpar::CosmoPar, settings::Settings; kwargs...) = begin zi = linear_interpolation(chis, Vector(zs_chi), extrapolation_bc=Line()) Dzs, Dzi, fs8zi = get_growth(cpar, settings; kwargs...) - if settings.Pk_mode == :linear + if settings.pk_mode == :linear Pks = pk0 * (Dzs.^2)' Pki = linear_interpolation((logk, zs), log.(Pks); extrapolation_bc=Line()) - elseif settings.Pk_mode == :Halofit + elseif settings.pk_mode == :Halofit Pki = get_PKnonlin(cpar, zs, ks, pk0, Dzs; cosmo_type=cosmo_type) else @error("Pk mode not implemented") end - Cosmology(settings, cpar, chii, zi, ti, chii(3.0), + Cosmology(settings, cpar, chii, zi, ti, chi_LSS, chi_LSS, Dzi, fs8zi, pki, Pki) end diff --git a/test/runtests.jl b/test/runtests.jl index 7b734f8..ff33689 100755 --- a/test/runtests.jl +++ b/test/runtests.jl @@ -31,9 +31,9 @@ if test_main nz=1000, nz_t=1000, nz_pk=1000, nk=1000, tk_mode=:EmuPk) cosmo_emul_As = Cosmology(Ωm=0.27, Ωb=0.046, h=0.7, ns=1.0, As=2.097e-9, nz=1000, nz_t=1000, nz_pk=1000, nk=1000, tk_mode=:EmuPk) - cosmo_EisHu_nonlin = Cosmology(nk=1000, nz=1000, nz_t=1000, nz_pk=1000, tk_mode=:EisHu, Pk_mode=:Halofit) + cosmo_EisHu_nonlin = Cosmology(nk=1000, nz=1000, nz_t=1000, nz_pk=1000, tk_mode=:EisHu, pk_mode=:Halofit) cosmo_emul_nonlin = Cosmology(Ωm=(0.12+0.022)/0.75^2, Ωb=0.022/0.75^2, h=0.75, ns=1.0, σ8=0.81, - nk=1000, nz=1000, nz_t=1000, nz_pk=1000, tk_mode=:EmuPk, Pk_mode=:Halofit) + nk=1000, nz=1000, nz_t=1000, nz_pk=1000, tk_mode=:EmuPk, pk_mode=:Halofit) @testset "Main tests" begin @testset "CreateCosmo" begin diff --git a/test/test_output.npz b/test/test_output.npz index f81b4e78a7de5fa558e5cb3a4aa6e32df1da14e1..bcf2672700ab5d2e164025496fa90d838dd1ae46 100644 GIT binary patch delta 8163 zcmZXZ2{hE-`^ODq-?8b8?SxS?oq7u@AN|rX%$Qq#-C5k9YDan$ekhF-h zR4B>5RI+c`=Kq-){m=RJ?HuRa?!BM8JomZp`#i51PQwqU;l)WdEO-$n+P@tq6Lw~9 zuF1HBV`8HHiBb6ftfmaIlgK3>jlS4Kc0-3D z=?B5{a52%Gi!8~ELIe!@%4*or4no|m55;_FgPzA%bwgnk;ZaBQKC%+Ir&jem@R340 zy(!|!+X)EALzTDckwIW{x60rBdWlz~m8A#}yfYJ|w zFqku5lovuZ{opWlon3m;^na~P`G`y5pqblweD_6)YQt(uJhKzI?a2F#sH?58qA`?p ztx6{>44x6lMtv}7wO94ak6-Y8UrGCy&Eue+Nj%v)G6NrUO8hQcM`C6>#_^e4p-L+{5I`nMQ zo%Iv&j`ztbL+?(oZOYofJ5ddz(PrrxM=8*GVI;MGG6gO@JG-ge@)Vr^kxTX17YLme zip^RZPQ(2=Zr{+azHs2ih3`+~4nh(4Im%SG!!B@n+Q%`>N`fL;|BHhq2&nF|#{^il zp>&iOz2}E4s10(L&v}W!k=Z$6_e4BAbU2nN;5})omvE4E6R!6kD|ZKoy`nrS6{c*) z?NPVr0N0xm?njbG;CYm7$Ju+6;4EojTEsICR#Sdcg?DjCf42~O z|0o;Ms~WCQ&E`d%waH;>R3WrXQ{oBBG9T~BRlJ;B*S3)0J_T9EAmP1ksVLRuT#gU(-=jGr; zZggU!rS4+IB6!6)=WTg82h)$KB(s`sxEa+XxU0AU)PpDN7N1AK2{#3O>&R3nJJ4}l z_;(DfNNRq^5quTSmKbiH)k^@Id-~z8^&)`N-~G+V%xUQFIxSLxBOQgFko|`ob6sJ> zp}yFL(QR-+LQ+qr%na~I``K)b4(#`K_k8Un3tX&4jS>dJkSbfH+S1PqmyXzI3wMqE zW86Y2YXy$(!+e)y$0_j@t0*1pPUOZhR<43`t?m)s9tNn{<3E94*$kZm&|E9I6PG80=Qi(Zo@BlS+6+7(tZXlEc- z!cvd~I&kG_sL%ctNNug8R^wPF)E#y!QI?qjA(k4G^_M1LNGD;>z@>J`KN;&KS=|bi zU?E6MIR`MwQi#8S&D@vU6Dh?T>aFRnRiMs$g2 z2%j2WK|JJn$}L{?tI7WkCME2VM{(7cEIcT2V|nRZQaSXcQyxp@(EZ4cQeQwhu*OsO z@eqI7K}VcERDyf3-@pRXT20w=_JO0Y0`HX z1xuMs`pB8eQ~Z=MnWKi`J*+avv5^DdVK8mv+=QF`vJd@<^Iz*Xh+rd6q#Hh4@t=`5 zAKKGOTRi36E5|}^S!;%QHM^qgDJXIk<%FCI`Q(bKfWGiwz(X`v5fvE+&lOs?E6Qhp zI5sQ58j>jJw&zB#XW`?r6~mA>`I=8ZGU~j2;K6MLp!lp|9)Hvql%X zz^;+2&;4-?C@kh4_&!z*+|p4(*KBgYvBD-!UaSCWH&=gOdp{op4~dKSIcCESIi^>B zOIJbrQc0|pVjSp{RTRImbO-^CFJ7eEIlhq9U2wNO%N@3KE9V>Q*hBRB&p*8D%pvHu zd#0$5HZVKhd~j4z7E*nzw%HU30o#+7`V9|p@W1{dVjZ-HqM&fpkpKSbZogsK3Q9=$ ziI*di6*-bZp&KlR%SDf!-sP!VOuYH}Zcme8+&{ihmKI)(DgSE9VJaWxzB7p&Ak*)S z=N*FZfz;G%zT@B^LSaiYnT0CnUO4rM3QsjjrJ?P3)a|Q(c83)=vNLKmw|XUrlBa^& z=9t6~-^oF@v*#pHs6jug)B{ZI6XR;V_1B_!V;26;nktCZ>vu-mH4PNJ)(2OY2)Hm z79@S+^tYo*ICOL8*!)(9KG>J8bi9nU0PMC8CQP`$0lA?Iu$Ni@A!qFrc1D&!P5FBv za_IxG8oC_0`Q|M!zjmPg{=H-n>lM>J7<~aQE+wsQYA2lqg^QyuagE2|oN<-vmJT(gh0Og;o9Hh|DI1_l42qY zFs@x;l}h8XJ|*kESQew*`hR6f+OS{^$5aEm~dMCZ3*P+Y}I6`D2>cx zxF_PPWl#Xm^Gn5+SnMNFG1#XeA%m1ub{XnSiK8x}pNe#l5L*AL$>7so7UZ#&`Qzc! zQ?NC2_{N=;6JRE-PxyMg9i(D%%D1G{!q~A*>&FhJ0FQ0yi9hK#Ve@0OQqFom4E>%< z^->9g$%?49*{5eANBi61oD~7!w_3>X=T}dt3m)G4u7l(PU!SB^?vCCDqYkqdMYT47 zs!^*PvxN>wZ8;DfH>ChV8hQ!TEh`}Yr&3tH4J)KcPIC_MO_^5T8a&4HyXW6`!+Q7Q zWue8pWYFzBs!~xwx$knCtX8aSQ+lNvrr81#CVU4$IYeP)w7@95ieKlrYvmMhW^ZdN z%$bLXTVk6+MVXP=^Mg`3m)Ox*a$L)&Fg}!_H}d*yk}&F9FOzKXVI_*Hj}H9MA%&#) zMJ7XE5zxJc2R#*v3dmit$M>_g66zN-%YE;>4sj@XdjvsUZsc&5K&p4ATzBo^LrZQfBG0_2aQi~QwyIxqFvK_3_`skKDibH2->|*`S7Ju% z8oevf6C3?BLMs%kDBq6k7rX`oxlvUS9TeDl*R6b1H6CtR)jp{^7Y6o&lS5VACjlSE zZqQtTMf2&K#h+^n?S+`vrI$I|NHF^BW%0--1iqtR*@ptO;HgH9tkw$xaBwtsmcI}J zl_JBg?ht0!6li?NopbzOTl>2{Uqd(X!LSfEih6-q;nt8)7J5xR>q9GAfS)|VjN_y% zgptVj4@Ied4sT)I)DxbMOZCwF;^AN)+h^FJac1S`>uu0*%u+05bvFb?D)L$94S=Tq zfr8++Vc63mAxl`8fVXeU1tW9*!1${SzNigM$fo?V$DkP=nMsR#Y{}(7B6;RKdFyzQ z$I#W5br}Lk5X}W2UnPY4W{l1b#0rs7P5ZFZvtxqD_sznQv@#wfJEX{*d5Hy?6Sj?0 ztEYffq~jDnHgQZ~z`O6vy$~p0b@}RKBUpSv30AF7pmnq+@!*fkaCrXN_DyU_kP$HT zGfBq*GS2wbigaxT{od5nna{i7Q#4bc)z=-+-*LF^V=W1)_=T_j(0AAficJ^2hm{SW zTAeRH!A}K7yX#*&aIS&7<;N_7{RANL_I{!}I}UhWyM8vd9y2vji;Mqe-evk3-*xLZ zuJIr1N1Yw`PDzjOG|b48e{|x3b-AMDrPWmq9>kwx#iT|xdt&Kvrz;$zO2oW|0MG3S zqyI!E3U;Pem12;p}HIA#Shg(kj->mc*Z9XoFaa|AlL>><#k2AC0!tK zB`9@9b^x}jalB&{7=~$sx(=43m>}?7{Z{y67S5Aay35B>vEBT1Q_W^UEf;KreCF8E zwHUFMOm80K*eWEg#w~z++6!5>-V;Pkla4=)bg^?GgrAX85)(uUCoghH^7EjLIu#a! zPngkz$`<$d`3VSS?(brX?uUc250h^7_JDt7%q>o%2G~s!>^f>z0f|at5w>wxL1Rr; zKvZ}N6lc`_$QnHkwTp(GIUVkh;!M2eQMeBxbkCIA>hFf|?fE`_e)iD#!6`9Um1GG{ zgXZrJs2Rihjp8+mN7bPZNBuSRUK(2MyOa!lgu!*TMp$$X2m6~vNmeVyO(U)ey{?<@ zGJX7JuWt?Ozvn_>JD*O!0_TX^_EYYk_tf1J`P$bXd&`wGH!)XZZ@ED4Y#Z$@e|s$} zj`o%_$!}fr@jq`l#U}bP{oyF=jlQ(}=A4RKjveLVQy6V8;{QG`C-X8H?~KI5<=D9L z;c}+(V}KAtiF$NY*Ti>+uq);&7e3uFl6UkVsymj7e8gPE9*MDYVXl(CcKfv*|8bSV zbp>h53#cnzp2YX@a;zy)vRw1xZ{fg`GxYt35>9AQemw~wE2N(vS()4oy(X!XK7NBx z-$khUis|kVbH)9yG2JceiElH#I}c~tBz@a4-3{NpWtCI{I}+I5EqXYd4{69qAO4Q% zu7#_z<8w@R2?sjuTsozYt&xOR3#PlWdPCdRC@CO`#bp^52k&~7&V{h02~nH?T4!*tih*bO&`>26}?NJ9*! zyRk!gpZM4?-JLG7+M<62LJEDpuhj^JPAX|81Jhm02+tj7G2N}qPcs|AboX(1)3J}3 z?q-Dz<-Wmm*L7@qXM2T%H#BeXY}UtgSHSV*W?oEpd#7rjjA6PfRzG4rjOlI@ML$ai z)7>AtSxR1Dx{C`Yo$tkTm;CsigU5e#*Fdw18;kNzk;|`~+%{zY7WS|_iX~u@n)c0~ z)en}GM;>3mzf^9O9J zHJr6DNnJI%KI$AMsf9`%+Mn)mAO|1&SUF5mwVrr*cHIy_{G1D^N|>bPgd{p;&?I%_ zkjSrdn54c|JbsykN$O$0;@4*|Ni7qSSQy15wc(_}Yc@<$M-Te0DaRyr!dH%P7n9V$ zBejl_Ft<4PGgdaj_Gpj#U%BF zrkjibCaI}{D>i(^Bz1Vvc<$$aB=x0~bw4HOnWwIy0jc}Jz3lQjD__G-EORs#S{<{W zkXWu!adf`Z+Km`ot1`x(|8M44xcd2W)JOSV&dBXXO$lo-M&L)Tf0>>WJM{~w?WIML z2gboLiF42qHw#l-y#4dpR2Xczr$u^>N4__$z|x5ub>s$heVi6V9+wDh6>4G#*XJQ& zHY7aEDA^Y6|^w4Nj zT*;oH^<)%mc_wwolomR3ZeSbrwi+7EQ#;&pUjbz-y5$P4dkgysVPQtWSSAoH$aH3q z88HRLIh`VOgI{vsGm>5}OnOH2QuFHJdP8=Bpldl4-ALVK|0o|6Urv5YQho?(eLKps zn=_$#XU>AXa}s>DxIBel69pPJr>7e}VS)JGVk6C2>0@Afcu}xg+YPqfJ`!#zYzLhx z%sDGxbAqkXyZ9L%El|eyX5`@sprUMjfUGJ2Z?edt)`2YlBuW(97in~r5%{h#OQ$eA zRZiq5a>{{+P2K{vn%6b1zaCB%DCAuD_z4sixF>Y|+92V(#l;vb`&-m06mBUPfC>rG zT}9#}@S)x?bF=6q)ab{BU(osktD3&IS+Fpn?CgRcO*1U0wui+lL!1LK%hzbi)bpV6 zDJy5?SwR$1DZEO`M+8;u3OC6O7AB*)!z(r)HWNg@-J;p}6L`>pBOO(i=UGu!)VsF^ zVV|L0lFy8pFbG__7MC8{_5sTk&GjXl8sO(Uy@$Cv?|`W==gBVd5Xk7~84wtW1Xop^ z=${^ez_iJKV!UPtWU%j+PqDIt(KTZ)hQ8Q9-7D%1^~2_HuQKvtUl?g4TzOd-ej-8- zb_;RybTO?1(p`3A9(gHXlkFkeIP-z!oOk-2BC4slxJ!QL_EA%=&p#5GRKJ? ztpE5=j6o@U$#_dBJE|GU3-6Dg(-gK>M-5Fb8n<#ATrAHV;#R(25mMf5YWInSX9Ll) z?Mm+;(U>?sBpQ^W$oWr#dcjlKk>dMCqVV0TWgb#^QIMN|QV=g9L)^RlEVtQ7OZ5m@POGl#B|I1t`zVH6NO|{FbudHRH0_kUX zlqiSlJaskNE^p+)R#5w&O95?H)ZW+PJP0GDYqi=wKW>$HUiXBgU~`*PU8(8spfpLMeXE z7{@LcBq`ohr@}ZXplv}DkNV6~v%D~lT|T?Y_#noy2ZyqnIWUgh@yfe?xOXLLsWEan zfpP5j96ncvDx?v=$5}-!jAK)IyZHrZ9J?xERaG^{v88l~XZB$ndvSz|`5VTurzO`f zWn&!Mg;&}3Eyl4Q%{e>8U>ut}KHJe2jAP49U2_ zK+k^qd4*#X!(K);RNGAl<%TJ~U7%46=ds_$I8ANjy3f(!3q~EH64K`sh)W}Y$FAMIDN%bGY*3g4f#j9}B?j{*J8(tgmP zw_N(OmR>>o5h8hv2kZO!j|!tN>zpaQi3ZM-eL1lv_c;YdQ|z2Rqwsc4gHc$RGh!5s z=Jn}cdFmDZVSJS^zx-9%{Cav5o%qf!_oMM!iv9xQBN~5Yd=$T+K<|6Kl9igcC`xHs zSRQ0y0q6}h44bUUfvw7UQIaur(4r2#$t!>FrsOf&g5NL7(`zE49LJM2S+N>{B}NKN zXRUfmjA?1)lre4G5@T9AH7#30^sniZ4^L(Er7=5(p`5A9hKiwTF&0os)nOE7shae{ zWmOzwh{z6;Df$eTV{wd)Tz-r*q>sMhPu?VbwaI)q9AmjOAYz2Cmc?SUPFT=$#Xc|21ZZBj5A@BLa`FP2Q9Owopf?LGmdf8v`mz4xnu$h zPLkeg7#WAht!I=B@i@lmCm+S*fKgX~#~Cn6BY2!4 zqom4;V|ZQET{Uz1`omdqQjE6StT=|3(UP9z->f)>Vg=b&x_2u&E4U8&;RWWAsV8=0Bke2KuXR_mn^zG9kbo!Yl zFS6rQ=vB1XlD2bg`u`s2emF@#zT{AjD>y z3wXr&cZ0tgY{-C=GV!)s|TDdR(Pp5-&ghYy{=ZdNjx2gmSv_fx;rXcaI2R?uEcnpw=~HlWEn z-I2EOE<4h^?(qSJg|F}~8$Djg(ju9N$K35BFHV%PkKeq1ccj7h*L6Jk{yOSCKAidg E0qh2u{Qv*} delta 8163 zcmZXZ2|U!#|Hm!smU862lADnGT<_^1DJc{xSI&xD>7dqCAw{-YgtimiQ98(#N|G~2 zy5#P1MQ;AHyZS#KzrOo;?Bg*zGoP9F=RNOv&3m3ZkVPEGB1+oXun~n>=pSu+t8M#h zggsLTEG+b|C}r`a1f`$Tj(l%p76&C%LQ9l3OQ^T{mHLR+E-^d4_sK}IZ@tiRK*t0HBmp~Ha_ z3aHuBEj)IUERs7q$Wh+^6-Jz9TSh}!knxdgL)7qDnBq8Z^IW|Lt`m3XP2A~(p}|;{ zb81;Y%bw^^JarR5W7Ki;U=GwjF>W3nNdqpj;n+cR7R>D43V(i(1nX$KBi7VJ!($zm zXG2Hqf+2`!?QIvnolsUY-F+-s7iiyPKMMC4fJ)YVe`N`P6yD6<7*1o>-|8JO z(elPy#z&M<)8`Hk6IS z!+Ci!3Q+nJ^wNJBQlFPQb?Af?GVj}RdiyFNls=Sqy+w;16{?Hqw0->po@|eIw}t$M zV?WeG&GvkSetXL)-qF`!qc0alvg_kEQf`A(NNHtg^m**fs)sjbPXWB-lxu5^ z0ccTAd@*zo;#YTCJ1pA=;1_aD^@4*3-1GZ%Zs}P&2zt=Q|LeUjlqIrytkMyNAXhVe z1$Is7m>#_Rh({ZO8q~J9jHyDmBI&-Tusql)r^K7~iu~i{O&x;YDfgcWQ8p00$vxf! zu}h&72#*(uk*E5gpHK7MgQJ5mW;Y;Nl{X4U-snCH=beTsV&Yd?E&=T|pV`be$bt3- zbM9lkya?6r5|TVRBZzXJsIWg}7elwnL1PY`(r8X5rVsSx(WUNFoQDz=QM+(b%QBJ* z3LZK1=>b&@1xSwkqODOUqe&g9*OMNqC{ZWi4Br_ARMb4ECiqSQZT_umb>b2)dK(bt z>Peo1hw)a5k8P&mW*T?;pQSzE=*xC6_hU0`S^}@WhZewEVtL{Fl1GsJOG_tSGzSP< z#%ZKBD)=q;G&1l#1Mj_fHrdF=!|QHGE6tt=u&~MLx)x~X4_pr2s@wS2!THVsr5m?M zVCDI(LAu@sJgX(1MBlT9JKh9a!|YW+JiQb|iVZ>j_sOz(c@0Q(4XV>%k^9HF1@r8s z1V#+gh7O#koGuloeBpE^7jgJd4n(&>@fn*`Ik^td(+u*GUDN}jhriy8E%^cVI;Z}) zX^uee+=&*ev2kD>%#42>IuE5S;^#zciKzZ%zP@lMH*%+WW{RjULGA@Pbw8AZk?ur> z{d!w5pX&=))yQQZZYGlJfivXLd=U zSIObic8*KYQ9BPsXOm8Fj04+A#r^h=rH4}Y86)jXb3y|c1Q}erG)SLT4{pQkkGy|k_MC+ z5sN!}6d+gBt>(O_edkhw0F zYP_2$Qxbt|e!R6c)w+7Escumw z3FVEwc_+;I&tNh_UUw+&M$!eoQWqCY=aI=}m`*8_$z{Y*krHPc#YfiPKvQ6P;PO~q ziOb@LDq@w-OUix_Ifw@WNnSi`dFNH_EkAjYiQ^FNo$$earM?@DQVDxcvf<0x)ir!)^| z*jK4Fdl1pfdoc=0nLOxa;!7c4O97<%GHB1wO(KX8XzxFrDusA+W|m!AERUu(nR&AW zDNoRT5DEGcD|n=Q*oPz zD|iIb6N_Cgm`%U{b@jW##(C%wyW%tbfru(YrU&m@@}iE`qPp1%0pz~z{T=RE5hNq~ zjr#hqH1hdio_8ci9!V8Mq$z1BqxoKw_21+*(AB`~J(?-nXx*>UspBtn5D!6N_n4&~ z+Q?x(?9B#VCM@Z+dE#rR0@Tg5$`=-Sq6jrvIm(S z89uQ4Z9C;*V{S4v z4d%8hsL*Iw&?l-v`GOejF#0P?awhq*%m-SWa~QjF>xt)8j183S;!0k8xcR!MN*h?` zeJfOX_61Z@rt=Nkx`D;yxsIL7511pG$5r+Y!`{PH2UmRw{* zF&ZzI-BRX8rtDr?`&$k<~>%Fhf z6fZ;dll+}?<~a9RR+Dph*WzVJ=0)Ua!AD8-&DUgV`#wPww=Lm}MmQVt?yCIe^?eMU zo~k-w5ceB|lN>fK7VLoJ^WO~a9;t`mxxvc@mG^-Z+MTjXEguGk&r2K?rop+~Iv*p& zlQ8f;vtxB*G)Rkm<+qkT2$RN*T{7YUki;X`_JPgL4Tk()EbsBP29}Ch?Zsgd5Rm?& zyz`MJ)Z46WPB?4;fs1YK2K>;2I`z3Tzkg|g$+gc;GA9*b%g7;%r&SXFF3|cNg$qJU zbfq!;?Y#5GYRWCwL*!O2uE~(|ogi{8;!k*RKOE67+9I`N2zr5hElYF^d_*Tw{Bx$E z@v`BuT}xTfvkh|FAD`w#q05T(qoerIqE}CSZ=4ZAZ8inECti!8{-@jdoj# zQK>vSa5ljGyNnWg_Na5wsDdi$Dh+EDU8|19Z>CnMzs0Fgbc&93*)LV(BfO|iLs$`A zH||c36qZ639Xr0a>hq(pcR^tlwevtcQ0*NnJ`G$kA#CTL_re=QdiiWuJ&dj?8gDbZ z1J}2Q z-M8BfU+s83*9Fa8qVt z?Z_<0UVe3PsDgsel1X`Ubhy1Tyjr3EOwJ^1)-Cz?WEBg_TNxmB9uZONT!N}*9v2GtKGm>GWf9sn z6xmmtwFJGCN=XzE5k#Jp{?c3Vf@BnL`r_59y#lB$d`g|;As^CFSeM3~!iLuMEJwSV z#-QhR?S~b;ePBAvrWC)q8|uufu4O!L1e{&&9Lp(!2l^$=!xS30&;2rQRLTOgAhN~7 zuqdeMI-%Te90YKF=K(`&7buex%BA|P19zU6tKG@^;Qm!|&we!rIbcaIYuX~&ZBf=% z_NRT47lJ~^hes4}_Wr4c@4-be@C)^|j@>Q_Vi_4T8`^|m=<7Dyi@d^KlAnp{phG{K2B zujd`e@Z&@4AA|+XaxFoJ|IGc`m?MDLsyU0Ybnv=}6ui=@5E4M9Ix@RPdHGP*$&kv4 zVpg=-G5+U}&?qFXnCU5Q>IL=$=kUyp-C(4ZVe@8X16cH!Jae%shrMgdDQB-& z!6l0>gGf7Zt?$J|*b%}jaG})?ByQ0*Xa+gJjP23euS-^gLTF=U>0vdHm2`KraZIku1#4#qMN{5eJWykbDni zA)ry`E)fHG{#h3aM|?J83mgvGv5j&o-rwNVg`X#`!F|J_trGq{`V}W*d&}|92^DR2`LL_&cRa)xue4AbY@%m>A%OzfjgI&=Zw_dh1|I8 z?!ryQQ#I%@<;49!1Fgf{ueRe*arkO?q8^8e6-T7kG!7LVuZ+37I8+|Wd^)M{pHL}K z&(B&gpyuJiB*BFXxu(>kg_`oe1)s-981FAFr=(5!@i>fZ5lPH!z`Dy)X|(M()?M}Z ztf>I3yTe{h_qoKzASL8O;uzLl!tb)Y60EzVdl?%vlR43B=wffNXns^t7dy|LDukx; z9((seid z$*oh122@e;Hk*4pvF;kCR+e4Ex+~|T^}Q48u7cIvwHH`-XJ#Jn%fq@WcQL9q0_$#v zjSQ&=>uy>@+oc4oyL_e>_7q{=eSC~GD}r_RfsVs@wvaT?;7!xrjdeG^-E2z^*4>H6 zBvgZSm*)D5B7NUsH$L0W3VLAOjdGL=I4i6SWnAmS&tTn+IyS(66YK7f_r|u3Sa+S0 zxa&i(?mDe$-B^fqm!F?~TAcXA)-{fo^q#W`p*@z|e zfX%(C*$W?F7ir&t9`Q~{4VFIQhb8r_>c6l=QsX{5KjOxcTHd6u;f*CV z(9{}fVM&d2d(ZkXmkW`?oZZ@$7a@BSvzQPpsqY32^Q)H%BAY3vOWt%z73>MH=faX2 zm%m&7$pb#r*2A~BH<=CDe@?h(k0tfUz*!juEUDF-_&$1LN%i1~Fweu1I{Mj3uMkV> zdLe-lnX@$bvO>Q^Q8^33A9{*5V@Z`d$TwPnC6&`L|0)NT)Y%n>15aT|y%Oegl?2~r}8{S8*N*HkoIxmHpvsQzgN-yDk%jsDg#N-g1{orl%0<9n-$i0Ii=M`*@YUZfEH$>GST0CKxP$}dq9MW1eH9=JR#g#xsM zdoddSfb$+#cR-1b5r z?*Yqs%N;mFug|wVWvZ`<4xGsE+F_uCIvmZc%(!2`T1sR}m^cwRrp^y}vazE4vdjt} z*{@I~)sd8qx}ig_&@b%HE09~ebA#=}djQ!524A&efW|&ZuoMpin(&|2wXTPu!{o~R z6^X+@3tbzjA>jj212aOOs_YQZ)^62ln{a-5v%VMLRehcwIrkHToD~ zO(=AboF%QI0-(maBDDLY4|qCkw;kv40AjT6V&QDN9Wb}`S0s6vJqS0}(yX^yLy)k2 z)Vkvqa5Zyu&(kleApA|D{lgVG7&6E;E+6Iq7WTxXkLJ}D+Fizws@iMuzbu4O@SJ&; zP_|VvFE1w2;;j^p7l9@>-ZhQ6%O;mR$+miyl^#_Rw_?BI3nw7$EZW%dw#bCk^fjY@ z(>;lQvW}DfYL1lPS5gJlFBW)6LBOlhg@j%Ebp~TeV^_#UQW{_T8_*8q1=5*|%Fj){ zbmo#-IigNyE))Jy19avxc zdYuG`TDtXRyp2|L>u(vEI@e8j{SgFJ*Mf291<)Q zH1xzAdu^uK!I^XTE*F^?q?wF4cJ)vdzYgZuC9~f?nPQHucF0lU4d&SXyLd>Am}95$ zrS5-*_nu<9#%h`Qh(U@aPe%3yE6mx7_ zOXuY0m}6(})c^bdbL?8}>NZo%vHgEt#rptr?A2#gY6+NQ-|qb~REjya_qdnsWz4bH z*vqsUVU8`*TR8C)bL{q4!Du_?*ekfjhdnUIzHrZ^)c8LfJFK~J2;U(>$(u8iyuMp#~@5G3e274YOzMjPicQK zv(iH};rMO4=@i4FN}c->rWh6SCIVu6HPHlb?Hfx>F$8WWS=VBUQ8%*}#xcd{zhP;- z9#af8v$$eQOfk|VbD#8Kiow!H{U+aj3uyOm*IQO*0E?lSxBEtKfI~SoTQJ2)i5KZp z!xV!Sl@}p|DMlHqd36n@7-h>AzuHcx82!8=r%N%#;JPawRE{Zzw);fWc1$s7Hd4oX>~D9{@(X0ku}zhviD5=1Eclgi9w(|FRDv1jRj=rCGRCY2c#^gkh}!+f}J z^P~o|FZ+}QqlpgAslnX1$#Y7H*_1S8#4J=zX)_BmQ>M&<>9i5!DSyL)N#;{$rWc+n zn$~ADF^KQPLO(jcr5MdHAENVD=0m4vlo)+8%h=~rXGJKDGYf;v%mAZ-j$x@Qxv*6( zv(n6|BWCp&O}p>8Sy6}S2EUzEWYiox;W$EF$&PE5%rURP4AyEm$DEc$>o5%{o7ZC&CgxW%3g=c3m_r=@Vm`*0;an1d z`6A~E2_}rumrmXpWmrq)ClHwC(vcI@lt4g?L8ws#f&ruJx*HR76pLvQ>N^5KhEYR@ zCDbVbVKuWSVkMX}OEI{_=s~~P(^s1MiIu=yYWhu-;W?^28$p`UiGI7JTCx3|Fo|tp zLi(ME*|h_AWsXcB67-oR6C#1R`qV%o0ho39M1nE1G)N?vFiR`g3CyTFk*B$aVSgMu zL5A6Oot?mpGWw;5`iq^w3@oZD$AbOgdyhY7EFCqDgTP$7{!M#~nI0?WSnycgz@t6P zUPR7?UhmHxTFtCNxQaQO7bk%kg7nJ{^%5t6#CUyrhR#@1>MSQggHc7#E$MHr%lLKM zp!NV``KU2m3+qR{#YJF7B0XtjHu7-)ZPepln24V5&>PK|2HfZV`+~o?7harxL{s&6 z{?-Ta5GOLa-2q+#GivE4Df7WK zd<$z%jQbgmcNUe*M_|Uh=fNME^or+yE9j>s-7TvaKA_7xBak*PS_q_RgZ+I>51(JO o;Plgi>ny1xA`Z8vMFbJ%d;D7T_l>mq|Avk~|KC8p!B1H8e`lI9ZvX%Q From 74d52648025fa60dc672f2959ba93d748eec2fef Mon Sep 17 00:00:00 2001 From: jaimerzp Date: Wed, 17 Apr 2024 17:51:32 +0100 Subject: [PATCH 04/10] z_max --- src/core.jl | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/core.jl b/src/core.jl index 30fb9c2..30bc74d 100755 --- a/src/core.jl +++ b/src/core.jl @@ -73,7 +73,8 @@ Settings(;kwargs...) = begin nk = get(kwargs, :nk, 500) nℓ = get(kwargs, :nℓ, 300) - xs = LinRange(0, log(1+1100), nz) + z_max = get(kwargs, :z_max, 1100) + xs = LinRange(0, log(1+z_max), nz) zs = @.(exp(xs) - 1) zs_chi = 10 .^ Vector(LinRange(-3, log10(1100), nz_chi)) zs_t = range(0.00001, stop=3.0, length=nz_t) From eff815ba30d8f851de43557090ceb500ff41a586 Mon Sep 17 00:00:00 2001 From: jaimerzp Date: Thu, 18 Apr 2024 09:39:19 +0100 Subject: [PATCH 05/10] z_max=3.0 --- src/core.jl | 2 +- test/runtests.jl | 11 ++++++----- 2 files changed, 7 insertions(+), 6 deletions(-) diff --git a/src/core.jl b/src/core.jl index 30bc74d..b4f4068 100755 --- a/src/core.jl +++ b/src/core.jl @@ -73,7 +73,7 @@ Settings(;kwargs...) = begin nk = get(kwargs, :nk, 500) nℓ = get(kwargs, :nℓ, 300) - z_max = get(kwargs, :z_max, 1100) + z_max = get(kwargs, :z_max, 3) xs = LinRange(0, log(1+z_max), nz) zs = @.(exp(xs) - 1) zs_chi = 10 .^ Vector(LinRange(-3, log10(1100), nz_chi)) diff --git a/test/runtests.jl b/test/runtests.jl index ff33689..209dc16 100755 --- a/test/runtests.jl +++ b/test/runtests.jl @@ -26,14 +26,15 @@ end test_output = Dict{String}{Vector}() if test_main - cosmo_EisHu = Cosmology(nz=1000, nz_t=1000, nz_pk=1000, nk=1000, tk_mode=:EisHu) + cosmo_EisHu = Cosmology(z_max=1100, nz=1000, nz_t=1000, nz_pk=1000, nk=1000, tk_mode=:EisHu) cosmo_emul = Cosmology(Ωm=(0.12+0.022)/0.75^2, Ωb=0.022/0.75^2, h=0.75, ns=1.0, σ8=0.81, - nz=1000, nz_t=1000, nz_pk=1000, nk=1000, tk_mode=:EmuPk) + z_max=1100, nz=1000, nz_t=1000, nz_pk=1000, nk=1000, tk_mode=:EmuPk) cosmo_emul_As = Cosmology(Ωm=0.27, Ωb=0.046, h=0.7, ns=1.0, As=2.097e-9, - nz=1000, nz_t=1000, nz_pk=1000, nk=1000, tk_mode=:EmuPk) - cosmo_EisHu_nonlin = Cosmology(nk=1000, nz=1000, nz_t=1000, nz_pk=1000, tk_mode=:EisHu, pk_mode=:Halofit) + z_max=1100, nz=1000, nz_t=1000, nz_pk=1000, nk=1000, tk_mode=:EmuPk) + cosmo_EisHu_nonlin = Cosmology(nk=1000, nz=1000, nz_t=1000, nz_pk=1000, + z_max=1100, tk_mode=:EisHu, pk_mode=:Halofit) cosmo_emul_nonlin = Cosmology(Ωm=(0.12+0.022)/0.75^2, Ωb=0.022/0.75^2, h=0.75, ns=1.0, σ8=0.81, - nk=1000, nz=1000, nz_t=1000, nz_pk=1000, tk_mode=:EmuPk, pk_mode=:Halofit) + z_max=1100, nk=1000, nz=1000, nz_t=1000, nz_pk=1000, tk_mode=:EmuPk, pk_mode=:Halofit) @testset "Main tests" begin @testset "CreateCosmo" begin From 7fdf6ce6f38ddaa322be6769c3c8c251e61a2cff Mon Sep 17 00:00:00 2001 From: jaimerzp Date: Thu, 18 Apr 2024 09:52:56 +0100 Subject: [PATCH 06/10] bug --- test/runtests.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/test/runtests.jl b/test/runtests.jl index 209dc16..16d83e5 100755 --- a/test/runtests.jl +++ b/test/runtests.jl @@ -617,7 +617,7 @@ if test_main function Cl_kk(p::T)::Array{T,1} where T<:Real cosmo = LimberJack.Cosmology(Ωm=p, tk_mode=:EisHu, Pk_mode=:Halofit, - nz=700, nz_t=700, nz_pk=700) + z_max=1100.0, nz=700, nz_t=700, nz_pk=700) z = range(0., stop=2., length=256) nz = @. exp(-0.5*((z-0.5)/0.05)^2) tk = CMBLensingTracer(cosmo) From 60457fc229ead706f8336b37d3ca3bff2bd469f8 Mon Sep 17 00:00:00 2001 From: jaimerzp Date: Thu, 18 Apr 2024 14:59:23 +0100 Subject: [PATCH 07/10] less zs --- src/core.jl | 37 ++++++++++++------------------------- src/tracers.jl | 11 +++++------ test/test_output.npz | Bin 69358 -> 69358 bytes 3 files changed, 17 insertions(+), 31 deletions(-) diff --git a/src/core.jl b/src/core.jl index b4f4068..2a38c33 100755 --- a/src/core.jl +++ b/src/core.jl @@ -6,8 +6,6 @@ Constructor of settings structure constructor. Kwargs: - `nz::Int=300` : number of nodes in the general redshift array. -- `nz_chi::Int=1000` : number of nodes in the redshift array used to compute matter power spectrum grid. -- `nz_t::Int=350` : number of nodes in the general redshift array. - `nk::Int=500`: number of nodes in the k-scale array used to compute matter power spectrum grid. - `nℓ::Int=300`: number of nodes in the multipoles array. - `using_As::Bool=false`: `True` if using the `As` parameter. @@ -20,14 +18,11 @@ Returns: ``` mutable struct Settings nz::Int - nz_chi::Int - nz_t::Int nk::Int nℓ::Int xs zs - zs_t ks ℓs logk @@ -44,15 +39,11 @@ end """ mutable struct Settings nz::Int - nz_chi::Int - nz_t::Int nk::Int nℓ::Int xs zs - zs_chi - zs_t ks ℓs logk @@ -68,16 +59,12 @@ end Settings(;kwargs...) = begin nz = get(kwargs, :nz, 300) - nz_chi = get(kwargs, :nz_chi, 1000) - nz_t = get(kwargs, :nz_t, 350) nk = get(kwargs, :nk, 500) nℓ = get(kwargs, :nℓ, 300) - z_max = get(kwargs, :z_max, 3) + z_max = get(kwargs, :z_max, 3.0) xs = LinRange(0, log(1+z_max), nz) zs = @.(exp(xs) - 1) - zs_chi = 10 .^ Vector(LinRange(-3, log10(1100), nz_chi)) - zs_t = range(0.00001, stop=3.0, length=nz_t) logk = range(log(0.0001), stop=log(100.0), length=nk) ks = exp.(logk) dlogk = log(ks[2]/ks[1]) @@ -89,8 +76,8 @@ Settings(;kwargs...) = begin tk_mode = get(kwargs, :tk_mode, :EisHu) Dz_mode = get(kwargs, :Dz_mode, :RK2) pk_mode = get(kwargs, :pk_mode, :linear) - Settings(nz, nz_chi, nz_t, nk, nℓ, - xs, zs, zs_chi, zs_t, ks, ℓs, logk, dlogk, + Settings(nz, nk, nℓ, + xs, zs, ks, ℓs, logk, dlogk, using_As, cosmo_type, tk_mode, Dz_mode, pk_mode) end @@ -244,28 +231,28 @@ end Cosmology(cpar::CosmoPar, settings::Settings; kwargs...) = begin # Load settings cosmo_type = settings.cosmo_type - zs_chi, nz_chi = settings.zs_chi, settings.nz_chi zs, nz = settings.zs, settings.nz logk, nk = settings.logk, settings.nk ks = settings.ks dlogk = settings.dlogk pk0, pki = lin_Pk0(cpar, settings; kwargs...) # Compute redshift-distance relation - chis = zeros(cosmo_type, nz_chi) - ts = zeros(cosmo_type, nz_chi) - for i in 1:nz_chi - zz = zs_chi[i] + chis = zeros(cosmo_type, nz) + ts = zeros(cosmo_type, nz) + for i in 1:nz + zz = zs[i] chis[i] = quadgk(z -> 1.0/Ez(cpar, z), 0.0, zz, rtol=1E-5)[1] chis[i] *= CLIGHT_HMPC / cpar.h ts[i] = quadgk(z -> 1.0/((1+z)*Ez(cpar, z)), 0.0, zz, rtol=1E-5)[1] ts[i] *= cpar.h*(18.356164383561644*10^9) end # Distance to LSS - chi_LSS = chis[end] + chi_LSS = quadgk(z -> 1.0/Ez(cpar, z), 0.0, 1100., rtol=1E-5)[1] + chi_LSS *= CLIGHT_HMPC / cpar.h # OPT: tolerances, interpolation method - chii = linear_interpolation(zs_chi, Vector(chis), extrapolation_bc=Line()) - ti = linear_interpolation(zs_chi, Vector(ts), extrapolation_bc=Line()) - zi = linear_interpolation(chis, Vector(zs_chi), extrapolation_bc=Line()) + chii = linear_interpolation(zs, Vector(chis), extrapolation_bc=Line()) + ti = linear_interpolation(zs, Vector(ts), extrapolation_bc=Line()) + zi = linear_interpolation(chis, Vector(zs), extrapolation_bc=Line()) Dzs, Dzi, fs8zi = get_growth(cpar, settings; kwargs...) if settings.pk_mode == :linear diff --git a/src/tracers.jl b/src/tracers.jl index f17a41e..fffdbb8 100755 --- a/src/tracers.jl +++ b/src/tracers.jl @@ -27,9 +27,8 @@ Kwargs: Returns: - `NumberCountsTracer::NumberCountsTracer` : Number counts tracer structure. """ -NumberCountsTracer(cosmo::Cosmology, z_n, nz; b=1.0) = begin +NumberCountsTracer(cosmo::Cosmology, z_n, nz; b=1.0, res=1000) = begin nz_int = linear_interpolation(z_n, nz, extrapolation_bc=0) - res = cosmo.settings.nz_t z_w = range(0.00001, stop=z_n[end], length=res) nz_w = nz_int(z_w) nz_norm = integrate(z_w, nz_w, SimpsonEven()) @@ -62,10 +61,10 @@ struct WeakLensingTracer <: Tracer F::Function end -WeakLensingTracer(cosmo::Cosmology, z_n, nz; IA_params = [0.0, 0.0], m=0.0, kwargs...) = begin +WeakLensingTracer(cosmo::Cosmology, z_n, nz; + IA_params = [0.0, 0.0], m=0.0, res=350, kwargs...) = begin nz_int = linear_interpolation(z_n, nz, extrapolation_bc=0) cosmo_type = cosmo.settings.cosmo_type - res =cosmo.settings.nz_t z_w = range(0.00001, stop=z_n[end], length=res) dz_w = (z_w[end]-z_w[1])/res nz_w = nz_int(z_w) @@ -127,9 +126,9 @@ Arguments: Returns: - `CMBLensingTracer::CMBLensingTracer` : CMB lensing tracer structure. """ -CMBLensingTracer(cosmo::Cosmology) = begin +CMBLensingTracer(cosmo::Cosmology; res=350) = begin # chi array - chis = range(0.0, stop=cosmo.chi_max, length=cosmo.settings.nz_t) + chis = range(0.0, stop=cosmo.chi_max, length=res) zs = cosmo.z_of_chi(chis) # Prefactor H0 = cosmo.cpar.h/CLIGHT_HMPC diff --git a/test/test_output.npz b/test/test_output.npz index bcf2672700ab5d2e164025496fa90d838dd1ae46..552a811a146456fecee84482e383ee25b19a108b 100644 GIT binary patch delta 45567 zcmZU)2T&D3vjz%+NRXU!&N;_3I|mR%RFaZGL~>9hBgz305CH`aL82l-2_iXXkSI}d z&RHa*g7~=izN%OM|F&vtw|2I>w!fKgrl+@iCk=lm4WC(88xNlb=YNiTuUO+LJoS-i zTpXPLIq9+FznQS>S9CGtQHv&My=7RF;A#(yj=A}JGfo^>vK{8H9i^l#W3#9L1LJyORifI>fpg-r{ zZR}v>K|4^tTpqt9fTpNcF*ZFHMUQe|g3jqA(IS6hv`+aXFlZv{F2Z6R5pE z$OuS^A+**_k;{*Ta9NgY&?OkU(WVh)^CdzG>%8J>Yy|AK-tE&24}_M4W1i;Y9`M{* z>j9SW&O@kQ@Xtu%=pGDh|7O|iuLlJhuRk^HP=$ZYX{5>MIN*cNfyvVB+|c*{7*6N&qT~v(E%px>q{mW;kRY~rnKjv6zyy$QQR_IR?Q$`zPbg= z*mHVsx9mY)mj+QOtrPeXE`GB0I)`^cS3?)a@X)fM*bb9YBD93}I`>up1^W93pK7gj z8uW~WJy8o06S^kjJkGL+9sTttxm#>C7a9#r(71Pa(Jjm}(U+w7(X*~G;ScNiG3XQD zo^|XsJ~S@t)%!l0Txg@5cXb)n4CrVv_U9HNBxur(;EYLzzfg@J`Q_!3UofjgFp*P# z1s46ge?z{g2S$65soq{phd1BW3DD|(fH=H?8(C)Fun0fBx_II#tTN=`O2#~c4DD}T zie2}D91EHu6x+_w>D}mM!u>l|@RO;K5_7Bp?8bJJ?NwfbulRG`=vv9cU!Fg{^oIz* z?%eW$4RJ;Yf6_d^;z$HzM?+g4Ait6Qs=9rJmi7PY6MMbj57^=wYV1vXTTHdrr~Xp> zWw?+TFMRvfCcIqGk;$mB2fJ4H2Fmk}VM>S|ekFSpTs!2q^ zU{WG2HSuJ4`>`*X)ALxUUT-nGH4_Ro)d+l8NuNRyF%=6NWgR#8KGVACE2%Y1Q7v8d zbT@$RdzuDcb=6^F=MN9|SUE^2+t>6Qmlq1LacQ;8(LiFGQk$bNd}z@AbH38~;6K4F zJ-^(9OZYz{f6;UT2m88_5j%Iq1f#n&dsr+s2k{KWa~57K!>wkT?;Ty6kd~R*ljh?d z6w>H2#^pYNq6TWUq|eTwkn9_AV?2Dcwm8dJ@fZ=>_XSJwq7(&MfjUTWDV_$ctUr0u zlFNwha*o}oIA%i+^7G!i8Nh`ecTVrFqTxq#h`Qb9RTn~sTl$CJ|0RGyU&=mlTOs5_ z?|9Ahva(>^)8o#b;n-9#Q!*;+(KyL(+GN) z)$(1=f}9XuHxD%E|8XJe?rgtaEdkouc5M`?&Bw3@BA$toT15a0T=5v45Bzs*`#()jY}1Zpd%;#Esnh*dX({Dx-%p zwm~w&(sYRKKKN2|`mQ|l4)_l}#^*men0+tS%Ktg3aIm^(G}tC)T_w%mi8@7kN}!wQ z%K+_bGq7#(b9@g+4gAg4Q4g-512OSs@708-{tHz$YMVl=l>%$2J@{C|Dj@L_N(w&p^d zYU<7xOArP;7jl++$X|mL;om=~@}i&nBqVWQ1Q0DgN!&BxXkBLGFj^Vr#Kd$BmRFmo|o6Bfz9A} ztTb5fAmtq#1h~^Ae+}h?fjc$N;u{x(04?QWKPIX}f{w-q}=KnUr_@;x!`2X|?w)^YdQWEYP z|JN@hC)Pb!bSNX%o97raH)u4b0IAT{i%2+!<0W*ZzgXB8a$NMZYPkg5(a&%qw@Zp&<{>SNHR`P%64I zQcEBT%G5TNcU%vI1QSlWIiDWG%*B%2nJjB~Yf-33MaF;_hP@owaQ>za{jf*8g8vJ7 z^wRblEFa%PC1S_J-IW2}#^)2IjffLR<766qMAER^V4g=s728w+x_PqmWAfho7vzg| zsJV-^j|uDYj|^L2s*Ba&e~YQ-=~dLaVES!7-D`i};-NX;`At8(V0xzU-l@)=73h*K z=v8~U1nL;B#~fI9K(6fwc<)~c3<%KYGY+kUoejNIn3~V<`gVNi%@3JSEyZr8;DYIN z%k+9v!F1j?yOvQGOb=9P5xsT6^pMX<_&EQ=bnLCl5u6LAPv;t#T)SYp zGtT&c*ag!&#{1@WFPP4=_=W$+|BvZd0W9ACV)8?i8V@TbVkgC6elvsg!pVAnIJ32P z;beWf%RWJV;bfKBCIagJpOba{0j&_WUE~2qh*mRePFDj)-p{5Fkjb^}gVYz`CAX~Tv+f;9qxn{fXLNzSzTEjUOJ{r8v4ZAjf#LR9=- z2a>z~ARc|B3z4pj#Y|c~$kZxEc&SAXcHag#(q8&7{Z-vbBAE~LB!`n#`@{`6E5MJ4iDhsrPlk{oNCCcu>{ilD89{O3qGL?@K)a>`+ z{&8`oM$dgnDs=b0b^yj21`hx=JB0_3ZY$+g#qS4DHp8q=u*C*qwe5`~V{M^%dR~{9 zvmG=XBF?k9X%D|e@-OhoIY8q4m!;(bj?fFcl|Rel1hczm81aOiq3*9bF#%;4SoXqP z>&abLh<__H==YO{ke_OeiN4?w3@S`*6!3I|@5pDWJTMH8;XMw6FG%ZSsL!CRW$5V+ zZ~bG2Y9b!+x8pA%-z^WA)GgQ=TkHvIjI+M@V7%aD4LFp!>kS*NQ>WIXd|-*0o2((D zFQl4@#S#Tt-=o9|}`XN$3GvoNQmr(46O@?wr zIDGmMlBov0fSgjFdjx*y&~L6WeK7V0PEEYb zuX`R1iB|)*|FFftr~VwXUQIDjNJ@jL^rz5Gg0KG>ExO7i z!^?KvK6UIVP@lMM;SXghEUGqq!AFt?D{A#+14z=L#z5>A1@${9RzAXiOz|EHTw%FK zmiiuMEG}e^3S~gz7LFBz`V1&0f#>kb-~)Vbz%>8t#|P*ut?DEhkqHeK1>-)7WWm!8 z59`IjEa-$GxW-laBpcpS`|PX5p97B{7aYGB%YlTq1ub+UKEh4QI*St3Pq6O6V(;s} zpWvamgWq6nE>wDG$n^AC9@MX0owZlZhh|s((;pKTK+<1-XyQf+;QqPjQc>n-Xw~+Z zI@hxhs%cuHEv|ooTQq%hhU7(%(q+t|-oFT9yvUoy^*4&(wF+^XW20hNaT?tEzM~k5 zl8qQ@D3rjmd-953b4#EvVX`baODUAnTRWzCQwq5Su(0Y<8El&D$lUTSgE&>$jqkV0 zAQPA3zXGdrn2wiMPB~f*|L)lStkkZ63Fw9R@a76A#3ktJC|?Ov1x5mIeXhj7{I%;T z+}u^LRpvx|;%ycDB~YDiPF4-??e$Bh1y{p?`&Z~fepSO^V!9%H*BYo*bWBRMPy<_K z2bkaAsfFIAJj59twQ%?|XWo%w9V~AJemVJd@Y#fPdgQ7-g{Uhk-JZTdYWVmJAs@eq6nz*gp zK26Zp=(gKwUK2dSES9RzHbIfU+60#EoArA?)X1FDxeZSYU84A6&svOF0hJB?z z5vrrj7k)i;CNn_`j7zF|VJ_POkG-Q#@N8P(?YoTsC@`@ta6J}`J#TD*&0@~G#@j9M zgB1OIEJG_CHd=lwq23DVOwD&CAGN|iYl~EgNv-hPM(Q=I)>g#ZBO# zZG(ey#nTL^HpsaiS?O)l27iIGr-~76kTge5!KSzkx(rqRJQ;0+8qYr09i6s8JO#5Z zj3vj{i+*JFkyQT*b4_J(-R!?YY~APr+pDkefpx(PRKZu+m~!P^Pw!V4SZaULz4H|k z(619AR{e9^;fpu@l5Oqn z5XY9_dE`<%>?30%{{F8WhV|dV_@pp*z+ENl-F3MR2#sI2ROxj=Qq- z7t#TL(|hlkyz794qdj!GRUJ^?T*YZ;Jm6-T_r*@ODCRI-$t5OX4;3ozR>8 zwR*5vC;W=`Hl0`RgvqoL+9Rf&5G6Mx{J^afextwr;#Ei|oL9vt+vp^B!jJ>Q#4m-N z(7Bjpu&AXI?z>a`)gA4Gsv}2F;x;;=J8sdhCug1TqdD&!0eKgUHWQ9jrGKVf7hD#&MJ?jf1rthsM%jjSL4E4Xid%_YaI%~G z%|=ca4u#~KJnr=yWvwSrAO+< z-4Nbl-DI)phB%^2`c;p*A)ixb5Tk!L#60uS1TR9rJoN?SJ!>s!as?CXXK^Oa7GW8IK?*52UCLN}x{ zS$NC3*$q|fY?kp4x}k6NiDlhsH;hIbKhePNfw2f)(ub=(Fegv{-4IO=2C6_&%qU9_ z9Pe!W@P(%bPKdw=S zr(FHjhYJNu4Dmi6d!P)1!<1%054?f9YXoxhugDI1s{tN( zueqvU*Ms{faTZ*bO+fQHC@E0Tu7grzpjH0G~^hz0#L&p!?*?jS!JOFs+m1wvy8ihHX1zh29K+hcEvb zxM~f8rIYz(`Ab89_Sgj*QauE^x^Myw9(@OY3DoL7aSsE*9_mB?(P6OLK(W&uHUcbS zV|wG{M*)^}Uu|w>6i~YQxxI@W0~=AjlEWy>IN16eryIC74zk3XzGl9e06C$7IC656 zU`24LOLl4!JT(4CoEIdBBz_`k}9F9^^Uwu_?!pE&#&m z8za4U7r;HH)F+f_3%~@2Wd^gf0MdpqR4jsvKp?ktbH#NLq_xY_j2ABgzU;HEy~9P& zMxJycFTVtaM7o+j`YwS0xeQKJ^%6+E6nA#~YYE5!^PVaBWl)vNowDk+44%!u;-V>C z20If~o)3O5gYVt1`!Q2OE8t<90i%!Y3ZOHYXOw!k0z`xdY_5*40K6L+Wuz3VAa|?s zlDOt7u+j{3bA7f7Ze1U-?krgawV(NduCJ|vH=HFED@<#E)KRY@?baGlAS;M`gjoZk z%8NH_iq=5>tO%*1kYL!!FmRQ$+XN=%DQr!ln;^MS+`=Mf6Ih>M z=7~tZZUU!|TXEkPH^KMnf*T+IZGzA!h-qvdO>bGry>JJKgAI7)9XVDO+)2%JQrB-Wh zd%gvzvc7zoBHsqBfd^D39NVA_LomZTBfbqtILgmmz&04&aZ4q+xeXK+wLX3^-UbrH zS?`T)wn4&;KW0}RZ3A->_TCQPZBQTbnlLhG8>qXk=Ng1>gJfX>4YufQ@S}h^=s12G zunzv17*5^>vh9&1<>}jiLd?51A!8fRP!|YbGPl8F23YerYa3u@EuYrdT->zZW9T-z zXh?6OOSp0IqLAp4oA!mEFy>8-_JtyM1@*cc7m9lt{W-3`*#-is8ROb7x4{!lIk7u| z+kmka^V7s<8>F_XjhMS`1CDl@Z42vdP#}S;;b^!G#J+j2xoK?!?_ie~p32(*W#iE8 zA+n7Dkw?9ceOR_ZSP;tmDakg-+*o9M`g;p_{9=FRy}AYB_%!TXhPJ@!8zE(<`Yk{d zxc>M4$1O0B%;RGay#=z%zMBBwEkHK!CoFVt3$$Jn1$5Ug-p>?ZI2GIi?1$%Bcw5&G9&duSNA8N{I-5YVV~o^K zWD^vNj~B`iZi0Eme75z)4R9Y-fyLHsfNNSN8ZdqX;HY_5Z@6rLEoNH#2(=BMOy2N< zm+@lIid_||x){VdKTYWMn%BVw6S{N>yAH6e^gWMk*MXgrKJGC_b{%Ae9Xx$GFa4cy7mW~Rbl18xI{ZWv?z-0>Zr9N1Bh8z(N@I z?4QsQD1Y*JdSzu1=td&c|K2SE?&r(o1}2Lj#dR)VhjI}NT!|I#?^*!!U9?nZfeXOM zE6w<+`~rBubQNpzV;<1RydlW=Fb|Yi%rRnjjOKx5K4p?W$vjA?NoGI?%z>HG;@$M+ zS+HT<%Po3i7Fheu6CM`MfJBq*%f!?(plYkgHq?6>+`PXk95OKlcy|9j=T(^k{^OT( zywfLvy`1`aEY2i2t8560vz`DamZ86r+QvaxPP9S?-#G9qTrUz29|JdEW9IySZjFMM z9^70!I-}sv>-_v5B_qIVurtn`egvGkn-F$B9R_MXQxC5$eFrNk4iT0&z5`~BsN2&e zL*SB$I^}1UAu!XwQoI&62&A7&r{^9HfF|0AI`Rhtpof~4&*5Sbvlnh5@qg6^ay!n1 zgWY?<6S=0763cHu6*C_*rJ~aVew2rr-B#}ggAFD0yr2un`E2K}t8@Ya%I_vkumk*Z zh$O;l&sI3G|EvNd<1vxl z$(3O0DmJN6zXC9k_{Oi%l>=V!ZC|0aQc&2FE@jtR0(>G0%%A5M17bysD;DpHKs;*; zUM%VhAbRSC$7b*uoY?BU-4DnIx+8LmiA}lS-gcH{7X2rns5E+>9-ISIg|l^c@v<>M zlrKAh_-!V*%KYQYE#w0@iecE1`TZV{|F(&=d-V<+7J?x3wKNdW@pG<_GX+ctMT|JG zB?04X)*vyZM1V&MGKW3lf#J;tf97s1Al9xFFHnvJIg|q?+8xnAUaMfhdh9hA{qj)z zkxnFVv95TCISdC@G?#C&w_(D-jrz;=r?nvfC(~|ubRrlO@+;fDV15pWb}Ry~VgiA8 z^0(p%`~cu;Y&a^d=?`WEZtU7M`vUuqqvLgtyusMJuavX5J%M4vP4drX?tt#a_@B90 zH{d@7@_H9tfsDuW^2oXqfCk=XHc}42lH;@X+ZtQYMToVpz&w2be*dBWNZ@%NTzQgJ z8mn**;IgxQTU#^(bc|~oxpT(g9gcC5vVbAb+9+K(E;j(KnTZ@l&-H+Nc0<@ZN*!RX zv{Jzzstv@^g9p(j*MV7<4U1XuwF}x}-jkfFf`HZYeadPGK5bu-WpqV>SB=L)YLkkf zEGl;r6Hq7z#(!;j*vv_Tzhbcl^baI~?)idY>3uPv+?J&u5+rmng^JYa(({9L)w`)B zbUfgEj4@@{gA*k15|kcmvVr^?ujO2onZP6JqU^mQI?zh|Z`3V`8qig!1^!&4062l{ zGK1`-AoXY~c>F36AkNPVYY4s!4h1MM>!u!;E=KiW)T5+-sPjQZrk=~EsFx0MrH0fe zD6+UrHipj!DD*KcQrWbN`f;WpO7V0XWf5wCD2c73o(S;97h;xCbjIo`lBNr&XxdWT z=({s0X0m!7{FZT)0_Qk2&CD>W&A}|#`t=}cE4XgmDySD_zZPu0(cFn5RmY&K7#rG9 zqpf#W5YuLqM7c=q*;yURumiX7rDzo@#(365jk^qG^*H!@Kw}ZgNbAyv>4AI{3Hlq~ z4T_Jb`cgfuZTd`9zM0r1m6&&^@Mnit7|4@R3`M4)vrpntG>x(bjh|ysgzq7pBv}*+ z(Yv2HFdByX6N|k+6c&t1!5rLk$3McLG7b);WSF0zo;^;ex~}bol8+I}2se6+GS`dY zE|hjbb-Qr=dB0+ZicU?D3y!u%&40FUAS1U#d86a%WNsRxK42ZSgi;MqVkBxyHpx0D z3X1isEUsE8haYv9#$DA>F>C_kEe{~7WFMdA2(5&A)>1{wj3bM}Fnr(Noo+~`@ol#qc?kZfRe7C18Z$ zc~8oc%AvcnYIy1_MumlS!>Dj-OvOC^_HldkkP7GP5w6+sZz_d+_Q{j?+EtKBNs>8( zW|i>Yv5(z)Yg9tTuE3U(GL`6_HWKNrFDfTlIA&6@pH$8d_^|4WA5^B@c=k7br>ant z{?ZqcNl;-Tl;&6Y8?8d6q8VX-6t3dL7pb_O;@1O}v6b-jX>m)Hex|+`*mgq|>*tzP?fiNw?-Ri0-Ciw~ zzq9*3c=y#*jCrqYQ&OT-ruSYF^tH;WEG~tr8eWr7AoP3869Isz!}CFNGe1CO_K)5>)gKduPI4J+4g4ekD7>{IqDk{}xU{Z+YO z$SRtOph>xAWA%Ia-Ad(zt^FN=lOpA8sUyj28K0Eh3;)Jos3$U%IqL>_B8QTc`Cs-3 zzsrtQW(wjVbX=;^ zZydyB#gx?p8Zh5oOn8-xd{3iexY(3k8$)vvRA`mS4+m|KMN(zw!aVE8fAN*6KZRoY z9{yFL%E0^9`}R=jjr(6>k#E~dj0Pe-PE5;6daboj0tlv+)WM_A&4J&Q!lV{288h@K z=}IJAUFvF45_1KI(s(sW>Xy5U>1V}C%Vn%9qobdcI*Bk+pU&Acl<+(V&?m8pN+hUY zc{-_RrDGQA%6hp_B~wCbu^-L>O01~bHS0c4rLSJQM5fLzO3A-ba`g8fDDe!{f&yz( zCC4YXXy5(i^C+$V z(O4NXVph_|h?U3VTTv-hz6tXT&mmN*Thi2f7L23hn(@qYhT&K-GXB-ow$LAn*DW0p zP_m@h$R2U>S$R_N>7VAG!H)(N)w#s)T;b|aly;%tfBLaOQF|53GIPC5kvAgqadTgu zA|I9f?C6&a#rk%skF(NSH?GVEYMa&(pM-I;opD5C%{6Y#;9xMK~ zWOptLvr`oNe9za|-a-*iC6Y0n-#{_Jx#;7`+I7X^gxi_nf>5!3Q{b-kzKmjz{HkLP ztFU6BksL{HD5qlE@oQJ#cXWzd|FqDA{G^IH%Lmx`hD(YfTp4k)6u%T+P!fwRj_xR^ zb?ZF;Sigi($W2=yr$tRDi1^X7h_MVPY$p#1*t~30=%jw;+r(C*Ku=XA+#FY=kfF-? z?AN7ig-A(fC;8-Lg^(labo;?*h1yu6p#$X*1+RR+cH6cm3bcL&`zm8@3IR8+=7#Xv zD%35y{la2Q6;wu*`#&YN+EMo|GZlk!87S3*IP^wAO?1CIhD zwiNfqA%j9JA^R==A7ly-`-`Y&h@hFp*F#Dq(g_|q?t2wd?6P_+_t%EF(&wyeBGz3 zw?}%=^7PN6@v_|v`2CQo?M^epiim3(;rtFj&@0{Qch5c~%A(Qs^pPOBZ*G0F4nIBRNIVGwM_xF{&BqvH zo_;cyBWO&aYfHIaXdS@4hxl zIdh~KGb$H>FBeH;6kJ?-B71IV98O}eC7Zjmnq}ZUCo9AP~BwKTsZ8z|>O%|Uy zS?>BkrL2&Oz`O>?mmQwCy!Mpgo$Q@~kvR8?7})|;wK=8G3)$aCAKZ<9c*{z#wS1@} zbdqf-A6d(KVIliq&)}wenXc@k!I zGAbCWw0AT&EM%%AuOCoR>B@|zm)jbqp=B0b65=0SmzK$7Vu*fA!z(lL>%9ZlC3+dv zYpV--6htym|BxY8`BeJdIi1MOkRQ@&HlFPXXLHiu$#vHlbOxp8OYSdcmbXajkBTdj zNtH>{{DRt1t=ZD&i47t9_KDKNDpx3?*1g(N zx^gQu{)VltbhMhT%u_$K^h{km-{rSb(pkQ7UWN_a(wXIM%){$6(#d3LWnzp3(uSGI z2f0QkQcLdk0sZgSr6#`u=G}t{DU-*Q7Hv8`QuIMoVnfAsQo3_b@o?x0rM4?|w5z?| zVWd1F4W=C&qoo=(9`(kv21(K3j{9p4Kb9hi?GK%cw3bR|uVOXRFpzqdEZDbvt|nz* zAf7&6DkBvtBK>{t9DPv#ra~|D{R2hI{CZ!I4GN;H#TvbiHSD5EXW z`Z5=fwNpt#(J8gSK}Sf!WOhw?qm5DG5q;jbI8h=A;}Me+;0yjqyMdm#{`f& z+<7o2UMxLuPTbZmK0j;joybukzFS=O2IJEr_K^!AFV7j?70>)l~d19JSy_d_kB)W9`D z(WzK8D^_pCjN!egA$n=KwlYd|qZ|{?vT5fpI{kJ;gG1R-G?*{$4OYTP^x@L_9c$#8 zsO>Oq!yRX7(X;AO#JiSLl+kvzwL^qL)LvTWhQyPw%hE8J!JAaW$(`C;{Tw8)g&45q{UnaI7t8nbWqP9pCy0XIut zR~d;|C1j-&YF!g4xa%Boms(mR8)4dYpyU)uTN0AV^&=BmIObKNU_BM?WOJc^LbWce zHRqZuY%wC-R*4_kz1t$}FOoJP*j^+Y#<2NJ@pGDRX@`h1m}z8!kP!h_WlMG^D#!q&gbA`-{(g^!^ub-n1m zP;kmM#-Q9eq2e;zRK3t1p)h^fgWr>tLZ19FUk_8WgqV7aN0xhHgzA!>9)Bqb5c(c{ zrnu_vBqTRfZ)w|VC{$HdIvf_GCUm@cf`0i;La3kFbLeq5i_p%t@W+RVgczaPu4#OY z{X;=vErm|Z_XWWwT!Q8xoj$>tLllaop;|CSWo7PsG+Qt->CMzyP^{o35K?lkj}hEU z;?|z?auO`jPbd-BG!#6(uEB3wrYfl6pkk16UtEw|rk+~+CX--aDSyBpJOV+vCNTmT z$9;jca^tt;EwchEB+9;#pD0xvLeqPKmh1Vlo-%O2^S^9yD3 zy!tV)&Oi4PZykdg=1;Tto2?FL;*YktPyBrRsQ*iGMJZGSOf%CyYg~KA@dI+#m7zbmk?yxJf%4#&WNeau+f; zp!X_Mx$m7u^)6+HaBDeXBZROOyK=RZh`5$3+!JB_~@ zMZ-*L~tQw z*USY1J-Hlm=mJY_SaSK8WaJ8JX>fVVe`}04li(u$pc_tdPS16YIhA1D#o}^ZyMDV? zqIQF`$<>qHZ(xX%^AqcU*wb3h0>+WUSEN~-|G+=jdHE<#{M!T`X(%sFemOl3ojOa- z105W_Qb`TY8oCO8%gYj+4+QA{*)h{|8hYYemuKT}GLIbo4m@7t@O0RkeWoN)#+nF&F`TaN5SS&!tBi}l!rT^On;s}g;p_U zwU%AB;LwZhxb7GmDxGB4PP~cj^T|+}|GgZxLaO*0zsc8Zl(!SmL1-_wn8r%?d#V;} zl9`j2&Dhl09Nr4;jbsY5t(SE!@lCqQd>dYAuN8+6AcunOcV$#=U2u&y1lczdMVuoh78zUp4nW^K2n6>y)CVr6Ypyf5p| z$okV1kF;I_msK)g%W+U}jm5&muN}+R&(cCIsY)|m#!~ioe@1&Og(ao&A!8QyISXrE zNq}IR1Iyu!4Imq!%cAKJEl#{I%i`C=_<2g0m4#N~W7buSB_7KoyJsHLg&WMEN5_Z# z-vG1ltGCpn4=R|q`}beOyQVSo@!B-T>IXBIo5x>CGjU`NHnO_Eo~q0I&|492v`Lov zodteIjVBB9BcpfBX^EGZV|e|i3)R+{0-L>#2vhr+ULBwQK?jyGy=%83a|upn+BhN% zyU!WKGzu#MF(#~bOw^urZ^A=vGTGTQU)j1R#e@@IBEO$Q&y>G^v~cct&iIoi@+Hp8 zCB}&t9Bz^%U5rwGKk^IzePKKq4A5FIi({O3o6h>T=*Nf>{3%BF*ox7L8+{l%qt2K} zQ_`jADa5F3H>NS!M8SAfOw_Gk>WJa@SFLi25Wpe6Q?`eyR3b|aRz^gZ!DHEDIv=u=+}->JQ6P4BIDPCgQ-K_9Wy zGi;VEOfS=bNfEkzLP4JeN+((R4(X~)PMU-QCg>9G5|NcCG}76oHPP;|XVM+~^zB@H z8b)W2Os`OryVBXYJQy3W(x*e19ct5}73gpW_qr5kIq5RJ*egkyN$5E5ELmpn@6cYP zHvh-*pBd1bg?PC2*h!PeO7w_}Y34U$fc+$y$7QSqEIzA@4>=2hsuYs37srHPyQCN2 zGGpVv@^dXlUDN(Il{Vef5|gO$S3j(N7K$*4dCy5ML%r9ceeIza6#(iIo$u`Tp#QQ3 z!>!I^X!+n}{?q2aFq-)wCeHv5{gC?yXa65Uw3q*6A`Kn|dif%2;ng)7^c~aFGc^<= z`nv_;o#}2i^xs_^8GIBMnlihV<1&^9{eH8A$t{Bq&Ehojby@qOh~;B34GsAKFIt@| z4Cmn}C%PiULF3wcdUU7Hc!%Q@AzD1n@%Lu-Da0DCxz|Phg7Qc5^RWU;&?e<>8U3iFTl$v%fV&78##S$*M`e_w4) zEgnId&&e4%<2u%`<+-VIIp!{m-tyO=PSJ$y4uu&xrwUNZ#7L*hMg-m`y=SX;mkBn~ z2j;by5W^qa$pf^se-NW{KW4AO|0-fB)iSUAzf_4uywxjMLQ5y5D?RCK32{NdI>N|D zq$vzUJg5ylT?+z_u(hgTbA^mw6H~PO4f`9fE`(X<#{V-qoRv+Geunc3C=TtNM zH(^~k?v=RoJy>FJqlNUMGy_F8(0)*`ld_HieN1|5 zY1WSx-JhX4Vk*mwZe^Q#GPJ{iChRpmvr6Pa=YMFWq0|vTJG%FFYUc=}+m6w`@pxir zAA{fOn9Pepm?RQn3fc@2^cq*;SG0$)KTCzYPi8o%x#FQ? z{Jp|f(P-EfeDlNazc1h~p*u=zZGP|`<&4vhvpcSE!Hf>A(QgGkbN}{2v_5<%Gn4LP zpavuP43j$WF7~v8O(jdlJTN{ZF2ci%28vk-_~g^!!{H8;H%$1iLKx$Ho*_4^kChwd z|7BP}(!FkVXDcM(yY8fwqZ_hGUcRcB;)^t`gvZL-2O>Fv()fu%FOWWZK2hA1mq?v_ zL$GgYB(g~TTAk588d2i9jH338MT-6VPoe|k5UUr82r4Q88IY)t4oH5B(8f;{t9?vH zIzmRMjQ6e6qgA{2$eHE~6A9N0qx%*mgSo5t*Wa_@7IvX?VQb2)j4 zMz`dc?dN<%J1WkyptJyS!EKzKss4-{CeE$LG!-It2OE3%9bYhrDgPa+RJ$U?b+yT? zccTcQlxLi0b1g=cbh}D*c8if=&nyEgj}k;In<4wtaS3vrFWYS2uN0}>-qK+HTZ;U% zkWA+eEJI#5y|N0zD@Vf68-H3u%aLJ2n}HRg3S?qR!1hmM1tM>tjBiX)iCli*qZNm( zM64|GI=|3is*twQ9-rZ)Dul@};ZTIN8i}1pdLO*6MzTWj^ZmGMkU#nHkp($52o2zR ziz8Hv#F0}HOMk9KMqmB?sU}f}oNxH!2A9?$EPX-x-{tC&xua3krJ8yqlZuss9Myo> zg++`RH8mjPuhYY`)f$l@(v9fB_C{oPGP&&shNlT(c&d>ZZ`Xu0@V@)1{HY1S%zXYb zu+fBM5YUw05pG6eUt9K`yEG#^Om>l3h0RF7vN)0AZZlF>SY}oy*@E~4Jqfe-XhG(8 z*riF!Tac=m#;%Ox7KDRQXq!i&6`91hkLvepMHD$!R=sLl5vCTOLE^uy$V;q#Am%k_ zL;NVvTXKPI$koOPhuP*fBq-{Jz9-&Sq{N-L>zDdhgo}ZuDlPOYLh;-vRipDOV(nqN8+xiEq%#+^tuVP|1aQYXSQl9u>%suP*kx}MF<+J(e>fYVQw zUC3EivZLtxE`fRA7h>Bwv-N_z8?hFY%BHpLM&^0$rN!lRBex7k)il<-k(hD> z^+T`+AtGNwKX>jyNW&Bz3)B2t{swvByz5raN}U8(LxhIn0R?M{1<53l1iSy6q6GyiwO zmViFQiJZRdPE#KeeT@)z1h*eq!Tr+mOsyYbY$ibohV&y9V=PiF?fnQpTe~9<{s6+H zS`+?*X8=KCc9t63uMHp%P3C(N>;{l;;r=-mVFO5*qq8{krvZc>?^yn%a{#d%ujlO9 z7(nkRa5YbPw zvefGzM4t4s9f5;8G4DQm@31yCsJZ%%{sFp%MB+h;WHmmYK&8GS1W!BONz{ z+&|=&zf?Yi%%vuHbc_rkIv69K(BmP5?H7aA1N!gC;*YPVathy(OTqzCxF+9`KO!{@ zEq>pTh)VUmH_6|T(I=J_4z=GA`Px_Z3e(?_=$@*~^Rw>=`I@*S&5`vZ^${dmSIEB8b_6+yV(EP!I)c2m`kDCX;|Ox253fh2a|9t0aDKJBK7u&$ z{JGUfG>VYF2y4w08bt_RJJvqZ#*88lSe-bfTt<;)cWIsD$Wi3ELge<)=TW4Ajpro0 zZxoSIV%C53a}@bVLQO7CF@^|98*7|Nj3G?CC#9o$W5^}`pAH}0$B?)5F)Lg@^-1Kg<2V^ZYCDO%eJt9z5Hg9_Xk-Z#gp1QC)pJji{cbFBip5BW&fJRgX)j zk(|>;M0$7{X|Lher8=5MY_o$B*Jx)Dr%l<*rE)U}k2)&E(`W{9%yd4L^qoPhV~cK6 zCCwlOuDuW)b=Z zjI(oPz%1f-*DKQh-7EtB5Vy-V%_6)+wkPBZvq(+CQrQOX9J1qKA79NihqyhQ48^F; zAzZuPe3fnH5SJfp*GWU>5DF?W{f+E7LK%D2tG+=lbW(Ae9BAiH@4kK+5seeZe1 zO^5T5W!O9t^^{NjYvMej?>WZC|7jkPW7IeauDVd*nH9F%IgiY!IpAAP%p>w$q9N@Y z^N8o7;w7O!^T^9@x1YQuUO*la6nOk#T0q*&Mh47<7Z67`KEV#<1q7qjBI+Y@a{(E{ z73hmJUqGyAoW(a>77(tSvLw@I3rNB?_ep!y0`l)P$v`Z90U17btd1^RK)et)&Ygw@ zB&z5eyHWoFA~6|O{dH~ui8Q)vD!#jb)J35FVsRD`)EoPh1ByiixyQ|+&$)=a3)quv zm0U!eMF=&-RTnYH{AT&x1pP$>ugSsm*F`B^fvKFyU5`a1kn$qz8oY>1FXDa^jax+S zJhAsn%vwa2!dB?d$`%ni?!#)cuZzfEZx!brHd-BxFoHUPRRUAHF&# zSVEd|pW@!9UqWsiaQijzFCiR}B=>j~mXMKDuerGEm?b27`Kj4I<0VAlUe<(-;}Rm= zL3HP<-x9LZyxPX`Y6*Gh66O(|yo5;TWs9EYEg}2mvnzJBOGpXBtAy^JB_#BiK}m3W z3DG&hJ`UerLPWltr2aWuLX^YKTI@)dFQ#O%*1dz=x?eUFxrRIfotnZRWvR`e!z@EM7*g z{~;du*1U|ouj|s58(Kz!(v`j@FJ81OR3d-&ml1ll){{GUD+tw5XG{n63KDbkFP$vU z3i9`wUQ(LO3SyOg&V+ky1(EPT5;6BKaz`vloWmMySCCB#hMN-JD@dM_^VPJ_6+{f% z-;JNRf?WJ>e{TPA1<^8h6aQNI|FCuE@lT;56P00eM%rzHqdFVZ#T(ZR=a$-H!OQ zn-VP`>;Lx}3#A33B0Tpy?`{F_6LO4=-ww9G`jw593eGJc9QNdLkXH-rn%wL<>fZu8 zpRE;kjA{X`<(h@xQd^*;_0a*x|lLg*%RIps?tA#_*XoaCvwx!O^b`X12&hR^4d>)>@OBkkSSdt`9Cf z%wcwXr2YhJF|*;17G!4iP^uuV%|{2kz3MrS9eJ@E~m4Y|Wo`&^2o4CI;G}tzGZ?J=PAm z{jfuRl)D4qEc6`}>3}!QosWOXcYwH#f`ht72V8mV#Gho`0ZMkJ#cYQ=7!dr)+{fk= zv!QitROyutIGNpHN`!ZS$*pv*l!Oi#`eXTl>v;$4UU|;l^PUggkm0j@Q*rlz-8@phGRE7JF z?OpKpy_8K7*#+GPV-0u=yMQblBs{IV;9T)l>Tkx$E@+WirKftS3+kzn`C-9b@PmG^ znB!p=_}kR~Jd@c4v#&O{HNNkH84Bm@uF5VD@fDU#YU+Z@U5&>U4tIf@P{|#~#och< zm3y|Dw;TNL3mnlB?S}IwUsOF&>V}SAWhOkkx}kRO;~rPjZiukof3L)$8yNdH>|v;$ z>4ro-Z=FZJ-C%HqPjBs=Zg@Gf?|^r5H=HzzU^Kk!hVRGL<1j9)D|&C7aVW6guO-TXap>f=nxW62&+U}2Q7(|RCP;pR4nJv~6$ zFIik;-UHipPIhR!^g!DsX-0LXcMpUvsi@<-*#mFYLRim6_rTS~CwDLS$e8Sm5$*-uFOO}O z$oIm&UroGc)qBBtrAKwk-d<>0-jk+p*9+x*H|}S;_rj?Z#=8veE4>h9x1z&2tQU?R zQLm|t?}d7g3e#QLy^yn{w=Ju%7leWo0|b8bLQ!d@?#0$#fNydN4WqrV)^AEkU{N2O zvt6KX$k_*d$KT$%&eI2JalRJ${C(g*`Q-WV<~}gW5xy@b)d#b;W7nA~_JQ&K9QDUn}AIKTmB_FcrgFm}#p9R|Y!AmPG(>#|x5D?3} z+3C>-w{%R$1Q^VYX{=GE*ZM%-_HlXO%|3XzA~)>)?LP1-y4yA!-3OuxEMK-H^})hr z_PTqY^uad$!Y7ws^g%6gE9>>!K46^r<)<@T)CX+>Laq|!eXvisK*Ic2A2j?r<{t2; z4_0L*n7nH519!@v-nRZeFzGivAu`bi@rg}W4zl(`@ezBOpr!qwLNj(Q;O>XMHR9Ca z_5JX+JS0X$q#vpu9g5Z8)(_uT>b~}s?}z0JZBOJ=`avPA_2sB~KLdm|o7~x^+YcKW zrB|36_ruTP;kEvU`oSS!dCm*Feu!Pb=hN-n4`*0R2SiWx!-6Y!$IZ|8!_sH=cY?1l z_n)Md7X z^XFdDQSu{$c_8y~X#TYyHYYxGxbcJ8LF)TG&+Gdke4B)JUu!>fn4cFB@9T$byN`l08F&H z%s$8tfYI14zX_EAST9+%UrKENbgxN$HfQJzfZejTUpI^fVEM|2zRwQ~083WWoqn4E z2znH3A>lLtD;}zSuyPxKck5|iLp=u|MsY%T+8yq01TK0L@)wMnGLq2 z1ibvtY(RT>y#MzA1h!?TNVE(<(zh@l>z)B%&r_WV8ySG&kDeQgW(FXk#i@OC$sjB| z`1IhGm4om*TZU!t+CdO}zN<1|;~?miT+=BO8wAm~Seq%CLD<@&&A(lF5JHSzp0g$g zq3LIi&n>M%@VYF;PGRAT2~o61N_N3h%X^)<*{6(uc8u@RNg3sBvfI z=d*)wLL#ed;^H94U7=9KuQLapa53H&JP7N;5-n~-4#GJ(uO%NI48j%m4wI>rK_K#y z^<^^$A+Fe?&-T?I1h}qL34K2ZT)Py--xUvnwIM^bu>0E}yvTq3PPCS}eW$S4s%a4H zkE9HScMU?s`s#qY7kB}+7w7F9)bcZ((2&yAuxE>J7@ahN6jHR~e zO{UBarBCV~vK)dF`EH_N4nvT#p)bAo_z+ap>2}SYX6_f&t&zRJ>>xhF%i4Dc*2L^8 z3%oT1>lVH-d=~W8NX!r%POU#AlRN}J)~#ti{B#K5SGiE+%OSX-*UMd!KLmAd z3a=(V4FOSjWSw}$5HM)&O{2!uLr`?(qxP-FAu#j*vinoV5EwZXPtOhv!J-qvmI{+Y z5PQbh&3@r9Eavke!)8qqN{7zdp8WjMJd;K@y)}qK}Ms?MsgTJg!)?U zC=5ezqfXct>M(ewZX0FUISdKQoH=Fp48!N$7k^po8)m?Vm%?_}&4=O3jS~NNhlfGw z&Pt`JW5Xc6dfRPzk74+=d??+~dl)*~=@;%^9fluzJ!~a`!(hARyXHjrFhu)SQAMMN z;i%%1;RA`oFcP*wFf4r-y55~(Eqy)=jS2T2u)ZCJ+b$L}ibcb4`togFhcf0kS4G;w zelmvPfOf;FcMZeLS4#OZ&^8Q~zt`TD>K_L7d)L<6j}JrsfqwIQ3r65fg`{;E#|XUo z&b^tHdj!M}dnd@Q9|1_wZL<&_fwcpN+WjO(z%5Yl^ar^SaJ5+QaYl6nwg%oSQ&b-T z-jWC5M|Y2aNhB#0wRZ%(u01_g$~ZUzKl94HCv8VS`x;}fnDYq4W^4E#avydmdYD}E0+5fITVh7N4s-r)Rz(HJ^03~qG|*pFDWK47S)fydZC=Q;|Wrl(dEpBn|&O(Ba%FO4$qfWJf}Zj1ug z(%Z#$A)}zV{I_z{-BGxw6YccvA@jIHW{zx+N5Sl1se($@D3F$2WsbR{KsAk~h8K*2 zDtpE0{E|_)_Uroec;zT`y?Qhw|9cd41sKDhoSH|$yllT^Z1*U{Ie*Ek93BN_9qqQo zGow(~^QcZ?$rz-6T%2sZVhk3G1o#B38H2n{k7Ygyj)4#5buz2i7)07V8BmcKgEFJj zOONju1GTe@`(t4YJQf)wRO}i9rK&v=3k=2}>f)%El-U?eL`cRSvKoVSMqzy2ZO1X- z3U28wJ23{SzP~Hk&W=IR;JEt~|kc(KxJF{zU#E=QzBz z&ItLzGY*LjRfQY^BvRq}$ILI|na~!0`!Bc_lLWtrx$OOt57ShJy<@W~40?l#A zE`K|tqQ~5J7u`8wG7eLsL!9vz<1oPvWpDaD!+*QH4dg0_nN-_8HeC?R-H>bm;-h3 zvR1;tICxlaN4rii2hR4>eaJciT2?l-K9JI{}|&5h!qKUEj=^= z>krwJV|Ejues1B6EaTV&81St;eP<$dc8X71Qdr?8=m+$0R_6F zGYMZOpxtJ1Xw8oaIKF1lsuc|raAL{ODWYuxo^Sv0j^W(LJh11(35 z8lP|NVda>F+!HMAlvR`9^4rzHZQUdY?Hl`)v}qE=LyY)p#V3L5*IMJ{a+9Fou6drS zItj-)>q3sGO~Q`Y`rDy8lMuogyrjfv5*EFTiCuDV5+<5AZUNg#c-P1qbjpdj)w#}S zPIjAw>#SFWYdt4{RdKZzr_UsuFLgIj^`8W}=QoV)Lnq-y(z++1_aJ zvUUnOeS*1aHcmmf;GY{DTc=>NZSREQ_9-~|V9eS|c?uN$ci*{3PJw{w**8^MQ=t6F zzjB4)6tvj}{nXe$1u6ld#b>OiV83--YRZu*2URg-5mPGI}HcybSIVwPD2vs(+OG8 zX;`h=vB6P#8YI%}RybQ|Zg99j2kor1Zy<6VvcT+#^o*%rpZ^PjS$WT$l!xbBe3te5Ya00--;(x2B=s zPUYj(cc!7R)6!ilb{hP%zwYo%o`wxy|A8@pMkZQb{*%?ZIi#S)kM8Imws?uQTi~ z9xt4MyLLGNHOpq;VUqJ&b-Q ztK*@~z;Vx@FWNh2VCp@Q@3m(JqBpg@O5HaDl;sh5wdONG6)^q8VLt<6Bhs}BE;CRQ z8B0InF#|r5(tdFa?-}L?f$J-OU!8%t18%~+H)r7Q7cF<4+cQviEalVr=ozrwmMol} z$ZRlS!l+B1fu!QX`eiR>K=V4M#*Vi$u=kZ?k3-Q6Fp_VEJuI66$q^0B-#=%-YAT3@ z_Yd=Ufr|w??K9BvMj+U`e+Ha(yB~QnJ_BJ}gLN1U3ua;J^%cbxOJ_mV#X>@rdluTR z1-2YsKMUtolD9q(o`uHiEGKFuWz}J9H5Br#d^DML|bEHMP&%y(q(@(3tW?@Z;9pmll zE3?pH@TXcQU=}>;rdck8%|hwxV!6zyS$Jb*exfOU7J5VrU#?D@h1!GbSJJa*fkW~5 z3FkcK@h@_Z-YcAi<#i_&e}0(-AC1{n+&^aF%Jt(cI`y-#LE{U{h1Oa4I3~6_t9KSI zY`VUqWpox!ge=&xjD>|6G{H3=dizf(rN9pJ?;$c7YeV!eSTQ&s=!RCO786^!*7x4k zs3j6vW2^HXHWMD-`U`lEcN21&QErX-14P@#1-6g+MhN@t{lc%+Ckda^FFuUr%n}~S zeH8I+tWah3Qmpm*B3N|cGUe^m61eq=%c}SS2ec6dbiuWp;97rjPtBVZAjC70TyO5dJ18yA(WyK;JlI-LJDH0-9pwtrlq^&=}mmVllGZ4>E$|TrHn?> z4s8JAViRKz4gpwoKA=A+Ul2}&t<`i^6N0PH{6tE`H$g=7>Q~%rh2hp=Q>$Z3MVPa& z>pthA&EU0Xgj00U7Oq)B@Q;vkhfHnX|6{iNSVo%`|@MI=vkpYj++p zOP7T&ziv4^wwHqyPr98RipWD|(Z2A@E%HDq`Mvdax&nx}4Y*4&6ruNU3gr%CpAx)E zE|gc0-vKU_ECE)0%5Z#<=f`V{RY3S?+c3wZ3iG>$$HQ*ps=yWX7UUOD!R6T5s4{LE zWXLHHkYCK4vSv#l>(r-R7L8uOYRYM@nh+~#GM z8YHGxy}6yN4!hm1^f9JfG=N|2<%VF|PTUo9+R%N<){DC_T8rn|o_!R3CC{ zcPgf48i4Z7-2Sf~hA^6X)gVyQ2%akI)QdziA= zPInD7wug=4cSq}F55q3!IE8$3E0_mi|F&R$pN_}<6G(i!mvsS0$kb&6#?e!^E&L6B zk|IW`Auq@4N$0XBqM-|$Ne#d4Ha8+VNpGX?rUh#Um|4Adz5b+(l6F_8ChdOAkQFhu zz2CIGBD<)qg9$(>+`j~8RN;FzNp2$p$W)pI0 z^&;Q8I97-GX_I`n$5rGctjVp18%~@0Sd%P$MIP^d9UyD_wBEkYG&UrcIIaH%0XxZz zvdE#XPs*fZyXkdpn{A|$uk*VfIs&9J=N2{kjTIyr$n9`fW)|)n#bkZaW`0I_|6KX) z?z(?Rx|->WzHbY;f8~fO%lgt)?|Ss|oS(gXy!C}0G5FDfv%h;Up%k)jrQL>Y#O9@+ zzKWao{|DQNd3-crV!P9fkKX1t+?#%o$7R1y9I>k>%Y%*8gzqDMqVKfs?H0(|L_fv%e3TjF#-In;t7g32wu%1x zyvPfGb$)ueY?+AI5I5aUC?n&w(;_;5YEI&J$|%|J!l9J)_z;=wD?M^g_z&6cY4BxX zb{VFq`@MIn1_4bq=$2@XTR04Usv__SkNpTlev$vMonMqm!J)gxvLM$ERMV-B-*B4#u8 zo5O6)h}nt|vvne7i-?%RZ0tjb*&ZQgOGV6PgqZCzVm5(ObC~S{Vzvv2*(MRQ=~&ER zwm*p3DiE{rn$BT1{@OXr=8u@o2r=8kdCWGAn9Uq9TODGy0mN*?<2lT>A2C}aVz$jU z<}llyjycR`f|zan$Q))HLd^C8F`L-UIn1_l#T;hisX@%9hM27pFbg^1Z!6whI{F~n>h^O&s=F`Fi0wlTzPd5GDz ze4fK>I}x*8Ldbvzf*%l|iNM$A@?m`xQin*?GugL%xxK+M*N zm@N`9n;2p?7sPBzh}kR=v$Y~-`-_;(9x+=QVm5okY(a?G91*jr&11Ibh}i%!TQ_1h zH^gjxh}oEdh}kwHX5*X3Z0?BJMi8^z zM$A@*m`w{YTOMLIWyEZfh}rlMvn3*C^F+)RhM27cG24nybC_)lVz$nC%qD@Djf$8} zA2C}gVzxxYY-5Po>=Cmq&z-|;ClIr3M$Be`n2nvq`hS>B5HVXFVm2|vY-~k{*-{a+ zokh&%fS4^4F`F!6wnD^gMu^!A5VLtAW;=|S%^fkDC}Or#h}jlQ%waYk#BAR4n9U3^ zn+#&M=ZM*=5VLtAX7fbMri7Ty6ET|`VzzgP+14OtTZ@?OE@HMR#BAn>*`($%+jYcj z28h`d5VMsbW@B-m!)%ujvq>XnJByf26fql#n9UwBTQy=fw(2>|MnTNRjhIb;9<%L2 z%vOz)jj(?ZOah?wovJZ9@a%(e_M z+jGQhzKGeH5wjgc%+`UJO$ad?+vz#X){mHt9Wk38Vzy$$Y(|LLBoMRJ&SSO+#B6Ga z+2j$kjU#3|gP2VjG20Acwj9K4w-B?LAZAlQ%(fFTn+0OF5yWiQ5VMucW4619*(MOP zB_U=DM9g*vF}t3u3XfS9cgF`E)%whqK>>}%#Qn=xXxF2rmn5VJKQW;2<`Yy@I9 zU&L$$h}p6cvjre#n?TIwikQs~F`F=Awiv{0M-a0eK+N_UG202mY*!Goh0bF(Dq=Pd z#BA)MbC_)%(fLVTP0#P9mH&ch}kY6W*b7xmWh}x7%`g}Vz$K*bC`{Sm`!XRvpqn}){dA> z88KTkVz$?a*+dYtg(GI$kC<%-VzxxYY?lzTNg`%jikM9XF)jjyNj5u zA2C}lVz!-#*~$>JT}8}hjhIaVF`Fu4wgZUSOc1jPBW5$2$83&>*;XNDTa1`395EXY zVzxrWY;lO$nh>)EB4*=3%qEMNtqL)l5n{H7h}l{ZvklK-x0HgAZANL%(e_M z+ZV)a$%xrF5VQS3%qD`EO%E}f1Y$O6#B7m>+1e1ZU7W{k@rcbVm4#MY_Acstw+q3fS7F!VzvgvY?lzTwIF6&H;>sqBW80$%$ACn?H6LUGljy?M+ggqZCJVzzaN*?uEt(?`t4hnVd?Vzz0- zY^M;jy+O?8g_vy^F+4dl2yN{U7 z12Nkn#B3Uf+0+rUQ4zD9f5dEVh}l@_bC@j#G21j^ws^#B zWr*2|5wmq7X7fPI=7E?kZyvK1BW7zx%%+2wEe|o9D`K|Ih}onOv+*EiQ$oz5VL)r$85I{vn@i*b`dd~5@NPDh}k$1v$Y{+GeOL@;MW{x(?QI3 z6)~GEVzzaN+0qfST}90HX&$o)B4#^{)gGD5VM6NW(%FiY}1I@&Ld{yLCm)3&m3msN6dByG1~>iY?~3Yg&<~Q zL(Jxen9UP08$V*UuZYu&w+8HGY_v2 zi$0xC&))7reAQgeWwlxDKZ!rlBCno?2&ErYWq6$U{Wx}j;L&mFe;qwUgv%1L7uoBG ziQ1_Oz4|PI78o54&R+j%o}2O{G_xR`V*x|$i9}N?aMOJr{ygeu%}W2J>nNVVwSeAn zfa8;A_aAaWfKQrC^(XQ=|K#(!;R15E#*4xXA?89!kupnST?Q%fCRuSWvpVXFx|$5a zB%Tzl$TO8@R!8;Bl%^kGR!6;+c5QVnvpVYd>y~d%nblD#Cg!}o%<8B&b|zIFliy2W z=7aGhul=GydS##S>v54M#j_q5KiDcl-p*pKJj&uG|4vEpX{7(BK(N{6uvKm$*Pgnl zUa~wP^$)b!c%OSRj0tUY&U{$;545pO_v8Kt+Wa2*`2PcKpf;tDNeKVnT%TXG4Sm;c z#(w!O9d04H8qCYA!c&_~x%Qz=FxjCi>dn;xZOhaOPuy;W6rqoKzs1{OPfUy7 zv7!#}zp?ybgA z!%6<-!sLqtTiZ>_BIFgfG8Rtn&1A@(nRDA4w~&;BnxRgPVx&RDk2n=xak8(iIHk5* zg4DYf7#=sqi~>_Lmn);nkS^gL%J2T!PKMiCpAgtBN6tv!bT;}UPo@op9%lS}q)6s} zJ-{QMwS)9pQN=sFSe1&Y(3{!vzdN!*%APr?ep0?~EvkH&cw>3Wac}rxB2%x*_oL$;f{|-ilmNoS*L#h* zwi2!Xu~Y{6)+HgQvQMgxJZ@@k&etVg-1yz~K3s(uFN*0BTE|Cxj$-d51jcCso8Rso z;{5fW=C6~S`a^`CJE?L6%eAjukwAzXmPwLapF~t;f4JKh7DVj%K7M4|S9>Dj{Tik6 zQmOyk_>^9ACxv`??xdK0VE^mI`qqS2L<8l~VMBuZbp419k2Fy?>%a8kj#--J#Cct!gD*1! z*Me#Vw8Mt=N&f1V;>qdISzNZ3O%)m7(KUVvH0@s$T>daDS2ClU^J#O8=3|uqQ z5qGX+2Cm)O^!QIkdnn10Vz6q>^fl5dOvKxY8MtO(yx(srGjOePx$(1i%)qt2&(D9e z%=nOf6NLs{%)qr2?s1{ve*)M1uEnrwGXvMI^vnL#UnHjhi z^Mjcg#td9D?HY>nU>cPQ!~!&UMt9ysZTylg!J9+Bw3Qi_$-++_4Cu3AakZnJ^OrEB&TYIymj)Do+DGH zE>i8g@I0P_{)xBww=Gkqs^z|ZbzhB#KG^*=t@1G+-LGPSgDX>}R;s(SY+%aN;MR~s z%a}4%u{V5`#y>Lk(+1Aavdf$3X2cSGN2W}bn00#W%ap0Ew`Es4Gi7Su&5+0SOqqJH zF~GGiV2DgPR&_3xDO0r@1CJ&$W$ND)hEoMorlx-Binn9R)K5wjcwc&wWwg(*|*h$X#EOqp6b@rKyR zl&Q8a^RHGgWoo$0vGjk{@p&V{WRFsRVjfRyxpJCNa30?sm9a}3t%-s|tKW3T8~^*d$z1mui#jUu?vM*< z%yU(=acExibNgW;JSkXe>5;>Ptg~s)57q0o+?Pt7h-eRV!@{`m)`gO%rLp_F zFF9Ig(nd+Ek2^~>d1pzqrVN4SF4s;G)w;>zuB-E+H9&J$se-q zo*$w}Au*Y-l|s=Z{||m$wf8|JOU}W@Y@r?Gg$Ajksz$rXhTEAU`7vkrtN%*L<2%Co7{qoa;WZgKY7g3CO-JL2mgbmHGIR5UG}|&iB=rlPt6M zoo;cOhOr3i&0B@qfv1+uHB77VKU#KyWy`ye=kJfIM%wHZv|-y%H1?(2e|Oa;-dx(7 zXB#6*sNB-y-jhB?6L-Ixpx~GL@0+CA$<}`Hkb?y$RpVpxyJ>%{iBj$x85PTHh^~Dv z?8=9B5!|lJeGlnxCSn3tj~J>h{Lj>$d=k_hvgCt{YF3&l57B)q&);LQ-&)$Gyc1`K zDHs2$6N`c8;;1b_yt3f1h?yiuFu~D`mgW#ywpQBfNZ)nJ=&&&tjyUhKEf(D3? zJsY;GISmmmZ{!!*-TX~BQBSXyKmLL+5PcmIx*+^NuSe*%pjC)a@ln-Bx32Ij4_Xs@ zc+^@Kk~V}1tKHC1&?0^xE%Rbc5F(Cn+~Zm>-29)>|0+Bb4k=(Z-@K#n^g5RZMC)_i zbsu%(h<#$FB_&<~L_PbsX4+jlqJ42wn)Z6h|BRk-%()}PvE+zqmh=jPe82W} z=NLq|SzfRdHFhG})>QZKHtqV)=mF(wSs~Av&G{ow{aGhvO?Z{PJ5uIkNObI8yIEUa zme5)iWZQq3i-2jScBOlW@(hYtpZKS;vrr(#U;f~+CgH3vLM)q%#N)xW$* zFVhc^F^uzMuJ!h)&043(2fR6_Il7$4#`skVan7p;j;sy3%RJitqlLUAwjH-P~1G zRB0(6GEnZM=<<4WOGxJ+Ar-iT<#o+4vDjnf&uWea!thgsTEmqrVzc~;CrTaX{yoQ; z!4;_7E3>U$hHR@ksX7+6(mT)m0dZL%c2^T$91&%_XtaKh42c^#O z*XNzD8YnDU#|%A78Lbpwm!0~GEzD(cJ{Nr4M@A`U_qbyX3YrBN0lWf%OvMe)x z9b8nEpvfn9s@76jeEZKZ=o5 z19hn}@9Ub>x7Mla?;qBu)<$Rbi~Zb3o!+bOUSncGtvh`9fG>mXFje-bqW8NE$EbVd zHyxDQ=uTZ=#{SZwz>|7%xANGjjhCtH$p*B{J=duayqL3{E0`*uzDP*&Lj+aiO_veJ zU<`HJ>x}eLi9~9)(eKky_tU8?j*aXO{hm=hLvKEk(tb_V+{R;m`60WU1)R90+_bb{Ebtv-d z%I3|J)SZswB4OKDY2UwZQ61RDPK%?*CSQ={q7CKS#8rmy(nO?^x!AM#X)8N5ziF-z zp}n_LX(s97v=b64YwMqHry1XM74PR$qB-cNTF}I3G?mY57$Ym>)M*Zd7i_7dHf`)y zeq>ds0j+EI?`5}g_tE+&mZbOTbIB-$5R@6ye?c%74JN*();_3nH^Q!=xTyHlJ+yOuxJf8}ZpE!|Y) z%(1_3Y4rwO+zp-|XqRN-hZ1)c)B5~2cfZpvr|o|9uxe-hcN)k2H;R1mHMC2=S3Jo+ z@rSnLs^zIAqs_FB(Z37&?{(6Ie$0+8FzKgBR&L1@${nVOGW4{pbmb>$ddnNHk<%>1 zhmR9dg0hQ=H$RuiCE70~ZvNVR?ch-^VsZ5g-h);=1V`xo+|RI%kf}-F&7f~2*!JD7 zEz%bu-n4mmc9n|}b-DIr*Ip^2q%yw$F}oa5(m6bvy;+Iyb5m365}^=Bj%n*qBT1sW zSJz>$mIkv}v1aBA@{Ja;de^IEzeV(jMZLA96eAb)%(5jC-`=gXdhb94Ellj7-FGGmR|zC*DxM^4BCUQJC%F@MElU@#VwUv0 zGVklr!%ktuU!pG^r z6~iP(0U`Hs4TnhoN5apwyi_!^n2>1tN-3HCO2qxVY0B$bL6{p1CD;<*31u3^!PfaF zv3{YUZZcafG1)YE-fcw#K_A(wxyrJU*sRZ{Ff!Fdd@MYqe5<#W(2My-orF$ePnz_? z^a(LW%2&roY=b@p5~?lcyO|gs2Y1D1Ivo zFp9oT7)Y=}m!+#@uFWC{%o6!~y=n=3@8D7O^XGu9F8WHF>zwe;p`pI=+e+vP8!ga& zz6#FW*z$({WHl5|<-c87%?ted;w2aC;Dd)g&Bww&tq0fgO0lhB0>HuYVIAX^o)CnU zJqWrzvZ^`Rwa78riQ))9`EUWlqH}7-RF0U zHm;^Xc2`U6xh*u15q!vqJ8cXQ@3rIe;jbh-Pd4JFu&Tkgs5iDLYU&{TyeLaBNCRSm zR*vlM*$K8mrYj@NcEOc5ee$tiwLtFe+sZy89kBB#QAlsx4L4Wq7guBKfn!pcjb@_y zpqZ(kf78hTo~|?>Fx4`IHbcV`w&F%`%p>6M#}#`)s*EwDJUy}(dR!@oZ2uU;V(67O zDBTBBuePXIoaI>IL25wpcCFrc|2?;BBpH8D$`&uzW zG;ZD8Y5&Ixc6&Tf6;iZ@j>)=kh0E4pmCR!PsMs2my;HXTUTy=YwqC40s%-`Hob#=Hb9Yj_D7lB7GkZ-jb5*{g?B#>dbw@11(m)E z-~G(3@}$_M!_56>6l8m%*=&JvG-z6W3$tOu_umKCmfOIX##7U%6dP#w+_pH)%LeW| z;cUv>WdlaxPuqu=G56b2GxI-L!_vp=pB%Yp4UZEUElqOP@TjFXtEkQj9%=MNhg`7& z2DM)0iI^4G5&T`>H zdk(@ILB-r7_6H#2KSD;3eB$=+JCAI4NudKG8{=wy7C{i+*HD*kJ>3xUPQDah*Qc z+rPJTaAlq>T5VQ`o2z!i*(%!ptGjfdUcb_|ra%j(>MFh7?AZl5=_cRp|L%la4do9P zU(2w)gfGlCaaVi=e3ktl|p_+s8`--{>8lWeX|b<2a(P zvO)#iQw(mcQP=@CW2-(XyD5ThYLu8ZqgWp9f0Je`qsc+JsGxh)$L-)A+mu08N`q93 zdyefcNpN-O{=oWI9IOndBS(|OpqQ(EDE0Oh_&xqQY+b4dq}5D`tG8|fSJiJtKI%en zDel>6;=BMnHaHzJ{d7I_g$VcRGMDBH-OLPg`N+#mVCn3sj#&-qnjEYQg``#B5U9U1 z=GjWvC$Wh0fD0#VZ29S`C9ZgGE3|UHqKR#|rN&JhIc)vM^7o z51AvOQ$!tM^MF_H6QkQj3gztM7`pHN@vCpP=Nhft~>ci*$I zlPIRQu_#uw65*yo3~kM(CgNzV?Wxtf8wvBR)^8d|>WQyHd`o5{GB*-Fp@v+UmU>efAYvK_x z$i_4L@JAw{e{uPqox9=)>h*yBX!sbdpl`=~jArc^eG=?IoMrc-y56uQdTyMNE|az*;u!If z$KMYL~hc{fDcI{Cl)=#T`bFY;p*cpKX&da3;hT=BI=a(bIh@c$SU9G|*MD|_t zV^M>GghBvseeo08=gmDzjJ zvN^OLY#mEguVm3m_rDd~vf>HNpQ_h8x+a-cJjQ;ZGvpyn!ota8=JI{o{dK=%G_r2f z?p}IsyTdM+W_Z{4x%>rx+Cjx9%2ms+(h{mgbUUBOPuaUU&$TCM57&w$pZcvS67 zwHD1mnGu!Us78ByDsH$lnM&(rnW1RCRHP}Ys;t;~ayxA@ouL<0yOmb4iA-q86{gv% zSDW)O6+d85{DY0AXeTqN_-g?zMc&F*Fm-~u+^gDn zuzry0scOReG_;%gM0$6Cz*IAJKd)4vhhYP?rSVU;z3fk_R*Y3_U3)q8@sLwszUyad z@T()Hd(OY7no5qG`H`MW4G^$fRm%H}x>fV}Mn>{b8dYf<5uWojfx6g?Md@2q4E0tx z&(Rnnf?6oTw?O1k5cNmk<*Icr{HVq&CY<}8UZN)cF+1~pvnTa!LCfgPG&d?OIQ}Bf zIcI8&@A*D>_6o!Cb&qRgqQjArwN4EIqbH(p-wHC~S@da`AI9g8M)TC?Gn?@gcv zaj#({@ZVIT7M%&!-nwo(^}UMGduE;{^^NWz)pnl{wcyzcg)J&;sjAIAqx|Q&sJXYE zztFN{qxSBry54q$g?hr|dB&9CD5d=CKGuZeeU$bmenQD2trVAjk1tm>8YtF7%P;Kw zSw&g&^_n{GqB2VBqSS+?R)v(km&Us$Ro_rJDj1RCy9=IC4mCx;E7_Aq`NSP_@bap7 zN=YVvaDZAArQuWbuSm^MN^5*{Z1O^Xib6Fbf8oB%6vB=~<|kx{h*9K{hegY!&6<#@xL>9oZ;#*_&_Y1XNKAN5ohkss4SapzqTg`~Ng`<>63&-@nPeR`xLqW{hnXW?$@C zk}L^X%MvNnhp22dwvb4fGRj&ikzHjkv?$psib$yJlC3P^=RM>5Jiou*_nvd^x#zyl zdFGivUiVxE(;XI)0&8M3GQIKnCLSq_1a==z@Nq07vZU{q*#2+^ljRL{SqfyF|6PUj z&h%lx$5)M8L2isk>8peQYdc02tH}o*Z_b!mG&0k-GGye9=RY({rZ9G;%?O-##xfAG zZl4mRkPKpzm(+#Vatt})hKGvcyBSSwyb)n_yiCT^8Iy`%&sOQxV}2Dl?K%4Jom5*0 zc9ael9xy718lWF;J&@CF@qrFGV!C(md<(thT@E9DseyhP5z%-UQb{K&eXh@Pct}4` zl#=m$C7YiA{_H!~QWjmX%X2q#J%KLGx^&o}l1Y!*`C@bzPcS|6WbAUE&PjSiC~MBK zmFY$=JR6hqaMPAPlK+eLT4_K1gxsZD^;fpuB{{C?q1!^F3nsOG$*fYNZ?s5ASH6eS zue)1@PH9Qf52D>kp5FxM=Y9rd9V%F-*=!o0=<=DTeGko@Z|)kUUCz@|FVY^Mg^V4Z zDD!?#OB7tV*Xs9*cJ$_>b|c4H+9#Rn$`edn8STQ3x|I4ad9+o*u^B7-bXt^#$@!+MRnDB4)BE+Qen|W+P`Xt*#`ntSmsBcJOeV!hpd}nq}P$W9H+p zOH{9bj)-aXNovHWi>J-@f1_@m@#I^-{efyTEp85mzNV&tI_rY;YNojuzlPsvyi7f?x2HkLuV+*RrgL`8AZw#o`a{n=h;2A_}O-lP`-4PY4;exN~qnq;7cM(ndOh1WJ ziG=}YZNGx3p;WZbaM2>=xon(8&FB~<5#6*IaI2s4#c4G6Oiw4J(ain6_TgqqQ(C(l ziuVa+e(IJJNa7*o^fTw-o(Pee!1+80B?yOr03oIE>e3B#uQ&owSZ@6}RZzhU@To0$sxyVtJ4}?Xg z`ioGSIuj4p^=^;@1Y3@*n$D5=cR@>Lt-q5UoR(k6TgphdU)m>Ieq8;sEubQ$VTU1 z27T6dAx~)+LB&=q$4^Bu{QCX4Ssc{VCZ zfZSG?$S}UMLShY=3LAW#BE|KIze(FWNK!PH$f`ZrMS2uFyCLwliNw6IC-C~FVI_%T z>l-~$bC0yp@L$c118F2pN)u{)F`m?RT%qUj(@;_)2;BDix*zGIQBK}S+!0c0s$ODq znl-8SD#At7~86iM3XJeU$dmLiSM*S_83%SYOT^lG&h{vjII z=S|%HJ4U>pYabeVu#ZU$&sZysu5Kl!FYi`J18(2lUG+n zTJ9P#*z&>+-}jNkZ(vbD`qe<lph5EfxhfwOuEP(JmJ zW2m=;C6~qEw6Yq)C1yC`I<=56epT9kwdOX#Fy8W9DD?`VeC0#O)6#IlNl8B^SA;(y z*5{|Fbc!pX4w^1#F=9a&1C=7bDe4j;Kr@kJ4p>6OgtNV2JdE(ZCBeVFM2xWX(wY$1 zv5r5Lol`Hv9G$_loUa`#So(_RxxV02Ip2=|Vk;RZ-PwRYc>SLQ=3X)Wk6-`p;J_@r z`;EZgFlsXX37-h+=u!lJu{mtezsMgSAF|SN$<-A<%!3z77O}t|h3N#BCTQcCF5V17 zel$M3?ZbowQy%X`7RJRd2;tA@9`;%}wxX5VH88;J>>AgyNe(>~!06Q?R6hNtcc(>5 z+3-R||Kd|E&G=?JyzT=n9m5j$;7d2P2qI&6v5r_Rt3JWZqq4zTQXzV@S02Z-c=Xs` zQi>00t&ETV7F{#c;$7@||IiGt^>*&x&V^V-EmF)?p<`EfX`R{VJkkZ((DXCS7D%a> z)^tN0V+PF~>eqa0v#O@E{#NsK$b9^R@@mcM;qG*YfP78!=iNdYs%e_lPfkJjM`AT2 zPYX5a76xhBDJ@(qXL@Q1i9oEWo;I2n!B&kKHU^p|uM1^NGjN(eSu<06ilCZO3t+Xu zZeh)U?tT5>g=L)NhEBku<`@pG?GZE?*Ml>GF`t&N|9OGKQB^WA%nICWYvytg{0=U_ zk1>;*pMq0!e1(h9xqw^08t*(g?u#2S8df5X^P86+%-Q|LBe^Aei1o0 zjleYtMG<<<_TX~voj=`?xQUHsmOjBWPGjjfkj!(zB6`#4HxS1#s=Zl>L;XA*`fJNm9$?{5f3-zfgM=d3rTjnAR#fq^YXuluZf z_^<(HS50R{WGWUzkFGIy9|B|U+C|>Ioz0I?zIpt4R^tMCr67!kBMhUPuI;S<(bbM7 ziu@D0_Mry7%==K<(cu9)MAj)p+nR-5_lf-`TNH_={XwC7oBh!ZwT3Etl|-uE`9h(>p5 zU@z$kAtdTGkioyo*q`$?F5z*we#Q+AD@mVkPvxUD*1AvNH#+?_B%Y!08hVZz-bH(T z(%p?UrWtg^aVTEnypVwu_z4rL@q7N>M1{DZ#-e2&ZUMQdZe7)3t=K!Ht}=!KTT8d8 z>ucBFMI5YBKY4fdovVGGy3XHcD^9%E)K}wO-zda{sh^D!yq7xVt$xaKWOPmbkUH63 z!bW&lS6%nxt6}0DHTA!b`oFa3OR0Ml;SuWBcc}B&fq8#PO{sk}RI(gqLO-eTM-Wjr zgqqddu~LUpyo=SQYYciUezDcOYAr6&&c&#K#~mBb=lrMkS7Na8_(LbPA+f4mAtlCY z`cJP>YaVK;IZ@gL^E>6$4sJ?d%ai%kx+_JZO|9ot01tQ80R~YHuDD zYWuNV_3GKHpFHC-Ro8hVsZ96cIMplKXt7f}&Zs^zt~MtpyQ=D(hVkZRo2s^3Wq;Zj z#j8qbs*4y{LsiY*?1*E!@~ghE6P9J%pHq>~+bbvXp#ri_NB?4+O5}-f z37#HVl^UMhc39X>l|AQu7FiL~sEn!~$^5^2P-*2Eb9qCJs8a|lzT~+A6sEM3ae0b` zsztQ8+ofDYDgR8mAnNCX3R+Pp3ZJn-^=gW(&H8Dht}7Q;B@HN}o*o}kQz_VmYKrwZ z>$GE8xrBV&)&@JGe7A>rUBfc!o$_#Jw4z$;W91^HHqUJD9OWyCEms{%f?5rHBPfYr<$5?sZ+t`Vu0C!%A|)k{H1gVYRtwI9l+roDe6eF1vD`Yn^*_2B@hH~{Q@qfC z_~5nrK7HgqqR)22YPl)}fqa2IskC|yF{IXNJdo>-Ak%IQ*oc@Rg3Dfoj2^)w_6MX( z6Ly0UKfe1Nx81W7!L%wCo|5ICRPt|Y*rQzYLFq~lm>501>K&QP3`8*%~t7sFH(?a(n!p!q3Y$ZrrY z7qL|oqb#vexeP`Bl5wL>HbPPFd%+}iPH4+P=8jCMc?HJDq%>p6K83sGPcaepO$x(? z&t;|O@)gK;4m%dcUQ;*)WB&XT%e3Zvw0Y)Tz zMQ7*+?}j(L4H8x)-(m@Oce{a}M3UfHOEN#IsSx>;pr$cMA z(Os}^dA@2MnWwN@QE^L&FUN^Ta*o~V%&D#&yotO zBNxOPuESwzw!!ABH+f-|wXDyx#3^WN-NpLrQ}3ZrgBUxHv?^%vMV#Xi-W=$y&vSa! zZn02(reCCggD-TZ1o`)4ybW|Cq(|b(9}09-l=c@d28Y%^SFM}B#{<3CC3E)nwF$^& zouEmkPHGqAb!Hsy@VO_Dad>K;;;RgZIXE%BXFdu-p3I8w$NE5oS;x&^j#)x_rq!hv zOo)*5G82Wk_i_--T25x}%%*%k!m?32+A1%j|FNBxS0_$l=BsFVqEeEc>DfK<@=*f=RPZ0L$$4P$cQ>)SkNb6gX+SJ{+yXXQRvLdk5hGwurb9iA|8fF1}|iIHqBO0Wl`_>U*1{A7T~ z3Tirq-Y9}6!YqztO!0!-B$f8<5ucPB#4YblAa}`qRor(pP4uzcYX_S_kGgcZV|+}T zZ+T%C<$BDtHcVf5%4K&@N+44En(wQB zktGU$8VKiUl#OF4UU{3CCtFwcGyGIlysXA;H@H@lzwCcX(PHYpHnOz?ihnyK$g)iq zT>=}SVA(ya{$-U{8!~pu*+nV3BQgg?5B$GYua(5V4IsfsTRmU-=eWkAZg4lYIB}$xyK(XccVm&i zp6^(78-^U)Yu~3+SiGBD+mU+RZ(O9DH@Qe%-XwC&X7)+)u$WuXe%X}g2+WR2Yq0)p zVu9dpVnBO2-#`%8Y1;YWm7BF7E+$)#PB(BOW{ZPB+>n8Hq}kdakRh;&%>;p{9G{$r zkM-AN*xeuyoFf6=GG?!WKt`My@BlCSlsHJ1?Ya|$;D~_szwEf3AUbCT^fU!h&;Gg- z#H|dV2PX$eHi8EP;}`*bHQD=kwnGJ)L9(y$fCwBN(A$w~KJ-6xpkE`Kj~B$L2R40c zz_A3{Fm7d!oxls?lmNS&7sTLbfX<5S1>Wt`16MnM>~OQkxv2r1FqUuoguvwsb|D`~ zo3j?U$-tiB198&SZILY#s35VuMEI!FDEoRMNuTD0&T>$D-d8ranVl` z+fK+;2{Ti+9`9D2eiQ@AaMl5QC-$7!HUI#9iA@sU1`l8(0cFqmWnveIgG~MpTc5D+ delta 45567 zcmZU4byQVB*FGsAf`D{)ch{XW7f~q%DG{YXr35JvkP9Lp9aoT01O%j%l&)L-fQ+V=h^ejIeWG~5Nv%QVA9dVC!ogppX=W97^5_;j?`#e z9Gw4Ybl9t>jM!Bo9ZdP4(ofBxc{pQeN$+j33j33A$TSDGq2sl`E}D*~!Q--fJ`w;4GIM_8m&d_|mci*Oy&;f#&zko6v@i6Mt73fX=>qqZ zUG4kE^zGr%m)d|I^HP?_m^Hq>m9BiFY3};nnnlp?O3697i)4~LHErWTS+&tqgm2Un7m~<(Ys@+?E2}P7&Q6EW$C}3 z9OynUp`7=b89f&@{^^g=HMF>#TT7uNA=>=}75F3h2zo!Mx%F9QA9j*eFeuv3!aD;Q z*U9!<;7*82O6pTAoDJJfT0TgG(KWue>eM{po3*@iYb!4}dy{`e^7;#y|AX||tMO;> z(Z{gq3hu|Sl++)4IIsTznpQ3mZxQOiQZCo1)p9gEx8!oBAe4etJ0QeljSGsdlGgpU zqlS(rzpq)w62Jm|=QIlY14J(7G2t!T#s8Xs^58ITiXW>j9L54mIuwmlkYyK$≈=-+C} z6GRi%WNj*#3t`asv^plyuLaOOqgGj(-?-7Nya7hP@tDvqYPNAcv7~5^8)~xu>i^3Q@G$~^TZ&}9A1C26yCRU7ltvj z$Z=GwKmxQ>%}Iwew6V9eZM)9}FB6vbGb>z!NnA}Lqsw?O+EG*L+T8Acf?IM%*nms? zKlz0`nrjamS^}x8T-8hT*k|ECjY0xXvL$FP`(cT-ejPFeRne`8 z?!ggNBtZJ;7-r4o89j18hvFLp<>oOm3uaU8!H9~Mph&egYcSblf5rDY5@z?+mH z(I3LqnY6pAe$ekOwS7sCD{Nm~Uf>tOz}KXLQ}1hCAwOk;^Zo-{Xq0CrGI+xrhPp_) zS2AlqgjDK*UA0PgV9@+iPSll!_!tgq`g0YS75b=Fja~{G+BfqkZE`|s3~$C+Fcs|8E)3Y- z!-fA5K#4L#+5Z9Itq#VpzZsvZ;KvvZ7?6g( z^XH|Az z>(C5GG_Vcelo}*q1mwxXZ;!WEQ3q8lsS)~bEfKt)r;p+SmP`bMQ+jliT zNc_`ICoB6I%$P6EnaZ$)x90^5m8JDA!;m+_YxX}>p*QxBN8o=Ur`GT`!SeDxRlJ(d zIPbE>Rgd+1?a{o=QlHI2_3FCXRvoSLiwCx{-#~g>&cc9V#D76v{o!@X#oEV+bvP%- z<{Rr^Z}KHzYPh=<)GwHRn^)(?zXW_V`^OiPPcN9BZnXPHTYm{!r3koJ6Ba{V!_}C5 z%kPkTBMh3I7ehZ^U0$Q0a@bMVO@*oY0@XL-g0#{ z0~bu^joLAdxL~@!a+ApI3#JEtNhHAeAEv*(T|SI+!Stydz56#Vm~M|V)+c(w^!Bly znVk!!^UN3UZT^2u$MR$G{}&VNWBq=tsPH36KQ3M|+6$4DaJD8VcOkO)i9QbBx)50u zEVgq7|DVVTf4}q=+a_#<5ntcXbl1HJ<%JqXxutJGBUHu93RX23py!|+LZuF;AKHi? zT-AUC=FyMu(`iCHzGp|w0$R`xarMZ$aT`ija`ag^-hoHM)Sy098}6}{Mlm?+z>}za zJ|^0`uu$X6czV-an0UYc(w@66443LP*Amx*wzQRkC%bwWNU{Cl2B_ADSnu;IcHssP z`@R{|X=4bbD;kuGZry`|p897ll14C0O4jZ%?|t}hH6wA2%NQon?`0YZJ%A;?s?)Tn zhtR@W@>b}56Bu3m`kcw%6k3bf^C*^@!8<-e&+}fGL)QxZ8*XeCaMdNGqqNHc7EG8A z<@sVPVTQ88VY0jxOy?^I={U85YEn5a-c8o|w*@7+qZ@2bg}3`Sr4aBW(MTrfT8+ z6fRMBN($sV!JT5;Q3E$;SRUYt$;Z$?gAo>M?z+v-;QYp}V}P3r{Nm)d?=9>KUvjEw zzFc>OQgT{x9Yt&nCdM7+;*na-#6}(TRXjm!Ro6o zUdaZ+j62FMlLSGq*}RG(Z#)QYistMHRs_Q#Nwd%IKfHnQLIrt9SO^?hAWhNt4~2b3 zRrqv%VKBEK0z=as91iJdb;qg`BB0mwWG!DwBox9x9)r;+Xw^(h-0&(I-cd*#{=*sr zvoRt!{TpIn1R+DS>*H9cQQW?>e+>&+Ekt=a+pzF8r4tpe?_21CNQQ7K$HBiQ0e=;6 z;^FjuJWX(CJRHqW8x~G}2R%9n4Ygeopz&vn-x--sB8)FxWDk~3f}cy<@gN1D;}rU%S?O zWkTM9XbVlgEGS=Pb>ugi1tY^4hz!F%!(AtJ$y$|cD4xrUj{27kXG|7x2di`7VTk9M z`^#Kt+QwOHrIH7ipWoc|x||PlC^p05hV$XyKLMDppT5A^8`UK_ZeL-*7j`pSwE{S; zdo9xVY9X9wrev=4E`*p>^18E|YlZM!iwU!LuLzFu6&dGzFM^s3uW+u*7sGKi)N*Z3 zF^roiy+y)Y0%xV)y`YRLf#qnt(^|Yz=*&3(eZ`{`UKg$dG3%wUkWJnv%lsR}?)@#l zHu4QdEqT5u*DQlijIO=#YAl0PNwhbN<;vlE96|xrFXb4RAz0bN!&L!`NM}xm5-K3k zAECuN8u6cUQ^gsvu6(_tsebYM8535)k>l z8V>F%YJ#i*@~U=^wW>70RZ5Rxw&xAdo$^KYUTy=V zic=u~(+yDKln6PeXoPlBs^JpX8{wZ7gtg7B5xNim1)Z6V&?fn>nDR&?R7rmEi1AVr z?`)yLEhgmoc$T&~iFu`DJ|*d@D5g@cu><91F9rkEU;iJ}FzD z1XY`1+h6i{A*W`TEgk(Ho7fDkm`%GZo10;cOpbNvMl)o1ZC=Mf(*nIj29qdJEzn86 z&ez4d1tP1hP$8@Zq8S*5&5Bx}E*{(O-H{e3eQPRg=T8ePDhl(&JY;W$ij)o3VydmM z^~kE*{&6ehOrH^D2yKP>FR2RU^IM@mSu3`yyA@_jaObsewL*$sfgXjcZICL9Fr8DZ z4f=gpdHwoM8)P4q&oXmtgGTFALiu5B@Fb_*#UrN;y7JT>)VH-k;S~{)(1kYmEq>T~ z;Jgj~VhhIfzhn9iM{8);=4HRb2c|Ys-|l{g6Ib*VDK2Dn2Qx*_R^WH|m?DhL@Z)z_ z=UgDFS@9h?o$foMz6I)a-JbS&%fDx)ecDd zY4d%nK?md_ZubB5r~|&{>@yI2-T_~D?G0K4cfbmrG8N5t9q^{L{o+bi2aLnKJnEOd zP*L>!yK!=@9T3IO>AOAH0pFAoSDVduKyoA27a==5(s+dqjs}t7gP}!;&b;1$a{U|2uPN=xZ z_qM{R6OP?1{Y>ZG2{G0jdby~ePPkk$93>Rf2|pX#sLdyL!VOe@m3dYtoW6-)_PMYV zrgD=$t*YvTcl=>V%@`s!aU- zPUxLiK2-Ck6G~IK+_*u|1r7feUQ4^u1uJW}uJ%)RVPKLn+tDy{7p(ty=2^hq1&4yJ zNSO$CL4pVUEb~%b7yOo31xj5o=du!TQ|W>u_vuLd)w>{bLx0jEoh}$*5Z}#vuM1{m z#G;!_x?qo|&07WQE*M(RP9A301@WSEtt*|opxQx+Ppd~4-0){!E%51r9b3A8JYQkD zU`-w@99jLhU!U9hh-^x$b-7o<4UO#If=1#esz z^__3)`mZfpylm6|?PedyH~yarcTiDt#&RpVVZ!%>YnFFHf&Cb1?({?wxNj<@dzmH! zlt1RX>`{~p6pXF3C5)vM#Jg5RkqK>j9k+q;qF8UL@ zX#?0_XzS(JYXrtm>*9?9nn7Z~?U#*|t)QZUY}n{S8xR^gvt#6K2c2~`DE#~m49L@d zmrE|z1ya7_^C%Yl04}7JMSj9P;LD^H^K4cxP?-?O7l`Tu&#pCbJE->qCB~rzX}ke& zTb>CwpmG3=MK+J>ISqpDrq!q{t|8!%)EeeFG6de%rO&qq4}&fSOzy z0rID5Yss;r;H}Lq_CXY83@G|eyJ1$wfcvr;Lu%AdP_A5^At^f!svnCn$WDxdS7|78 zw(kUZFVcqhhJ6y~`}JB5HcSF?(m~M>vnlXXkt5>s-zi{?F_{ximm>Y-P1A)iW^&pzhTJc0Bdys#SeoyAZlgx>gtC% z@VwynHfCWC5QrOTFbm9s1;d?{1;=^7kI!K^QZx@%tH?w54(35nZO5FP+yYP!b?(gc zTmZR9dbU#K0w9NG7AGeQVEIOK#)RA=_-W-9zv#XQp44g2QI#x$vVZNcbg^fR!xIL;^PuXIFDZ;A6o)GD^w+96w6=`{YXsY)-qs1o&-6(Tm~!e zX2iORm%&wkllN*X%OK(G(fI=73K-(##Q$)61<3VEwYXqbK$uP4!rFj*!)wqVNZLR_fyH}~% zbZdZ}8s`!Qu7NcoZ?wacHSk$@OW!GK4KUl+bq9Z211eX<7{86Lf%Zq6K}Tn6U|@8i zMvLvTJf_?#9?s^~dtY0979z%kqYW)RpP7An)=YIh& zqRlU#&wqiEu$v!4nbrY`r=h}2*>zxJ9&hZcvksJa{@6y?u7m8pjLtm&b?`wEkseK0 z2a&^F$K(a;pkCJIv3|=sAfAs~$Q)Y-asgLV@HW=LgtNhp-Pt;z{*_=keRTtfQf6pA zVBY}A7)r7D1+fjVZL)OY1~!1zcfELWtqqXv$}?PSv;j2VJ1dx4Zvam^^;}A)4e;_w z<+ELk*C1LWwvIy4H|07P+CuKdv(K)Yz7^|70U37TPmd6@hJa}XM`?l_d z;lx8OF5L@9ToDBfcP|`s?j>;PM{NLqbaJ)fn++gv&-?s=-v(%iCRsFpz5zl@{|G&H z+yM5WWi*bK8{qT20ov7Y1H|}gOn9naL>!1ry;Rx&W%SrJf8h-bKsz6#zh>S5K@;~+ z14%bP=8gI%fv4*rAeI07tL1f2HsM9&JFpJsZqDEKs9gtBHf0YSKCc7b1L6?V=yk9> z{50~8=Q^ko3-eYmTL+`7lx>_h*1_T-dk~@EI`Fo@Ra_xk2k$erBO7;q0a2p}Y|_9l zkSM2Q<4}V61yoyTU6ioD0RL@kI|7$qU}Y($tX=yTAP~Syd?Wk|WJ#$`poxEh!scz0 zlle8!mF_c_RkH@%3g30=#jOF&7FU^5hcz%gl6pD$`Wm=B8^Vt<=eP^b`1LW+@H1&?mi6U7p)0#Zd^b^4(dFvG?fl$5pt z_LH-e*q^QdHLX7^844?)atPZiK(GSLr)ITFx|e|@>HHlKzYKnRVgv?EmH}rM_S}wp z84OF)zQ6Ko2~bnr|6KKT30MlE*w=D9ra9RKrl71QVLJMG}dB@~@X&w+* zbv=WhYfYjyZ4?H(# z0ol;*LT&miNUGSykT~6&1v(GB*WyWMfx9`|Jzw7$P{q%6Fts=h_}`#Ax~802Zql_x+rNVbe}V)vj^3XwW8n6&HW2~e7&w_ps(unO3U(_oODg;uBcP^KuFqV11YA`&m=Z1- z1}m3uaeSd025q{(qi|mi0lVqeRn>(-AnNz_N4&-$cz*>u%TO`^PU`ydR+tCCT@oK9 zmN)%i_DIZR=AaKae$>2q)2a^y^h_WByja90%7@QaLwkUm$kWf8F5TeTJ7sHnvmc-a zvwo2BQM(HeQL6kYSM3C4<+(Q_KnEz?(3w+1wS!(79dQ=)cMweegi2JS4P4#~Hv3@I z3SOJj^u;)|04%bZD-hHSIw*BJeKVUto=WGdoX$p2x*;vseAoa&e@0%L@ULp6a zx(;Z^ZI*b;)`H@i+UKCQ8en3+E|O{KRDl*5nYE1V3h=8>z4lIWIq=fhUKfd>23uX|x1LEPqZC5oFLz}tI-MfRM@Ak_4EJ|9~m_;QWz(YoS0An@)= z?<==BaER|G@U0UIUhDPAC@IH+tlp&Fw(rrPYm|aZe=HI_COG4HsuO;pGhWO)9EJc~ z57W7Y7ECZ8*e4(QSQ7}&>5F8Zj|YI0#XKhCkLBnBjO|9ba$+p3Vbr z3fJsFIelZcw$u}FJ%64yrrHMB^73?GFfXkD_rQm&WH$>C>p%1(Siub3C=dO!x9|W& zL=1tkJgu#SKyn|JGe1Rvc}0H~%LxTwx8LrIDK3x&D82DH zvso#?=vPG3Y9#@jj_M2wEk(ig?N*_XSAsxICSmY1Js&7KeaSpT#|_qFB0jx#;{fEp z+}uypSb@SGIk~YiBe-ciy8O3@7CbsRWx^(119C1ZJgut~;M!lduXG$_;8BQ?=|3_O zK+pbrKIb(d(9^-glybV^0l&)SWte=9GXB)Ss!Q|-#kYs|B#HVMRgL-Nz4v7w#b@gs zrQEoKI;Qp|FnYOxnykJ*IwiJ>qK+Dz!S!E64Lp*b`eQtY%H7br_snn#6&}T@x70L- zs`E!H-kTmmS(3rW187s0o5IPiM}v!*QTdxTlQ19Z}-fKEEk1-q>6h^W?hTx=C995G3QqH7a5OqR9$)2sj zgR*KhC&JleM`fzqmKd34MByma)b(i6pe7_LY^ld7P@L_VD>|u{QSX0v-x0wjL`AhJ zD@=#spm^M90!@j2E1$V-UYc+|P~Hls8gII^qwJoW>%~mBrcCsP#!3Feg7Ug(yMI;6 zG)9>uaHe63y{?5$EJ zyw zk2;H1_IV+X>O2fl*3WM8mi-Z+oHeugcAO96t2{DO+3&LDq3k`jc5un|nX(gJpiQr} zoiei-56kF@m9k1_`PTzcQ{`mSpz}{{hRTa`9#lL0ca_&3Pow7g)RjjD+=AJyt}9Q! zRiepIthQYYcSC3o6MrG<{bKC8DwN~@e1J)f9*lzwNBvW5R`RT4Q(yVifHK}m3$ z#)8JMTq$Va?-R-6LZ!#_xP%YVvX!1Hq?%$Pex@l!EI54bA52u*BS5Q+XT>Uohcs=O zt%NG+%5-|XG!Iac&dW8AoAy>}9XLb3!gEuSpi<4pFK|#w8l&Bs3$#{RPPX>X)-_Rz zh@VURyRNUquiEEd(tBIUpNxknQu>CH+ZDXW9V|e}+Mcob#(f#3D`RGrzhy*~WD_PZ zJ%+|SNFN_s?;Tb&O76m@dz>Ca1*bqixf#Zv>vd+T=< z_v2W#PIc52^ELn z^U@n*Txqll`7W&N%iLrN!>oZ`ja7IGH$mHS&y^E-oZLk{lhG~tEMfn$`}GSLd2-6& z(Iw=kJRcOkX3f$k&n@-vRYy>Z{0-$3pLY5x`D2Gy%yKb>^0asx-D!B4@^5A>)+Upa zRVyi9zNX51=9LEgq{HE51UTwd0{r@5AfTV9KB;RAMyUY;=XMSI{c za`|~`**7cP`0`imysrCwIhK2im1^Y!n{pm|&LrwY^KwcpOU-iPBXZl)uXi1hF1fV9 z{LJ4)4RRaOJDtVPO6ADf9r_z_a^%!g3$31qq{ul<;C^gMiIuCt+*3zm@`B{7oGbNi zhI-45#(oI=R^%*aepUM4S7~dxv@VL5D^HE&OgJ|ocE7JzdJrua?;FGf9L2g$&py)-tx8lEo)(GUew{SExT-! zgkr6km(5XYztwhgMAqt-20o_Axlyx}J6`7L_+An=JSA;n#{t-Ob_nFc@pEmrEt z;yV*oN%>xvjmEwh+%1%qeWX)#Gr>YoHf=?q)jEw;_Vv4@lE;mdvRqx-m_)a50@)cs zNra^8SjILa>-I;3b(u!9I+ERIGcs3ny4ZhG4#>nV*9#ujw#d-YUGg{iUM`~)@Y_8@ zK2JtK{eVb_`J+s_&6L4tZH$bFtg();*lQUsn}2!xYaTNBGKB=#b6Xjg2Qh5{9wsvH zwVElP=j+I92s|-M9aE7RFT$+;%~q9>d3<6$Fm=Eu!$0To#CDBQrsM1VXWxa$FOs-f zP{0z7jC;bX*J3Ms()*=#S_W)O(xUeVUt7eEO53;$XW6QDO500oQu)Z&O0zyOj7yg* zl+IyD^|L%klm2z&a_%vTmF8|ONO4^Yls+ArjSLdSW6NZyxfC+&zC_o=7QGLdHX zsQ3~>rz0&(>!yDt6D_TLOusLwDJAW%u=YZoj7NIL<#&Yp9<4N0`dfAKzr@lH8#khV z3jL9~okLck6|gB~d(Ds`33pa%V18oT$Dm(Iws1m&t)@u|W>wXticMD9e=U zfmj+TgLhK(boXDUlZHrP&j0-qLyvn&t-eW><8O4367Jg1EvPq@G9GkfoU+xCYC5Pl ztMNfg)&F9`Vz81@!$bcp%?r7t)S1gG11G4Zj_l27=%_A9B{E(~7|=VGye}!%#htn; zDJh>jHGKF}QWFWGantRR6j6)VlP<22yyMxasc_}1B#)NwXGiOg7|9ABs{JFYXi0eW8;-MbeUn(NWy@ zAM!}<tv3<9)ME`eEUS(1di53k^mA*2I#3gpytX^X>3GG~o&>@R| z;^gjbt}Z5<;^Fqs4qxa^iQ_Ik+%!XXi_;s`=A84^il@AKfA$xH|3#d;tD!1*^u72c zbKlK@kKy7bwE+^#_kG0uUfezs{$nqW@0};>7X3inUUaQ2f?ZQwYtYA#DNj*6sck^d zL|jnZk)Gw7dm)23BrTy=7a|e2aHxzzQcuJd6s_%>$X3PP`2eXc=TR}G>2}QT<~A|! z8Ug8|t7T#nj-FbWQ`1Z_H%-w*_OEZnXb=9}bp9707OOsxMtR#+Y^|;9&i;E#vB){H zj;M{hVv#O2XK6Ymz7mK;_zlW)4E=3a2 z$2CFKes7LNQs#1hC9Exp#H9NuB;6YliGESKvebW33-L%gINzp7Bwr>^>nla7$kv>d zO!}t?5wcoLNZqEkx5(N#oy9Jftw>zJOD;K@Se$M78zQs4);aH8r9@(hYzmCq zIYczj6fR#SC`6hauIA0V{1fI{-d3|{`XwB_;hKL_V@#MJ?q1N_OKrl3u203E?Ux82 zCQq3%kfjS}&WCF>Yex&Gk^OZlfBI6m_mcD!AI46Y2cse65!7%`*em$Ng9?2^c$HgI zKO09%xD>?`J95GWx37%TCrklo31(MHo@q2%4jK-}L= zLSH9p!YBI+h0ZyS&W9>K2$8CdBvA>53Ssg_>a*~l3vrsT_7N#q3#qg>q`EZR6`BM? z`c#~V5JvcpuT}gFA)!5wfu^-dIw8BB8r`yS0-;YO>_SeGdxGUrxha^!89@SZ$9K^$ zx&)anUuD3wmJ6nw29er_X9zA_F{_BLixIRREDCu3!B>!?L{v@G(oQf^JgdyB%uw)h z$&mkM@O8mwf^O^NZQ_EqMao09&CG)0mE^?1$;22zckkzkeg_8v*GWV^{OOq!fNYjp zIOraM2BX*c?D>@fN9z?U%FUSqZb7dfojS$}40Mu2I*MWh+H>-J)$iL0XpZt-$x}8I z=-Fiua;Q-eaI^au{T3r8;6Lv9VA7aTV2kzd3nrXP0$KUu~5XU2BntwJU z?Cy#mrh|W^ENnVW{u@7#*h_2ql*Z5DW8_Ss8NqKig{M2M;mO~MXTQA3Y|UTi6+fr= zL5Kh0-dCNMI!gSMFUMHgp78Qd;f36|Y)-|myLUghnDLD7R@kw~m-JOWO%KC=Ke&eY z{PeSbrn@xoQBtLAvybHRQPE2J;7-T!4F~`J5QgdW=UaJ<^9e3E@NKJm5-#4m$0xu; zLP#>8!sj>W(afzY%6ENvwA}m=Jztn-B=@@y_E$9z;{O;P_bmV#&wT}A@Q53@-dpnK4sVT*|ZSP zeO~EGsstJyFZDb$Z#fQ6pLJx6*!UVZMuH4;v*=`y+sP#R)$~9ew^dd2<%P!2+~mqM z%ztmia({^uoJCRiaR0JFTlzTIa6cAHOxD-a;Wm?q+}oX&=Pv5P?Zw^XT;!{yiGIOBj8qAioa`HIinV-*K>VYY(Yl*GE%5d>Wg^ zj!CfSoW0K-!hXGTJ9+B(8T)tnBsRCt_u1t?NJLGltFYr~E*6#%3bGq{q|C`WQn7#A zKgP?z`OU^F->pqm(*(zDrPGX=KMq+p2*e+veOB6{n=7x@>}_m`ZdjV77M15!KH{omghHn&lR08 z@htmmc%+h+zATTK8g2F?tXVGEKJp1@*JKGJ+a1+zl4SWPDg2=34g-tAGS$)<87|A@ z>1Y4%N6XCmXGOc~3q8yP+wZ2C8%vqPFp9=D`&##!X6(c2&7M?H%e~Am9 zeXYYRug#HqXHACrZ_^)>8*n`>{=u>rGD!$~+qz_tsj`txgl%AHAn=zx> znO@@hb@Rn_efmy?*j9ooiu5{Etj}wt*y%5GrT=2&BA}me)LTBrS)*>UxzE&`YOW z)%!|FzJ%^26{&c6*E_mlNj-l8re-X8C&>r$}MMoT8rp=N$uor6EqPZw-{*UWFE1+pk zzv3QX$L~K=WTNRa$@^ghg2?J{#f{W~chlPQ1$zO&Uzzvy?GZ7!VtXj>l(zO?^|>EU zRhRy^ls3iK6qD=bo}Tdzr4O0VbHGKI z-J1$6(>kCqWQLEn*2WujcyS1|)3!c-rr(Du?DCxB?o*I8G`sG0Vi(+Xz9qB!I1ygC zQ{{1*k`A@X(Dd9IHgIF}caoiO06eG8bAFo!M*b2ADAXi`A6Sk^ zR>%*LmQSnBS2h>^>jZ%M!}9+tl~Bt)Il~g0+9~d%Kw#Wke^7Ef|467I7zpQXn*Lq! z2d~UL^=q~5zzYJ4o55V@|0)RH5w5St=9|92h<+2)tcst9N0;1U_MfdnK|fyl5Q-hx z8YyJ(fc_AYWZs=9$^HWaKhfQKMT>`4seR0v9ZiJZeHSr06hMwnmD_IfeR~bPB*XFU zW)D3&xHrVBXNDEcdX@U%RVybtx9`hBJwGqnl7HdZu#^CLnNB6I+gJ!qEMRZvy?${J zrm^|hUyMc&-M?qIt^j$^p+C}Y5F0b2MRE?!ADED%&n)n%YwXTpR+!t}ruE-&gM7FD zi^39Yh&?d=PWT;~gVZrLxp=q`9#R=7`3kCICmn4lli^BUa8uw67Ha2M)a{5w!cXd5 z6+ip~px)}iM$`OrNDx4zFj9YU@+(D)saA^042HH<#gRhHih~BwoX8Fl50PY6Z9=Y4j+ z48)50(SwknOynhQisNp076Jy9hmEkGkyjNsrj3c&h*G2I9DiyK66CYo7L%2WoJcei z|IW`tu25txT9o7?6$S^~ot0k@L(}OQiH5JpvfWN}Qd;I?_NuqjJTe?<= zK)mEDYYs&S8=0p#%T5tO%kFM}*R>dFn$}VeJSs+@^BvC07bVDwxxv4Uza_|BVYhM) zzfy$thpK=f{x<|Q6T1}>^bIk5$jehpQigE<5xKJ$UWU*GS+MX@lp{&r3+wl><%rU> zcV;9lrUJpERc>iatU%f&8+wkIE0Kmbgg`L060v``8h+nQuFAr=gM`6ZmFG* zWJ4ViknwZh=6XHicf^2AZmUNYKlw?lVz?U+(G&~qSC1Ny4js&0=IjRKRU!W6vb6>z zIJI(0Tc{C9=ri_Nb7(}|@f_6yzBVE{8{K3qJBNJxE2Z0kh^S67di3Hp#&;_#6#`vp@R@Td>j)@FI79o&c9>S@POxSzj?ah(2rngUkg@5^dr96+qywt`w^pC?=y9K z`;h}NT@m)3euVbfG}$J_07CwNnYKz|0KqJ}Q(q6(9Y99&^JVm12aqQkUfLXQ2N0Vp z-yd(44j`72_wQB@4AabFE-)$KUB6r&R(raG~ zB86|S&IczABCt`L*r<9C88=nUVV@jC@}hdI_WlkcMpNpsRV+hD<7t^$I2ghpE4S9Y&5Gr~hlT8AhDR(vri2h7qUe{7K`_!^jiqmVEy9VZ`%8GT2)kM&9Vq)i;xj zAb6tO31NaGNc(g&Fww+}AguAt>;euWh&xNf%t81FqGNOIyyeRX(rM)r71lF?Tw+Px zH{HH4kbh4sNHL12*l$%YAAzN}08{&V*kTI7m;}hneNV6e|qyzO65kGdwdT06*86XG` zdWiXnki8rVlSuuExXr1q{A>7$R1K4SYK567j|nSy`JyE3R?m((Fw{RCOGA zIAlYE!L=DjHgO8cI|Ij&SlXbB_^ff{vX+>`lkek5O#P6#`06+^UqzI7Mm&LdIck6D z6_`M-F`oR6*O)*;7k_XOgm(xR8LHGxzm zKYe0zbrNCM;?j^7pG1mx-CGIoP9g=KCm8+iXOjq#w{?0_>?9KK;su*y$t2Py@T*j2 zXc8gjpL~LIIEfs2+*KQ*nL=9VBkrWgP9ard+@3b~rV#n?DJCV)DMaA6-NB{ADa19u z>3K)x6!J{v-IKKODP-?{hP%_BDWsI+b%_-7Gy>X(doH1-5qS+W?5N2!qRfoheV*bw zjZjgakJ){kMmYFxErNz=gguOFn|N*-S#i_n?8coze1jP#GB{@tyK={Lm+Lb~g5T6* zS?d|(#TG3-N#G2!tV_E+o;ic)_F1f^x6L4q8>mscl^I0!@W*F%f>}iK=x8s9dKT%9 zz^blt&myxj*WT&M%p#$fGF+SL8?%Vv(L9>Xa26@#cbb&=jy`&ms@a@*ad#T-ddVUfOD(MbH7(;dg${BEx>4qpR0uk@%h7 zD6ZeL$d!tZZ-Oq*A$D;Mnj=FgkR!?ey!3Rg2}?tPP1ywA+Ob`DS{u) zA)X8}Ub_x+$o#pbgDT5C%Xviq=I<(F*LfscUyHFWU>-T=9NT4mJC9(C(%yt; z%p;-n+E))t=Ml}3hY3cl^GH@s#odbGd1MHu$(nb09+`gTY8!PlkDTP#S?yn1K;(94 zt##=ZkdSM5^YwfSh_IdJR{{A2gE>< z$e3h(s6glfqCdK<8Je_!sIMdl?dC2Z99nu7Ce;gwe7|F4OV(mZ^ zv$25Sl2X_l{arxvMBNq4$QF^R_g4s-Sr!rN;g|IsqKn8`$+~wWSVYXb$y81;cNURH zD(iS=W{XH7_GV?h(;|{))Tzbew}`yy$k2|BUPLD4*6#jJT|{OIT<@3^Eh4*H^02yb z5g``V{Vg!Ch)}1bg}j|#MCi^q=zs4mBK0{lsmAzA$jA4x&GpxokocKSNfGWPMA;

tr?!rI+mM`H^cf17L5hmE$|?8s77H`3q-hI`+0j~3lPCj zJfzYBH_tBIe|S#|7=O}b-28H=1#G-S-fnej0k0#IHMhN6AWVaeV<@l%BCk2Gu!wGf z1KII^zb3Z8H}m(w+p=1~E$5X=SaA!?cIC7T*0ey%c?R9Gw*^)_mue}RZGrX|<=%>d zt)MBIyzq`#E3k|14;+$jg@FqjN-gQF@cK$r(zpGsAibF3b78x6D@f_etA)F@g0Dph z%gCixm{od{VHwg2EH;}r7C&r-^37R6N@=Zd;O+0Oh}>2vJ=9DYD{qC(45=!M=2rNU z7`(H1xD}EZEI~Us+Cax*q&8|<8=P_H7oS|;2Hd8a(}#Do!K|;pR5U@T>d*^{R zXkpO&+N15-AXB*M{LJY#DCZUVZ5q%9oI5l(6h*Yb(U_`}JD;|}LCu>vF_~@foL!ES z^=lir8HTD^RJVai`#y`p&NfI(^G@3_)dpKmTJqfIYX=sA>ci}7+95OM&odj@cKDDh zU71U3hw=A0LEClP;i$yVefQ1Vfl+$OM1aes9Ud&&`^}!w4r&dBe)%`sA<(PvqC#vt zn74;>#3r|cbM5b2+&S&=$M@kH=hAja))jgBt-c*T`7Nv3IM5D4T29yQvUY%7(F5QX z?f|hRF%u^wI$(2Vc*c)y9U#88;*5fN2Q&>EiG&$;fWM`hC5LSX0|p#-J#jtR0Xvl6 zO#Jfg0D;W@{mNk-uw2^s`R(`)2uPS+#**FvHH|VxE}uJqJaa1I_pc5xd(T|;(AEK= z<6nrFu?{%C`%~luS0~&b0}lI@oj}^hB>a-;gmZh}>(W&^fyLz}YrJ+RJPmk|&3>p8 zw!FLLXX(i31fL&X{H13*A@jgnZwTy!^2-~0AK&eSdwL}&_!B!}xv1H!)B8@?r+qp6 z^S4gO_SR9_Sf9}Fz<^a>uYRA~SK{hCE-S4H{ zU{}V^uCcxwe%_JmN>u8GFr(3CLCtRPGJ87eY|;%csUJ2M9_xk(LVS!tJ>3n%dOM1q z`FF#s<(m?gMs!1{s{J*uC*AOuzCWNUqZ{-lMec9^(hZDU!H#=X-7wZ^vqrG98yqe~ z9`c^-1|^-b2aUWvAaeM`#4eE@xYiKv@^EtxP@>ZGXQ(}pucv3|wzmfiODApW5BC7~ zC*I3@oO>XRUdwo%e7*-bzln!G0X;LcmGZPN>j>?*?(JbJ&2zugOewDSrYAN4{r zbvJe2>s~nFbTa={elNUiZ?Rqes~4=g3dB8Hd%-46(dOG|FA(v=vg_FTAi6?t=Rv+c z=v?Xb{l=0$kdn_G&JpQ@W9EJd1Co6ZD)^*TOtudSbu+#9ZtsIDPd5nq)B0eHVbQZK zdsiPcwAi!`@9l#lk}Wba#(mIbUOI8myblh3UTf-qv=7RsPwq=|?gMB0F@qk@J}4C= ziq|vxpvpy|%Ir!XX#Y_>b?asyeC$k$%em7BXY6>a`(pZ_vS1=b^@*9hBtp4E$D*+fty!lO8eliy++f4-^>lVyYRbjeIMMAzt4En z-Us~lf))e)ePEk-s$6EG4|JP&g@XxIHd|hMP23z~VbLD8vroMg%6P9~wHP#QA z0kZnHSq6aX#icv>Tmzu1@bS#B-~jN+UmjbxascZ4qNfek4Zwocrl)Vn48WJ)%yr}n z1K<;`v2RLw09pcBpUJ5XK(K*Dyd^_>0O+kF9&6?kR(J@= zI5O|=VqIh9J^)2}R^P(T4uI&sjbdLe55U%#@g;q{Sn;jS_Z(P(o5Q`X8_)*v^0l~48YQ7uXcZ(832x?DZ^QgK`=71;#3eA zgv0rPmX^y0VaL`@me)lG;hS8O-c*w&WmieBWFV zp*0BZGx-diBE3O4>RjzKbzl%KNv%@aXf+6@Rwz9-vLA%wpJ=YvTnAzB@saVoGlNiE zS}!qkaS$x>lZLll9R#zkMIOh3nH>wI-V2W$1mytTlb;_9!b#sdqk~Te!E+?7T`YAF zG&hcAnY|r^Ja${A7i5ag{g ze7a!$5R?hDEmznw1hRkgl+2Zez*#tnBYq6uah0-61&4FdPzBG#P?J zoV7d-7DEsls1y@@jQIqPi><{jL$IY>vT6F%5bRcN3f*{N2*MWmj2rq7LDLDFqMNsd z;9Ke`*)QQku)lnq!g_y*dA^@9y5q?ZJbd}X((&aGyj)>${%*z)EGt}W@HJ-$7W}BJ z87&-wkb>ju#LI?&F=XCLJzO;etb+8{cN&I3J7BQ3q+TMk_LJ z4(lyV zE%%Ioae83bRf7?zYTq37nQ>?Y98wlqvf3~kO7u3VIE}z@`w|&fj}cgq=eOmN_Xt$_ zjSp1#jDRZ7qEj5#N5CZcdFHm;BXBj9%rcK2fijht!0V4jprbHK`b**nBuT9GXL~&Y zMZ%oVs97U$)I~YSEq?@d3P!5L7mq-8i9o`yiV<+Tbh(1TQ9A;>gS7?QT1VjfU!n8n zy(8fM;0N!`(Gf7SW4rNnb_9aAuKUH#H468*_*-d$qj0I~`Fr;jqp+x=)ignD6eN@4 zvnr)WVXJm|8>hl3L>)W)TTz+$_XPP{7W7d#!@;xq=I&A0cfPyy>;6$t*duSiVLA$& zPT7omg!L%!MEaz8IE+H9r)}tCw^8OE?dw06f$)gtHNxXM}fEe2+z&#Q3$oVV^c89?4bX?Cfm#?2(3sR zqH&BtFr(4!o746XIz#%3t^prp4A6gzVb{!jnlzoDq{dD=`7LSi>H@ z+%y5NFRq_x-ZlZ3#>+RapiaQn%11OU^$Ez{{AshN&ICLn^I0AkPC({z`uneD6Y#20 z#(lzO0?N;wBX%;5Pk_ST1zqkZC!o*y)!D>z%#QiPyPJF`;KlS)`l=ff@Ux@u)ZVZO zFmF13>-@b5IAprvT4MYJIINPit4f-HI%NxK_Vfv`-4gIa_QM1e_HDj;^z#Jx3oCj* zESUfv=Mcl1U(6k%Y6~x~n}7n*gKze>O+eeyk7|rdeG_1I*8f%7*aWn&zvORbnS|aP zx{fgSB!md4x`NOoJc`NmwqH33(YE93gV#+$o=0rrcbQ3e5&mZ_|JF&c7kTcbsWJ)1 z35#+Vs*|uwg?-?a_9WPB4(9LBp9Dkc3}KN&ldx5W!KQO$5;9M9)Shvi1nUG=M*Ji9 zNpKJOZd`J95==^q-cMYf1RW01Rmy>r;B=MS`ef)Nj8ZGYU)-GpqwhXht#OkOuet2U zn#4(%84Rr0|7sFWxONp>em@CXamr6q^Oys@_VTL#HVNLs>B{^Slb~!j+NxB;99ZJ~ zE~}PFC~5Eh9n&)j*{c~FP5z8bg6gS~56fmJp}jzL^*+uiuv$22?<+6`Ez2jKWh|e9 zjTE)|PSGj2)lkJPB0U9XbMmFM<)4^39G|ataKWmB|`;>*<<;#li;{t{t8N zM&N`(%AV;dFetY^e2RS<4h#*j-RGZ%zMmbB3zkkp?GLRAj1(TuD}B zY?+3Mi-ge|rD@Qot#;}qreXaW<1><)(;y!+>us`c8qO9xJne5X4bGE79_bd-@YurT zMANZpfYO}f0w<;+R#(tz=jmw%tlCb1?S*MjwQ}c+^Ph&pzmJ#K-I|6B1?4-dBc>s8 zMT4Q?{b~4X;K_gO$uuNpvE_VtISpGiqg@6vreSlJ-b#s&(@>nLkhTBIG-UdpwLV|Q z+|b%Xd{3&TL0zNuZE3?a+?e=skEL@O+9KW?!yvPrj?^OWNyaox)$rxLUN{4x+b*2# z<(+}c)If#xOJ=}vl4vy*nSr$P33snc&VXC(&;41OXFv$5Il2{RKrXNF?Mm7VI7FwN zP~SBJ-@C0g9^X3yKiJZCqX`ja&Lj&v^#!UA8X@@SK6I zJWH#x8RuspY_U|=u-^>4IG?~Ob8`l6Y=5@K^v(=euOn#zF*88Qp0awKFav^W)>rFa z%z(%4C2x4%%s}G0`yC3|Gf=N^WYvj+8F0Sw=3V0V8PKwluxqcJffVMVEQ$IVc-?GK zW!gT&yqKkV|5pDDeCc*^%^jbCcRelIjFAPiQ2RYoc0JE5T%4F#X}EY6WQx)SF07h` zxO=PuPsL}UBD;He*+%B`?#gnrZ<~e6nxU6!lv!XgQ!~GyHVfB$oR+-ZGYjV};wlFX zW?@Uz(9?}(vtZnG?3#tmEVMLodj&hqf+(%k^OMIccwFPZG~hi8v6j~uQKG)H5T<(K zllJvl2z^pC?tXg~lxJQNQPHy?9bZ&a@n{yZ*%~I6Br!W!e^7M)>si=(R8!|#)+`)s z@4xsde-@n2!sqegSr~YlEV|{_ESwV1yJlHC%RKhW@C3Eaf5|FL}<`d$!DNLQYCAi0kf}ue>f-CxNA)r!1D^9`Z&n59-MWmwBlAUHo!o7-?3}5(y;l2>EEi|GT_zZkkWm8BlCK*6VG_W zCOEA-bNpoSX4n|};^7V6Ef7}cE&XM`9Pn!Jn5O8-Gk<^9|8{Cx9^M8|o;&tJ0lF(c zX+Jr-748nQ_`Q(W2Ia1I8~j_g!HPEvitoKp1PXuSA!WvP5G=o>AIUIMf_e_wK+SDC zz<-j{)oH~}xHr5gIFel%9&dKvx@J-tWXAq3^chzHao-E+yBAQw(@!fRcQFl=Zo<8V zas+HWSAJH}7@%!LvcONBgj?!%emk3(w=Yy41zx3tU%IV9%MMk@Z5-vv=~9I)rmII{ zpQyoLn%ODFVi$Gzs2aWh32hhLHc~g4U8(`;sxd0B$234(X;(-`izf36v2Y7ly%vnz z`o(A5wj263q+E@f)`p>zC)S9I>o6~j52^{6>;)6`+r@*?x}dnU%SB;qAGF;3Wc4yt z4`=~)yX!so!^?Y(>~R|UAd}ZrvsBUm4mF=;D`yVCT@Bv1bzcs`^+dT#-b#lc(#Yo_=La*Wu%BA;=%_h_S`_!6SZ)a| z&5KN^d{)qubt317pf#9BIgIzNI|3>8Ym4o6*g$Qy>F^q`h1Q1=$slwTuKkc#YsfS_ z3X1I~GNu)5q5hp!TH|3$mE$QANIR*ORB%6<;$7WC9{wdFAwD%o z78GO&`8^&d2l8SUT~TA9@5#Dt)27N?;wQG+X74jD`jRAdzNAh3bepd5Uz=hD=|c1g z?~^MR)0;;r=UW38Gw6kXdk(%kAw+lmJas(!EgxN?OPclmt%Y=T@%;m{xx*y;lG@V( z%z>|Zd3>p1Z6dv-KHd%=DkU|3(}Q+bz9B373WJZe-zJ?V-x+Iigp#oy!yTcDhGgtH zx35d;b;zAfk&&Z6&B*cq7O&k^CS*&kQ_+vFdgMm|;qdmIhPz1VvYU=Q8WTmGbK_XBOq^&5ZG>)|I`{lo4BrT=)L zE9YBlLhkw6tHdTrz4T1f&XU%;wd$jx4bjp6O0RT_0g-v~U>H}Q43YHR!`?Wk?;mWp zULSTdvE3=g;1PE}jn>cPC8b@HtL;k2vwN6l9I@r3&@fSRO{kh=$tp^7O|K_K!p=sH zO0<&uws3u67wjQNl#6N?Iy2qdqVPYd?PH|B-U@4}z**+rbXdtgT1fwOZO>wt2~K*{ z)#bnLJMz(uDi@1w;T5F2{7qfB@2L>IY{;mCqQsmRai4Q*WtI!k*Ppt{wMKxCzDx0g zM%x=!dSdo3>l>4Uq&I6zJxf3fc}(I(=)k(J+AL0`E^>Q?s2U(8DUC9*`yvJeU;8->FzWmeTUSJwUGVad2z5y(hTA8w^R9^>(t1$(2&@1gQ0TnHl4^?R^75)MvTzi zqWMty91CG<{GJdwQ~S@nTqO3}hCKJTS2^sL@9NIX>)d!>lQox_*D2UER{fv6&Vs18 z>;K8?WRcn@`s4rPbso&Q9TKuH;FzqKHf_xzeKNc5@Mr%-bKmUKE(aDgPQ_%~T)(jH zy>~|T|GZsO8v}yoqA|TTnLi8hzV5dD!6YM3lFyLv=$^@wc}yqR1a6&*Q57RP@030f z3}PX6T1{8}iT-EnLnLlQ{x@N>cxmt>!S>UBql|h}rxRv)w_= zb__9_6JoYj#BAs1F&h;z8{f7$%*G!uhuJa_v#}hV!)%F&*|HF`86jqSh?wmXVzykw zYy*hd;t;d_K+N_IF&kGLVm1@RY~6_2o+4(8N6eOvm@NV^+Zn`cT8P{zL8!_8Ly*bQQiI~j@Fro%$9pZ^?Y+l82|2r=7T#BA&ebC}H% zG21)DY)=uht@573Y?6rC&LC!^B4+!8nC;X&X4{0A?J;6D8^mmz5wr2%o5O4*Vm2;| zIn4G9F&n@39A=9}%ytnmTLWUY_lVhA5wkr&%(inLvu!}k#;!Mq*-R0$9YD;s7cpBf zVm4{SY%_@2o*`y?fSBzDVz%Rm*?8sWFxz^>Y-016Z3r=27h*PV#B833+1wGc^&nOM9j7XG211?Y>yDLjU#5ej+m_h zG22JPY#erTnC(1bwu^|_&LL*=Ma*^&G205nY`^C*n;c>`Kg4Xdh}r%iW)nfoR)d($ z7%|&J#B6JhVvA!hSM%vOw;tqd`n24c2lh}lleV>Z?WbC@k2G21G{Z1)kf^&n=u zhL}wnFm`xEe+ZM!ZQ;6Aa zA!a*^m@OVLn;&AfX~b-Hh}q5~X4{LH%@8qL1!A@i#B3q+n9Uk7TLxmb6^Pl+AZB}r zn2nB@%?2@>HDWd$#B8;Q*ZOg*nV7hL}whG21P~Y-Nbq_z|;F5wi`*}frWGe^vpjF^o}Vh*!K zB4!Il%w~j`?FC}C8pLcZh}l@5&S5q>Vzyg|*_05oU7W{kUWnP=BWBY=%+`RI%^We? ze#C6Dh}i%!+e*Z2m5AB65wpEQ%=QT}n=@iIj+2RnhWg%uuL(C?JnC&@Ywm8IWyAZS8LCj`8kJ+3Nv-u!qyN8(VK4La&#B3Fa z+0G(nRwnC%B*wsgd7#)#Q2A!ZXs%*Nq2huMM=v-QnmwvUL}Y7nytA!c(&%w~<4 z?KNUHZ^Ueth}qsEW>Z4U#*3Kk6JoY}#B2+*<}jQ4JZ7sx%=Qd1+aAPhw-K}1BWCkJ z%+`XK?KNVyLd0yV5VJWWW|KzDwhb}cCd6#rh}lBsG20YkHUcr*3dC%ph}mu;W@|;v zrh%Bv8Zlc8VzyI=*`6b2(?-nZgP4s6G22VTY&+*M+fBr59}%J&4)75VK_>W^+Q! z=8Tw49WmQO#B9BY*<=v29Y@R-iI}YzFY>Vm33xY?X-FP9kQj zK+IN$m~Ag&wtU2F7Z9`6AZFW$n9Tz*+s=8+rh=HQ7%^KTVzw5}#BAH=G21D`Y$1r*vJkTgB4#^>m@Nu1TQ*`g9>i=0 zh}nb@v$Y~-V?)eVjF`R*{l(>aUo`- zB4#^-n9Un8+g`+M1&G-O5VM)hW42MmY$p-32_a^SM$Bf8n9T|?+egG~L5SIk5VI8` zW~)TZwiPj(J7TsR#B9ch*=*)9n>Av#IK*tmh}pIxX46B=b{8>Q8)CNah}jGgvza4i zdxDtl6=Jqkh}qZ?vqj8fwlTzP9EjQWA!Zvu%$9?gjYQ05jhM|7G22?iY_|}znIdMp zftc+EVzw5N=H2_a?+Ld+J8 znC$^#wu^|_ED*E#&ttYO#B6nl*-{a+JwVKsh?va_FrPQ+~0h}qN-vmHjv7KE6s5;0pZVzxBIY-Nbq zoDj1CVzykwY##HNO%O3#1!A_dh}kR=vxOmMdw`hDA2Hiv#B9=t*_08p-ABx}5i#3& z#B4^0*|s5O_IAXR;#B3>u+2RqiokPr4f|xBEF~?^O(&YF#*3J37%|&7 z#B9CunC%E+whqK>Lx|bJ5VJi+%odB7tq?I=8e%pl#BAn>*>n)IIUr^;M$D#!m@OMI zo83HSJB^s_AYwK>#B8pJ*=7*49YDeJCt@}^#B8C6+0G$m3q#B%j+o5^F&lH8$bT`LBVsl&#B4he zv*{sbyE>2A0ui%WBW62~n2l}s9A={+W*b4wCXAR(12Nk%#B4+VU^Yd>Yzt%NFk2~N zw&nAf?IvP2Q^agmh}i}Zvo#@Ri$=_*kC=`1(i~=qG23>;Y?l$Uxg%yvM9ij& zm@R%Dvxy>R(?rZ>jhHP9G230lZ1srQQW3K$B4(3E%=R5ITNYwAL&R*m5wnFLX6r!A zwsjt}EvTKtY>|lB{1LM`BW8Ptn9U0@+cLy#TM@IJK+INwn5_~q8-bWj6ET}UVzxNM zY&+<O?w(mx$Sd5wqDNX4689=Ld;f$n2ilFTL@yd%ZS+&5wn>hX6r)CCW)BM8!_9% zdCV4xm@N`98yzv5EMm5gh}onNvr!PUSs-R}L(DdYn9Un8n+9UG6NuS55VQS%zwQ6` z$M{UW5WbMJNscP4Y;Jp(8C)UGMJ$lM@obFR&JKJ+DaxXN5fk#{DHzUCJnmQ?*G z@#o5uE$Sg+FC0{^F4?3TeSd&3?0w>u6f;C&fXrpCq}eLh`!E zWO~gov-X<$oBR~9$7GMhcMD==3c2aSlkEn~>Zlnt)hUGWBXUDomdO@obyV-vEp$_6 zb=07wE5d)6)lnzjxBPg?td2@GKD@M-SsgWS*W-%g+YBgV`h_J=mKUp&=U$%<=yBdg zZcKY%_&{2M3{M*b@3h5aM#B&5PYD$Rs3qmm6y&cG*=U0OD43rR-QnSGO{7$ zj53$}U_u)M`f^SlCbYS?!8pp~|3DkD_%qj-gz*2(_0j2)cq)D~7J_vxwqJjd45r7A zi!V~i0MjUW+51D;keF$7^n*nnjFgBS39c=KJ|kJq!zYWu_B@#kL8TzbTl{I`&=1hu zug|v0rUI_^Ocb*>RDx!G%*pG{)ll$SY9^AW2DW=yv|Tw=2bB#Mv`=R;*Y!_z{N%hP z-Uy=gg*8V*o8Wo4!0&T{E%5X=XSl$ z&1SCM6!s}J=Nal{#_vRjYc%&k%XH$4JC*&gW`V*(m+}FK9z4RY{d*Af53eE;T8H3{ z0>xxG%Lo)N99nWrVw3^Yr>A5`_K!i(8x`f$z;R&Mg(z+Cf@Dw+ z7Oy=oNnX(XY_HTHMVfd-8~WI7AcL+tM;u%#Lk@d&bIx>cBuo8&6~~QjCew6{X8CAx zq^o;wM=a{ss6Jfac!dx z+4sRindg}vDG_#Y)lF9eQpc!R^T?Pn*{zj-EjaKHxpUc@#b!S&Np8KUcsmYTazu@1 zM|CEzBkA@TBqN30$l1Qoqdtr?q{sT7*V~xkT26P9JXqv|{v_{A$Z@KaoY2xce*+FfB)0Gq$k}mL@eD=#nVyz z)ii%R;Su!pqr|Gm#6J39Lig<;;#Xs-?Us_G#AN!Pi@(3g{&V2_efx}BNM5?*cFXa% zqrCg9h(@Y-(kt8j#Nof`-A9&eA>Q-^e7Uz{mNw=|`W$l1`QOn)Y~P&!@6r2`xaZc8 z&Nm*6_4oXHSKMhJn`n_kE^cjPNxk^TDrVr?adm0id}iR<+r}H!@MMxaSG~sdFEen> z$87J)wLzS8LF0yHPCkq1^6SK$YM6m*8kQT4-ZKN&#E;hLo7b(S>rzAzwEfT~I%LMRR{tk(Euw#y={aWL+9M@9Jx6BX+QRo?ISS0c zwX9&fomgayZ7!4BQK=7 z)NGKyMApsfjm9zq*R*`?P0lj|*RqqNRePC%Yp;FEPW)iLCgLu4J^RcIT(j&SIq>&C zfom0N_7yvsfop-je^1J@X@pBmcypTIR`g^#?<($=RLSz)L4cNr9rhZnvK6lcoR=s)i4ZcLfV{j6A$ z#+0e&zuZU{Van9uTmE4MiLInhf(GBusvh!c^1aa4OrXf6ulQ#=VVu;>l{ekWl&LHG zwqLo(l&PsQwcE;{a?uT43<4#YGIi@48~gegLAn6-bkK z0j5liCRF9}nKG6AN?p)>rcAxGF0c0Oc0-cj*|;HxDO1gLw{r7Tn31y27Ji^HW$M${ zd-_Y5GS%)GH>Ke}GF81sU!a{SQ!@^`a|SSF>a)#KEP70t`nPkX#~G$frJqpZvS7+o zPX()W%1oIWDYR6riYZfjrl}Jh|3{{NTC3d>a`v75b`@p4hQz0jvWj-Qs?-f9WKkEd zx@s#PoF)Ikd9qI1itzrrxt}xH@PFSosVm=OQAhj5DIY=_vfQ>u(lWc>%i9vo&5{Na zcDBTO?Y|caVJ{J9-^BU4c0IA^gPbSJ683-k3HX5hT#2jx8KR$aOy%XUp=>~vFzltw zj@Cb33dRD-#Mk$;M3n6&gT}~8;)}R*z4_f&|J%<2!5aCokQ*OXhVXqn$ao{U(aVWe zK+fj!TV8J|Ay22Qe|+snIk|p{cTmZ#nhdJZxOtUX>3v!;Rrps*EBQfW*`c)6y<~BT z%6*+>LuBF3>ps`Ej*~0PYN|CBvd|wrN@*_}T}c1hxX>YXH5Z*#s!(BF2_Jo6c(>`s zF(LYeY~ht_U6#;49P*((^A={%C!Ch*IH?NJyR8B^1cLbJM~~Nj)beMi$6fhSpnUl! z`FG7CRaWsfl1EA7)(hQca$%Ul_H@nfWXG2sFPJyqAzR3#(U~C)XMG0ng9&U=~^C2diZC z_l}Yu|0re`()?{m`X>jKjEu-_4x9&x+zYcq4Q{)MDCwt5J!01r<67*N)-T3rnNkp+O-DR>*l9>1u zdsJIx;XhN$-`+72!jZ@PYIl=b$*X)b>wXvS#^)a{WHo?(`9eN5LaREgsg0l#16|9` zMo4r2)BKwKMxzkd&myT`@?z!+k~)3Po=e#j6^Yp!qFCs()#{Dk(*HQzBj4R$JZC!G z9!Wj7TBgHYXvN|B`d^28V9z(XmXNYSdzB09xgTHurV@t7YuD}LIYiLvwXcghs}gX2 z2ggXyQex(KfZPq)s(=0+z9Yd~LQZ~hRUy;P{?!Q@AX@yN^lo(=B0dSenl``rmndJO z|5w)K4Po)ygmHaA*gyNj*B{OjBKFlmMYbfFTXN8fS-R8t=?GX8`oCG8S%Vf~emJt2 zC0>lsFrN9iV7U38(N}EhX$Z;vPxI#1ll;yP2-$<=mn@w)LbNJDu-N-L@%-ZAJIPTt z#7fO&@AOt}{%7=i=lEno?27DFF<6lrHpMvB%} z$NxJSY(8C{A-vz*R2T<_UO$+6M106xA)r>DKoA1@%6sgB2pP>60+NP~L}S~wCc!4n ze@3@H-MJ+so!NZs-6DQ5St}y!jpOSo$Nj{|tD9;O42%@XIbf{I`aOK=gNxoX7aI(D`i)EC$qZ!Y0)z} z17zUh2shK-Q8J}XPx^t}EIFX(Q9Y%=MsHS)jC1JXre}?@n? zPo`mqy-NsLv@5XS@o)smT#6sJM(Z@m^QF<`(h)E6{a4-zCABl;o2GRXuPRsaS(3s? zisUg;ayZJct6|?EGL-dFSO>cv*>GU>mBw^B8B%%9m!EwbS${xC)wErL)D|-mA9Ywl zj^@7NFb`%URoDWTFLE4$i>rS~sxN5z--AS|^!JZ*@2W|qldT~GrH&<-OT|-X|07wuK?5ez`D&Ze7{; zHL6O*t4=++t+r96la)MJ8{Mr^k;v#-%Qd2MQzrCLuNezPo&DG)V|6adEtzhe&7(pT zZAR~(t{bZ;-1`D24@5{%tP76bv8mro*%TgpCW6>bky}yFdV#rmN;OSF;C8(RrD$2t zF~w{>3jcx<{afA#C`vWvE05J!F({@deoqE)J5c;Ashh6#xlwe>K7Vw&be2-BoOojI zz-3BK`S|RT`anv<+X#E-+E5D18l$@BqW35+tDbJIN_|9G5U^4rX6FkEKabav-;S>- zjTd+cPvLBeOUXF7itRJSOn92M%KbYZ-~G1=d$LP}vp<@QPiOr_xjI?^?Ke2ld5*a zdPQYFji|XV!xfhpm{ad3iFffcIBco1@Y5k|)p6>(+&-aJ5f5tAj|6G{oU>G3M|y?j z+RM}re@q4T=v<|?^6mVbyeOEu$M1C5{*U3*t&jK=^#@|92R6-|QS0^aCPv)jN~*v| z2k3ZENA28whwZd|3-$Zapvzs7J=86CQ?HrQhNwAuyBEfbPg12+&8$cnR$6~(ptJEV zF4}nBbI*+mg0y$~dN0OC1C0eU}=`LkS8qK~ho>7skphjEd+ibYwkyD^|QM2yETxO8fa>E}t_7yY4M~L!=G>NVg zqo*97(C(_cp43-QrrA(_<}WHtqp?@=*t+>-(hOtOejaVirm^$Hune8aqvf67xo*As zR~nz5_v9|EQrg$4q$rcxpEQ?KM+QqDRMSrXwYi_-Tu)o4k}BFi)J$`$$=oe;x09y% zXrDcwQ9rGJ`=NcmGKOhB4AGLBjq;PUDi3a#<5Mhz+(g~(Tk`D0x{I}x{3AR>u>SD{ zVRnK<_#35&Ci5l4P45#^g|L!H^r`hU08xVe{&0rpUI`+x>9Ag1$p+%rHla~J16d-O zT}kye=T<^9!)Vn-aV27a-_WH}j6&$@Zp*qBK@vyS^0DM;suK)Z3()6!uSHP%1W#m& z=@M(B?}wG^=@ThMY!OtdF=0X(Pl_BcB?65!W#n}&iTrHq(zui(gj-R7X4%JMg!4-I zVXkN=g59a;=f`cXgubDij8eP@;a}5}6Z_#b(PDY+{m|WWgm%RSaf{@O#67MHcjZ_5 z5>IlYf_g;*7{s>Y!5=r&UMFn4@0BHG1`!gck3Q{Yy-hea37xwAY3kP8%MBidFJRokwDzvIP`{B>N!!uLw^+~h_JI5bld(PHjfkk!1MVt3O z?wU6c0X!kfrjt#?RRItFwC+}75^Ozl2xbB6fcpV^$9jnG!t3g7Li&joPA3L;#~`u3 zNB4Ox0hkIu^po*g2xR0|cR{i+SXFO3YWaI9I1`62 zzTdtA)VnPV<3Fzg4Z?DiKXeTo8;-SO6zhsX%CpubMT6@g;!3{p-?tJF_O;$1JwXbX z)x0aqGB&_tzu2q$hh)GwN}j@VU^9FdDyUl5wgp07`~J==lZU_X$u3rFD~$Uuza8{z z8wA~YTS#Fpnceh-#&(?C0eKnw8q!3Rfx1R=jg>G3R;^#fWg|rctuR+c+yeuE#P|DE zf{RH|JNQc8ZGkFtiL4(hLsbpLM#?n=gVdpu#cESz_b$lrg{2=%G@m2n!#39ypiPnA)Z`G8slhKE@ zW@8NZMFtSl$6#;P88HBk%Er*#IzzAzi58aoW(1|{A1-ytG=}9@UwMl?HUY03cFmOA z2Vk;1KCATNL9mc}__*HD6bd&Q7Wx_-g0Q<>hRQq4;B8b){Q6~w!QDMKjOaTIAC5ir z|NhAwEI*qiu-~$Pg!-z+eH$&|p`gU_E&i6w;Dm`AQR4NMAeO_KJgjI1zw6CSv@csh zdA)U6)K@DIw)gw_klz}7x2p#4+HDO1(_in*p0ftZxu^HSncs)QZtc;TF13a#zbu=U zK5J;WUM4Qhbp%#1++Mw0eguTvB^=yEkAN~{_#9;3O8kAS;K=;@vl~Ly&p3|&qwvD0 zhyt_W6-9YtF=j{lWv%O?o-*%`S-;QlwuTso*6=b-Yw$aI&}kE&HE>5wy{s*?0t2&{ z%KaCuAU%7yt8c3nY}xWEV|lG5Jf8W|RN`w1+C8<>KQ~wc?;73*8?ITvjIGqB*er9X z4_EO@=r{~r_p9KA;9&-w%-9qWrCdG^wx-ZN+nkkg_8=Tv@H(mF z<^kZkU(F`@zy!E^e%5YFGX{G0gBJt&M(|2s=i91kL-6FS{<5df04Q!tTa~#CAaTQz z5MK#>$Zxq>6|cG<(jPv`)3ntC_CIG^lUMG8Rdy4n1Kb&V;q{B(a-o%bASq3tY`K;; z@T^?D!z@<|p3zwqD)wrE^np8zpEvD-5PLzzjw|Y*u5|vhk+K?8h`j&9J5GnUX&sAZ z^GP_j*)r;f8bC<5eYoK=8vGH_Z<1%HK$1M7QUuw%A^2u#GWcYtp z>0x}`22cFz(;g99;bp_nk;et{u-0p5qfhx3U_9D+hh1YcGs02|~;7B*kl@A`~wA}yW$qgwY zi2@&TIG{D^r9`SO8{DD$-?Lu93a_u-dYQkH1+E9ko$$FcMI_4i^*bFLC;AUnUlVm5 zA!d6cD}PlC60A=|gBD!sC$gMX^{RjM5U&KsoVvw23B}g1#qWN!5%zNuYWki9j&Znv$ zB}BlMsYOSweota!V^EsN;9;I%wdDw9xZ z5MSls^qQF1(%VTWq!O}H4B3zj+hk%nN6pKUiDv{o;MS6>6$u2(khmDTRvaPk8!_Cx z<33T7Cs0)_7fpEjd{N*^3n%y$SM_k!hZ0^ZKR4LL1Q8khTV9_{x=wWX@H;*q@+Y*O z>5MKN@*$=u`7&37E)rU1lGi>to+Etw*3nmmo+hLptXR~~@F40K@w6{O!mi98dY|0b zbIgfo?vOUDx^s;9uHvqicl`+A{WkP;kgO$P(4iSQQdGIDaFH=_=vG0;@MeAD zQyqt_Qqx{yY3AstjK3C}vE|J+&D;6V~{mrFk{YgHjuKbbUU28a>2_x>FY zt5G0!W$<=dFOnq~(wlFZNZ;8&oZH&t>nA2b+$%ctCwhD>G4pXyX@+$r;b_5g^xmc= z1h>PxqZbH4f|t)Ird^DO7~U`Jnm)xwI4CS0vn^#I{O<328{axXb3getZ}s_MTH;;V zguPPzw0mSmZ}L(BoxXq%p>8IFq%!XeSPeZdapu&{D*owwH;y(0sC$STfp<(d^t#b4c#D zrUkS;v+}HfeMV>ZTBl_*=Pifjx^UJ~1cj7cg{n;B045|pN-tA0K<@TjC1?mm6 z6J-K4J#~FYsb8G5kMI8MvZ`1>%iPo2Ogx^TcJ8E#zG)t$TF!1(xe?h-1&T$3?@TjQ zb?<4;lKpkmm~(qyeUz)D@@1-NnzWQsjkR)Xm!A4UJx+g^<#XW^Rm0uSrYYqewdr_* zLCdmN)O|6>`x&7_NmM`1O*@>D@H+}p1df`fW>ZR7t zr&tHqQ}5LsI_YgJM*aND#65ztoXSBP?-#isNPQheyj8W}q=vZ`{IK(3p+d_RJ#@kx$WGKb^Qd?K| z>@|u`6L(mE;bqFE{rcQH9(z&V*O%r;9&)4bN~`RC-pnb{Gj!4)+SKbG757B>FfvxYgbp5@WrUJgDMyH7jFY0CH*UuaGA@}NQo3v0#VGxv zRrcGxkr9U09~J%el5xlFfW&o)QpRTXYKg8z4x@bh`N`1rM-0=I(`b*fWJYD6-QjwH zXvUaE8@%>0i!qxx)i1byiLov$qKBcoGJXp@?MF7TY#AM{Cr>@#GH0yCzU!Mr>NA3} zo{z<*9AOmQ5*JjdCo?|T59??aDKqMSTyQjPm1bz&t(Ypoi89{Iyt>u#m7n3CdVR(s zZ=J5J|GB=-ZI0g0`WJYze~j)czEjSZG)SLTp@tlB_&|So=+uoamwI})q^TYLR26** z$G0EK!WGf&$7I)H2QumVLvE`#Slpq9v{mgkY>A~G@}swFD~8goKKlzeN%+$X4u92< z`|3s?+;8}o-^7l-57nE@YBi&)IEOD*uIbWUv-oZl3ml?LATE|QMG)!bm(gao46>uLTH1#h5{qH+`=hpT8cRi+_<&_Ws)uD<4!YX;3TSH!ivK80o;--ge89!-M{ zo#d*q2%(8Fuk@z%`_Ld2kHS99oToXx9;d2?e&u_5h0F(#H#r%tnemB5cEm zow|GC2B>6@x>Llr9aPLr!l3}&w^Y5+Rcsb-Ikn_+^@lOV9O{Xi&v->s?o&$x%7_WU z@l^TxN1rauU8nNDCKy|J`%?9qSz~ST>dsVOYEr-%cPnbqZ!2Wdnh~|l+c6u@r9oY= z_p&x2lc{D1-LCj_BvljLq42;$f?5&hcmBf}A*v4no{xU9p}JRw7i}`VTqx$Lv>i@W|NuHN>ylBi`RqkiJY&_!&xCotUDiglLV7o&8F2 z4_s34@am*+zYl%5=j|H`EB5=0UJ|O5V$tr>F`S!4NiaXEF`=GHDHx<0-ux3uDXBb6 zd-mcog|t!eax^df>5_Rphuf+#-ioX6;UY$`Q<+rd)MF0R8%g6Wo$*=9S--=)5Bah_2 z)wzUPBx%=Pk-pl?8YgKX3qtF!^^rPA8=r@in@R9Vx&6{LFG!6`eX^g_o|1mR<^;Z5(JBtF|2?$ar7Qm|OnquWKhNElatxw?1j#K_U{)?wEBUqt&MbnC>2FT@lZqQX#H zJMnd7*?+CCUJ|dw%fBdonokT73$aR%eMnp_Tnn5uiX|?cB=X!{xk8+Kotkzo&4ZY( zg%`Sa2q@(4QnKzLR|04eX1a&K6i}bRe3__LO4G*3rZoJ z)yH^TYz-qk)_5}T3g%02D6kk#^gKrx>pW{@TxCkwzws+yctL|Oey5q~MkW%hN}SYK zezF88&$zsJz-~fL<&PurLK}Gd$Qf#L&NM!(<|`{x(s%$rWC&Yr+-Sz{ynZ2kXJsXx zQVO8{cxK}%zLKbVSSsFj+RR?CFAQJ(=}Ssp+9iCcm+x?tzXM)9ao}UXX%qafiRg4K zeRcemdedW0F*tm!Kk#Ss0Ss?98)IA3Ac&V0W40u%E#WHaOV0a}$8pEszL=+l_u#6J zvyM-8{&|Pn9c2`U_bb7P)W2xa+?$S@1D<97jf}_5-J|QpNL2kvWST;Am{Bab`Q_sEE)xfiFKJYeeWy@`!ISnpeOZyL+} z8z~vJ*^eC$)0>*R*NDxSnQ@SDEyuoQ*|em3I%Q%jLl~N0;uEpwhsEGQ_gPp}r;(?h zy$5zCV~H;0cN%NEQ^~BTLkDXy^0CCqor29v@N)Lll*iVD3e7(^-h)Mo${l+TTgPm? z1Rjgd|HMR#j)e2UKVe*&)OPQSsKeaLTdLhBSAx0Lih0+t@gK(JpxHM|$}J2lcjSDh zZEql^^I)>(?S4lLc{uge8?qT@$Y3h%*?>AmrR#j*Y9bb6-nFaTc0v;4MK(qKd&Y+` zeidH&>g8Ybv})=8O8F7=SJn@QqUY`CjkAa+3bilL&U?Lle(F3yb2$P}r8JY#Pg=W8 zO4F~Quf*%G#gurU?{#T6mRDP{(2)jxk0T7U(GTZU;&Emq^h8j$U)3=g^hG21U!5XB zbazR%r}oQ#s3(Ay`k>W!RQm&>wTfRS3KO#A_jaKgl_8XFpZ$G~r+Xk`!0 zs|Bv3swFo@6S_Q5y0QQGMeQw7iR$OH$zEEh8J3N2lnfE&-R@J(lgEOg&WIujj?4Tg z+z`yzj(=W78=cK#_tMzohQ@9|v`(m|({gN{M7GH28hg$UprN4LI(4%{bl@LQ@wSLwEC6?Fj8P;q@gpxsW?mYUR zkJ4eHo@~eEGfK}_QbZeW=_pkjc}q|95|ws;w6Cco!IWUucTV3t%BM6R`MNwKepWHM zpMPRc{(vH6zV3TKbG_nj>_CL3e6gZ50jQp6yr($s)8^CuHB2$c{HJ8#5pTsc?z7)E zv`;Jg4a7A9q}Q^q%1RfqJ=Y3Gb81xYljK)#$ZFi zu;M}F#pv&J`FsUr0>Vfg@Y{`?zkiqcabi&+4k2M=sq{@j?j_5%hH();0&FC&U z>gbE4ZlFqTb`qukXIG}&l?*J^Pd7$R%0S{FZ|EgCE?D7M%8<3(&!(1Z_MJ!NjC9R1 zg?3`(Y%MdGp<91(3vplngIisdO}dwN;A_oy*JFxq${HzsURk(x zUsgG>x4hr_21~YVv8lNo>n^}E{ zUb?M}_Pr-DWmmOjy3v0t?#E$e!ap-R@yEqv;tr3vV}@9((o$dWTE|4cOLutdobH%v zmJZesxuklhSek1{N3VPkkT%ZGv9Up2ljigEZ7dyel=i9{OCdX|)(q zQ^E>V`pR-lUl4&?TJQL&oW5%l@V>x%E&NwI;pZy(lO)j8(UN?4PP%h-ezjX3y-YN{**B~3CsQ} z|4gId11wN5$6M&~bC{d^v9QU-M=$|RJ>}OfVKCkz;r6Z}HyFY{Z?iGW1SXJ1mO8#c zh4mP+W^Vn4!@mEm8}AL^fjJx<%jY)yDRpXh)>)Rjd#4oIyR5|9pj=9_F}nOPa9_%g z64LECa9wI&?2ufoxU&Vc149!hkh96w+!=vRR$=$Ea{WIdEBtY#8b?IbEo{CtQ6DX9GY#JMlsC&sP8Rr<=-v0Kd z1lI-AaGu2uiCJs65$?`X3G6Ws%GKyQ61vY=)6eFtLnNd|D?L2kIY^{sHI%(W>q_9x z$WhFKP^WBH{^T~&nRoMfpxij%>ZNj$po3D61 zFe|)e-?NP2u?R>^=3P-|n~EeL<^Z3aTN;NIMYbkpg~5I{t`OPJ)fF|i<@Ma3X30lZ zDeQqOt#p|)Epe4It!kCZ)&%Px?n_=;rLzYG#!}al>1&*UwCMQ)&YZY4MYj6E>m8d( z>j#p%w#T2yg&&jFIop{gca#7`fi1_Z*5PbNg4Q|FxfR*!0<)mY?dwSPI^mHg=K$i4 zEynREeI`dJ@|1ud2JI5I+uchXtJL^5YD0=ntI+TRS*s) zTbEAHvpohxAP|nr;62j7K?pqi0rucz0}06UK;Ud6u&*Yd%d_n& z*bEX#;DM0XIxFFzwqe}L9uUO~VV3|<%nM<%HDG5& z;4klX^x)ABFgu*=u}^A%5r*y9ju1S20X*FSImn(1o@4-i?SOF70Uk^MuzcIWfmCaD zo`Env2q)1XeVPO98@_F5K}s~oikE-e3Zyx6L>gNnHst5{x7`D|%Ro9m!e-8ad$tK}@5v^^ff51XZ3)P4<1)A*Iu0W6x#X&!H*LFhU_UY*Whk3W^w0jri0DB(DcLHX2 bZ36((mjH^$Hh3T#2`qc|mkG!ffoT692&A$R From 19c73122c4cb5314813b9fe92db074de82906cd9 Mon Sep 17 00:00:00 2001 From: jaimerzp Date: Thu, 18 Apr 2024 15:44:17 +0100 Subject: [PATCH 08/10] better res --- test/runtests.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/test/runtests.jl b/test/runtests.jl index 16d83e5..bf90fb2 100755 --- a/test/runtests.jl +++ b/test/runtests.jl @@ -707,7 +707,7 @@ if test_main cosmo.settings.cosmo_type = typeof(p) z = Vector(range(0., stop=2., length=1000)) .- p nz = Vector(@. exp(-0.5*((z-0.5)/0.05)^2)) - tg = NumberCountsTracer(cosmo, z, nz; b=1) + tg = NumberCountsTracer(cosmo, z, nz; b=1, res=300) ℓs = [10.0, 30.0, 100.0, 300.0] Cℓ_gg = angularCℓs(cosmo, tg, tg, ℓs) return Cℓ_gg From ffaaae724365b3aa6c77dc18f4f24fd88e093262 Mon Sep 17 00:00:00 2001 From: jaimerzp Date: Thu, 18 Apr 2024 15:55:57 +0100 Subject: [PATCH 09/10] better ers --- test/runtests.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/test/runtests.jl b/test/runtests.jl index bf90fb2..e74e4da 100755 --- a/test/runtests.jl +++ b/test/runtests.jl @@ -707,7 +707,7 @@ if test_main cosmo.settings.cosmo_type = typeof(p) z = Vector(range(0., stop=2., length=1000)) .- p nz = Vector(@. exp(-0.5*((z-0.5)/0.05)^2)) - tg = NumberCountsTracer(cosmo, z, nz; b=1, res=300) + tg = NumberCountsTracer(cosmo, z, nz; b=1, res=350) ℓs = [10.0, 30.0, 100.0, 300.0] Cℓ_gg = angularCℓs(cosmo, tg, tg, ℓs) return Cℓ_gg From bdab8de33413dd370ae7744337e7b749b94a7551 Mon Sep 17 00:00:00 2001 From: jaimerzp Date: Thu, 18 Apr 2024 16:22:00 +0100 Subject: [PATCH 10/10] bump --- Project.toml | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/Project.toml b/Project.toml index c9024f6..16e65e9 100755 --- a/Project.toml +++ b/Project.toml @@ -1,7 +1,7 @@ name = "LimberJack" uuid = "6b86205d-155a-4b14-b82d-b6a149ea78f2" authors = ["Jaime Ruiz-Zapatero ", "David Alonso", "Andrina Nicola", "Arrykrishna Mootoovaloo", "James Sullivan", "Marco Bonic"] -version = "0.1.4" +version = "0.1.5" [deps] Artifacts = "56f22d72-fd6d-98f1-02f0-08ddc0907c33"