From 66db010d0e24f5f252e97b5183d796ab761c8c4a Mon Sep 17 00:00:00 2001 From: Panadestein Date: Fri, 4 Dec 2020 13:59:10 +0100 Subject: [PATCH] Use flattened array --- el_nuc_el.irp.f | 89 +- jastrow | Bin 44064 -> 48560 bytes jastrow_provider.irp.f | 12 + orders.irp.f | 20 +- orders_inp | 552 +++++++++++- tt | 1800 ++++++++++++++++++++++++++++++++++++++++ 6 files changed, 2421 insertions(+), 52 deletions(-) create mode 100644 tt diff --git a/el_nuc_el.irp.f b/el_nuc_el.irp.f index b15b10b..512d62f 100644 --- a/el_nuc_el.irp.f +++ b/el_nuc_el.irp.f @@ -3,55 +3,62 @@ BEGIN_PROVIDER [double precision, factor_een] BEGIN_DOC ! Electron-electron nucleus contribution to Jastrow factor END_DOC - integer :: i, j, alpha, p, k, l, lmax + integer :: i, j, alpha, p, k, l, lmax, cidx double precision :: x, y, z, t, c_inv, u, a, b, a2, b2, c, t0 + PROVIDE cord_vect factor_een = 0.0d0 + cidx = 0 do alpha = 1, nnuc do j = 1, nelec b = rescale_een_n(j, alpha) do i = 1, nelec - u = rescale_een_e(i, j) - a = rescale_een_n(i, alpha) - a2 = a * a - b2 = b * b - c = rescale_een_n(i, alpha) * rescale_een_n(j, alpha) - c_inv = 1.0d0 / c - do p = 2, ncord - x = 1.0d0 - do k = 0, p - 1 - if ( k /= 0 ) then - lmax = p - k - else - lmax = p - k - 2 - end if - t = x - do l = 1, rshift(p - k, 1) - t = t * c - end do - ! We have suppressed this from the following loop: - ! if ( iand(p - k - l, 1) == 0 ) then - ! - ! Start from l=0 when p-k is even - ! Start from l=1 when p-k is odd - if (iand(p - k, 1) == 0) then - y = 1.0d0 - z = 1.0d0 - else - y = a - z = b - endif - do l = iand(p - k, 1), lmax, 2 - ! factor_een = factor_een + cord_vect(l, k, p, alpha) * (y + z) * t - factor_een = factor_een + cord_vect(l + k + p + alpha) * (y + z) * t - t = t * c_inv - y = y * a2 - z = z * b2 - end do - x = x * u - end do - end do + u = rescale_een_e(i, j) + a = rescale_een_n(i, alpha) + a2 = a * a + b2 = b * b + c = rescale_een_n(i, alpha) * rescale_een_n(j, alpha) + c_inv = 1.0d0 / c + do p = 2, ncord + x = 1.0d0 + do k = 0, p - 1 + if ( k /= 0 ) then + lmax = p - k + else + lmax = p - k - 2 + end if + t = x + do l = 1, rshift(p - k, 1) + t = t * c + end do + ! We have suppressed this from the following loop: + ! if ( iand(p - k - l, 1) == 0 ) then + ! + ! Start from l=0 when p-k is even + ! Start from l=1 when p-k is odd + if (iand(p - k, 1) == 0) then + y = 1.0d0 + z = 1.0d0 + else + y = a + z = b + endif + do l = iand(p - k, 1), lmax, 2 + ! factor_een = factor_een + cord_vect(l, k, p, alpha) * (y + z) * t + ! cidx = l + (ncord + 1) * k + (ncord + 1) * (ncord + 1) * p + & + ! (ncord + 1) * (ncord + 1) * ncord * alpha + ! here I try to use the flattened version of the array + cidx = l + 6 * k + 6 * 6 * p + 6 * 6 * 5 * alpha + print *, cidx + factor_een = factor_een + cord_vect(cidx) * (y + z) * t + t = t * c_inv + y = y * a2 + z = z * b2 + end do + x = x * u + end do + end do end do end do end do diff --git a/jastrow b/jastrow index 7b89a67c35ff14ae2343a9c1c66a4526c41adb2e..4ed8e351004c3fba7f31df83aab3af5b39461b8f 100755 GIT binary patch literal 48560 zcmeIb4R}<=6*qo2;bIh$4VEY<(G-L76$B)hDA5E1++fg%5kb)%@>NFta|!{-H*$<#_)Rwk7(J2h zA!sr`UjpRHt12NHCMkL%Aj>UN`LS?-;?q!TNXT;CU0Zv}91V?BnPxeYs7lAe-(D*D zG}QTaef2Wmh=F@jBwa&iIr?T zPOqt+J$B6LH5FMk)wNCYv*zcF%^EvqL}T5E(ZcS&P=tMI(gl|ok8FMK%yl>SsXTG= zf=wxBh6Y}_;pxfLn>tV)d3cduRVV1o)Acd`zquxY>DV}oY`%E(if-sHq2hS_&v(P0 z4tfUuVz1I}0PGy6qp|Y&T zb8SPl7h>#OWLDJ8^EBi*Xl3z;<#n~?WnR?WP+dFM!QAMj;EJvkPX)A&CfC(h);cJ* zX4SRc%DK=mRymu2syhmzPQe;mh0R}H1slrYye=Y9w7jOS(W#iSwz0CZ!c#rQSxu$a zQ&;b;uB&ZyQdQ+n5*W)G<~ACY^Xt(>Xr1z_J>^wbd*+l?*BDg3rmozWS2?e|{yJlx zOnE%z^UFMQs%y(?s;_6mA}p4xhMiYdU2Du~sH`+V)J??*nN7kxVq-;CIZ7}!im6{U zl~v3u^Hz~Cl7w=}IH|a3{CS?yBSwueG#bH1L#C#V&NU`Yyl`si)Pf5xzmrs62GTxJKGA4!hL#*HP5yMwPGUXFZ=TS3zjGntwHYmZEQM zrxe4`be-P%M@i@0=$s#SIOtH^zCsSVTDpk5Jq|hyx37Z^x>}-%yngMH1xCRWkt?HJ z;-o861!1s*&UsK@BikiLem&QdlGDy#!ysL(qeKdA)h-dz`KK?dU1FqZACXeh&R@x^ zmp)8SZ|AS1_m-&moz*TAl%tkff>7n4yCsYGt#{D%nu8hh9dzgX((IuBf{f}~+%6H; z%e95R+S(;XdYUEzE^U_y(vOp<_-${O3DS?3sQ6vhE)%5plc@Mz-YyfQ_m`;nUC}NR zq@N&B@!QcZ6QrLg(IdY{JaEJVM?COfdEh_plfJiFzD%9xKL)vT&`$iiPtFI%bQSZ$J34UD4TRpfi?O3U}L)v|L>NojGweDo4DWgz)<0rHm7j9EHH41{?CHS+yho07upDf=VYt`51 zTB{BwTCVlhTR(aGLck+RKvL(9Iqs7SBju@IbKW~BYc!p4iPe&K3pfqS_o+A8^5;!k zg1(hj4^pTcI<7aL0Hf=QHFHL|PhLbdj5Tv2IQt*vh3A4-=~qyC$;H;qUuK}g`O~gN z=?cr=W6)qx$D+)BAdl*3-Qr&Ku%LWe92k^x1?qpX8EJO1U}0t|8^iMLF7{PqCRx4# znW0P35Gy$=qY-|XDKbNUJA^jz?GEjdP%yJFTTT=j$ky}iUdgVBa{qy8p*Mb%g`Q&} zztpF~8YUHNO1ExaOT{n0qF|=>cW4+|IXo2xhVDe6wgNQgpy#0j%(RtO2FEcKj#Xyt zrgB5Ia;=c1HOu{tv}lx4u78Y0!*rn%N1>Ns4$c4Z2WinxS=e?1@AHKlPIMG{{w`T4 zs0uxWLVMk*Ch2?K)6KaMN_Yp-p)FwA>+WN&WLN6qtKLpug|2}h+0Qs${Mi#Bp!Hkt zP=9hMv5nqQ{sTn2Pj(g zapIIfK>|FNb;;!*2fXj51}^<+?|66Wf{jT8V^$TrHeC?tv8vz&p<1Byt$xx!MSewBEY@PWMF#Igw`!LdG;2#QwQ%o z&z+jNV{)MHtHr)Oqznil9>bUX>$}HlSqCfi_S+0&LhH_(MnW%`@D`Fq{!0wY{nR=u(Yy^w*l)G0L1!Gd8r?nr zU}tAR?guyi(`s4gvRV%HX!`#n4sG`HfHm(|+oO$}WGPW1V0nuz6sRCEo*RR0nY ztVV?Be68=zUvfj&pInJbqr#RYBlJf02Md|UppKw<7BCv&Ppd!`jNipbLYdcn0Ojce zkq=GWizW!N`8CPI%F{=Z6i=QEe8|5eM3(l5Jp^s3|_6< z7M|rdi;=<*UbGGrMj`5FM{IZNXGkIv(HLd0Jx`9i-12>Bt)Qp`EO+3Z;gBU}gH@j= zT0V0KH8}82s6XrVnF&m?s(WFB-~1<95oP*;O3N)@dVewm&7Y7VC_(cVIx}dd0_CC# zeBZI8FKeX~5Clj<1|+wh&1qy40&>u-2bogBr-_DP{9AQ0lFJa&7#~`V=d`?WOCo_Y#|8!&CRL!sh8@IHSFfG~ux<%_RUQ&hf|(zlI+eRE7$rL8XxFFlfHWph4FVa~gCJ z2+E*yBIv~tv;?R!C=F(KSq^O^0O6v>5pps5*jXJ}}`pe_k0-jP;xJ zSoeXB@q7!(!^}p2sIfi}c{j7Woz35Gt5PJDPurg>SBA*vKJs+?jHd0N@ z+>aq;m4__f_9EX0)(d13gW3JmAsh%$_ZmV599u^bI`YMQsX{!`S)n>*+|LxC3FKqd(eCfs{`t{ z>k$>-LZ{@T)Nq@xb=ytNe)B%C(QOTIM#OC|qoBiW*~`Hd>$Wz~l-m{#jCR>~-w2n@ z28W0%uQM$9vL0lbam5ZO^94nWKVtUygQv1UKm>XSD>^vrJfwuTgXYKxIw67{3p8p7 zV#y(gVDVlt1pfg+1uGCR8B>Xhr-i_2%;SPWV;&>sG-ehEFy`G7(1YeRNGW5=Bk24HIss^uF^ff=82AgRP9S*? z1lh)HYmPN$NQ^NrAkEP54WMldrXLs}QCCuj}^MyvhiQ08`c$IA6L z`t!XvA*ktj5Vj5iA^9%yZERW&!fD{x^=A>kLgypU1A6tk|D@22o5&Y*E~QZ}vFkrmazAN$ z!D@jesKN`KQFW5#@1x_wY;v3cq<)iIIQwSI{Q8H~7eq zgsupc4nc@iX*V#cc^lor>f)*k%YgQ$89^x$hUk!mg`sStI1?%8wW1Y>*jOxS;0k_ zTLF9E4`@Br#zn^bjDluY=+|G;@oCVi5G|X7NxMJr0^iy~_ftFg=Qj;zx}O?Ubra0q zTHspC8xybPugOd*c!dJ3!2Ox+K=gXrx4obx=qmO-k+}?1RJ)u2j6liP47AYqZ05tj zLt~+U<@;LhaQL=|KHJ6gcHGA#g>C~rA@|1{pM+YspmubzAe6DncZoZQZD?oQQ|=_m zZ+=P;C?49$(YUBH0jVLIt>9|T_UM2sk&XIr3uR&8Wj#;spn1Qh2F+gq6}()11kLM_ z<-jy2k}WA!=Wq{rZha`}0;p=fK4- zt9(mf_(@jz>q6@Xpf!TY8^Rv*O@JwZ;*p(2)cRMbr?~tRRR1{WRJZ8B2gO4=id4_4(t&|IXm-Bt&ZmmsO^oDEdyO*3-$xLa?87_E8kriG^zx!$&hyjA2{D@Spl@9Sir z92obx6&SNMbSf?Imu990u7G`c&(o6fZPt+Yh4LW^^C!P0>@&9z1jh8Sa@V^TJ%Fs@ z@=vWm1KsyHvXuUad{%IkP!IaE2=rIkXpL3WA4#o0Q0e#0?$&$hS7m!HV%maXMXnF5 z@~@QksY2@^ThRpTL>Up60%JmNKIGk?7RrG3HhM3|$;gU#nFM%Z>2?y-EO znHiz`QNo|~;0>&f89wti7KZOLAR#fU4LRZQ2#i{e?`svNK4X~L6S|yLmZw)yg3o+X z)m(}hTf}*ci9eCn@`WpGfoKhR52UZbUF3T${Di$8EW$AHo5z85if?UrD>&u6VFj*t zp+RMxzPwLa@W2Oh8bXjv4nbtlybomAwx(Yw5|lU6E)LG1*#eqq8lSm{*+KJWO$nMe z0c9rv^#+B8lWo1QMp1QV9ZfrQbiVs@QSmC2P!&G|lyOkp!K;e@0-CD$ab^e2WttK+ z9|5W=rvHhC2jI?0_zeS4)hk*hdpJ0R;c&YTvohOdRwCops^aQ z08MH10<(kWbD9z~{{~cPv|OR#iM*Gah+1TO_|1mv^n_9cOio6Z0b@b#zY&)cPKljA z-p&u2Llsrt58~cECVfIw+MZ^`{!iYa`7zV1+W+a~f+E7uV6 z!Dey|Ggkq`XwZb+>+eD43iM;4OP3v~52#Dm`@GOUcZv?4gVMZkzxJ%7Ya#zVl%3N!jLN;w;-6N-Ye3$J=uYJaeI<9fww0oks3(8-wUVo0dDW#-1bTIcE<@q z+ufn5EUmXY?vr(h?GB%0SUE^Y%-#frLZ3ni?`hJJ7hWpD<051u=l9IPJ_M}c#Vvv{ z$anZodn;%T(V@O}gc7#n${S-HchSeuj{68{`k+5cYsWPUnex7_M3}RCWeRD5>@tWC zJ)sKi1Pt?8kD6CP<;=jk9mepY-&P~^#P!?7+G3ws3N&<-5Sq@i@Y@gq>rQJ3-??p?y4fpm@_(7nbqy0HUXuCUfC!p=mI>-z^E(T^ZP+|WbY*jhz7cr~~ z-CQWGN2Sj!11a<@DuD}{fWrl234-S30Llec1WgC3Trd=yqt%&d(XYSQ~#b)ej{1hr1%} zkKPe&{KgNXjh`X3Mf(HV#&dru+z8`Cw}Bj4FFW>U)E;zp0rZZ#bysT&>^k#ws6TlR z{LfVa>Hf5<&>=A=w*cV5QtN-<6Wlu>g~=G$Bw)f#oS?^|1SaS^SY(*l41gQB^7LU$ z^VSHrfr0Z~*zh@KK z;tNd~@?|mhd2nyDjxvjLv8TNu-+ljj3#>mB=5 zyGY*Gpm|$=`4em^PoGH%+ZJBrH|GdhP0vR9TWA+gMB-NWt}G<`_(Li66nte{y3pqh74B(oJLnzvT+v}hh* zI(g8%LT=E!#m>kg4>sy+1jss6BZMD=6b%j!Ld)p;&jJuJ`bHC|L6?iQ`7+M<*iV(Y z!=1Ti0>rsnZv{tq8=3@LHGP@y`xJJpW>ws|PTwrB^p6Pal_|R1RuJ@(6)NH-0&c&5 zkqfqFm&rygivDYdBbri#U;rFEUlMvdstd+(YH#5m8Xyv$UY9@ zLmxp13-w1{_*r_>m-W@vEJHW_USM?7V**QGN#JnPpr=8AU$*2>jjXdpO0DyQW(ktK zU~C>x<)*2M3OAhx6q(aiCfrmcDOV~AGJm)##!Uw)k#71sbjM6{D%!+$(?wOWZu0yq z+D*HWrp2eRw06^XkQsg({VkLMGWIr5e{080L$N7%6q+euVF%0HPu&82+eAF%^KxI0 zTjXZYrj@_;>W!$`$DyUNtc6<(<|1NwU zG?zmHFWiSNS8F)R3XFSzR51Jz3<#lUUb7u&=s)kUt8S7pLTo_QY(NRSF=toAHfGL0qZ@Mr(hL&@OY1O^2btj!P$KjlvhDuuPrj#8D7)k* zFxab8aUAPnJnYOtGj0WkKkqrXo)7Ns0>pE@hhQVNz1N|X5xBGy2ePgM8=ltnv#_y@ zM|Rm|Nb>~Rfh{=NHV*7~#8!+Q<-pAKseu9)p4SYp0zKZLfGMthe14OG%6d^Q4v;ZD zcupIi?;)BMa-Tkv3N6|$^E$fAwl3`(}-&@j66XB-qP45UdeHn4n6R~b`MDR9$^0acwY%gJk!irMiQu9t1ZG$n2&YRe%@y73goe zGl}(>>Q6Hc)1P`|cI{4|j)tH1yvczxcmu33Zy!ndGx6AXhU|Z;lf4q#$(9_(d+cv% zqY<33>&yOj2WknL=^&{7b}CT6*_Rw0O#_2wPnAUtju9GMX(oUMOU=EonET%;g4@01 z7||X-+aJ#s#Z%K?6b?%r^f(~;z~r^QQfl^XHr7fy%BMH+Uc=wpLno1_lKsA*iQSqShuVD9(r!`2nxB< z+<-_*&#eH&*l%VeLMt4NN)h}2LTa16K9T$X>p+ZJpK#59^$8Xt@5A5z>wh`Dl?nN)dn9!2fIRTlC$m7M^_D2EbWGD`X;UP0leFn{0U`b+@`>7vb z{1udT;GM|LwKyJjCD%z&$HS6YfzLdS>@;c|IfCYBq@~g6XD|!rLNz65W@*%yK7yIP z^y3A{+B#dXehNTU{Atev!x8waUlneAvka_Ux4eLQunGw>Ef#JIWPgrE37w=0O#=+~mPT)xMLl@_b{SY4Mu$ZP zjXof-^c@6FqhAC;8NDom{v(3k1~kg(OAKXnCe;ok&qI~A(e;#B(A>zxEf^623Jk0a}M8y}L zNxe$uT~46mi&+M->P+fT_yoOd_Z2W37h}xxfPOQNRsZan)PJH}?BUA4fhI-}&!kRH zK*do(XXyqpT5blH2s(qwfqss(^h^35n1$!=ni4dZYSfo?FEf4Vw*X*w_9$eFRtcK>kfG3{Bh+_6r8f5Xc^C;1h}xhG#^iSpWRJ<0E{z?NXRVJOlbev{ zm`q@4efRSSWa_i4_oBjRW9kHj#>^w;G$s=S7?UoIISVOe!D$h6Tm(HDXp}L(*T$@Y zS}_R79A^F9Rq~ zXP$^K<$`R!M)EN89{^Dih7TS#$h;Q_bc9)>Rxq!vg#&`-MabiIm=(-9V4>^P3g&D| z>#~9op&}}noX-k;=01!j5lni@V3LaB(#-TeloT|-1&5#n&0dkr-pur+?`Lse)^h|w za}5Ah@teN|hRuC!hIHB7zp6_xVO$73j>E`7a|X(>vR&&{JVCs8+KoPZpNbcuQ__&) ziBZ|E*m1ct&RS%M(nLAsR3u}a5*q@}S|c3t$#k%CEwWpNfNe~RX}J&jm?IMG6&V~n z2BrJb!z9>a2<})5g7Vzc5p;b7{SDBlQL~fn$hAlS%3#zShbjwJ3?wy>eDB2=HLzZ+ zMX=3fgw9-zQA13vMR?wkqoxzk9yOn`wAxM{D{S1Y*CLl^LdadJ&?Z1};tcKi8KR9e zdH`4)MjsLxG`bgb2iQk~)95!rP)0LmDfF!f`Ujvo@y0*yXXN`YcSB2cKl9BgV*G0hF#du1 zlkX!3Ru~O1(Vtcbfs6+Hju)aq9wHo~!S&z~(V$h3%_~VBW?l#o6%DxFxJ>4qPoSfL z`3=lh>y4Eiurp|G=h6*VWWDh^cz%Apu?FR0C&Y(96QP6a4St0udc852HTuk7flGvl zA87>kg^&iu)?dOb+-GS@&|IuhUzU%VzVrqF%#LR;h4VQ8GDI8;3?X7-oDi|NONeMg z9JmhQUA%j*))DviLL1+rH083`0P@|-!d>IRAp*!-GJxE|bkrJRx}b0nO(SNH9(XAk1Wdb+W(a4Uf|T$y z-W!adgCgi&^t!0AbQ+t80i^g9F_sQNkUf@iE{q*Z-~JdJRVsqUDRLJs4t|clBuG1o1(=}JQaArhlw3U`vL1-Yx>d`XA-S~wu#tu*diuPAy~E3|Rvbip@@MylGot2@5XkbG-PRsl z1dW*jvw~)&sunLmye5q^@D9fcpp4qjgOx98EJd2bqv(}p{Px5j&RBY%1@xO|Kmv-q zganHmk0MAi&ngEKi?<`e;@qg5kNG2k?}@aLfcNg%p#GOc6?moOPa=u;cDb07WaZRnN_H2J=FLbmMs@G_jxvRu~gw}sDPJagd!Ozn6EM5f`%m1wK5;POT7qAc#g&Qq8x1JwMflsG%XMR@X5w5=`vjltUg+tc)zIf|< zw92BhAl?_o8emC2Fz8U6X{pD~5DDFjGjlMn=Wf2~KfW(oK6bem@qG&L-JjWv8iCv` ztnqCaviYr$Z+-42_Z@@4fwvJ#xobbYIPkFOKzMQJB>2Sg&&-@3=zY9ZzR_}Ruo8;s&ATWGnv2Ta#kn5eoA^Qt_Ul!v`iVw@|MaaV> z5$HP{eIxlq+P)oLE6z;C?@T;na`)+*$=c z?eU9O+VG2~4b#nmh$_H~Ge<%rvoFIBo@VuPw-#U$hl~Z87P1#)R^e5(LKoWxNNG>#>;RC^Hhcq1J8F>)T}qT|v3EO+fz{%j8Sg))x*1dRgLPvN!DgX)6mGAodD2mbnl@Rs@Nal)PEt5BH^ zT_X%yl*u4w9}#F8h;1K6JIB+)&_{S26g~-J@STv3U28&rdJg+Dd_T0%pKMO%rz+m*+*ejk} zVcqvD)CX}_NpVe2gY*XWIQ)5=Q5kj!{{lk#J`l1`7Q&kz0--nFTJ|3 z`u6u?mNKE%-G1|Wq~!Z3^}s6NvYaOyD-qb{;y|hVIt)Gr^8iWeQ!hh+^Ybe{eEem4lG1}8cp1!}%sujx zFPoVdumR(XW9)y$_g;$d?U$vO<3lfz(Tw+z9p8JI&Q6CJVlgZg&&2V)msv%T@4ajW z?)ct|3wAiZ_p%bX(cgP1c_!ZXUXDR6K{GTFmEhuXDm8Q~IMHW!X(D>HI_QG|aNRbH z;~E4{#+&HECd0pZN$s4?aV|SLGvzPY)^lh7gg~7SAFm05HfiGOTdSO z6TpQj_axdf%ya`p-T(4?FDJ;nUIgm?SAOqh)ziX=MPR=P~iJ0k}^V5 z@ZEw4nhz8!CQPd7L&)dL`U33)0S5q*Xb!fT&PF4s@4c)5nQ|T9dwGlJ{QjHoz1*!t zDm(6EMX;l6xv=9ZCRtk!YJ(oLnHC#cV}F5x7TMyyPGs>LIy6B#tCxhkSXs^7fXhfsiXz6pX`I%qMZ;zFAPAQ z{NBqzVEo?8F%#sdl;;&=#&abzEyox3^mxYedkhuC-9psi&-)Yr!?^VlI3;L4kCYtf zj{y^Y;$U9E!G(Qs&@h-k22g{!Jl#i{{Du>ExIdFD|02+PT6%~3izi_T&Lty{*I{PsKW zOkesh1ju^(Ji+=A09Emue*;FTEszx_0KLGp9fLv)-;5!-J1Tb%vt zYE_Ric@C0YOkQ#YtM+B(GE)rXdG4rwUm`Nt@2?bC`XmBp_a6v?@I}ycN6^d&`px*5 zK7SkAftx40o)CS02&%N-8M!nMeIC6OV_to<>9>!g&z~+S?6!kYitq}S)_vYDWXe8& zH!}40^h?RCYjMztVS2LdqEQiqW0#y&%;b~@w zcK8&EJL1mI&S77O{5IEc6J`h||onJ@L`y*&I(5TpRGr8Fg<vMizwBXX>K%Rdr95JoP_in)@MZWz7r3JnpFR=ny3qZcOc*r5z z5qbq*R<%|gO7xyQYMWXY6#Mp~(-!(ZFX-%x?@!>{dN?`x1*>t~`Z!#|`qscMLdp~h zA${aJ;(;R`IO2gL9ysEGBOW;7f&X0|aN&cLMN>;A=8p7uE9cdpepMO%t){wbM^rb| zj{wmqt*XrM)_KcnGA7HMj5+uxlIj|?K809Kc4F<$i|aG~@j6+8c(Oh94RzO4A7vV1tHL!=^c{sVI3%i{G*mX0m(^6-g6u|?jmlru zsz!WwDghq+Q&@48u=8twCcm=QQz<=jcsby!)gBbEIVxowjOtI7H6Hf6${O23rx$4+ zhG1=^!$sSVULIE&mW%Y|?pfr(ah;fWHosAn4Mhjvk6Sh1OCP`kk_5xp0Js_H4S(qD z>QkE7*!8u+}T+V*bfJKI!Fh68*mxmmw-IS$tU0zPFUvyPTzzFOn{34A0~Vi@&Pxz z2D^`f92_v74w&;sXXi4&>3}Z*Rsn7XT!Gi&QW5WR@D^?fA?`EV0OxNz(icG>!%v7`KhiGdmok5hPEPo0RHx1zk&NCM&VK1?lUs+HwJ$>_rX`- zGg49~ru3Wa?saWavvKYTXAeCib0DY&1HerDEk|G2PQv(<)Y}sZQu^JRSdfy@(qmFe zM~{RTj!DS?sURhFd`gnls~_0f@K=mFvsk7uCABIA4!}iuPv9>fedz?&E96m`8mzMiVH49Tlw)9MRHDxeN9iP$`MGSYRdr)z5q|efvl5u@YsJ9c7k-9u#OF&c|rRg<4@Kxsk=Jhp`^sJ9&ziqP}k89 z{>{XL7F`GH9S+*a7+OAP<3Q6sC0{9MIWaW;i$!A`H2Qvlk~fiEJYGLwTPAwq`kpqw zg!1_>AST#0ck9!MKXwPbSKKQlGFR zZYL3=>`kHL2h<1hPJuXHtYZ&oAAq)x1X0KN>?aeW`=VXP%EQ$`dHI+ZZeQEkIr_iS zIni_-Mj-v#26+oG&s<3TaIubuL0bV@f4h#u>b(6hRS3WAgS^Ff0|DoZ#Kn4((1-2@ zO)b4eY#0Dq8)yY=mu~y#J+b|B5*w8LVJ6DAV{R1xHg&hJrloH=EZB9!^}Zr5U;9s{QsMS>QVN~Pa^S> z{bUhnA<)Jjwx8V;+fJgN(Qj)}{vg&y@#dgHaKr;gJaEJVM?7%E14le? z!~;h>@c)7bmaC1sSqgHCgI6;4QF-Y!-@bOH2gjuOMQdmNB$e(h$#TOMIF12|Yw=tO z_1@F62V@%WcZmzP>f$;qMG+L%8%k`38kx>qlICe#A>X?l2ZB z{kZ4Ci~A(J@Lo&QrHbK(Uy@ZZgU_#cac`6two$}|cc{d*0eb|zG(AbBwcZ^f*Dw^> zuK$NlX%Dgm>rX1ZSQYG4X{{GtBNFzPsC=9z6=|KXs`)<-R;c!ht#~Xq8F5WiaGHX1 z6ud^kTNS)V!ABK*PQmpGzN_G83Vy3#PagBdb)15yDmYfbi3(0raE^l4D0r)aQcC1^ zYrDkNSpz(5jd|3sSam*g<&)8f3}QMz&ZBa{sPiMw`(=4e7tgyTDT?2I1Re9}SeJ9= zS01i>^mAuh-$;2!VW3)m^wZ)9PgK5e$fWbm!vTv)7hEzVW7NnIqeqO4Fn5C`A{s5z z1JJ>fXe1bG4fC{fztv=)p{GN!cyfk=?lLf&{s{$fX?L8d=sVg48Jfw9|LR3NIn%qLYx()z$j?x8yaOh# zx!uUQu^YP9bGnkF*Fm>;!~bA6^k+b)A2usF+TMRC`tlgMwyRCi^?F9@Qy)X$06Fpe zw!IrV|5=lG{3haTJUZ8n@#tF5R5ebuU0TmQN>APHBfF8)-8d=&zYD$Q%bQdKIzT*q zE(Je>ZEW1**T7tH;!m^SKekJJyHn9)Xi4EnJczp1ljTCZ!jDgOLNC4IUC#t(vDB^+Q7=0(bBR_k+9o?0LY;>-ik-+|70cf{2DVmIQ zpvUunAKZwF@t7~gT&SdFDtg9jNf6K8fu5=8%O^^@IJ*J#d_k8q=Wlg4^rfH=bTTUZ zStWnuP)X47>XmNf>;#?dH7lmwli(=QGnD`HmHd-Hr@i$VlA!%RO3^#icob)pz;dpV zv-DBPuiLRw(d%Q{t5wlMG40i^2ZgKwnjIa{esoy57Bt z-mL7>aodggg7!8meHP1N2G1mr-tn*`hLI9B&L(bj5BoD7wDS)%1m+ z)1R9MN{WsHOB6jnrXBB9^p2U5Uz|Mwx=hjAuaI;de^v;(q&a`zRs5-QCACszex~WF zUail5MQ=V=GK%w9Kuv7?Q2*tlC4Cmg39l0bUE20ayi8>sA3 zVhj)H^rxZvt?th^2sxT70)MCY+f{$o{&`XG!>*x{SDcvy`i7#HD7#KkqCQk|YGTIk zK1ClKGwx2rLS5(?T{NzrxuG0J9FXVunuD@T+Y9uE>8Z`oX;h`$2-qcUQ<_F*$813b)LC3b+gNAJQa9htkF}}G~Xz% zn^#{`>8-55Hey^3&m7e4DQjpbyUtTt>%~7-KBu8L;@14cXa1_L*u=8Csyc6fCf~_M)88$l7a4||FrdS*9O*LaQDH5G~kZoLVvd3Cx{yw~NK zJEyL}3#WVd(>Mp+7w_}cc&=@z_EzTSPnvk))Y7R17kF|Ysit~1o)0uOiSB7cN|x8v z;@w0K-kGefof}g|*Vd)>T92)6WV>BIo}f$Y_OPv#FJ0n=wpvH!vcGs5a$+i$ea|i% zdAZQziOCXg^+mL)oL>)7#YN-K^Nb!bI#)Q;uGwQhJE#>B9f~x1>A4Dw_}W;<&7D_Q z3-_QI(JR3FUp33h$8^hB%=uP&>vFGH8PCZaTj>Ov?ZI)Hwd0}*y3LrH|D z#%Kk#|M_Ag>x+#Dwi9*6oTi$ZNZs?w_;RYlXCtz8XOT9~E32-J&ek*SF4EKNelH>o-xiJx1#K}5bVvtBU411RJ|Z_xYn4}(HQ;6H z`0eE>pI^pkQ&xlbW1}4{y0~Z~N1ZAdh~ZT|Cc5g_&cx0xp6x~21VgQ?wxVvH-QmQ8 zOYIEZQm|5!L`7&I?SZas?i3!|y5ezZq-(`0XxC>~Cwwn;wrRR*5zaYgy3P)fCUaD# z%XLv1!Syghp=)O3p>LP*Dguh5PIO9pYOH9i^Hjmx7<$^9G7!t)C8H}FIikJDarz+A zomB9RX>ohP;`AD+Mra%zZ5(at9QTf{gSTnFEx#%@ea6}R>UNso6O0e<{dpIW;ZtK zh!bm+;~avAJmpt=%B!x%1Xj&3q@wAV!*sQX?OGAL7e)rBa&wp7EM{&9Ka4{n^1F_0 zyD`};a=~fG3!TP0_pL-f;Wuhbmk59X&ZS)U zRZ78KhlLn2hB2b?x_RER*?`^#39B^5HCscyF+$vg&#i45!5f*1tm+CyDi~js#e4D* z%&RJEtTIMaTvrQu5_%hC&NaB*#LanB!h^hq%9=73P*^?gkw%C)VT8AGK5*V|Bd<<$ zj1iSp>JGcA0%dhdmXNpFvV_JVxNKf^IclwgEc}cR6-Z$SgT@Hl@z2AJyxN}lp9Va4 z$gMAItBFgWGty9Pi>h47YbStkzws#i>GJwKl7^X1X5e~1Q+?+-YW+=B-rlEC=?eh42dB&Hb5I(#tBjuDp`UblwVh>511-@}UZ0oJ zaIqE)9$spjQ6})q$0a=5sM++nDh+cKr`BKdX?QjAc~(-V_4z9en^k%0Pd&B%ifmjD zgg({f^*JpKI~1d?UzgYRFH+^lD?WXmOGABLl;ycq99w<~X!Hs76PG^UrJz0+tFvR< z--irm`3zM~L-CCcT^KM%yus&V&hmQyPs9D{d9)Ux%j@)gG3E7nGYzrr8g<3$-;P|? zul;YIORJZ~4zkf4S8V;ukmoF~&v9tTZ_+rg*!rJ~DX-7DX~^%cI4|9Pn$`gXTd<;U zk*3doEDJCdX>XZ|q_v#akb$jKySzR(xBR~>zZGS5-)0^9eBBCFK0`TH_b*+q_Txt= zLlvMQp@1Qh%Lq-DoSC`lDGvqkS8|{jqxIzzHnKUvK)Zlv{tYV}lD5%TF?xAaSoY3yO4jDPB{tfDR x3CHpI_|tal_`&X_b$1juj7Mc^4m*eA()DT*kJ-d4Kj-g~gZ~tS^NL9t{{zVt918#d literal 44064 zcmeHw33yaR*7ofr*U2DB2L&ezk&FfnvIt0!C}}UXBkW_>m?~gDr4tzg#Z*8qm3i+JICl}9D#J2 zpvn9q36LwFc`1o7UC~nkS#E{Ocfx*(PeZLCA;q zI^U{qw#+viup>j#HMEzbZ&q9*^H+ot;Ur~5!!IRWm!msaZj~yxN|n=RnfkH6g7|`Rk?aDOE5~c7%bkaE)brncwA-ht`nvPWETzAOs(Rzp zn(8^DMxI(zF{Gxtu6eM2 zU*l*St-QXjyv&bE8msH(+UQMw%Bbi{@l`;G#z)_^a#dn(S0nxt#>Yvc|bhM&*JAI2(R0 zzrt5O?+Ra4S#^y;57yL|8?}|SHXj1x;s#+~CEF?9GyL!)tQ z1Y~OJh(crHgmb4&pXQn38;yRSBK}JiF)syClQEG$lAflZDIimf?=D2iUI>zmokjS$ z$Nh;lZTs3C2G-DXjd|dc1fL%;EDK?YojYuMP|N@&iP3{!)$c6nSMsw=wf}NGKy?;)~z3} zjjnyfys0+2T0)5;vut!UwDruk(K9V7+RsLJ+vp88x>{O_+yyqeJzlri=ttY+++?G3 z&eu=7jowES0hik7$Jpr0Z1i8*=qqgWY#V*0jee|+ei%LMfx{j+?1BGRJ@7C0@dv%F zpQn3+=BDcm!`s&3Pl;{uwmzS}UYHas*adQIz;66@_b&p*^gQN8cgFAxc!z0DZP6_v zy^U#3Y0(am-pn+ovgk^YUdJ@2u;?<8Uc)q}u4ubRKfyGo=xB>bKg2YruxNuw-_JCs zuIOx$zK3Z}-O;HceLK^fvZF;JeKXUX(xSse`Ua*sbw_hW`YNV5wMDZ;x{+y4X;DL@ zFK3!lS+w&f5C+U)I*aL@B7G6loYJCOM0y(2oXVmdB0Y&|PGQlNNYmGme?X2mRM0WS zFrK4jkvl<%^nw$IcSP>~Hsm#{Rm6ndWf>g$Mk7(6m2cH-Sz2%h*nLf6xtMy@KHX3JOA-&j<$o zWx)}Q>nAz*_NkRIDv|c}@oeN&2pJ)1v^S#%9 zo{J6@>@yRkE4;z90r{ei#d+Bv5ASH(=3cy9P(CgV4H!BZ^*`N$G&@naI4_Iq1aDw( zX`mu6-5cna7nzBMSi{*JX7?p~U~gpok7yIXJrW9L7Up6?6gq*e7udUoT@U46gK3d3 ze~^WqVxgear@}i(D%hOkz5aPBe(|NAOSHcugV4&+X)rKy8w$01(3}I_Dm4t4X(_D? zPGTw?tIXI-_Dr0;XjGABYP;j2hT_JL`iyO((nyHXcl^|bpcG9Q9uKjU}_ z7TgE{ZQuA$31-|uY?FU@u+J{YY1`%>9W;xGb@w&H2dohRrJh__TGnm!!IALZ4;m zj^~FOVxGXd;Mh-*U*+DryX4z+-rj7+?5e^PONC@WCYXW;Eh}S~zJv&1GUX z6KcE$&2dPzbu{-5o1=j-5QoiEfKl=F@aIV)-@W8W5Ih&a6$@Joqc!3x4)*z*a7v$H zAb4h!28aCD)pXL>Y^JA|1_!)~^bGj#d~Yx(2)s1pe=FS^>Um6QFy}Ogg_rk}?E1(% zC_}_E*Y?Sw;`h@hg$AEE+4c71P|b!cPiS1sJ7~ynDBHEg8(I|3LV$ac$2Z17y_xokbW$++s7AtD}oy?K-L z)Z&}qxqs4q#lZoQs|@35;-#S>Ex^aO?Y?#x^nwYGBUuuhVR+q-ZSba=mm}$QZTGgW zM`s+f4&8nBfmqB__}(?Id0RKQysbZ`HUAafJp;9b&20SU<86!ivm!Tw6Mep~CSo}t z6Yr1ic2b8grMGsdIRgs|n65qcB*Lq0NPe>e^pjR^Kx4!<*gmm&#e z9$g4fo|8h|+hXp;;|1A#h2%lza{$A)!Qn2i+qO|y@B_=t0*c3%Ne+* zCuE7)VC|==-hlZ_YVhqFkwVrRFn29Ll2zRe8-nKFm<*c-0BE^4kn^3Se4r^|^Ivha z6DSu=;Cq{Vfg$%0gw4MIuu@2FJCoDMW(4H0ITd8ejh-wThVk!@c(0AFKulwNXnWVa z^P1|);1cFp~q27V;Tw#7(apC@R7VMIxUbhP-N@{ zb~3T2Tc8f7U?`meg0MPl?!D4N)8bSYP+@M^d;ttJ(3*NYl)-6x85AlFWkkngn5;UU z)KJFY`G&EEtFO|)m(QIF7`gmwOohZi;f(ef(u7B`G?xIJIL9GB`VxM@AP*{(2Kgb` zX3#8=L4z(JW;bX!2+E-2-V$5j-A$ z)Oa5I`)b8C3UVH~r8DDaAZE`Q3&s3k$kleP7BY&u>BW&gZ6OIdx z^#H|Kf0o0BV?6*m#`7GK2bt#sB#iYFkawBPJC{I50rL&0;9mSHOo+VlFIX8iKLUjh zJN-WZ5gVx{X6^@(@|H)ufgL4*_q?mgBnGqlu^+vGj4-%LzJ1-h3O$`qe~vfM#2_yW zyKg;#JsqQDl`z3wd>~TF8UtoCxWr(`)+IaqZAc5AK+f&V!uml|!saa+4Gak}GmtYM zpfoVxbf&^)F#xL!nwh|`vYC`5Qr*Ie9FYp`hjD-RNc8P%)wqvbd=L#j2PrwfM6O0? z51WtsrG9(JE~ncXKvQnJ?u0~_ z9lKAsY#ca5TzQ#cDKO*)rWsePfHF%^#P}m-jX!uQ0|Z2%^&Zf}=7~rNZ->pSI65GX z?rU-k!P{Vh7%Qbv21D?D2=c7zPiiP*)o(Eb`3A?BG8`~+%vTtK#8gbV6KRg7X92At zcsuf=8}SRq{2di08WR*08WSL9H>L;#81qdL=wb5$q?9q|#L-!CbO2Dg)~fg>Q6~oe zbyO#mk!3e#PlMB#lFme9HX%)8z5=vl282u*Px7U;YmO!@RIs!N9=SsmdIm6hKBl0^ zI$%}pV!b}L3MyhDF)mop4~>A8dr2^3T_dnf{;9z}xx_?VI0l$p5u^j7)j_idb36Q_ zp$O!FdgiZ4QUArmsKDS?+9#hKekT%Z`MI17W(ElfT`vIRGht0<*fCl&ZBQWWXyR6)R z&NZ)Qu2@e0*2*ncx!zE5-WV(d&plhb=uX`Ld8}+6UjHuu%^O&p*FsvzoTidHMY1?= zH1oJ^Lm5ThK$zvSO9M4|Ud)eUz=gh8oVS3j;|(v)+YXq9AD|6X8xQk~ay>1s$dzBh zRPzn=wBn#2qRVsM2JQ(wU+jKtCx3!w2Q%G|4af&{ZTGmI=fc(R4c6qPd!7?56}mTX z84$hR3heN-hFztBNAgyHifUI9pjx2hYXMpucruR%6Cx);0dL?-`34-=5qa|qrgvbS zn;y9t`1ryfu6Yz{jYI8NWqwLfzDLYD@1mXYE|&(EVuCRVLZzp~yoGDsi(`zj1Gjj? z>o^-=Y`#P7X*!ml80BwKPS|`#QxW%p3SQ1mVe=j&Wfa`tZo3}RFr$S#1Vh^w*Y<_) zvs_-d3@&#_q8T&}a--VP+2kyLix095Er zGYUK1ZP!4I*1T}@btjg%UiS{%TH<=%8${oR_AirtQfSPl-q6VHkt1k9aC%-==u+4> z_N_XUE`QfM@ExH%zX#&o_bXC^=01YZ$X?#U7u<{QLsn_|$KFsQ-S;B0l>Xj<8@=JR zLOtluBG6xDqczr2e4o`a(4n4%eA^lP5S7d-_f z;(K1+wVzn4iO$GhDJ@|Br52?F;phaIqsN2A7bwlkM}X<{2BzobMi!w&aL61M2r?sJ zJ|uADK|<<~sYpe~iN!z*wZVd@ixsm{omtg<{OUsP$nEI>Fq9T8o`Q#8<8(po=r zMZF+;2fhtbCP*cLm!gkY>xB{wlb|^Y^vQwequar$_HRO0xzM1p&cIlc1;2eyylF5< zW_$*bVRKE5YPx{=Fi?i0Rn*qS!5KCm08KPaz+BGku=z(#37h{3l$`|BKPWVsVd;g5 zP1QYQAnnl6`HrZl_(PPCoyyz>)LHRn&{V}QGCOQ;)ReHf0jR2YokFABJ`;Y!095sg zRw*b2hcFy&4}2@yPRtTLsRu?%u53QMA6;Q-v=gXPqwSz6joxN<*nCq{!scs0l}1|> z8lAwk{z<4swnxxB1vzp;$pt1SqoWiz_h*R9{B$e-;N?;<^NGv+PGGY9E-OE5ZdO$J zJ`pz8Inrwq($+LA_7BG1hs=a&SnVIgIC-|Hl3}wE7!3^;Q~<-W0J#@gxurTcPOkKe~!^kVv9Vc15*@_oYYAaBTgyRYLKJp+6F+MRp$ovRk4f~|^;<4V{b`tw% zaO~-_ zEr3IVX6$|*1Q~1e#k+(qU3N6yqb^VtE@u3G4m8vkyif`t5vealo7oH1g?3)Fg*x;kSVVVRH(A^4pnlGz+Nea5+FNzg@t( zgx^jW`E!1QXQbawrRAZFz9?b&?cH*x-@ew1GW5l1nUz1+&1Fp{E5Yw^yd!L3xB?d z5|%%MbDaKs@WVuZ4iwt%jhx5QdH`dAYD6Cv1G5FFu)h;q9X9(#44YS>l-8p%VB+i| z-d|A(T+j?0E*M1+HZKNHF7U?D9H7btrvSBFa6PLNE_ikL&$-}!s3?Zwd#D#K$g{iP zrZT4s{=75M1t$q@_eQ3$w06M*kQqH)^$U)tN4`t2Ut5gz6Z))b|K&gH@%41joB}o5I<7qhFLiT)hP-$KF=qh70>53X$+#aLNV`Ey)yS#(wpxBQ%Lat}HzcGB<@Ne*iJ~Z|tWJT@)f7sjtRCrHp zIK<177o4k+4cT>+9nZvf8*BZ;g-);_1ipmz&uEbGg&RAp2O%D71&$Z2QsxnGm2xB) zVBjN64>IoqaIO)?1pAx?GQYB1CS)wmqg315N)$<~5vYRt>Vd`LVt5$lUHt(X8vAOH zzz6Gum1>MpR%pz9q=F$?0WqK^$g3D2ryNeyM3eHpyc6F3(M0{|oA-$x%KrZ)68 zAkDe%3qY%nA?!n@jM;g@oV{XO5WQvRPzY&Pg;oJZ@A`=`+X&{oGFy0QAR1){P;|qc z1o6(-V&GBzuWF@QXJFbk#k{1^MXttAW~z-r z@wxSU4EE7z#_iw;jy(ol6CZpt_QBTDU_mMXwkAFmd{Qnb zLz-WhzuksU*Z0CD_ySzY&5}!A$O?H}*l+9S4W+$70h3+(x&M}n%B}z-A`=d(6zAnv zwPQ0P8_}$o`z*Pop~2Zi@~4ExsB<{b2pcf?JbmGKqXgSsC9UD~D)&sYH(VCE?j0Pt z@VI+TJWx{je)H@2>dT5!u?KQyODP+tN)F;HDkTZwNPi9(U6tPB+;iczZD8f{SJUIT1ce z-hY;YBRKX=s0qV!o{(_B{X%Vf2O14N^sEy08`Xz$K(hMK0rY&)!Bh`|8~WvS(~X`7 zFY+tU-*jgZ`x2@@U3iH8l$O`EJ8d`xej2u&11EnYtT127r~FHBvf)D6|5PV?&e_S9 z9Kw6-Z+&iuxp*sQf2&6=Ve=XgRDZh>XwdYNqoetRuvw+DsKF$m!DZ$g&|s-K4Uqd2 zb%NWyq=#sapY4x-5XG}{`bs$9NI>+<-O|_g4q1pCTZcSBB+(%;y04Yq?w5zbE&61o z<$iQmt8>EML-bGS_k`{_Wgy(X`pvl8XPoSG`-B6fh2J-KI6LiKPPeP~33@I|2)WG6 zmr9!1fEfEN7{FXgMB!M({vM>ZTk8|KKQjQtg!Kv63|OCFA@UCV9V~br9TxjBY|TGH z4!(1MpjvYk(4hHea$vt=4*4^lVG27m#YkawGXzPjZXP0_(BBKPIfdy#=Gg!V--M4x zHIroCnFRWq@Ot;+XJJg_**7==nXg?$(@+r>$$tivlcD(Vi;o02)e$y3z>>}^_ha9~ z_jM*ymTkaK!BjZZOl)wk6v#R+O-Jh@CmixJhXUuAIUX>2C0MQbF!J;dTz}ksp=92(SoMp@`o3W8Cs~W7D=+5yOzhJp zxUv~#ZLVxY(&1O!J=avs`Thj!p+mM=r~HS zt53qIuVK-F_ddOXY(<)jj)ws)nV7R6Q?4L)NrpAONedOcoQM87MHQM480{&I9s$Lq z(M!N$GrB}%(CC>0%PAtT8{HQKWpsKR&55I*&2$(&LnsKNi>Y=fV=}6=jBYvKX>{hc zM5C7=O{1Rzw2Zz0GNVVJZ$usg**00SKPJSJVw}f(2VwJEJoIeZ=u)C58*lybn<1;`4MU9mzx{Pf^%mY4Gj5?nSq>L060nU08?Re zC4dY%mjgr4IVnldd9q8;X^(u5V^v{uI0~`uJ)3mA_!}1C(X*_0fw=S`Vr6szDC`GF zf7J{9kytE*f&3glkrzpm$8-X)U%?FJB8m#Sf?=1AKp!vah(!G-yh067j9)^9} zk)DFr(FjiWvQH@YIs?%)uL{@R2o4d5{)9Rah+bftfyg?Wa<8CpFy2GV8inxac_838 z#}-%v$IeAcv`W~V9!IO==qR8GfoM4_5P@hVl);#M8-lDc`O>-0F?rQ1iDU90(j1c~ zu(W>r`4eR7vnk(Aa~kuipwO6?iP?==3<8WvmB!qIlrrW&;^@6`v>IrFF@Mm;d<3;% zOi#Npu_;bt9(_5{m?6rT87!@hxmU=PQS}zdpwFiKkb|gt07B@2T;xRu;T<$`5->SG zCd8j$!3;P*N%)D!kNq$RM>n{h+1j7obB8uIU`Buj`y3<(_Es)}V-R7EL=q8(c_PA; z3$i(w=|N@=KthD!3MOCX9ZR4ij9$U~@Dgmnc`bByK8IMr{3lpMua+yAKOpIxcHOs% zP?2zeEyxP6&rWs*(IaFqX+T=EbRg%q%)+@OO$nP^MZ;eh!<-!(1-6)@gj2C(MVl_lx$b* zOWl@aEfRnnn^SH?(&-du2)OD+;gEC43Wv1ql_B6vro~$1+kOet@-&gb!82Q6IVA+v z7=kdIcgHiJhs`RWnfz(jOgTKM3=>wd*S*HU= zuH1}KLrm>veTy^)K!29jqvm6lR%?+9gpGUkTEqkKk#(xjXMp0w0OH;RC?<_gNi_OW zkwK%c5LnK10=vh#zPK_JxoJ zh7OYRCA09Jr72-EqS3&RC^G{&+W;^-{*I}zxeS0+2F*rb2oV#KgovBEgot*&-ypn4 zm#H{_^n<@d8{eSg!XO;l37f}4D0SQOpk9xh$R&2d`PRij^PUpvviXNxLv*4|oi5t~ znsS*lfE0#>yY2^v2q0T!0NKTK!Wv?gpl}d9Ma&vK@KP-Zn09Z;7S3u#N_aYKHpJ1k zI64hz!dN<)O~jX(8YqLY^gaYx0mM7eIhIDfkT{l}LYe{O13+smCXR4NUw7%Wj-5uk!=kmTB_(%@PR3l!JQKn0g6x$aSEa2<_CQ(XHs*N*2!jXDjk z46nuY0SX9Nt0mW$Dh;mxh|5}~xqc_P=Eb>w9p}1TbJa?&lsMP-ajw@j*F?$nUYu(N z2Ab0TETBTKUrVl0Dh<8daqV|%S>LZ0Hg?8kJs9V@9jK7?hUB_gr6KFYI9E^2^_b*J z(_G?Poplve?2lSIxV&|`6tc0u40ki6_f2JL?s!;&_bDX<8A~#efxz7hq~O@S_&t6& zIs#6vh;AOtSPI>-rMML(d1a1!h@-(D1v&6c-c6`va;Q0r>Q4{&?%<*LP~O_mI2X5#XOV4oXd?Q2F-`(w6P`!|7GSK26T~Obt0aVN zTMq4EJcH{c3+{kZgH%-QO|8kx+E?t(>a+VOS#K^g8&E1uD7t{jd6=J4^ciEM2V=Uz zsp4Bu4YMI{QC_a+8G3=Tgox2q9(O)=MAHs%llgPjJC)2+tH!(YpJ(R5)O84r|^xiqj3AR`&)Op+>5`5P%^y7 zUvCr#-W~Yb);-0A?_4*zB$SqpQ@LK(I=JNxv`FuZN!1Ah1-9)Qh8xf zCw{$&AKbCLkelYOB4Z~nV!kPlKhwMyaUT=MZNGB2P2qgh`hu&ZaG$$vAQC05&$~*` zT%OkgZQN}yaz<%|hv8yBT1zs` ziBH7`g6P}WKx8-H^AOkb{9axh%$RWwuPScxdy9kSgX4&I_|FZFosX756ulr=@Dj0a z{FeqZ9zstCXGO5j0S_>2e+JDPN$c=8VpBU^euwA#z+15x#9buCH9rp0@gN1q<_f9J z4}p+@-HBj9y%65~00=+e^&*%tLDbQFF9=tl0|fhoL{-hd2jL8a-k{lA)Y!ZvvgdoQ zCtOH^5#~7V6X5$7gd-rt0Qhi#>HEPTS@(mVLk->B4+a}!Gm+x^}iG%?|RZ_6{u?)P4VT5uN|zxi+y&WijQoanP> z773c@)#~5}20-=u7{xURI$ic_gBcl6FWLK&xZKi-9$W#5VRIKSwr8-<^Jl|GOrGj~ z?_`w6`_N|OLlf(MZym{lOdmkP`w;KI> z&ELh*cA%+}Zz1^tL(U)wn-c+8DJ0wP_ojpF0@-%IcP)2V|DF52yNeR+_z+3jaoZYU z#~LPC+h|TYfgw*YEjBtLJ-M*3wmdh948E>}iCHUrD1AQ&)DGyaaddec{UcCeF4o6D zX`uC90n1P!-c>rF5SI6AQLVM#Ip$0(@396En;pXeBR8+c@}8L5?C2rP*c%za()zt0 zoyEG}`yslZH0W1}1|1L?H0T>*c7wKppbUC8j=mO0w*XZJ{T^t7L8DQ*FlaASg+ZN1 z+6=nl45vY>S0x%$h*ErCDPw7EP@#}1-&dAOhBaBFg$jHSA4yY%1_R3by(56}e(yzR z$x$ii&BXnWjK4F}Yx|2T-V_YM4la`30T%4j4-AO81ppYvZ8PAMu=xumGHiY~R{4p8 zc@;`vzFaF+Syv z=Yy=H%QxGrsZ_#t#@+uA{eBMFIesrBN7%d;Y0>0?oa>ksHY+qGY%bJjAm=J(26D;; z7?MYT?;ikEC1|DsqtsT&N)mupv$7=LZ5a*Qp5wZ%Jm--yH2d1WCE3qhjVznVPa)aG zcEfyGsQse|I37dyi8=U-~Ux$Id>5}-2(lOzz@u61VQ*BZ2ID;KaP5UCiMB6 zS&)PH`=>;oe;2vd8PA#6n@l{Uzv1ub^QWSW8pPvKihaJ8rFEYV3YoIcFGq&nuH;{c zBy*N3G#W5y-i%1WF!mpCVtmGpJdD4_5H_okNrR%*V(nUFL{9t;1En7p!g9SH_x019 z!r+4`C_)pT<8cYZ`$0 zi8NjNF`(5BH_&lSFevf?$hPGS&$I9^m5iN_2A21ID-fvUec#QvME?ufCUNuevXhQ?)%OJhraK7^?tVRLM5)t`ka4H z$}_VxH1_FXaKr^Afj2!fN&;Vdrh5WE%UYbUs>C5>g8qlgU|I>4y|r%7z&~>eO_g*zusR~lRHV~7u7`hYxr~3c8C6%;R92jk3m#v2 zeSKqvk=uknwtPlzejP}IPa7CdSJqUP4<>b>R3xdK${I?lz#rIF+IC_o$7ko+c9l?3 zWn+`2l1%!l>l*l{peu}|qRYVLn_pS(kE=w|z~tO>x^k1$4R^W4ZN#OMn0rcYof8Td z*sYO&9^0}esrLNa*_LSWhqUSEfeEeISlLuwR#RySvYJA=p1qyvY-7LqAJe9n0z7)IEdjxURmbXD-KJ zISq=ecqA6v4*bBQvDiqg5A**Ti*+EK|5z->zbBKvG8UUk`s1;G}0L%qE@J=i?8*ssnSZo>LN00}&>Qm$+ z?kvTv<|4odKKZo+^6w_B1U&F%EcOy$_FnKa4fq4#3jBqV(O3%4!;--dxB&1bz>NoD zv5ll-Nf-g#4mhlb>><5XM{!*>&2TNqa`o-`i}ZH14e9;xGaGHkb*xeR3%-sSdH5NL zpN$W}7vM88vnFI_PjVl5WqONo_Hk#Pa$4RApc+_q7?dDM)O`Dk6k(RQ$M`kWap3JOqnd#mmv%%JmpJyJ9#fGpc`-8| zrjER=b)X1e&4OU)W5hK0w+*L zE`Eww9^>y0;Ph|GBsh57ktydW2bY3e0iJoQW3hu61FYlLl!=+yOH#*W=7!S7W#->v zdNPOIc!Vc&bn7pmeR^tHT93@ppurd5u^h<1W-0#8K9A{Tcf{L-wrvD0&p~?+v|$cf zCun0p(>`F?H289~gT^>B(nh1N27^X<6WAq_wHtk(>PzY?)>DJ>ML26P0XQFyu@=x4 zfrd91!`c}OvKPN9bR2^|v+|`_>_mv;!#bvd_7Z6ONf34XhILFx zY-_8IHHWH$@|Ho~BE;Vj|CY|Frt1&_>DNxkt3ez;m-^vh9chT&D?mHes^gG4?>Iyi z!Y|VyZvjpx9xG~)eU^U#dNF8fDJj}>5oj%-d7x7_?TP+$r?WpzWP`FlJcIJ>m~*Z= zWcjC^<;UsrKcM^;d{O$b`~HObZ*rC&ugj0Ycy7lg?V)4G!_M-e{(jJZ5zAxl+e*J9 zcBH=K(qUP&%e}~-${WE-zvTT@#f>el)cq+*+*fzDyqDi;>TFqy<=mR0209 zFs_tv%#Y`oAJ6y>TW>tVE8}}nU9#aL#?38|vkY?3Ee}73J#g3qhdprE1BX3u*aL?> zaM%NfJ#g3qhdprE1BX3u*aQC`dO-iXr~2PH<(3d1{@ba1beex#*`D6#NOM2ip3i^V zl}}FzK3y_3}AR0z?11tlSsjqw~N1DONAi+^fQ)_b`?! z1Y1wy$yfR6^NTS}>Bl`PKHOL0!~2i+M-`JFNRXk589SBz+)L(ztse2(MtT-Jb?8{XJaY{&}wO;u2B+{R$itvdkAzxMVe;w$5 zSKe9iC~gMhnV{eW3RWpNU%?v{yi>tHEBKUxFDUq?f}be(je{L3$~aT&yPe*F8r%-3}BW!zRo)8n7C6F)U4 z%NP5AN=BFb#P8IWK5_ZB%7NjF1Klh{} zy6}6m^WncF&F|07Clx1SBT4CLM!FslwnC|h2727v$dGb`q2og${TGHFM~U=wLyx~i z`Y#PVt`q4!3_bo6>4-L}eGoUgoDxqp@Rx+txU%P_8ox3+lF|{K zZ?|)5IePYp=P@j&8DGWX6p3*j;=ehr`%kW7fF?=#bGo7P`?$++8%48$;L-CD&Iu-! z!+rV`%z9fK?KnZv;p&9vLPgi}vCd>rcX5_beo)hoakZd#ZI}C%96i5lIe!)WUEB9% z(3AOXcQ^DOyP@Ymv1D=vgYH86cB=N(?%BKVK$631>+bZ6XKF64BLZ$2#CQjD(g{8}tHbJ4!eIQA&y^uqWpx>Fj& zPl9~5fPNI>!_Jv9t=FUPD*oBOm-JZ@7~czim9T*UnExnes#;g;b#Q^8+ZYvkA?Vbn zE`rH8hDQB5#z*80fbxMAg3V7n%GxC0?Gkas-PF5H*ZVsuxJb$WLeX;- zUF(wpM2{WCDnO2)XBb82NvimE40N!fXXQ(}K0kD3H*zimJy)_O{H|4UW)GKSt-~TM zf3l>DZ{$F40iFKr=V;%@6@RW8&|1$^>7SIY^LU5iZ&xo8;yX6vzbfdmPW$gK75`LD zWO%f_KY&hs(r=e^B*ilX9YfoqNCJ4eq0bg{^au5#bDb(QU(vUWmIQHz1L%#Qv)*Nn zdLIHkSCzGc7X*KN-JWQ=dFXzdA$Fvy`5? z-_8L&w`;ZH%tTg_aUBr#b}jIzpl28@%0K5SmJUVlRO4No4M5&&ir#>$OMJxH62MQp zk#FLCFIhW|1pNddA@SE#!H@W$#)XbE7k4A41$6o?*U^s46}?C~xJXsAM$6BY1nvKq z6}?4`BXNcbEbl5gt5kf`?e0QDu-;*g_WHG=uW+>27|>+#bi}I{6n&|~KOd8x@+;%OZL(hd`zuA)?~|JT6X?{x z^JGcU<2M8I3+c1-C0(~ywxYLBll0e=ppzB7;e1KZapQDBS4=iAL-Frai$+Pl z{Sp|zRdl`Z*X?q(qIXV^jCy{&UD3A;lXR7B+zWaZ^x5a=Co4dwo*l|w-4CA^a+DG_ zuv7Cp#!DI$V|%rzb`)nPQSMiYo~7(MQOP<9bk;lAA!n?jdmZgtq3F2|f3_%khhyDx zAL!Jl;TTz=IKvC{Q7y;Oj$0Ldn4{kB6un64S)}yujfHx$cIQ85NTJaD2jbs#O=YU?{O})z$f6LL5U)%eHl0VYn=S_;f($UYqRB{>| zao{_}U+nPD(HY78kT2-)|1_zO9@nRX&Uww*56=Rfepuy@U!~}Ij`q4<$d|J1zn3Zg zS&lfrNzwBieuxS=h%;ycK02QNpyfNFtyRbyFgrLUs7wss+k*iyKuRP9e-wd-43U(t+PJiP0r$X)MeCC#uE#Ad4f&KB;D z`R2jap~DRe9cI{=Y(;Tvsq5V;ai_?}%BE`n)^G%rMnLOS%eoC>=hlavNtZ%$m5xvFJ z4By1kbH{l~edkV?aQ^t|zUiKErQ^w%(25m_geHGkqu=lqUpU1xx#S$=S=IYcwQusd z#lAYRLtNeDo6}rf<2UBiR45X-_2#nX)#*xczsfhas=m<=r~CM;sS4c}_uguJS2kAr zD~pOIPB?ek^l6?czR{3WQ#}X!>P^j}dm8bQ<@I&A!{@`@!s@!YjxxHoF16SBEOq1C z*lZHs*rQ9V_OPs!mlJWZtj<=s>@U8?(T<8`-?PfbFX8!ojx2GjF0M`Gf(D2xEg5%? zZ^Y0Mg~FLu%|2^SzD`KADFUr|V-O?0&grsvUy6lRw zhK4fqkoj>PDMJ@R8PVSQn+HT#Z3>+biJB7i(eCC2O4jEL0G7{mMpbi7O}y^fGTzy= zxoBvC?iMFjH>63l&P&<9{ISjPZ%jVQn z8orWq;pB?yI$v`WZa7voR?=?WkU~pcU66<=HdE_pErEU|f{dRn*s7 zy+eGg(`M?Hf|Z&iDnbKk4|HvFyKragicf^`Zs%0cs?VxU_+IL4(R9@!)Utzjjok4j zvsI?cbx|3C&l$bq-!i+5A2GmfV$tucnW>_w-Zu{p#kkQ9lW|Z66B*ptj&TJ=j@`_7 zk5RGD(WusZ!PzZdjnFu8^dBnj*(WtypQmmzX-U?NBKo?s#pn$7hOV8+)(z}|Q*`-c z(Y$LTS`}NW=oVyLqo(mRE;XxwEyAyMzMVlublrjHG&Sk);55qiA>2ov@+*Ah^RB?u zP+b!rks}Y$Y2vnPMXdf5A4$s1UAmr_0+aS?%OP?3T?eAom~0lgaAJxfD&A)2)-@018$rd8>Iy~jj4K(!m)oJtn^)E}&lp;UbL35`Q=S#5PWYORMX{2D4MkirlKjiGpTuEiU%ZExUz2jw|87k;okB_4eaMniqB zKxGSFjT3U;4qNC}d3|0+!#pzC9=$)RNuz++%jS(o+Djm1pU%&*^B`pakRdhp=6jSKC6y1)$>-npIw(=h0A~?*of?=xh6P z84^58r}fw8el#4Zn6>_zPs1ya&$D?tt(=(YHv`=N496!&JPz zWh$Q5a$Z6PzrAQP_4(SB|7Q8^D69K6>(J-FR;lv2%CWkC>3X#vKSUY2NY}5=iEY1* z((pHFB(Q$#@)~}E9D8}|JXwn@vdStc5nWEhy&&7m>+>2NEs|e<-;ycrktaT1BS~3W ze?#pv@4H^|W$S^fisdS(!9gIbf0!mHsLQKuZY!6(eC*KbyE04sPQO7Wa@BVSj^j!A v(RS