From c1fdfd1d8ffcb32979fad2842ed120d64befcc56 Mon Sep 17 00:00:00 2001 From: aichhorn Date: Tue, 7 Apr 2020 20:50:37 +0200 Subject: [PATCH] Added Tutorial fro basis rotations: Sr2MgOsO6 w/o SOC --- .../images_scripts/Sr2MgOsO6_noSOC.py | 118 ++++++++++++++++-- .../images_scripts/Sr2MgOsO6_noSOC_Sigmas.png | Bin 0 -> 78192 bytes doc/tutorials/sr2mgoso6_nosoc.rst | 113 ++++++++++++++++- 3 files changed, 218 insertions(+), 13 deletions(-) create mode 100644 doc/tutorials/images_scripts/Sr2MgOsO6_noSOC_Sigmas.png diff --git a/doc/tutorials/images_scripts/Sr2MgOsO6_noSOC.py b/doc/tutorials/images_scripts/Sr2MgOsO6_noSOC.py index 890c78c1..232c087d 100644 --- a/doc/tutorials/images_scripts/Sr2MgOsO6_noSOC.py +++ b/doc/tutorials/images_scripts/Sr2MgOsO6_noSOC.py @@ -1,19 +1,121 @@ -# Conversion: +# Import the modules: +from triqs_dft_tools.sumk_dft import * +from pytriqs.gf import * +from pytriqs.archive import HDFArchive +from pytriqs.operators.util import * +from pytriqs.operators.util.U_matrix import * +from triqs_cthyb import * +import pytriqs.utility.mpi as mpi + +# Convert the input from triqs_dft_tools.converters.wien2k_converter import * Converter = Wien2kConverter(filename = "Sr2MgOsO6_noSOC") Converter.convert_dft_input() -import numpy -numpy.set_printoptions(precision=3, suppress=True) +# Init the SumK class: +filename = 'Sr2MgOsO6_noSOC.h5' +SK = SumkDFT(hdf_file=filename,use_dft_blocks=True) -# Set up SK class: -from triqs_dft_tools.sumk_dft import * -SK = SumkDFT(hdf_file='Sr2MgOsO6.h5',use_dft_blocks=True) +# Find diagonal local basis set: +SK.calculate_diagonalization_matrix(prop_to_be_diagonal='eal',calc_in_solver_blocks=True) -eal = SK.eff_atomic_levels() +########################### +# Now we pick the orbitals: +# BE CAREFUL: THIS NEEDS TO BE DONE PROPERLY +# AND IS DIFFERENT FORM CASE TO CASE! +indices_to_pick_sumk = [1,3,4] +SK.block_structure.pick_gf_struct_solver([{'up_1': [0],'up_2': [0],'up_3': [0],'down_1': [0],'down_2': [0],'down_3': [0]}]) +########################### -mat = SK.calculate_diagonalization_matrix(prop_to_be_diagonal='eal') +# Now we set up the U matrix, first in cubic (Wien2k) convention: +U = 2.0 +J = 0.2 +U_mat = U_matrix(l=2,U_int=U,J_hund=J,basis='other', T=SK.T[0].conjugate()) +# Now we transform the U Matrix: +U_trans = transform_U_matrix(U_mat, SK.block_structure.transformation[0]['up'].conjugate()) +# pick out the relevant orbitals: +U_red = subarray(U_trans,len(U_trans.shape)*[indices_to_pick_sumk]) +# Finally, set up the Hamiltonian: +h_int = h_int_slater(['up','down'], indices_to_pick_sumk, U_red, map_operator_structure=SK.block_structure.sumk_to_solver[0], off_diag=False, complex=False) +# Solver Init: +beta = 40.0 +S = Solver(beta=beta, gf_struct=SK.block_structure.gf_struct_solver_list[0]) + +# Solver parameters: +p = {} +# solver +p["random_seed"] = 123 * mpi.rank + 567 +p["length_cycle"] = 200 +p["n_warmup_cycles"] = 100000 +p["n_cycles"] = 1000000 +# tail fit +p["perform_tail_fit"] = True +p["fit_max_moment"] = 4 +p["fit_min_n"] = 30 +p["fit_max_n"] = 70 + +# double counting correction: +dc_type = 0 # FLL +# DMFT loops: +n_loops = 1 + +#for first iteration, add the outout group: +if mpi.is_master_node(): + with HDFArchive(filename) as ar: + if (not ar.is_group('dmft_output')): + ar.create_group('dmft_output') + +for iteration_number in range(1,n_loops+1): + + mpi.report("Iteration = %s"%iteration_number) + + SK.symm_deg_gf(S.Sigma_iw) # symmetrizing Sigma + SK.set_Sigma([ S.Sigma_iw ]) # put Sigma into the SumK class + chemical_potential = SK.calc_mu( precision = 0.01 ) # find the chemical potential for given density + S.G_iw << SK.extract_G_loc()[0] + + if (iteration_number==1): + # Put Hartree energy on Re Sigma + dm = S.G_iw.density() + SK.calc_dc(dm, U_interact = U, J_hund = J, orb = 0, use_dc_formula = dc_type) + S.Sigma_iw << SK.block_structure.convert_matrix(SK.dc_imp[0],space_from='sumk',space_to='solver')['up_1'][0,0] + + mpi.report("Orbital densities of local Green function :") + for s,gf in S.G_iw: + mpi.report("Orbital %s: %s"%(s,dm[s].real)) + mpi.report("Total charge of Gloc : %.6f"%S.G_iw.total_density().real) + + # Calculate new G0_iw to input into the solver: + S.G0_iw << S.Sigma_iw + inverse(S.G_iw) + S.G0_iw << inverse(S.G0_iw) + + # Solve the impurity problem: + S.solve(h_int=h_int, **p) + + # Solved. Now do post-solution stuff: + dm = S.G_iw.density() + mpi.report("Orbital densities of impurity Green function:") + for s,gf in S.G_iw: + mpi.report("Orbital %s: %s"%(s,dm[s].real)) + mpi.report("Total charge of impurity problem : %.6f"%S.G_iw.total_density().real) + + # Write the final Sigma and G to the hdf5 archive: + if mpi.is_master_node(): + with HDFArchive(filename) as ar: + ar['dmft_output']['iterations'] = iteration_number + ar['dmft_output']['G_0'] = S.G0_iw + ar['dmft_output']['G_tau'] = S.G_tau + ar['dmft_output']['G_iw'] = S.G_iw + ar['dmft_output']['Sigma_iw_%s'%iteration_number] = S.Sigma_iw + + # Set the new double counting: + dm = S.G_iw.density() # compute the density matrix of the impurity problem + SK.calc_dc(dm, U_interact = U, J_hund = J, orb = 0, use_dc_formula = dc_type) + + # Save stuff into the user_data group of hdf5 archive in case of rerun: + SK.save(['chemical_potential','dc_imp','dc_energ']) + diff --git a/doc/tutorials/images_scripts/Sr2MgOsO6_noSOC_Sigmas.png b/doc/tutorials/images_scripts/Sr2MgOsO6_noSOC_Sigmas.png new file mode 100644 index 0000000000000000000000000000000000000000..d6c97c731d72d5879e0444ba1e165e8441d37d08 GIT binary patch literal 78192 zcmeFZg;&&H@ISnOqO_oZpnw8`G=g+04H6Op(%sUXiiAiAilB6NOT#W9-MMr~ES*a& z^}DN|?-T#PbIxIohK620jCX^sqTrhXE>c=9YWC(X?nX{# z5IZ9m2U~j=TPx#7Ze~u-R`%~6aj}CxtdA^RTpR>BIR5WnV7GU&;0R6`#e+Z|LF6P~ zsCy)C&3fAEB8*PX=kdE1xn@p$qjQ5|%;d(cb*i=)?A{9HEzK=N4Be3P8u@hsCD+Dn zy~<$v0#T&=B{!65wWn_|p*VsEt5q80Q!wc0njWp9@A+=pH6`^alYZx1kUP*Y=-<1c(GBr` zZ%r+?1OB}Q1QA0p{=MSSK`{Tlkx4^t{Ci^!hur@6rgb0VztKxM|93!OlmE5hKgalA zxA|*_|DQz)r<|*c%i0+t1%A9ez2YLMp{13Jb&GIqdj6`aLEpZm_^|zf{o15OiiWFe zQza@LAS`j*@)PsS|f5;p?rG%_dgGMKx)R=l_=(~pK`IxxZ+g% zEh(wU-tREM>txpbDCPPr1=T_heI#dXWi>RSEXGzn-MF>sb5vju(E)4pKH2CHyDIBX z6v&{tKJA#QcQ2t|Ux_mqH@F-#D7sMmH&eq|PAtwY;j7anLB~HIRn^p7FHRjrF7{e0 zpjX@W{yH4c({%cE8-w$o(5s^ojx$O&Hnyo+=bVmS$8+{6g=0b(jKUxsn?`91ttX&SsDsK7pQy z1iegqVt}Zhzbpm=)32Y;UtblTGRkRdYoo4CU>tI|A+4wo|I?5s23Y?r)Ab$_FX;!y z>WI2tR=zUrciwWV7d`wjARE)((Q!1cW6!ou8Emvh4m%B^w>m%G7W6sYOGHqK`e;p8 z*?pzI{6h}Ab64u05qKFe4M{u^N$&J=N;zuam+KLpGEvkKG%3QKPIB-X-}5RjZ|_~N zq_g6Yk&%q!`5;Vkdh}k>8+&2knEn7HSUtf=* zPZ_AThX1H*cPjp$NjKp)(LZy>eaCrUcke>{(IdI;-d^3x%~yIR0&s21E(cA`IBL&j zV!2iON6gI3-;t=lQ#3U4_4O6>*h{6)%gi)Qk_kL|33){#ARw^6)W%w!l$iMKSItUG zPh2<55@d+W7IV?vQQ}0|9h!L zulA!lqWn)=Y5F~to>Pv&t~O!vF$af-N6TN8+1Dv~6BGx855wu5gXinv*n$R`M>Nc4tYg;i->yEBua!q&TA!B6FKL-JL}p%e|F0NH=`&5M~)cpJ9XJMLMyuu~`=vKHKv! z1OM}C|AIrrKtWlVqZGV!di$?2QJW>3KP-}8o0(-jBn$zodeqB3u%-z^#p?T=nN4nA zF#ZFl(&$P2R?cznJ2lrX?(PNI+0rqiV`Etx8yl|F#D6U~-a0Vl^={Qv_Dm#{cq8K-)fcE9L$T9xjT}HB<4{*4NI`$wFkuGv+f7oNoQ0QI=}NF|ce5 z@g|oE(+g4x!A&CM`Nj|T=9hK)kCYkW4CSUhv*)BfFc%T7i0WwUB0W;w<|P%RV(wPb zj2=l#w+#tPTrabb*|J$h;8sv|%PGh9nd|HtnFK!ir&U1aDv{LlWkSor?I@93&)0BY z`Lm0f*l)+p_cL;ED@sigl3L^xc$g3NC~5AiZ}qt4PW&jIdG&$%4bguAUntPZxlH(x zXYD7UHx|49mdT-WU>Sw?yX$VgZrP9DCRiKuQu)pWe5JP(G}A^N%jEkuhTzNoXVicj zT=){n;#ye60oZze1T~jsO99w3egt*BrF;tqJv3p(Cl0a3)p?!kDsqt&JK4#fE2y|5 z`=ap5AC^!5$O?%F`HyUpwb6uGHkn~J`o&Qkq8KkILkH%#VXJsj zF0sA0z130@1p5bG(QPyz{I^^)i1gPIbYT5|f6#~Zh;4dsl3k2@^UZ`mua|e`Qjhss zJUn5-yr$?2rSH)7OZeV$pks*ppZ0io!5QYg&SztMi|{7Nt;rS5`Hus-$(zS4_q1Vs z+)mu88b1?=qrB8Q)Qn1=!-dTeI)vyDdI!c4*y>Ahxe0djxz#DUu z7uwMKJa#3!fEAvYq)CS7ZMHvQ(cT#T>znHXItQ`;*6zzk%*9Y_2F#75y*&GtoVoqH z8ibkw>j7|VxS6eX5J=A)isoO{Lfpy7=C^4AiR7Oel{xfPSmU2 zcYmeM#&@6*SPDeh(z6m_4mwe9`;8j@=`?lT0SN1V%hN=1Ke*Q9-9S7A^3Ooo(i~1o zd!tRg&O0}+d}ia<13vwUmyYj(54c9+D{sx@{5Q$oJD(^iT@tuY>NAzX-utevkVZ`~rnb`#Y&0+@MuW<>= zqm!KQ-B`_2KbZNDOXEjMUbheTp2n(0E|-`ZsppUGz6C#T$eZAQ?o@iOmGk))PCg6S zNu6}#{ve#kIv84s>jR=N{~#wY+4G_IvRN*_ke-{6BT>@^Erj7QsPp|%IsU`S{wTHC z22Zv3;^I6!JYJ~voTDrbUpZG-R~D!=C<@WE<7D37z@k5KsMZ;=>{LrHXK+3RM;?{< z`?kT+@EI`1MRB{prnjK@b&8-&<3^Z>piJK&0#E9VFFuo70s3sL+K!^bz8=tsuKvvZQYv*-9zj7t zbpr!fOEfyyH*I{9COE3U*m!u9@BACzJ88cOF;>c0vtiUY9FeusgylMq7DVPX@o`!C32$hec2gUuX~=wRmOs&TpbHsC*N5|0$JJxQ!?bd!df2&8<6aBiehHeIdHwms z=z?5(?b)dByTW2HG;I%SEZkgPUS4}Je;Q0kEu4FKeo~N{DsB(EtV;Gd=#p$8KpQ;N z`AlNt;WIje^Uh2`X@)W4^P2e;gE5<&DZ6~;PIX1ry|H^rA<&;KuVPJeDr&I<&M`4K zIa=r^I0UVk2?5&iP6ly~g zR`=%_TI&FM#nYG8G|B2f;#5fAAew^$8N)+(#7f8P{xU{{q@ylcXxLT zd~>YpDw#3bo#T@$7{!XqlcOAh`f2WZ~%k)Trr}d`7Nt%f)!eqH=+5h zjY|(A3jsBCjobFy3lG2*PlP#mcr*a}t8-b80Zj0}AzA4uc}BT+Mz$B_W+@NW-q^mG zFG`aG82H>|py^}TWM+56t1)tmysE)G6fO6`7^#4UPW%%KL9w5hP1b*uC+UY~XW8VO z1e3Z422xIDe*UD?t89U8Co6joy#8{=(gw>>m6V>6gh1-qhhZYvjIPd|886)b#BEL4sk*-S zzVW$F+viA1iRYdW+D;oESOruO&^Q5zp)X)%ALn=D=mj1|-W*7xxbP93&yPUQRNU!y zaE3L_5B2F(#;Y2Ngh>$dBz}$MGcI71F_pxk;6H^E6}7r?owpzlh%g*JkvJYs(p>F@ zpk@%Fe4EaI6~<&5IxNmsV7_0UUp6y0FZK{FNCOmluH@wCXnmnM(B?TM5M*ewhuvBy zz@Jz3xOxg(#F}a^*RC@olDLF+==WSAhmeRyf_&nX`qpSR3HP2IvtvHAvrFU*(qw$9 zCL#yoNtP=P$!_Z8#O2oz|D<&@t**|iF$9bYp)cwZ*D4C6t|Y7(6I^@U6QEOsMST>` z=oX|LHk&=w!pgnV_HrJKiY%*Vhs=oG4p8_On&i&H%)Eor1#Kgx6w#utF1d)PsDXOo z&e5&9n3$N%hHD!1*dwNo3aj#(s>9O$PEa9nN;jscFKZJ0Xjjd##!F(Iya~AjiL(aI zS8z}#I(^e?0rVOG?0ZkK$8CelDxlU(h)h!4M`U9=7ruk z9g?BZ`&x$XBXCM(1lP4U)Xj-eoxMfP)`bxv6c<+LwGZ0p{A!`|G1F z7I4eQ5Qa~T$sR8~Jw1QM#}7vNA?X@V|9p17!rnEwS~2+3F(7)rMHz|iq1V8|y#wp! zcE^uI_bb87wv!2W7beC7=@$3AE{Du|)`2>l74oDsp znZ>bQa*Zk(o0%!x^+5=EEe4UBCmjM_y&kKnA5Bjs;-!9z@L~L6Wp;~^o?g<_^t7hy z43OF+O&%ucYGbL|=gZD>Qt#3@kT?t9Jz!iu+Qs+P>2dnQjp7=CA=!IJ0C?d>QT5#L{6*jV(}#qhE2cJox#O)O<5#se0*F0mYk4KplRT%L36n#vwKGQ z`3J`7@4L!t6J;i!DMze(>9=-v-Y6DNT=jj$QpdBZ>f6dwzCdpzhAjKV1;Ho4$aBidULyul;6#3V zSeJn3YXB>JXil}ZwnmPsDX;e!T$KYk+L+VIH@G`y7#HrmauiS5SksPyPeysp(Yjiv)fr?!G+#(8RAfF zNS}Mvb>*?Jc@$r9nr(ha;B|}i+4t%mk9kRp4>EWK{|Y!Q_c=o?O>e!rQ4wN1zq@q; z%R};TmKBRED+{w?x8EiP#k_H*+yzA9Eq(Ytjmi2)9$tSsz?0xC*woUM%;-2`UN;sV zo>@bjH-k%__sh6?#!l(BN^}9 ztLsE&N4U>sL@4(NS-n5)sG3S=>sO>@>oKWZ&rUJV}bt(8t^vSkT)r7|lMsjpaTY9xcTUU`L(p@1kL z2HT_L%8i{cx8>jeQ0=dG)8v=^1K zF3qK~yD6_3!)w@GA~wN8nj!AG)fs>}Th^pksi1Ow{sB^4Be^OL!8{Tj)ASmvv|RRi zx}QCa;@H^Hes$l&(C42pYw6K~Fezxh7P3dLAAb!g%mkC}E~s6q9rtnje&?E#21dnw|xg-oU_zN3WScm%9QzE_3)Ew=lmz_08Jfk@o@e0zW|# zwr9+#Z*4lc=Z(zfbqyX^+~@dhexEra{#l8jsWf!MK}+NDq4SgI&@Bbup{?@uQohaC z_Sn#0;Tj%{O@IDw(*~S}Gu{~w>UT`qzCTtQ`Al!^VVqED!xQI7^?Emr6)Sh3B7E6O zSrKZDS7ExOuioU#VJ*qSDUhn_8sXmrgYHjSuL5(t2 zaF1)S@Q{Y8t1zy>6~#KDY_fj!j(I!~Z}}H_pEncs)_b>KC)v@AOh!HfVPPduog0%0 z%7E>nKKLBayNFFiO{(}|o&_S`gZsieHA_d$ z>hdT?UX-NzxvjC}E>ZhT?EN)rK+yf*w5k_T{l5~LcY&~@@bm9SZvpuQHsuR6uzPSjY1EFoRGtWJ+(lB3Ve6s=!Sd8gBWV`q;Yus)*}< z&H^k&b7W3AXss-Ue*eYcVVzBJi2K0&26V#FkL_f!)%RcU+dg{lAW2@&;DQ+!KrYS=w3 zoYI;oNFqL$r;W~H8|!0ECO)M!?kGcSaeX?E!gT)Eo003K30n*7f0hL0$sehfae+GP zzISlxD*2iaosV@CG^J^fdg&U3h5!mHg`YZDhJl>;5*b&oE9<9Rz5i5AF;4Y$BCIv5 z>vf1t)?v@JWX+%F*J&*meC}H-zzfoNK$!Qm{M-5dWUrpHMQ~$_A^<`zr@NwV3x48% zj}pX$irpK!MP=f(3pAaj9BM1Zxk;x$U^?J=nVS@mTvGeLdFWv7ndC`tr*aS@9=x3Ao-t3E4EP zV>jP`E*}6pEtzHq^|-o@qn$qDMsH^=8kb0xns5p~9r+5^)HIb_8+b0MP!&K(3pn3G zHk(OZ`oetFq5z7wNb{$`%O378CRgdBSf7#(sIcCW@XOL1iO3hV#0Vk zFy~h%7#>!t%D;WtL15_j&zT6M(u0eUuhwQ0=+L^(!9?noy-z1?hmQ|*sjQ!lE-D5I zPN-zQ%rS|T#y4aG?*3vo^_}9(aM357IuXZdVQN?<^T{epxTHD_Zk{ObS9|7DjgXeN z$P@=jwc6E_^(jT*Ic^`E}4nk!{j0w5<_! z>MUFP+GL;On~qs~V|u&Rg3WWW{*te%XDe8x#DBpc;0C&<*}Q8yU5aA3P(bD1rV~^*so$Am*c(m(Q(18Kqz{mwFbDzcJZ3c47Vi01iMJ z7~442@8Why7)Wf@h~873GqyO);x|57wh3>8QZx@}e=9A+loh=pP9I_6w-;%C)mRCG zs(5O9Zt{WBAix_0M4#_FDjSO7{qK(f{%Sle|G453!yCuH)V!-Akrq}yqG|ODLh=~c zcY1^>%E!*c&N{(T?Of9PrGx>97nF)t!pp~YxWT;$6_>oi;po*Rkbkmt7C!jV zT=MwW$0RqkhpU&TVb!v2Iwnydk#{Oi@MCCUVce0U6JxP^19^CejSv8KH*^i6d!4jP z2`_sfAvhpN!`nzx1Uk*>>{L`T{(4Xu*ofUUm`g8NM|5h$54LGxf}9-?ZW zS2393u3JR911mD<<r0sn#yJ@)Suo6{uI;`?Gxp*d58;37e|gSYD}@Vv0>GH z9&-!=3iN)S5W}F`0kV|lkVM{g5<^fab%P!G#gZ_~PFoSdak8_4_|Y?GtVtP1|63 z!F<&h*Z$<3*>MaE4cR>rtjm%--@yBp+)eni>1`P0>!TWP7n7~vcj2=tJ(i!*kEk5R z$0aG$5!Px?3)z|tg%9j-g4M}BvHC=b>V#1o$ydo`DVndb*l{fx;`dy_Oizw6KD1CEb8B*)J}Cog$MKmW4#0BihruD?25rlf<0 zErSgtoTM+_M=LV4`-9ay_n+^Mk>fJiJmWj-Es2$PkLNxEdlyx*Ru^25rgkZ`S{8Rf z(kYb$`%aftp5oV=)-a`Cge-UTb@@Tx^{8_-eS#H6GDhmiA!b+!FrQj^pl|{%JXQE& zcOlF`Xp9qpy_F3s(hl!u$CJ<6&`s9ZZnxjG23xPo73Bx?_iG@JChtI^pYP1vFb!7} z2W0-CSBEdEZ_Mv<)f(16Lp^3}kI-CpS(8xALaC_zpq^GOWW9BL51WMTAwhv}x?W#K zUmdTLB42xNR>MxJcUy_xorJM!Hfq#|seAZV%{ObTv_~q=y4{q-gg$EuGe9Puo@x@s zv?nD5)L81`1 zf>NpaHJTzUa3Lprz*Et(Wru6{GGh{vuJ=J_S4$Zjba|Vi0IQqH;VVXNjt~Cfg~(P; zyI88FiT-MQMe~kX&F=L$bJV>Wu^V(R01RwwbdGs!*-+VNJgqCAJ~oxi*w^-*HW|?Q z@VQ3Q0F=5=7P52WsF>}Tz?BOKCpkJkUX0IvJ?UtCfo2saI3TfB!o&U0B9Do8)1e&P zY*I4S-`D7fzOB1g#WLAP8ZcvOHK@>Mj_#L`VdQ?zK zB0n}0HkRZ`O)|m)raR%cVA+uFJn^=ffv7r|1RZyYL_^z-d*x+_U>|1GaA9 zLjsOoNx`l?Yc=h4cQg@ef6@&+-v#^GTl^3$uh6~q0Jm0|vsCgmUA~*9xA3tZ1vGs# z15O%x;|;P-Bdyf!M&J}Na3-MUYnS&_wax`VyFRl9wWtfhov~Yxf{&2UOw}(dD)3y6 ztNofd3bbZY)2TYJhQR^@{VhG@!Q94&S0arCFJaf>wQib-o_+@ErgRV{pmBkg#zLp* z;yRHyi@@;&B3;wIxZN=S4K-BpJhW@Gt59g0&ZaYH3j`e(i=nV~42TvbSV;LRsuK)T zY_N(CbWY+mul>rDVGR=ntD^99Xihe#n(}u3;7C35uL5Re-krmntd_K#X$5qeuAXb< z#rOt8D~mmp`%#5!JT@j}_7;zFsmjb61P<-`qhsFf)ajH>`+VEw6;tZb9pZ;TSY&`P znHrX|XVFk&xJ2P}3t7-(jue7FtIFNYW9to*f!)e?`sCLn{2@B^UHIJK%HsM2e8<{g z$mUJR%PET6)`;+|hT6lvEoajXbEIM$^j21kohoj^9z>jr7DAghBY^b<9rY9U8L*5l zQSvZ9zcd+8?(rD##jEPXXXVg5&zSF*Pp9D|Z{-{nnk-O|i0kdxd%_d-yu|Wrl|-f-E`!LBN!%PpTMKi@0^ zsETRxE>oWzT!yA_lx8PybjSMp38wo&a6_G-lLmPRdwa&b{P%-ZXwDk_1#Bh zv$2AsDw+3q4wwT+mKo@* zbT_K5NT@C{m>**_`F;S7-Ga7SjLf5lL-r*E8;ln9F6pk8&F84B2|F(zg|cMd+_wOj zG`9-tZA@r*l|>aEc{Qei#a=3X*jzUIJ$WD`T8)o~LmEFjVZUE4$~9Gl6X zDL}sy6Ev|=D{cCh%eHaF$Qm+X7%hhc&Ad*Q*pyBdo0$MwX068uYCX!jOd@PCw7PZA zedM5}_Sz*H_aZHoAEnNgP5KIbcwt6)5h0YDxvdy_=P{~%?B?+=YO=vV~xPWpOZVr&!K00o6jzpoqjLDmM=YQ}1S z=QXtVg0vB|;dg?7+kQ=naVKq1C-_&T(6trM&uV91m&2rG zC_l%xkQivGgQhYl!36mnlq}8%Hf9zCY6Ewkz?AZhECD<=UIiTAz?@nl%6%%>QEr%^ z$2@D*zfmAyDGh7si}CG$4q~tc_>l{4xw-v~7JPtZ-x(-*Tn<)Cug=$%E;}kQg2>~y z;p))^$~z_F+}`;+Q(y2(KY+vvz#lNpiv>^5*10FJNN&n8YJYK}Gc2#CcLuTEr*f^n z?N3OyKnpX)h5`Q(Ojv){Sot!=UhD5Q))Z&##L!6-jGRqiJtN zDsof>D69^X@oev;5+9TzOOXTA)@!Vg8qi0%z_eG~l@7JCbam(n#h-p!$g{>!}DU zseN?lDkOj)68K^(1F+l7_Dgt;AM#zjK5R6{aCxrBM{S$;vbh-7cXd&4yM)ECiiVdS zB)9FI@%5__v>Xf60N6jdeg&_!CAi_*L^Fs1GdRD%SeYliuvMK@;yxV{wG~BtV_a2) zK)y%qtsKGJw;FUS`iPB~u6>y6Np zQLgdwW%$HB{PlqOl0eZzQNJK}Y}hD%9x6cHUoziYJ%5fP%Zy=-3h)d^zy(}(bTJ>| z%dFr@eOp-&(}V<69x0z5Wt>bMZR8QwUMzYiBV;Q#iQF|eDv02q**-T#a7XJs=;qtd z!IAp&-hyN9^5obYPw(kX$wzgoXoVOz?|zvti4lrx%@N3&k}%&Y5Tn)DiB%oBdkXwv zBK|3}te`u#xE$qSFuGUUQ#U1)rl~YD&cSia{E%nIQyV1JP|_+2bGwQ8lgFQfXGC{c=Q@6L zuckILj=~QDVWj}z#3Z)a59w5Df&@gc-HBYYph2!>EmgaCdmJPqKV2+Z^hFMl5RG@gN zcmfUjjJl#Lad$y~>$P+k-xw#-@KLHvqdC$B7n^4jwiT1TzMy%M z`|U@{-acqwOoE7(59*&haa*Jg09G{aq#C>iT|a>EO!k~k1Bal6z5x>*$_2?;m&pyN zb_B{V(771g-&qoz^Y!VfxUhctScD{D#eJ{etWh9#c$G0rQCXHiCM&W0V-z^0&>IEN zN!cjfYd)IRI|3* z7rH&}{d;o+n5%Me4@piXu*%+;ko<*Lgb6_3sXY}Wz!kkTbv8Jh2_TwcMvYtNLlWeV z{K1TUU2od(3mR6jyT?Vk)4)IN-1UN8``*m>eEr16mot-2l=~E zel_1K@Cr@x0>w2KfbFagg_GBorl2JuUatM>RH1AdS$2yB=`UQSJ>BCB_xKoa?~}BP z=V|%ehL=J|T_hmteNz7RphX#2rtKGn^2x^;&>GfP;jp{<^6SjK*Za(DZ=MeV=?X{N zW5E=M6Y@xajo~6|$__P8oj-nkWSU|9d1C0EP0xdajoCT(BRBXoFtM)5`!X=D>(k7R zMn93ALo0C61DU6DCg7B57o`vZGXd3sVga`YKs09bYP!#2ebEC9Vbk{a-N9?$@V$`b zQL)Jq6TJ?k+sQN<^V8P;Ty1{?1?wi!nq!G_)VYuK?I5$cq5%et9UhA_9zWT9(WY{J zz+=E2#hBqFg#eNFr5_e^ALbpqQW8TJ{8o-bkO29;&VI8g&Y5T$w9t$Ll5I_8KbDYW zztPvvHzqvQu(rzfRRLnQFvhx;s z!iyc~)v`HqS3kw==7(~Hm_Bw5sxOpctafLm>Sdbchdplq+=T*je@k&aA^GZP9|{t_ zmg}Od`t_1e6!e^i4V)7l@FZL$vUzd`XH9*Vsd4fsdv3=|C*x}9I;HVEvZSu=n*XYf zZj*rzFNTxk7+2v~T1ACeUGkDNEH~6xr$F`AuQxVH`~{X}xM5+jbz%IkU*kUYzD@g} zMV`~hbS}Z0d-vo-=DCBjeLMx<@}-Y3$`U&73_bCC`pscp0O6ywDNtZESdseWvo3pZ z%?$n3R1=*rp&&j7-sM{Z{Hh7-LVZa|Nz5QjaU9tfWZcIpERi2#4!kF03!!(2iSuxS zjrNYl@so9aeHM{P^h3lea%Owess&eS>@kkDYcHE>PjQja==L5+^f<=VC~pdQ=Y=%C z<0G$~!FQEg1Z^?B^`B&xVo2Eg-JPR*q@q5& z`~Z1Lo7So^c&uvaj5BdJztWAEm>9?<6g?>%yW&yi&(LCSbkiGl<+*II2okDrK zk|OO7-gGSTD!b4w)BH)ym4?{plxDtL(8X!2`SSfvX?*50EDy}DMuUx|dvSvUimGxr z@vMVBkY5ilwmCwui)>Q8fI@i%_aO33VFV|H5a(Z*Dw2}@uU|Ma*%wZ9z`kcLu(Iw8 z<(RcP>@@zOq7&k4w|%q6h+wA^3zRx1FI#J(^%F)dS&5kN;G=5s?kDdW=k_RKMcsOt zk=i=C$pQEa5KlK`$#^?o^A73VEwcZ4-PU)2uy~VD> zl@V0rl+9rG`e%d;Z}K=dZ&b)*cD~j_J$xr}LvDP_DLhofM?y{UR<~5PHT3i&ktMn~ zTF*$|K=N^bz{F3ZtmMl6m4O3%Ip zJX`#gTo&h5)4B2|#)t0`Z{}lsVA@L1pwYcw^$%$uXgZIf7kmCG1Yevz?(+|Bi-W}w z8gNSPrIFENIXSs-N}jO#ox0DfeF-Teo8H3I-#HBW!2zhY?pT(|GLu{2;^4>N;Ne1T zZuW+~mwwmh^COiuaAp)~5)}^yF>!RoYZcvY*!S$cjq-b8^1^-yqz!3LV20put(7`~ zCLK{pDbSXzc3r;=kM)bYq(lXENOa7Le+Z11FMN?8?8&v%9zn;-8vaD#`yRNy784iO0Zxo*ml@xjt#cJBI5`8i zU@HQ=25=zyJF9jJpk?B2ZoJ@L8$Q|Rp!Ja)&B>DeSL(KyF-y!X zOV385UrS2*Y>9mzpi8kc>*~+&>&R&N5YJV}cGnK!_Y4C0_qhNl9>$4$Y#Jx8FonJ35`iqggcy7q@yC1)m`|5$uon>nL0v7?7c!9Aj5xVKb8hJMU*hI0J;ww<@-# zD)9m=l0AdLUAym5N4xVWX%Zn^rk!_m*~$G;BPm+oOs;mLa*C+$DY#PFpB1OSQGqfV zOnnhfBbI!1dF~8OQ78KzYrktedh@o@+I1XOT%l2h6WS5b@WD8xs2PnnwV&nB;KL%l zZ>!?cRp;&E$tpfQIOntIF0DdcJ5>0;{^S>ZOKCIVoryW{e# zX`6w+7z+t4XGPkKYA6Z=G|n7ng5ci)3*P|FG(+ea8J%|d3i9@05M-BHh}9%q8FbH)`>HXq!=Sa!%>2GN?aCmx+Q0ph+ z!pE+k(M{g{^D^#%&P$6mG0GSj3Eg_iTrDuuWh8p03g1~D^Hpy#_VUfT+%>Pb-j=(| zsAM``s10G&DqSeVY#MPKk|K}wPl}0AW!rURHfkn&8Q6R89*x0J^fq()lf_EasL;?G z!M+^_o~I*4tkWz0_4W-SA8$!NaUQ#==+d8(IQs158Fc3433_f-g?2fK_RsK9b$Iuq zca7tFq)c0G4qjW!LUAd76~7OY0u$?x#B73;`J5dbl!V|@sIE8YH+mUWOghx+7HXBH zPi-)P@xa>NaZ1 zWdx~Y_fJn&(VF-HHBl2+l(q2v9G`6O7Z`CuJ%XpamNSh2T|0*EjHILR$k3sVi3`E{ zU1itdzcB-RP9C9M^k(XpRVx6bNvN-P&yuKA7H$6xt`3nA-A^$rtL-hCW*y0`@{6j-Sj7bUKzWoLh5b#1i4GIeS`0=B< z6&pZL3T|`A9TJijaQ=PZl*jM#NS#4MU|vYP11;Vd5ivZzz{8G(X)t*}@MOX3+l-&$ziK;@uS_)9UgMGgV@ z-41ER0IBRZ6CyZ-Mht6U-Q|rZSQNq-;$9B z^z&;aa^H-%ChNR;tjLice(0P;Jp_x#l6>)E^3At}1vU2?k!O3(^_T9V!a9emeI~2D zKg1`;=R0Wp8{`eHU$uSoE`f}a84^by&M!7#UR`y)flmEgbZ@XRzNf4=5*w*VYg5C8 zqfXu0x`7Yly z1tXQRvT_bUtQ;Sz6w&}|OrX=;fd~GaPF(j_gLfUT7alBpys2GtGVC@ZT>6BpqUs#? zp%e9=uX;55`~5_+9Nhvd7ratZ;cm;^-47*n5f4%iE58il{PNXBM~*Y_+K;pPZkDwp z>(>s4cUKUIIDcP3onYk-aAmr7BfsD&D=R%M?FVp5Ia%#6>i6&8%y}$tbppT=1%1m< z+PPk19N+0V5b`ss=TEMwao5qtj zG1C$4E|d0;@}mB^iRL{>a$5dc!*q?KNs9NUPoE04D}upgUMD9frK#J?@G)0MM@J2( zt9fwGm*Sc64Z@JXFJJD6UY(ktj<>!0OLRzm5&|soN<%x(elI=qFliG_zVu=YdP1GD zKe1;>U6)RmQP1=GYn0bJomE~k+aEpTw_t6H0%t@wih`Si*!eQ>g8JiTciqvqKD!3vYoV+lN9W0CK^7>j3AsX(t8Xg}ru} z1Y83udt20SeRYoifOjj=pgRBeN1L0eU{=gLJ#CFNS6b`&Fq=BQm1Y(}k_UbphmI>3 z-!#8n~$$cGZxmE5Y-&4#;>3W(Kwf(>pX;05dnrpl%yb*UPJCb!X}_pzd9;6 z>tQB35rK8@9WNPAH_W>N@ z_1Rf?StWh#`Hfde?E!36!TIQ^m$tWUGs+tZ?PBX)Iix}C&Cgz&DSuk}c6eTgSJPhV=y>k<(d_|I8JkV2 zJeW{)o!oUD;hpB@<{4MGp78O6iS0Zz36bm_35o!(KJ7|tQs>oP48zVZkC8jG)4Mj8 zln=pmtOnr)q1Q@EgmHRq6hiJhHx`1(^$ib}yEIM<@&ri}0>meI4bgS7>u00;KB%Pu zVPw}Er&%bXmZAOQ{+28GAj~!JZLX3v__6J)%gF01&aYAeLQjR;Ft6L6Jkixt>0ZlR zeBO9!F1wXwPK}0voEJ>;kRO7s%)78_(R#OSV_-ILoj$~IshyaNtUcz5LcFN2AoyU1 zk3m7eis2PrQ7tVmr6VZYvK12`HQ-Z1C`5hu!TtI?^DhxOkl!J3`mY`h;jAw*`LoG_XV8j>)m(F zkGCd!lSQU&luS+O0fc6N4+*e1q6VK9QK(y6e&IFO;MtqNx0^YDPr>yuOCc_K>^zC( zYnNm=Tc}C@QlHr9ZBX*;oSIJ??Oa<{z~`4>Vbvjf$V$E1TnQgGb%si#={b8(|5gRo z-g$*)rZ+oJNXeaxn8Jp-maHbWv6D6#(Gw#?*VVjsQ>Bhe(&!8Pq8ED-U}6%lUlXIr zz(au;CbjIr2G7GO_6ljK=?{&f{N^hxkbA^9E<5R})R#UFc8f*U+%NL+e(tTUDfvZt zeK}MxdO>JV;uU#9i9fYf(Hh70dp3`o1tIHJMXJW{g$vO(k*S^n3bDnC=Y4#ut* z8nFz@TXe&XxZY20xfhRd)=b=Ov45-gVt3Bnn9yFUbu0`FJxT+J`a;{A+mK4YYU%fu zYr;-@E>cdd3+yw9bV`H3-plXnqgSlE>Kq8-)<^F1s-pVm!n>j(PTMgg=iPUYa7{W1|RVl6b|$>}hpX*B@STshlr4@TS{GzilSA6~rDp>{71k7W z$K!7Jnu_gas*Bdo=5d^i5{lvbN}LMc2%A`_U!RFWDdGq0aDu}{bVJ^fCF*CQyy_PW z`Z(^yaRV7Y8E_HDzZTJ7YBA;#ZXAa8_98mJKbVRY$(|Mes5)tfwpuTdCGC(-K>r5GR;0oEG`62A7-ApA6&>%I*dpR%=i=wfu+*N+lfK_Lz5CcN zewT34k`rpw_FZ9()C=8QrtAP6-_go}mG(Dr$O>JOS!)O7ps(fU&4g`&$;v0Fa6Aj& z$%07jcn2nVR3a&U$uI(Wr`<*9YHqQe)XuwxE)c)*I*PyOvXLyv+*u3GfIZFE=xns* zM{V-{abT%fMYm z{_Evm!;0mKgvAnD!rfHil~Qqu9FlMMuAj-@1tpF}H4wgPtrsWwZJO(HW9IeiC6!HG z_+kYNJ-qe$?g4;P;Jcd&f+#p3%8}DG`!C)(^X;A++ytx5y-9xIHX1^cV9sfWg5@JI z#fB;0+*aD~Xdonb?Jk3=W-!Wp)XdDht7C}MT5p#QrPO;&vlI?7GAbYem$X!lqzG5x zq3L9D464S;WG^Mg;u|WQTSG;jgN7P({qo%CyoiWCVV}x!&Gy-hnD>t4+`q$~C&<|R zUMyxruV>CZpE=&MfqkV`{IGe6RWL9X0pwVX}BB?L6K~!Z91H2yV$Vg-=Xk)Nye0N4Cj?Q@A3`RHGd!2rgNF9we(^_?3q;0u7l}mWT)+?|rA;f^zMm!JcOa zje33Qi=*&Vc8%Ceg*{(!*jtsOn5ediB~ATq&NcVtDiPAs_`+Z_zIYYx#ZNN2GU)hp z@TL9~L-;d>@(&NOMb-XYCds3qp{(cFTnwR}+6ON7QG5Hu(rfh)^PZ6ild}-CP~A_*kGrVDsA(zko&-JU zSsfT`3?=pBe<=~?KB!(kHgMO38Kxb*1scj*ccBs=LvdU`ofa-V0kYKRG;aoW}eSQ2+(li;e$K3H-gVqnX5f)!Pm1>gye%^WbWP$s^K zn9h!0!rR=x+tsEnkFCvi#LJ&V`H9Jd@Mh}K-90JZzjB& zb)&C(eYXHl8u!G(#X``=2f|FJ0XD*4pxC`Rrn`CSz}nqn*0;gM6_f0<*DWvQcwWX+ zj%G|c{SO-?kn_h#M^<9Ei^*h{2~&<_t21G;hYx~lAP&C&6&~TwHF&`04>))Z*ZD~e zb@ZJECSn&^yw++8p7a~aIRzwieSYZHTH|6pfL)O{-bVNDp9_grJM)!!({tkIgB=n$ zSEhc8w3j$=^s(|V98Gt#FeHA>{@hc{=%356Q;>Pk6QhJNM%=Nam3L8Onnv^!U%8nn<4e@^JHN((_k0z_+-Z z_bFtSd;VS+&zgimWxx8smc|P4mnp}I;CBkBprF)`^JYX=aMHY&veuXn0~$RHSSmM% zSbzrtiP{5gAd{iM6KXt25I-ENRJ#%{&QQ)PNVe2%35?OaBa~Ws(qhPRWcQuI{M&SI zF~w!MlIh&oHkH|y!=rOeXK{IYxJI4%1L0%^HdY=73p*Ezv599fN@zSCzI<3ms?E-7 z?SlCqpp28(8kV>`xjp2F#7C6+;d!xVrydqc(XD-Sqv%iRZMolVDH@~Vj}nsG`;bM; zAdvw_?F;uch^jAgA_x~1MU|^&cboXe`VteIMzi%y806$7b)kOH6hCnog5>l?y3A6L zzlDC+mvFN)m?V9EZy&p9@p4sBbvRa$hdkWE8u`K(JN6hj{R6-gE^;?TkHeGy)IIyz z;U=2j<_E5lXw>+l?(>!ln8zU0OZu4yLIE1;x=-v6JKx7Qt^?iy(uIPSH%QeHxVfg& z@XwbiheN7#DT+{N{!73*E@|{1qnmXbTD_iBxMBNbG|ik?O?~#}t$dPB{#X5H^v(os z*H~zhSNo=|A9+a$+v7|AijWddn7nk~JO4J*RZ$oq0=L&+T;Nbbp7jZ$_Y`GY|8<|_Zge`n+aVOg!2n}u1dUxZ>wtzoOJb$p1 zH=zK2h2B3Px(wlIz(+*SUXu;8O-bClIvKL3*LwTM&miJdIRC0NL%C4QkO)g$*KW8PQ1Ygo{Vd?$I<$8G<);(uL8F)dDcEuyI=% z)XcO1Zhk6^NeO9^Bw~@b|UMwq1OWerAQCEsMG$_#LP*DVd&ZN;vWul z$NesK;E6{~b6M-Q*ble))-D|_s5dga_czedQPdFj3VUsh1qb^o#%gS00yUfm-1JYx z13oMBUTp8$<+-{aPrrc#JPCOPvXE=H6>jeMoQpTUL#=uiQ2 z0xWr6J^5EUC@w}j?TUhoDE&9VZ_+%J;$O1j|Ba9HFg_KXh0B>wIrFT6M8%-< zjhvm2L@co6UhK}SwO565BO z&_n17#!&-HD)Jk7A9sVyDfC@+CY7_zzEc7g52kk00soLUjlW^83ErH`2g&xHsVUu^uy+;i$jcv{4GRkNI4Qh5nDdh784k{K{k{teevOoM6yG`W@BkCpOz5WBs zQTWwoiAvj&xAg&W=M|X8r3wPipC!hp|A|E@E?(}gD2wW(!eI-m{Kwd;ex4#|H$y#T z2-ICOA38DYxgs-Crbg6bW68YR0UPR;mTD_ma{$@DQ7-oNsl@2^LhH#ZZ}%8GG0lxj z61Uyo5r>g&k}#jl(KpV`!xo{_0|Gi}SIjIc-lGO#tXkS)v+w4n?<5F(*$d9C9snXQ z0t5>$&;MAwUswX~w9FcWO5$LF43QOqx>CFh$vwi$Cws#-r2?a!X{LVtP?bR{We&7# zHBVv?Mn6(wlqpx$UBlD76b*1YI+*!>PTX$*{4^9yt`qaZdQ6^Mj@y_2axhAu7!~3^h+C1wlS2tv1YXkm2a%9{-B&1@`8v=HZ=2r8#B;baU~Y z&$>VU=8|51HHr=*+$cEX8d)lnKmqbMmWqBjxTpoYaP#BzdCb}<%vvS;*sZOZ!uk2d zd0k9U1mgskJ6t)f)PaA1Ku-34WS4t*wZ4vtLlwBrYXc2sY@OMug0@%7ni!MKRlWd{ zkGRUmQpiIQ&0@=agL$E(4rpexzhd0-51|!$SuHEyv}U{A{KAw58h%arN)94`j3%3Xi#dY`XmM~zcnof0zwZ+)M4df z_K~gZu2VYG&~x3}2aKDnB#X|3;~J+Mr(IwYTg`Xwajj??_(sfsBf^bw5(+H+*YGn$ zzjT?+Sd+YM_{!~0xr^ogjGf4qNT9>!o3B@Cx%^Q|@12nU=|0geEv?pw%I;M7^89{- zt7j4?b@Duy=b6Yyl7d_6f_Rip#%9@xPO9TCF1&Ajx^MT-Gi|>cA3&W2TnR4xiE0t&BE8M!QHzN}_#);MnM`x0UcH0){KsNv9cOB==o20RT6+KF zL7nYs-U+CK-VD{$)sp?)E7*)tv7GnYZzorS1#$!Rl}rKex);lHI@(DWu708;mo=(q zD;0#eFOGB@WR@rXsd&lJfFzoNIQgzww#?++ZMrxnzwMsyr$sL8kJxhy`*vNnFTee| zIvq9wPCjti`fX}{Jw{uE{PoVx_`91tu1&=<##1d%j)i@+nG$p%rRsExjM6I4YOgw% zrM!?2Zp-;E)!ct;w>KP+_Sn3$n=r!A-FO6I+?pFNar#-?F=KUb+vQk&PN^<1a5&7Y z+5PEK!ClcIW$|!W2Yo%=Ou<@CX0+^ZIP72Hx0C7E;hNM{^k*PkB7t4ukPlxn(!ibP6%^A@{h|mu8{8$W6~$|L@LSu|<*@!< z6*fa~PSVfR>`v3I*7V$B=-Abzlxu#s_J;}B@>pkT$%{$N_&qcrURPeQ zsh9=5iA}k1Te0u+$0v`Sai%@0x-_DHa78oB9#GPx3QNml)#2W+-FD1CUCigvarp^6 z_ZgaSw9gof2A*0{&c@nk@mwzcRehX>)7b+wlhG2v-Df{C zDG%FK6|X7QZtmwXKYV}w$_u)Y+b~GKaLm*!`<(o_ zXxUcT9<|L%O*O~!aM^e@9{stN{`PJA?n1^#2G{$A92DuQ*Y~cXq#*eK;)O<@dw%PA zq~NVR=e>`OQO5$8JH-|@^{}H~LcCrF8_nn|j&7r3$3Z8qS&D=5RJO?-Ym5;nm3eUw z_|p^c-?RuK{l!U0)BM&1UmV_ztvO-Iep}xR^JD z)?o!JQRZfuJ*)LIJX93EH%^P)QH4bvHf>k`&AkEw&Bu{Xt(hG{Nv{uR`cC2^Uh16G z@kj1=t-fPUI9ry?;CF6n+i3>THWkJXX`ar@e`Q(@rZG}CyqP6Oi+b~}YTdWPhGR&K zuAZjSmb?weT2%&Di}<{CXFm>hZAt&7{u@^WAuI@yuRg-(5)!3x2=#LtBFi35?05;q z^Z1%c0({|su(3Evva|8{#yFSeSstl6*B{1hpHEF0HTN_SH11xpxE%|SkRZ2{q=9?c zhT~C_4aOAFT@qtQ!;iY#~G{Bbf!w%bosM{qBW1rq-wCjmq!g4!!)4#ar|$*A@3UqQ`glkQWkq zYC{qvSK9Q*#~7DltXeA|Ygb!-lQiXN;-)SHhqEuV>A-fbm9ohJ--I`1@w#7!Id}V4F~gWp74SqhNAy>2jOUi>s&IwJF~p`*Pez4`LF{b1t4>W z8e|BHz%gb0(+b{C`Ywux;)}q;?^k^PtQPrhI6%}FTk_LxkYYkpv>*Q+xnH86OP120 zyxFc3VYOXT)N9dw4r{oj0H)$W3;Y~NHxAft?fnWIFf^5`5I=-pcPJ^gjU;`$)VAgp z$a1nrZ9ZAy@K&t6E|w>nMc0tv&Ah@ZnFNzC9Cp7F&D3wXkr`@uaCa zjt|G_c&epasLU9(`oFp0P)^B0KeRie#i8ZE8#}N-x|IMF+5;#wcwt72O15&6hHSha zY$xb%aJ*3W-A?zt#7*syYl5Y1b=J6IoI~%FOx65F%QM*OIhDoq>fdbw=D5r4ElJ*z z6@>HkgOfM{BNRROM@o`}O=mzJ0+5V6_*K$#Y!axgnhA^a)({z6L+a#RpaBFQ*STka zDNXB((8I*5;ET-*=Z5`uM=kK96Rn!YP14+VpQujcVgq{I>o+VFPx5e|9c%X5DDn{A z5|3(Pc;35vqh@~+di*u>Ys2+w{=hT~a@cZx-i>Dss)NhYYK6a%^j#%9+V&rl1=nPM zrpO6y=B&Bh&ru*aWS>jp#UmX;&EJUj+kul5k%Bt$PEa)*z#~QofQy^{M(Xn1c=L%Ol(hQWPp2 zKxJvJzDs!F?hvjmApVh?eX9eG=ue>?Lk&SjzflkHCt( znRdk>9ppPXR(Om3`(h7H`MuS~%G|ZOJDZ?8pWeqY6x0@@9{8um$ zZ$T_V=GUgxXsnkmm}rTJAN9Vm>aH0aE(IbhMKpC8)Qy|0U#PhcXD~)QWidX>dXPJD zczn@~ftzTtPv3LOQ*0>KDbZrqeDZ~YxnRyZ7-an)g+?So4KOQNitL|ro|}J^DJUFq zlzC+5{sAOA|6E+&vrGrF`|Xj+x~}n z5|8ymt!`g4`^~yT6icytQb7bv_S0y4i2U|ok&uS|$P^LbFD+wY0&>MtJV}7J-<=g zUuo9GSiumbz!52$7p3#6X_X4~^yU%T9jf%pi94a&*v>*xm&;WTc~WiQM@MxbWJQE2PCTYiUZTfty?Y@;IwNjv8p@jvay znCH&Fa(?CPYU z^j64cFYMiX!ttK1TB-6trEDyb!j=b&55w5ND+Md8M6u@dU5zGhY#7Wfb0hOO1`2E2 z&`EG533wW~RO|}t!z(QMi>doCC!+VmoW?H`|B}Bh?iiM5`y>6GMjuCg`IB;$3c?>D zh=ck5#{-Q+96_H$b5$rs^2^b$YOar3Lzy&}zxVL0xi@r&FyCUAbUMwe*Vw-aD{_l; z5t#LM?zqsrUfO00`JqgV_kapLN^nr5ZX-Cb-KRNPoFnhi&a(U-t8fzbnTxTG!l*t@ z$kfSs+>PW~H-GpJkVH3^>3FsanrEUftd$O{SFK^Imc171L4OO1Cr#}G60#@wJHOQn z96NRxm03Ag&&G{8KF2T5pagJ#1U<)QXJ?GSH!g!7)8q;a>40#{x%P z>tupk+$2J~qN4wVa*YU%xH#C)w{7xp&+*yCs%~|hqu;^$> zcu<>>)Yd!`S(;hU+nsuL->rh~BR%}tYeKymy_>%^ps3WL#2>V(eTdNy3Jk=Oeu{xw>YMe@o z@TUe6&4*5>ff$z3ENY6rR7Rb;YyRichh2KBW1!KOMlgk@FeOvoG1wcyN&-2xsD}d* z!h;`BH#MX~(x6X+4dnLK#$_V>FJrdp*OrH?j;6x|Zi4xkS>k$MROR{PJ@$A+w>--8 zW|T7)>Cdcp8AdZ+ybHvHuPW*$dj4BH>!1NWh_np+ctuZoHssVw;Q;0T9iHB^*sF3V zHumey(WRp3axmxNIpY0*oy!kCrZ+3aajm>ix89|CTY$R`2OyO6i4O|yT5GjC6AyGa zWop7=9qor>>hLeEmks52#L0_?yp(S8XrWN%WM}7j(UPz86Axgq^Za3=Wl*)zZ<3Wk zLuzbMGdvfl;z5i=xmFhx@C-F?g5c{;d5Se40xnTU@A2-i>a(;)=+rS6NtOc|!_*a? zizVh7@i8$ga_gUZt;0qv+#S|LmDql-eo|e-8JU!K&Lv6RPk6gZHrn%2e>|L1t{Gqm z{_0=X22@Fx%Vj^nHkMLU6q1p7o-P=&QdEDC=6GT~YiGOGgWYs{EDAV!-IxNmADo=H z!8_3DIXV7QQ;KzV`+egI8iOI63qlSx=3DafWF{#57spAjzI0Zd*Dk)NMT3cAxPe{UA8!D0E6x8>qt_pr%eU%-LGkN)&iQks<1=^#p?zCEeF|t z6DbSvmEYYH50F_vIhSHNe%1#C`H-zCiV=e%|HMW(Eiwg$V|-z7MNA6rHZe8E103q1 zeJr_xz6c46e#HmK)hUY!uL)IeIk#?jm{o3@I)K`*1dK~ct-U>dRiAmna9esqio&w% z2WQhfq`X=nD;DoxvKxJOyySE|*6j)*{HCb*5|FW%d#O_>n3%#!b-K&6JD=6r?%K&i zby(tl%vLdc_QS|@zBphrorwhhK~_qt<>w=02hi8(5_~q7#S$b&)d~d=%2I_R;#IBy z3C*jMnT9WwJhMGcdi%B`mbn-rwFwFQp0px^S?t;z^J?)RE2pFSjbtL9cekLu{+CU` zN+~tI7;7y4X2NrqUH69ZuX+Y|ofU^Xs!}p7r|hFjB?do-6-|E715<0DIr>xXznIB? z3?In@VGT6R(wi(&gYcdVu!<#dSGp9q!|u_Dxg+0v_Fb~7y&O%s?*OdW;T>1AI&ckE zyN=(iZ8zEd{;J%mxLk`R&Y$`91qkm)Z3{-O#C6&f z-Cx(=ax(zn@2Lg6S^53JwE(87at6L(XkuyO@5|jD-82(rV)=AllBvTdE+_UhR@b3^ zg6;m)%bttT#C==e^2no{qZisk3&3Ua{7T)?+^~KxhF^;k(O2tr_=xFZ?EW z?bi`B!KH)Ybdgc6a0b87pw&?uXMN4#G^Hx+nxJ`EL&oGKbZ#UjyLdP6usBO^->bUV zlj-1$^xn3C;JyaGE~J?4RwceJvNR>^jPPDuD7l4@d_n5+;m*}6`QI}@kJs&eZe?x# zvkcD(koHaNRyqO(;@MF>kWlge0$QI@Znsh8a{_|E(d~>sf0o8DzrU;RMU0Hr;CTj0 z!e-fKHf#%lS7Oj5_r0wxh2_%1mDc$DN?;YU#l$lKxoT}ebXE&>qB@LV{ivmGa%t8t zb!r%si;SQ($NhuV`eu8ylZd0^Xn(``s4>5e)&w&N@g)u5Amp^AJzyJ!)E17ahhme9 z)R5YLXgtu0Wp7 zbf!WJ@E!rjv{A&^;yqDDR@QMVGO0t-6{lWbB&}=8D7F>`!tha)UntqPOl!I)9>CX0 z^r8IXy%kn1^c9$dp?I%`^Lcd){UoP)mwNA~U1w%p-fn{^RXbz;1nX~o#*15EB!15Y z1g88&Jn1yNIoc(zSrWp`Hj*AHC+~ zP2Ae*Q?P_kP*RpUp6G{?@QHvLDez-(*94%MIgIz%fAL*VXIl_cQ_r7|UPG6*?*A(R z$XK;Aeu=vdPCf80(e^Be>Yo<>1&A>~!H8QR?bn@@?B5}yR33Hpdi(jsxbNF=55I*J43ucK$8uS3D4D~=MMSu>FRE2c zrEG3zmZ1ar8s1No{LB3;#-tV9m>6a=LE^)6+?h4*ED3fvY~uN5f2F3#;CVYw7mjx2 zn1wB3WZ|{Sm82pIu?VVbcX~pr9BZoF`nRtTT<5^6?rV2u^Bdlj(#sNu`*)|LHbGE% zIl2G*2M<|5ymm9Kz$PMc>8quF;%0VHn`Cc)vtMmJZ##+bDj4((>Hh3R;=-!Qlvjh7l)-i{VK`d|#>q3t_|q>P3FRf;H~hEQ?T1Z^+>D94PE?lS zVObPkd?zFS1iCBQf4)$7T>I`3DER0_G~oiU6gXP+n;Hs2)a6Kv>E)s zn)pLY`=3|H{pq2&pit)jS6rZc-ddlN(=&{-@C-E!Sxrl;X;x$Si;D~I(b3WRr0dPO z9#DVYqYfP2%Ox*3-b(Ucd5bbNo(8LPva_FMtO7dkDX8HZv2^yosN(SJ9u|qkq3H51 zeS7ARj>V3MeRSUOV%Fu}{oPsF`5b>mWtGPSM;P)%17Z-V4-nive?ydXgHJ)Sq}?T~ z#?-RWV#~iyehp8d6Pqt3Ctql?@2WIEM0c*~k_ZTGW=yk1J2vue6Z6HOaaw;fe-cBx z!Rk+Bk|of@P;Y;xQY@J_4$$liCPrbazhh$^B%q|&E3O#Y$2&v6lI7f11T~YGMGjc@KE7{S zH#G10m!PWS*Lk$G|3xi31V|nRIkenR4zszbO?7Q;GUJJYHv$4FiHV639Eex8wrZ-X zs+Ebbi~N7--^xMb2hd}q2ekJCMOz;#8x_8Wr@m94lxVjG28hlEQKFyU)?POJ`ZZ6p zzu^g|^L#nreI39UTeRgb6^waF6a|o*9G~7MSrW`X6NI70UkIFuz9LayIAtm496o;) zgc)Se^_6w&`D7;T{~`fy5pd?*vso{*M@n_RZ49Siwt#kbTa(2pi;jRlo{xVTXely+oRI`EPlG+wAznY=~Ce1b`$x@|1?!r6g5#Tlb} zIBKXGNyQNB;7&^cEVT(Fe+7mY#foO@9fwvG?5RN3D>^T)J9OW7lF^rM@jCP?RU10$ zQs9f%Yp3$c!;9=G0xl=WM0lPZvAd>y#ndI2k^?GfUlAQI!nGY zonZs$)?Q<|`VV=v@oGILe<#LahK$UGF!TKhX2g=gL2-PScX*t}cQ=s7 z-pRk1$SkwjlrDBwurVEodcl^P8MsrLUaxr&wKQ-@yvvWkn1GQgoqZ+)GFZXrtPa+@ z3H&sxQ^k(>98SWscx`^YOtwDWE%?%!iAI7Y=IpK=X|_7QJVng@-`) z#Z-ybBMrAJ%6x@l{lz9%B7T?H#6%(!4npyGwp7riH24*(FMuT+T-@(po6A@YkR^f- ziHiD$A&-A3r~a4El4*`~LOZxQ{B|7FP<<(2ZK|N%W=osx;@T{X8Ouy+ zQko5uDI%XMy~L<8Fuiz{0P+#CPGYSlBLNQuEo3rWU~g%F#{rfO#G|6zsj2JnH5LD79<8fJwT|QU?Tacfj{; zY}xeiT<({K*0|MZGV?O|(;)u+n3?3-B%XVF4VdYHT`;sw@O65Y^@nV&H`NB%Z5~JOMh*bagQ*`vNv?wVkxTXzU+}t%I=3ZgF z=1Q@)N)FpyVY#>|U-*et+)x7@ZVaJp@h3C3cRY&;I}T9}3w}J$a%7Rf9Vk;TF9N$P z9kc1OI{2}jZC-Wd>i(O}V3%kxmh^gqiS3BR8FdQ|d?;jCRxRb?nx$&TnW36*`aM}uxP5!Hu>sP7p z&~YmbBJhgLgtEIs^gi&qw0fw1_#~m}k<{<182NljJFEJOC*EF~e@6fyk3AyY=qF|K zO4A^5jtxDCFqAcr{;_%7$O_9#Xvq81{cnU4W@LT2iF$1VZF6Qf-Y?Ma7F-V=X{Y8Z zu(5Y?wm}15_&N{^LS~ZLK_znU8^?y%;-bBX)mpSPeh;T#&l=%^NWQL)V-*F-Cx!lR~Pa+5Y)N_ufu zYDQ{bzoI=g0Pk<#@{=5Fr-&V;Fgv=Axh4#H+Ak`}jM##>@W@}{$HENR{Cnd7JB2ZJ z9*#PNy5t&NT+zkPdNoh%JI4%d zxe7YlK5U%`HrP7HmBD$I#@3DUr1(B@QMN85=krFBm;Sgaf%!K3#zB=f#wYV7^%)(z zslDVjVvF>I+(9zqnIpzAY?v@7kea>?5;Yn%jefx}$)2K?{+r-r{-@!Qq0`=4HyfG3 zm;cQa_%Ks^!`6cF)rGUTdPI^(brj1ZtZ?Xl=B#ahLjQW9>An3@7wSZ4r!^-qdDP=sWCKM_uIi{Z@CHW3t<7aB7jOBa5)$*9ldl!*uHlFC*X$_- zyNVZO10xHsOqr|vqge#xmpz2+_xL4Z-JJj(`^jt}s43n2KQ6#vXlLhp+mkr2)03k{ zM^4oo*L#nO^GD?qu4@QB?&d|#vi6KMvIumVW5@<*6A{_`xdo5)tr5XYKLj3oJXPrx zcY~Z)y;=Ri{m8<=F=XfVV`t~F1lyPBLQ1*1olI!3z4HCAmXffONuLnh)qr< zsc<=z!#E&r#0O>~dugB)oveHwnXG!OTduQO7c%;pnlFPwVKy4`5t!zK3VhcxDK=iN z8fAZOnt4`xex!ZZU~__O(f`{ndOXKKe^!4Tvb=O+`?cSe>q%%)A7%INQDr`?;ovr$=$j$_g z7?AT7Y9a5!wKc{;BRo8wNRq|oabp=A)#HwC0;mYxlo)b@}|K?rIVwv^c=mh_COqO5|^->$ZBgN z-E5-Mi^&Wc3f}sLw0wa0m#o`m(}L&T%-QI&8p?RHx;)>gph{iAWK?J(Dus%Zu~_9xDnrNWh9wC&5FlPdDEC>2&oBq>8hNq zTC)gIeArZ5Z1P>5?Ox$d#7wqVJp9}|Ho&<}{3d?)n}Fl=u<8wDK-`!cy!MKQwm3Ch!-gi?s9!fnKT`0dtHI~x81#Qr6Q1iWN@ zXB!zvw%ti8>XScX<5c+hnxnOF80R?~)4#$(#u!}sNV|g@z6_L7l=vxwG{%<}KQV@R z8~C`9;T!x}5J>7qfcUQ+k=<^*`;Gbxj;QD@8Ew_zJOa}^sgmq?2v(I!S><<3*M^B4 zJ<9Ig?Uc}w@J1Gs9olD?#+h+3d#tx~w%nC${1QR)N6LGllYKSs$GKM8iuMdyx`iAal~tE$L0vx)`D(0(3??oE>7()n)M-l(yKc{n(K z@B=Q7qG*69dlz17o%|3RyV^m0p$&```Df_y(8&Jvq&Gxqaah^D0EuOI*`%KA$ zUla@hVWVi5Qo4QP4@g@(X(if#*Z2aMaFg5 zPgx;aF`j&K*LN!YA<9-R1thgUD0pP;Gl=^$INd3Ge@l)KE4{iwht0h+Lmwo_3+Ezd zCo}N&HNaFba@pZ8hC`Bv5FC>>EED?48=qBJWxb67=kOYomSS;1+)G{T?1!tWYWZjxP8k@nW1)q~+53qe#_xVdO(^2y|8H2~9=0xZ zzwy4ftSnIv`$#!DtP0#!-7At$n*r0iXMGGH?^Rqd^0omb4CI#gCi&sbm3VVT6160F*6 zg$D`JbPhfr#W=Zw3HKX}V+x3gjNN`YA$>%JnOLAT-TrAFyp#+c&kO_`%D`1?TpnL- z)$LYn{hkPj^X%MtO4N+9wp*h(Pi=GPUev!ny!X2K3(@>HvUh)XEw-A(PA+(O-l&WU zilf>Zzz5Rg0s3FL)&Wf8Lg3vnBr;=;3W`8KB|FZ@hj4S8LOFXnv0^G{YEx{Z<|?BB z-%y@t5py-eNRnQm(F3Y_`bC$q!Q&3cPty+z^G4!8_^22ZVEgv%kKD3<8^jBFX&h>^ zHcytSWb`3vI;IS{BZcOL%tpzVihOL2-r72s7txR8>%P*eGmS@zvg^t0W5NpGatbi2 zzP%zqyw^LjjE~ZO;@Z7Vl1+McEuA9(w?UAnsK8xuCh1S@6I>0Vd1+7^6Gk0AK1 zlwMS$V6l*AL5b45aPj^5=813aH$iPWg)}+R1u_PQTMwoa;X3l-;sVWfdvyC?mFnO| zA8`Trm|qycRy6AR$;ySJ{HFOjG8*%p;i;W^7t{&TW5WWouP6+-+g(aUByWs|E?4d5v(2=! zOs1cJQg$`|3z{F5E)t7-6nrP;|3pGJRrWqAog{Nf^o5T&2_3!t%_C=(eS2uBeAg{6 zQ+404t&z9D-+a&FM^tCJXKKZz_U0b0dyv-+hof)X)U66sjvaTq7X(S6qIOmdK$X6`OFjhQtZ9ur4|?wZ8NkP8fbi6fK8)0XaR0oP3B&L6XKPt@ z>SnFUCA`@_8+@p4JHjEdr<3;i6A(Cf5#56D?NY??x8ED)6%CaZ4&?%aQ3enwOr%3nNuGnLN?2k7IbXMNlSN17n8uK>_^3V;gPb*+F%Qz-CgzkXX=Hc5k`<-krlR}qZ^;u=Z zBV*C~YHT4L_`=`jpKp&Cyx>5E@g9|4Z?1pCzPF{&2&p*0Dd+g{DJc7iCjvZKkS{7v zcm}scwmoD2Hzqy+Uv;arS9K1=zYD7TJj7!Vih(6X6`|d=Ia+KQ`CXVP$OppenBw=! z-_H9=a?5Fs?o@jfrgN5lJxQG2W9u)oflgi+nYF5>FN}`Lh9#c*t4eU_JeasF;4Pl> z3af=!7LKOY%^Y)RESfygY+VTX-~oGkJYm$ws$sZpcd|>%`c|U&&Zr=?N7Q&-q#_}! z9(GPdM`^&xv!9ZGpF4T1!D^IylIQn);EP%`D0v z6H%*UOx-#lgORnQl3UhOfX=oEt8VyFYj%IDuhXp|TxWc%Mli1T{hK{_&>P@g2psSD z7t(8}MD6JEP-k8E2!H;eOt4chIvz!>sLYtjxfl={1B=eH$LIO#CpNm?JnfY^&~(<0 zX#TV4YL8(2u~G17R!CwNimUPh%>(NE^?+buzp3q|yXnghROT>T%s1|B{Y&~?9M;`;}r?Ux=qCWMdy#*oi2I-a6l_?y7mVB_REC0hgIMvsR3KXCa5 z!o6Q3Yoe8Dbs15tDigps1)8p#Jbs=DGC`+Gq6v-2R#FOK`u9eOby;<$JM!Or| zgEN>2f0&j=L#)llx>!ywr_6M0cc)a=E#7|`r-X077dFHXU%YiYU+S0-CO@ikQ%rsZ ziI%&GXSn-fi* zPU2SH7H{ImBD;bntbW9)EH{cVgB{y7HZGg4O38nUPJX|R;rC&RZJUQlLR}k7A695xdupN_c%KlNvs8pMC}G``dJ|eN z`I*JfyY4f5#ap;QOLDDsgS*% z9mA$D?A+lz6_be8r><<;wJKnhcYk?CM7=w z^RCt^!us>vgA~{>KKAf`n(V;ZaPnjqG-=g zTLGeiw5)N&?w?H+)*0tnVE7Dkr6uFmAUb1nT++SUJ4|`Y6TUtT|M^`^$;PL*qPQuPNCQh&$) z#t+Nx$rbh0$Q6B}kv2y4^TF*|`U|p+WDaxsnyYe4e-2pbIekJtq;;e$suX4-ZKOcs zKG*NOSaqHqm5MB*7acBkP!D5t`~Hz`YR@4{;7{)Qdfj2H0F7Zs6o<_N3?z|02-c^p zcSRP8$hIk3Ku+-%)a%rk7^h6eF#1JC?lwtVQ*56`O-AibFqw_-6_zIu_aEJy^^aI3 zsYt)ArrbWX{_ZjkQo#-cl= zMY_8i1f;u5B$w`3Sol7>-g`gi{QmfP?mfrj+4mDO^UO2P%xh+3t#)=-;?(~Z~a=~Qux-YeGGrMhcq3pAXJZ(#n0 znB|EE%}FA{1=>xXCgHgoz0`l4l2{c!*JU6H>*CfnQ+tGvNcif7$8EpuXLdAsvW>f# zLWLVo{O7M0uSJrUiiUF`yG^}h$$R&(Dfl(XGsqU1`^WE|h+!u+Wfh%KX0oI#Y?0IA zoV;~8qU@2BhnLVQcbh;{boP#tqceN{U`|#HBWs(Eq{nYG`D+iiQh}<`tW-9i7{$*t9iLY1 z3&HIn$y7D!8)lgc4Ii&Jqv$1hY;cS%8;d0R-6Us>$4@_vJWA7-IN+FWN^L+#jeC7M zq}EVloFLmH7BDgbTvnFKK~kE+q10a`?E8!+vhyolz6k)q?IUItvSG1t!#0~|mHBYa zj2zRR3osHVbS126Xi{C?{r#y)?ak`)fZ*3pr+TnCH}|JDgX3pbM(yN9^OOxEkfRRA z%V`Q;zzr8uW*%kDz-5M1@3;)dhC&Ty(n|}5cPdtKRj#m+Lz9PqE23@5S{!RUls;Di zZ$rQZ=EjLWCBKqa#JHd&0wMk=pj9E^b#{uCwrVdXqroTs%q??fRMEcXG{H@AmMUl$ zb_h$=BVk&U;7vfdu=gRv@#K` zaU+QwkO#>uoSoqVVPoYKi*w0c3%1y!ketFHvl&NI?}3SK#bBKH7u@FyC@aUF&GX>x z&pocCARA$e^EV5DVE;26Kd>~fS~+v~Zw13EXii7kMNi9-vY@1%DcILEm+gZ7W=>v# zK_(K{O>G&*M1DPkw}D#RAhCq^Jr=$a%dUyNhp&qU>O#v-gTaW(K#|K|UsaTRQIF^s zd?mwqW$XpOQQY(wYm@-zk5nSvt8|+#E!30RE zeqRDtj;^yVL54`2%-5!k@4Enb6`oT%UZh>W7hD;mKzh7{p}(1Y=yUJ4&#$c_%AvZt z##g4Dii(?^v^UbrtSm4Qv$PxzD4k#&1X%QJDw=`4+oKyXR*+Ip=>oeWE@7WjKEzN8 z?)3ZP;6xO98pxKNvg7qh4QVq!qE^PSV#EAey29$P7Lb+k?HfFv<}T>g(P*o2$xJ)1 zr{{!RrgzL90qau{@kH)~U!?i`>u%g+ArE8Gp7E*1WT?-p6>^^oi|{k^zZuZx5P|Mr zs*8v*xG0|XT{JjyVqH_2h-f6&YQU{pQJN4_9%yc2pt+0lFe`Ophde&e{Qk8+X|&0{ zIfI@S1Kt*bo;Fdoe3l8d0k5!Ah7$2XY5uNNhN&f8IfoPBU9rNO2_&fJ=@vfFfg+{6 zbHr{1uj!2{r6hRK)GekO2%j9E3~#&cEmepx3XJC%WrfHtZ0UpHp2kr3T?BRuB&PC~ z^)`K`5Y4v~G>RcJUqxQ35FKDy-D`^s(+O(eDK9Pr4$w8<;MaEr9$HMLC1`1zjwa&4 zhG?HT=ocL+6~3)2VSp{oqz~*2BObx%0@!2Kt{rDBHX>%{GFOqKK`fB|?nP)XM~7CH z{;kg#9!($`lkvPdXc{X9%{Lh7`xy&d$nCefdbGf^Z=EENo(IGy{KI?KQ&(7zC_BxQ zH&$-FQ!$;Z6nS>=Bt3aM_T`C$UGo`%j zNgJ_KvKUbavE?NY6Ay=y11FbYxDBGr{z6p<6WBe?%nc6FxP%RDV$1f>$^Xdgi!6?v{E~RL4j?!Rp@7Z{dmU6WEVE8~a zBHiWAj-iBC=URvn45hr@5qrjq{PYWyh~PE$v3n5x)z^L?h68;{saxc+nOW|87tX*%KV{}n1FKpM)^@A% zv%Y$TZm*`sBdG4Dc`4~v(!7}$EId^UHK~7}s)KldRpJnahC>pq)s5j^`l$oM9 z+?BLTspg_1(;}WH6wJCG^)Voj_Nq4~(BnS2xw-b&h*L9wYVH2@>%E*@1gQ3z4iH#o z0BT4M-vXFqAVH(opkk!FufXnW7v0&8>K}!F!giAwKttCi46VXvTfLlr?O= z3LuW;J-}Mvhbhs^1(|C0aS5vB0gAve8smYKcI1IG(B>7lit-0N{K>%cmwmOc)EF-M@P zxE7`7S*os{9)OT$Y!)>&YNmJq#F1%F3?t7jjcpP#Q`(v(+nd-Jm;h_tzn*~Cig$uhA50o`hKfKINuxeDjar?_@ ze_q!9j*DC4V{tTnyv+o@-M=4kvyZq&xw@MQ}7; zl<~aI5QBz4TQJ|{jiOY{)x1qFxp**VY&(iJPl?bzd^`IS$h;vau=s^Qjo%8?$Q75A zjH2pG5xfpdx%mAV-xMI4HWr$it`@E@ynNl8LSiD27n5m|L2YEnOCvPTdec}__m4{b z9UoPSQ&pN=48qah^kgsCPeTy)=FG^;4x1vD!%Jt}M`nn=-_)aSQ$5XH>-x~%eSZIm z`CrRs2ouuwGglF}=P3p*$>@elNlD2DfbCJ49B+*I9`%CCuv$*DHojCafa+IpcfWix z$4tm4#X*&A+}Q`E!P(60NJvr|1!eM644QpSUEahMI`7XI?kZfal=#9nBH>e`iiA@& z8t@|-L&l1%MoO;=3az7J*`1hCsU~v7%v}y-yb6N(Ri%!!S?kuVTM;0igDN8Gx)5MEgpObC9-8!pkYiNUb8vUf=BEwA z-Zv}f5l^K0#&XkIl6$E5G)}9?X}dX1IVQ1`Z$4SzRV!f59~^D+DhN=s%Q?9vx(bW& z|GZ*duJ;f`$3K>DL7|@RW{gXqU~Z2&k1p4lc!t;8qA}U&Lmuo?H8fxk>FAZ5)oA+qH0em$*Q~)jpYG%8xr~2ZZ!4WV3 zgQ{DQkKpU-YpV|EP7;?gOEE}~3|Y6l>77!@q@b9j?Qe6b$o+msJnwj9?O|^E@ND+V z(GyMk0hcz+>$e~mDE-HC)~}%rTFeT|W~$-m0v&-En)zvUI2}P4JDTjg5^Po=P@0C6@L3BGBPEr-cG0&063t*!lxCrGU87S60m+2FRurS^9lG z<%2lB(&aCORe4M(2GQ#gZl{Djdz0G0<+x5DxD@iD#s+qHd%4uq9Vk>_$ zH6xAL28brIy+cq z^0`WP0&O9SJ+5n?C(|lZKpVIUPIiXL@(0X;Kwka*~1TY$!~9EGn@ZGKImY~R(vpx=AA ze)wtNc*RM^SVtr%H|Zi}&TT5J2Z?~etVL;y$*(CL_=BlC_4LrWS_d!+blPUy)GKId&K8~udWb2{9gJh0e~zKb8K)ej$#aAyhSzMv6Z_W-+Hb*6 z@7qBgE!935Tb9hge*o50397 z&x`eRDMvXG>3FuXCsS!+_=fIYRXpvAPEV_&DA=!Aa*#aLi&}$yi>pPQ;=Q+kzaji5 zlkm{B{-haG$k}$)!ZB5vsIfe!6D=Fx0bBE3oq0n=KiU6U12nvD!8bIz=vFCOanI03(F??DOB56}+-`Z@R zcCu4SzI~P0|9f*+;^@o>L~o0F-_6;a+y!AwCQ@L^sa07fCA+LX;s?oJ9o)3Ev?vo( z_a|mTZFx#(w0V#Of}8-7%u^jW%amp6?zT{{imwLW?M{lfRT#i?(Mks)(W)!Aua(gt zKtk%t_WGG?d1D~^bSUsRf7oscntl|m81TrcVsw)@p%(e3fP5hZ=46~F0&cSZr@eN! z@d9~n!R1HG!Xk@A$m4R(emCfe0ah;5vC1SbKb@(g`o$m%+6OuJd})d(z0rBi$_An2 z2mbRK|TN)0t`aRNrMxaAwLEJ$hFPPEc|5?K2VF(%$v)LBB%7W{)yB4GAb%6 z9W>oA?g)B>;ebfBoZJWfa)=Xxza%ix@@I(D$cL|3A(n;MV17HUUM_8`l-L}EdB4>2 z%w;YG3O6n}whKosu2pqLd&xm<}=$ zNG~11WQpW8BPXXe_sE!V$dk3n(R;`}Jk|gnlRwM)c&|F%_zF?GjXVA-yANPw8z1bt z_IYI+N%x+mk1AYyz{%$gF^LK8Pt@v9xYi1ficZC-jnoH?W z5*l1Kp^_yOs=@Xax?eBNrk^XWa5_YJWTX2-Hu>eUTms2o(gDoCyL2+&9%3RFY@RLE z{Sp@V&|SMn*5%u`HI}f~rzakoiN5*WDeXXP*^fQ8z#lc^=x< zXT(JRVr2DCH=b7^c|d(d?Rs2Mxx0Fttc;93Y>d8RS=7)=N8tM3#B|FT>|^= zGiV;@otezxxPe3$V1r0xa9v+q@SNOSyu*qL?$$;Lh@hSz_)N|jpu?h2{{vo<`8r}s zf}%*bWOHTI1?-K9^b$oS@mH%B7@y|idF^mNvkSISt&kbf(qzznJh?cqV@w)AXWt?Z zLWVza*Z?fT6d8%6;hh5Uy>Cs@u2!yyV%Sa-%?Xz*v@ktW8Ra7T@EDL(9J<5gH2BK& zEA=MNWsRPzk@4`o)jCMiOuZ;TlvTE+%!MoMzL;3fA~H6rIcI8boPvUda~=N@GhY$6 zCkA+)6sGO@i;4=?=OTf#`yI5%rz3e1YJ#T#i1;+3ZSj32nMgsgxx84i$Ei!(Xc_eJ z$wUq5(RuBUFutY+M*0_HQB58ha}OLuJV0o+uilod4l!-eLNNS;Xww&D=tYegpBHtm z4B6ClP7R?TGml%5o-PC4-l_M(^A=BBECXgTIx~Z!=rf-L0au<6P(bV3XLEit5|52x zR)-uLpRyiimv1eZ1T#h|{4EtgtVPhQI9u>|5AF_wFru*=Ms z>*S*aUl|gYpwcL9?_U{iVu~Hzb2w_CSIS$`=rAFG7hjnDJdyk8VR=1*-j1< zxIWOsAEEth8#H}++z#JNJbtph^!g0_4L^%p!*c=n*BEn)MzaIpzB_7V?)bIjSYgf1 zljxLn)UDrTw9|Uo^8OIKlkcdKj2KA1;6P~;4{N?Bz-1=TMAF69J zzwQ2V1MXC~p*&Q~%_LB9rpI9?^6)qo|4?y4e?r$+sCMhyJbY$Agr0?IJQFWXLo@nl zgD_`SquY|)CVU?Mo6h#5hJp7%tk0w7{_>d^@L4H@W`daL_9lQ_FX>aYBb}=&=H-;o z!HCR;S`{}sp#>QUfDGY3+E+uI?K79_s=Jged6b^)99?u$;)l-&XX^)kcn?2S96k{? zY*^zeAh|ZPII-mD-ufG|;JS-s$W)H+qouc^*1MKuPaH>e9@SorqKR0Dss?*r6d9D7 zuVzgZ}If)}~y^Ug15*l*0PrQoB&(V|@lN~r5o~@);AOBvyC(=h$Y`>isxKgn1aCXSp zGfa^dDbBa%Qn)C`{=*&S}n=wX8u zKwYht+&W6oU^hlN3#?wCB5n|Sbn{bxznC;)m$BkB3!U?aq;=ua#9;+_((8FytCrbC z7EFc183RnC4?0hWOnjXbHOw#8cH98gT%-}j#9f;uP4+E~ zc%%&8M7&*)bKuM)^newgo!K)(#x&c?}Uwv|q_PQN0pGYN}L&R%0CO{u0XLUQ6;PMYo-MyA#hvTRL^dMCP2rTV(`4e5QRrq-67!Jl+I@sT6W(*)=CHA- ziFbbZ8sqFhHkyaLmY-feg}Ooj`)D*dYok_nA_w4l2aj(f39uf8oV}EwPHY9Y=XoHy z;sqCPxyWh%;93(e^jfn-p4p>9b-{dk@clx3w5A;94uP8iZr>H_p-zu%Gt|B*bM*jM}OPH_ql9@!z?d2lZ2g=aq>r(|Khn z7ctD{a}rrZR*f%gX?P6|!&?JByN1wrq&~{RA71ICc)s?Afkwp79Z|j`BLV7%Zo*(W zo1yhuStopoco<$j$X$D%AOp-XNsixGnC!)47HCBCdon}RIE(W{of7#Knt1KP=Cw*S znHgvE0hA5@ywRgH;s_~U@d?v-nm^h`U61Dvyo?QNJ6~kAXTc`gyb7>yYdwx-y==5Z zo#60Lunk5UyYDO~ykvB9Wp7|ppUpRnuYm>bFeJ;5{J$dqy8@E9R?KsSB-(r4yADI& zC~CL!+oFj(wO6dG7JGahM00D{X10Sy?6L&2zW1(f-sS>juL0x*Wl6sySi6!zYgG?vbGK5H zJ%5c}_p2j~`yMuJ*l$O8Lt+)s{=@ znVDv=#7cmd-(R=YJ62#z)R&J^1XM7jDXEl>XIkDZ(CYD$TG)p7XOpWd7SN0LWPEdm zt>>9d^ zK#BKjh@bS6&U;otf%P-y8#jkKN*NF?iz$lRMI0MWL=Wg0iKCXoIyB8>P=^`MY}6T%7ZetzJt4!=POh4PDk}i6*I7yhYMSmB3AisZQ117iUlL| zc_dg%Mc6jVC>F2*14H@=@5)O(jqLCGsvHiV`wKbV=kXHFnvCB6^k%NM?%nZ)>>5ru zJc;CB9`%Va;9AmKX`h;|w=O~Q>_|*s$bOZ9*GT4pl!(Q2@+$}In zB)yH`6}t1ruFFbiUoD~bjS&!~TD=zX0}RA*g-XxxcAW&_VN>f+Q;Y~=H0J8>^6r3{ zZuz#QCi-BCL?|dUwr5%nsx~BF*TxVsvY5*Ebh8hx&G4?<{ji_be6ouIp?G_Mcy3ns zQV?gy*9S<3Ck+*LYq}Z_35gW=4rJwt_nhLrf{U$PYXW}q;wOL@8Msioo@0&?(uy@)m)c~tlm zQjUvSQvREc>S>sQ@RT8vDSK5bUh)>-D?Q1~RBJPwL&>@{(*w1_HSTaY{@wZCnfi$b z^OefoDsYp&vra0NNLl53%w;F76}-pqp_TGf6^B{M|K&UxF>)XKg*P*9Gp1FV4crEk zq`SGE^c|?Tz-=@JhaY6D46BBU@pnv4qTXN-3?CpX;vu@gZGcq6Z>VPK*IgPk>nF1{ z*N<~I?=;1VXqe0Ym3T5TFghJ~2}%`bNhLkjEkKBQX$hc~mDxePP;|Rb2d050~A9$?ff|GsEV8&@! zV$CeNEh{v#*I-$Y)bZ=^{b37Jlp6CELJAx;Zv{}u5kMi|WpOMFyK+fZSQ{gNnRKK^ zRTCoYe+`)b^Xn7vQkSNtDJTLi1s=luywV)fFAYB`h69Zx*fZ2fD6?Gr?{q98!H;}4 z6bbQvRfh+1MO|Ga`auXh-xI>MnNi$UwJs54!kD5PXMhVI{}@-|V*g)Rg00Or_(Z$=fqT3q#l+K46o;p~zSuoTjin!oo*W)Sef5qI z$%GOXCIhuq%Yv^4q~;s9b!dM_YNgaV-TqD`(PwnMLbbf!`jKdQkPP+3HZVAnq(h$d z6BES{O^)7a9mZ&iY)nxWynNMCAZ%G>(_IvD7Cl?CllhO14F|^!rr#Og!x5P4I0ZiF zpYV!o^eFT+hPpDa_Hbpm0r=O%b^Tq=zsI2ysRAMiQVfHz&BVd!(J@#`Y$&@XiHl79 z5zTZs^INg4D1yG@L2CzUJ?kof*Y>veLDo^~08q1-q`)Lrm8YhEuZ<<{TEF&ssXRG3 z8JRlx#1aqrg?9)5?nu*gj?57g-7-T?kXozYYEY^mZyF`YV_Q0|l+e@W#rlUcz@fvZ z^skQLAp^us+KSTZsDA&Zs}vBa66-9(3-K1vIm+V(&H9~d$FyHX9u;kue)M9hr^~tesu; zK)u?i0DSf^VO3&Riqnk62KlOOo>o@aE-GZ^T^lE%r1$`c&s)=>K)FS_vhiy+*eLJOpIoAf9QR0@b@1At zolByq{d@TozsaA&Z7%mjxJrF{`DSyxfOX^i{JeW$V4`e5ATL(iQ#JCrRB2AzuP}lX zrgM_FO|8>%q;+q^RK}6cpB@S-#4`7GKGX0>PcjD{ZXa3hAv}rbpGR{s2HFSWeba~$f>?+L9FCvECL6`*G)%_(tCelD)|u^gj2rhmC)I&aoL_IvIoT^3 zz)$vfxnT|AcOM^pAqjhYxAsx{+o4py7b?xSyK&I5RNn88-J|6xe}pFd`n}TD&7@%e!JS#z+qxTk)8Fc<^K5r&DJNVvS`5b zmyG3|@eAN_{QSF^)Pg~mda(S-HGl&WQcxg7R+CXK(lP)w(ekyMt}2UE3)I+K+uA@m zs|=~{%IsW-HJi6e@h&>)`de4C&Rgm?Oc#WY{Ed|;G)=y_R?oM6QKs{{oYmGo~NjZ>}Nct->EoLf_1^OW7lf@s&c^#vSJ4&g3WmPVo1-0K+jeDLE zt{WS3(=C>`^J((qzL=It3q9PnMy5s_e`2{$i>MlT@80!hpnsq(0hl5X3cuMbFT?+J z%yp?LK2eIh5%jf<@Kbtz>wJQfjD|-A+L$igyfS{35f*c{R??2u2Vb2!YnvFW(mzBt zB=hrYyOlQhFBZV>o7k;pIo**=sgokW(L%n73gaPkCnqNe=_FjgRKY-@(Og9~*l>Nl8h+t&<=^XOS`^Nx`J$tUraLXT-&)NP~r#wF&-`;V<`I7%wp7I(v z0-#9S$I(#*P+SVjV=nP)BltU553*<@D6n%IgsAiZ)afhFb;W69Y$T1$!?VO&5AV7$;nvohO?bZCAuovFsZg8v!onb9VL%P}9fS!C_^&D-yKzjR0sbo0aYdhP~h3 z=;}VzcHbBNIH29&dW_zbnSD!j4!qRcJrk{0qWn$LRrJ$a2Ci9fGT68KCX)zo^^lTa zQrJO>DirlpjwJ`8A`U<&kl$Nd?mZRU*q&WoU7V^edWv~UE1`l{#kL!-fvSG5to8dOgB);{}ocdD7xAk#Zh|)2k51gu!d-WvbmTzchRqY$z#9%8kA1u zvi=Q8;IMc~Pmi6)etrBWEG$g7OY!a-3dP^U0vG$RpuFdq;ZL1l_gqh$%82F;QM;_F4z8=%_qjW`j>J?p78tTswZuITAulhWO^zGVhB-V z!olqH$p5xO2YLAZ^KIOKmHJ>3cHORsEm0zT3O+HQ$gQfGK&jN~YAy%~hq(%aYGH{A zL%ZwMpNGptrd-m+pBT649dd0kYO5Zbx;KU^+CA_KXwc`(4`_sAb?8@>j7qXNO zDCRRPJkTC}e|x5?%Bw-A&4+Jdd$wle3ey@i%MM^40ms05EAtr|8d|$+!JoT?#`TBv z?WQZ2U2{MOJw3j-UG3cK$u-fLkfz$S^sSB9`b%5=5HDGSAwULLnjy77wfl}+o=$2iXp z(ba_894kesUW_n>e+5f!DsLl7zQ4a8w7b6fYu?Bj6xOR5H=io~3eF6g^@OI=5a`%W zn!Kr2wWwBMI%n+(8_sZfNQ?#mg0PvYBE%V}Hr2a3^Fome6J@qzXgh?QRfz~OxYf;Lbg%0D=GWk{5=y{|9W zVRL-{>S8xCJp5xq0x7sx>5d#+kJjH98^7RpJ4pve?aP-hzyJJc1-R1lh3l)%eATk8 zIX_C5<$Ivat*U}m$f?J(7{-RQn&xHTPehrJe-Ta(WS;>v~)c(Jv+UE`e#@CaXAzK)bSLw> zolP`>F7?8Je!;=ijEun>V|gPP5=2_AFly8Bd>8eJvy&Ms0As_RskB&HA4*H)cZ&gs zufCxHbc4qu=du3{%Sw69$5-z+AMLVnljG}%;B=@SBP=Wq4^R4eX2e_G64HMCT`6uL zM|#p7O`_~q$VRTPEt=qC>uI1gb3S(+D9b;*v2V!BZ~AmN8xK zSz=M14W$5ug&au@qY1Rw=?-KGfxu}O%xlm6_ryPS8II10iEgigvm$ys%?$a_>!gk z{QZ9$>!JZ0dz-zf*=Fs|MctB46~Rmt((2U4weq`J^1olzFoNxTeV-+#*JhK?c&l)R z=)JgqFV5M9i&#VNJ1^}BJGQIh4~xgYCjKLNI^;;cTmjP|uY!dv#ilsISztV!%~oE} z6LVjEYIURanvWkx9)H5^TL)V3p>yg#0zT^pIPqjOfB2OX&2K$HK;{E`b5e6+QJ!Vf z%f+QXZDD>-R+mu~&+gr$Z38sHYUvt^BjFHjn_OfS1V)GDi#pe`h4ok3p2|7q0MdwP8tLu7IMGdnLYIAFn$U?tZN^)=n-D?do}_TjbY&Ih z$WQ;{Gu$;=X99kn^&Afs%&YbJMJXGE|i=1w!AASDwoabQ`I zXk$p)*wFv``(sQNvOI!fC=oj;C7nS~LhVNGCa2NVK{$875$z;|n*&Ok_ZnewQKS!qf5mH=h8VU9&@g5dx zaU00RV{LRln&c<)w&LSIYaeTW$zKe6{fYlmpV+^&WkT{>ydcc6h8K-ek>>aPiWTj! z@K%Gouh6STpB%OM3g;g~1kZ2%OQaMek$`=ha5rIJ@S<(wriN_1FmhXT#y89`F9u0 zU>B`9o3dj-9UZ5G55A4k7jNsCH*WKRH)+mHl^?av-QZI39>F&;;+yHi$)LiDx_Cm# zznf3jeUJ8@F1!RhH5e#wo;~tGo}>FLKhxj;BjOJNkveA#w{gVUy*FOiG&jIoxZChErQjb@ z8T^yFYT`_->ZM!gv*?mo|B6%z^0>-vlaS+dsy@UP{ll3Mm6;fDPU&xS!qT-c@lsW1x(aHLyWUwhW7ENB{e?me> zffSP5&EJ+>A-q}A&BCjv!DIprzRdLvIBl{|WXto$`~Nq7Bp6?3%z@PrUBG>%8Vc>E zoCzf$-rW`TY1gH+E9GGU;zaPTO|Bq`WSimy=cqMVf@fDc)r#b$BS-Rdi=h)|+a5u1 z@^>F6-2d(}1r5}%X}DSeeIlT)t^rV!ZvS5U;TU`)ckdJ(uEB!)1Dqc7|65D~YB2>v zL)M8VCDuowjL$P@X*gsITgIM|*}dX8Cw5(X$(P9YzuO8J@Gl<3WLZhk85ut6_;WRB zC>L-0yY~%~71@BKwpAPGp<`V;1~Zg<@mId@$pRTFuT@$SaIi5tO^8Yl>goPsb=^V6 ztjxsAMVu)j+~_m*!;1YM4jT^1tjMK~+2Fn+>(un;q=i1Jvf#Nx_U?1`sL(%hz$9Er{))Xh#|&Z7WwUs?ZsH= z$r4WB*E0g%Z25*H|D5hH&9cSy=b5)|hv}3|`f8GC5ZiY?Cio0=zL;=WOiV{Ky|S3N z_)}o<<>%+CH;KJV53;egZU)bBD=RBLKEC4RvndIxCgh^}XueZ?j8moJx^!VRpQTh+ zg!G06?iX(`;&JirH)Q$h82L8f4w=pEd=OqSL7eweKP?U7nWRjOb`vM5FG~p!rUS zy#5!*9B}o`rte(rFZKYFTVCYyi27%NPC_2`(2k`~PF4b6aA;sC>YSw4OF*G4gxBck z=zAw89qsMHf;fRm-LvE3k~D1W4eV%Z>4;asTS_ca|czUf+A%+7Q3! zJn~f=oE7n;vpXkDH$Qx=qA?;1#{KC#&rJwiKGOeOeB#d#Kbg?W^9<_Uw(puG);Rx6 ze3>y6e~r-$ZetmbdR+N%#iu6f6$OuS{`rCm*u8{k@oXMsMGQ*<-U+Auth7?JcKQTA z3(!7PhH}Oh&6v9xAEp|n|6>$}yUU@yziDlc|GlEU?PcZp*=?}%ryQjDo*3S&z~Oaz z@=9_+{ZOlw{%b-w>P||Z1o-dUKs4KRS$G5d#KOS^wR~NUUpr!5Prbqx7xxILFKS&+ z@1-eR*~|;qNL|FCy3gnKk4>VNbNMJkNe&%kObD zB5u)nllticp_z`*7P!>OQjgLkQQarr4MBZfIsyZOmUt&LKO0+J|JW>@S~p`S(~D7# zzvib{vUG|-A1IsLn<c(22+2{h7fq(0(kn_OLp`4*K)I9Wy9ztwr${qAp3rH}>F8m{WSR&wTV zvA{F+?%ncjkOfkGI`6jK4yxE=W=h7jg6=9nh2WENn16sD!9W85B1*~>z5!OlUSd#~ znInO5fL*R+`O>PTrzfy(r#A5S?+*k7{)9n05W<=c-P9x_TwVD;+d>JKG;YU?lhfnWX6JbiNLC=}Nz=NA$+Q1=Swejn!oIuB3-f7}a3~~jd9c2VpxO}OV&)Xn-a6iplJM;p z(3CoDHz9@>&UV8?XKtUA)#^Vbq$ND+LJ~GTs%wf$AIQR_Ifr3iyzs@z_8m;j=9!rU zr2;h;(50ts6Hxw}fCo>3dR5250wr)efU5hQAUDT~egHJUfn+P>3|q`i#V#;41LxWg z;2Xb*fw8o)Ar6uci2>6=x*g=5u&_>vieErL_wjVZ)fSURbu*|3?sB=`X)>C94@9Q8 zczFrP$a)c{3j?4!J5)k~ei9=ZAIM8d)4(&)kmZ^kAWkO^oYQ7Spk0j$!$Z@t+*{3* zUY9LZ1o_gUwUHJ$&-PO|0?;WG$R zpCdkb^2mCHkcN!>thcl018G$aq z)HptuS42!2-+>sI&sGl=Gx6}mudf>euaXPmcmg}RWoqitFUy;@SC=lA2n6=2;epO< zwN3rWw3$!CqUSZj3HWkbfXBJ(Pdn%vTNrV+1$fH^(jf}>N3vwffIC>JP@@a?sw_S@6Zz_l~k84(!?{`nTYQvQdty+veKsOQfKX+7*QzW;KVzQw^wCQ6g? zGbS25QtrpN-{X)%|L~)UGKZTYu2?NU$dp~2aR+%BEc?v6A4;B;TTH87s% zyI42?b>`QS-9o^x2>TN1`hiHe9FHjjSKryi1*_#Om)ku{KmaZto&%5`YYW`y?ny;d z0(HH;lE`xOr*lr{J9QxcO%Ef1)5`Gb@|*Bivr~Jol zwejC_zr2(A7!_!4=b#hqTc?vNgAf=!%ringU&1eYx=U>qO@{=#RpQTUc=amcF9N4Q z?1jZ-mJ_HKpqGzdP!Q5>Ru9tBfOd>+ZhlQL5)ef}2MN%~ zGVBXc==maOds+A4u%&&a3G&UEq|#IJ_FfBku(@lxKbkc!=B`012(0ziF}59dtBNTW z@zKJeB+NZ~7Iq|0%1J-Dl21yK(y^NRah}!)RidG}TP8y&)h@!$x>rn4?xf z2U5YGc)*XI(9-$?*ZV;?ebKkD_vuGN&zhf3l&BPxpT2ZN%(+_@K~2O5LkXTB=&&33 z{1HtJ@`dMaKNt3G_Wzoaeo?nL-RzG-cz0e5*0%ZL6Y+F{plJ_SE9f`^(bDSRFNlnY zP^mC|Fd%r%r=z3u)s(oz)xptb%Cus*NV|z|g8%Al>3m)wD_sz2?gCe+pTEBj=pJw) z2zpT5_`R0m)vs(>hC}@F@jwko1fC>q*0`v* zqN0buK|NHYEePDETV~S9e9$V(IV=sT=_)Hqh{MKcjY*b{mzNg^yhIk+#f%zNxu4sD z%v7JyP<+K)g;bEa#bz@6^UCP_U}B3xvQGG zx;1X2z99`J!OhZu>^muukp#$gPoU#MZ;D_FaQBYqs|qc31b^J>AZBhn`;i*CTcA;G zjrLgEt!KsQy5`T-7P#&Q1TGjr-l!f}Wfep5a~hfjx>6IMg{^8r^;Wv0(YU;>+@bI8 z(-M>eg^BCzeDl7JjS+^|73D(anZvvNc>kP$F4&E|L$9;3DgSn@WS^Y=3LZP76<#US zLtKu_bEkicWDYit-E&XSX-%E^K4hoy!V=t(!p_cIYXd23gQ-H7pl#hrT}O8sSWOTz z7BQHw%IbqhE)H5=0G7AS&CO}OcLTZqahz6KAF+N=Gd6{SrbqwR-FpXBoqb=TS8dvE zf*4RyPz<1epooB?BqLlf0Ft8wMF}DxIT}C(C0vmt2nu3Aa+YX7P|1R15Xm`%D6sySO%O+C&~A1!q!q9e6G?4RTXvYj4fq ziRXPTM^k%h+&gdctWpn7__;n#gU_kkyfHf?#7Ovw@tpr3E`$EXHlMgbqF`COJhrAZ z;>;JR$tF{r-Er|rADjrf##^746ctgfx&$^IY;pR4^x-v!C|#508a8$yulXh&!>Ds9 z;zps47U8c z0Xwf*lp7paTG?si{dEsr7iU}2okFjdhKeF!|LwL}*qG(=Vjgz;Y8IBi@}V%__44|f z+KI)JJyv>zPDUpiQYb|&&9sh*W9)&KmjU$H0?UEn;cK&z8Ias}3JVL@gt#+oOuAh0 z%qe`V%}4SQew_K~T-f1)J@MJg4G*5S@8xN#+GhD+5YP;{|wki%%ov?)p6{l)Aj}#V*N3^$2MXk5h&36_!%$iC^|Rd#iC_1o_57Tx9Bzp<%`O}9V1&$|ESy4`2XqtCFNda#-7 zD72=&(as8>KHTLO5+aKv>@~Ol0xFxoSFslXG1MV=*|s7oUxZ3Hhv3-A{YIti z^=W2u*RPNBmu+Nay_=PF1SP{hq{~PYCK(wS3C^=N&1f|~8+vG2K4z+z?5o<~e*VST z5fm})o3k$AJ$YmzcR}~kFs#{_sf}dWc;fLcOdd0f%Mm@5RQ+WzC;Z#nChmt1o1%pm zecrS)WN%?V5i(KcdjF6HyM1bPO!VT4B_m0c08xiiO*X2=Y3TuD9dv$ERakfjd7u(` ze7SB{V&T>{px11S2MNV4^xG-Dxf3BmJ>~Km4yWN|I2z~u#M}8axq6IQ_z3)^mTAca zXh3`+9`>L~Bo^napaAOd%~EjSTnGfnM3i^WA|rd6bF%a#czHbmCOS;YB;d_JL3SE& za?=qTK1ll|sh5g;w=;6+q%bi-P%d&?wqi6i&*dtbl2OaQdM535^6LAV13}kfvh+;H zDuj`OcI=+aj}zq^9W@MmJi@U9jN|>oxzI653OW%br2h6g&r0U4#48~-TC|fO9auKl z+1ZDOha>QYuJ=^B<7KXa9*Nphq1%Ks%M)`84hF@kIk!>HpmDyPfE$8(@!MYm_pPWa zICVf&R5TGO2^HTluMPayY;6l6ww9O6V_45e{BciH)>P~{B^lsy5Hx=GH zpVgxfORpO1y09t&Q7PAf|v`|?}l=T_Kdl{D}!EGafDet=Q5&`?Gxcm>X!-MM3Z7(7koV1hPDx`^>> zrhP;u$IQ&US?F^9*L``#7~0`A=a9EBvuEe2`w>_+zow^KJ3CLsYNU2oy;_G6DK72{ z1N2gMv1jeHZNp8+Z!Cta>$Hi|x~XhF68mE1vyP1oTarlEzIt;QjnhW{kE7>0a6)2)X8P51XMTVt zJQ{MrNDCgBkZ+V9w#Ar z`5|XYL)xafo#o*tR<=FJi`4Y_)KJ&(C%axe4+;~Nf&oHl1 zeEbNuR^>+xX3Ky*RARhGKZx_gcJuJO%8wr#8;gVqqQF-ssvXu9PeKR<1+h{-`EZoo zU9+>_3SetHLIq(+aJrQsPHqFgq}y!*pPmu!8q<_Y9fxuPW6EtYENiUuvH6DH$IkZO zqAetg6e!vKKk{>P#x8%q8g=j0=;uqfhNP1!OM->BeiC0eTr*g#Ydlx*pIfMbDdQI$ z68F5lJMpj)@>Xr!xDj#f6kK&UVHO5~#{Qu6f!x`_Jyj}7=6UlX7|g&kBjqk|A)g{* zS}XsfMtq{`JoteZ-7D?|NfC<)$_bRMl=1|<<2o*)Q#*rEymbT{J<&tk{~zyMIlgzf z?wD=P9vqpUoCZhYo!onDw9*pj)%{u&agQ~t6&tzhWv`XNENWjLOE{-dm9wvJJ8|_EHr7W3lZ@*Nk2!+fDA?na3h^dAWUYS)*Dt z&rUKlHYIbkHD=peJ&l7IZwG}}`Spp%BCa+r;#QrE*oRoDnetuyL-pT>E*@&`YtJue zsR{Pq749@RrsrfaZzkJchY10%EY_U<{ z)FWOGFyXuS_!trTu^~PI@)Ujf!aFSdC8jfVsGwxBE3rgh@9T}lppD%sr3HQ&7dHb5+Ni$x^Pn;T{L_uqd(Wd8j5 zlQz`%v9R!AP|!bad3lm>#{DE~HFXemRVO**7+Zx;DZq7KA zQC=w=^x)Ntegn3%rqNcF>ccM$^^@2V>VC{BICb}Z-Xz_&SS+2we$>~or~A_%lf#aJ zbGDR#wKblW*v{!%eXS#1ZJ|i)JmuEFxaNxvC_J3N0#(JR`e5Y-Aba7;ihxm!@)$-i zW~;SYm7$sYnZrcy!TE25ObJgGx~!V~u=UqURq0i+tj3BVkAF>Qrfq2{T5I0sdd@5+ z?K#NS|5C>8ZX&NkQ86j&YiHNltWeuIuHCygb{~fI%){jGl`E5B7T)SKd^9^ddw6`D zgsiz=Z9Yn|8tcg!cTJ<-xt`CyZCP=xvrO9P0RM(hPZ~V&(9ax@JZP*v*xO4lrJ+&Je*iC_KF`L=YLEPaNloUZ z_98!UZ=D*duN1j-T<5JR@Uu>1xC#CBFr>K03w629j(GMOe`T)RA8}OoL+z@Astf7H zYBq)1$zh^0-ve{(YYHDohTivCSe0ePS}hjy-)aBiXW|sOe?dBKbM+bT%=X$R>BHF4 zeYmC(kJdgeXK_I?>!XedJ53ipPC5LkJ?Ha&D?4Oke}qiAL0<3hkWT#nqp^=qeub|IXav4m)4Bmu z1k5L0^`TZ+lXP#-+IMz#f4ZAyz?$27akln~w6uhH%H8cFKF1~kxHx!w*o%igi|swo zP*Wz*(fQ-c5wn(e#uZGzkK}~_ges8kdaURv1rso$IAeU-56lh1M8Fx3K0!f2lnE;$ zV9%q=^}cgwC4oxdLn(%rE_G&#bcIGnE}Z}T>U@dgy*E|U7cLr-UGP}r4PHKiWgDV1tme`N1&oM>z^u)yIhy7^|5Cg1%o zcGDT7&vc(XJm}r#HVP-567LOiaJ?wEk`yIo7w)ZpnWeNTugWnexu5<(3spqK3 z(JB3V2dE_Os#I)|i;IYxd8Jm`rJp!}oAARg$!~M5?fc()Z0i$mtjlnw4;59LTb7h` zp=X=RgGIlmk5+unh|3z1E75)CczNn)*&iR`j^`l=&$pA458#CBWl2+;Lg7}A6UOoD zdcYBsMjV=%Vrt)((-~*Y7N;EZovAC!y6+U~4+;!?`il$mncv%sSF^Ip;9N+CxhK3g zH}3d*-H(0)R<^F8ll4h?MjN|5y-NO;yuV}Nhoa6)_7gZvO|m1GX=e2FuLA|P&VM9i zo&(stCr{mHJ&G_zVHLIHB)(!`5Ywi&kHDc@wm46LAN=peL7aIFZhJpGjyO@e0^8@4 z_m;27Jauoa_p}xIyui_*linXuxD+%s<;XM}*td1`7nBFYvwrdx6&;!GxRqd~GNNtQ zcA|7(&j0AUqoG#U1x`1qj!B$KW53t&_Q>-8ybArK+EIt^N~ZzL|G=?wrmmcvoM~jN zZhYaBt_+o5pX>VDP6hC(zek-Q1BR(J8zjy7=+k=T5l5$8S4>}@`(@!6*!stsH@BIK z>uW&4y05b%Uo)L2_ou(%3kb_S(Xlb@=D9U_>0#j(*3;LX#I60f(Dzm1D@w4A7 zY`EX;;v}n(3%+`=i-yoTWYpAx>MWu{Hy^mfh|ylFq@FY!xEJ(Ppk>9^Z)GitcEs!p zd-z4XoBM1By*Kb{{OIfR0yYHyWd#0cZuINEgZvBYPk;-pMZp-KJxULU{l{E`=|pf)A9wLdPo)CQKa10~%~?l~0#==;!St7S^U!unC$EsX3WoGP zyC@B{;fCG`x@jh57=eQX^p>yNbsD>ZQxi7xo=x8$a@`52EzF%7|KVrQcRA7gUm^Qn zlY6gCbH1^2e{-va1b30FtWD7|4oV$uXF0G~tQ~8X)+SHe|Gb3v<`*Vz+O@iGECOL9 zY};GK168YFabZew*V)h_5YSiZQ(pbcbMoq0TZ$)(-}aI_KRVj4u(I_HNq?G({Lwr8 zK1^s2gH6c|*OQhKm*~t@Mva04@95o&_uBVuDD14Sqx|-VOKlPU%s7L$5vX(Wr!Lra zveG~R5@7(ZPr!NRI-VUkigQo8FyM(%hQdib76@Mh4qE{TA}LQ*W7a%w%)QizO7&4 z+O;n&|CSZJOA3^c9Nnqw;J821QhDu``@)U~rBxTK%ko|hoX_HJS;~|$!Q+@6nP5Al zL;0JLHUCApi)Xd*+`;PrNSh8_+Xl9uPc4C>tE+q9ne3*Yt*!4~Dn#00pz6Pjj6%iO zio?Snn8&bM303Ro>zhxG?KKK_zNDz($F_8(-F`DZt2<2^?^=h1@9SEuUA5T5z#-zs zp>ylJMYuBg+4+V{}XPlr`WO}-aKGU~+GB9!{&NGF!k;d3WVi3I|MGuEeLRSLI6M0`7`V71Ms**o z9~Rcu&#^ukHw)dZtyQDoIO+s7KIo7g4rWwW8|Zi!z;=)l%464tckkXcx-QN`!+DsJ zl5*2oep&$Nes)&OvR&p$AF0H34;qz*YK=j1E!0nxk5S>L;4hgym^scOqJxy!Eg;}a zDhu#PEYprf%tl5=U=>#y4hh>2Hj17HJ5JPDmc)Q}{5@JR2U<8F3`{1+Fi2)yvK@%^ zdvf-1J4^kpPnE_UeRgpU>)sEp9D1!)EB@GZ+9C0!hJd>LoI#{XSRA97g#QKR`DjJ* zK4pEz$FDEU)Z`_Z)(L+21)9zd+;Ux})o%2ubV7~0KGYocfbHPs?kf|iTVnGDej&&MdfFzl;gw0|Vojr?>9kAY?f zo}X|`CUdCVNOuZ=KG|d{Z`2fehdd!P>_$rJKxsCa2$v0h{LXXz9##pbLr`P%?sdtY zIg{!-@Z*R4PmpSxnwD`ZGJ2%5dOwbAVCo_n3x)=r*X%O0A=|~2>%x0pxV9$;EtaQh z2rN!azSSFvN7qpCkYG6TVW)_n;l%msf7m>CnpDJ*I8K32H#ZbtpNOG;80^ZN@6U5y zEyqC^FKag>+KoDPH*mR0CN+DBzzrgP899!;HqX`dsi2T~Gdl(K3n}jLOP-x2+D5JK z-Lo;94}v%Z1n`xiIsY7lnu&;RKVBhC3W?lO)wSW;D9MquvoVd*2@wW-z%3jWeE>qk zR=Xd|C|F$|Ko@BI6)r!L=`g{Iwxv!OlgKASJ@z9aBQ&hpof#a?<7@E3_EukQ_XiBL0gPkB-_WpZDiG=$i4_#g2};`#HN1P%Z3oqst3W%(Z_rqr_ZE8jd2e8Ii3>Feu{T9-oX#_(o2H8r)kkum*+KUda>0K4mxISsV(SCWxSH;mIU%HZ&a~rL zue(!TbHO|JWvwe%%k3zGX9o*{l(2=OLa7qW>v}`No*6#4yTV5{`V=@iA(n)=e^;%F z#<|{3Aapo6Im^)ob$nT02y8Bu0Au82b+k>Xy<`MGrVpke8)J^6Uc9gx`FRqZEk-TU zWgkA=BDaYH=1H!)b4U6{d?HS~N3vxy99kcg{LA(X0*Bq~K)U%a`$idPy$5+0I)p%U z5Npdh$@!01D{^MbR$l-}@%Y&yOj@-NT%yi}O8^q?w;j~dvF~`KUuN}dSg+Zxjfo;A zE^a@NT61-T-jDX8;ob>P1=z8#GY^#h@@`XaF%n-KnAFuVURWD@Q?vDQN!CJhrj1d6 zweF`Mmu@M~k(%F6u=SFUce6#|7CH@Ecnpqp*?|J8%d*|K;%~-$%wrV@#IfioUG4Bt z{$*g;qx9h}2*ni_QZ=AW>ptn*>N+=~DY7Or`hQ9Ai zx`AxbG1~|O%!&y7{qm2gYajyfyJj3OEIX?(wX9 zA~ZA<-+lpzzP~9eT`d^^wjA%?@!dH3ZX#frC~!FA_a1JskDOS0*BN_c3Vr)EYs^G6 zb6t;flSHV9qq4f@+dYK~7j93h2O{Ni2IEsN!wvw3F3~3zZz@6vDnw|CUAv!jkcgQ^qV~CV(yX8;`|=lQUQj&IIF~&2F-E4b%0ejbe=K!k z{{(G`t#w)=3U!k&m9U~z!!}ilQDK7Q;COTZpol~dMA8Fo%}Q6MGp|5aBHc-p>w}HK z(_>w1o~b*YKY;$Z^ayg*Vxtd!^sS6V2@c~+bu&Ash=@oFos^`c@RHd2tl9 z(i&pq_#siz^f57>6>>bwzivA$@^w~L6a$?;lYNC$wR{rLM2@n(n5k+o+FX0 z7%o4~+@U^~TxZ%3b)_giEP2h>5{dYrpdhEPef<0b(Pv~)rvHL~LHUQ7xyG;qPShLE;RI-I6h6>R-!#>02zZwN`#ra zO}r4QFRu7i%ECPP8}M&Lx7p{s#rYvkiL-}+a!?1akU~{KpR%Zn=DdIC>6c%?kq8G( zM~Oh!Y~K~OoO7QJ{0$U7&9r_GpJ`ZM_uCEs(9g87v&)!QXav0kjV{Eb1c(D@TOqe> zdZds8d{nen8xZ%;N2)r@lGvifK;62Z4r>66-y~wn>f7|d9;vus@@}GxMfuj>D`1D> zUQl`m?1nEuGAg*mEFl8~{Yp1a!o|e}{j1l8e;h@IG8}3&PE`B#|4Fs6IS#_iC}e59 zYi(8Ftl7vz57m)Gv%4~xtJki5k8qV?*+I{+fiI)S;JJi^ z#ObiZk5J`o_;uJve3BKMOefrd?EqmY6)OdX1_w*gtM?1z^+{tO+HTQNtKiV{Jdl`G zVAkT?wadfkiTx*mq0Y`v#Dy~2@zv`hM?OTLWOR!oLdwGb)s$s>uvq%!$^8gFP6Ksa z_*9g5J&lmyK182k$GzQ8xW`jTYC|O_ z1<~$bdiwVT1@0i0M#w{)^Ge+#C0>7!G3Rp}thddI%=`MTez=aPZ!+rFuyYU%u3f*LGGww# zaWfo+ZEbDGEiDBhVs&CVz3>sm;O_vSouy%Nu%nAww%@I@IDI~7=mw-dDU{Dkum3Jf z=Kny(`~56;(6O+wJ??4*Gtlo?t{kToJUf(^7o~siZliCA+wVWW%R;(`{}<#g8B$gN z6}cgPs;8N50=~{_&8{}?$eC(x!=~QXYIt%FRu77#<#~I#x$nA8O4R?y3lI?zp)aF! zEpKt*2uNu84Fz4nC(Hy*@LCf^7 zRXB1nxPun*LM+qJxM9s<730U6&WcvsC)yoyjbB=?KCV?>oZ9$Kd!41U)X^U^^Y0!^ zA9cIQc;JxOU^-iRLxv^?HD}+!#6;*m+)s66c+XYQCnJAC_o2Z z&&v7|;~ZW>j85}bCmh_j-~&eV^VzTmyWhU8CQ?Z{rcxH&cpzV8aHMn z@~X!1q25LLbQa}zHJIBtJfPacyL|WU)?6C7M9Eh~WBdfd-W3!hnUhT!7ROOif$6)G zVqDFO5gR(v-s}qSX_)!4FwO(Tg+?r?A<-KaF*h_>AV^}k05b|H)fR(TBZY-#9 ztN*!?{$rRY@YZJFH40a4S{2kVrmbRXD8&L_ym$d8)b;kFHD+~5C&?i@s8~kD&O*#z zQQfwt04V)9!kj7~pK_jyvvRD4vW|`p4wd+ZkSm1Th^1qK{Ag)wTZ5X?bjzlq8_Tzm z_n5!e3DCXsXA3xbVe{r}E%8L4UkfG<2`xBRG~Ql+1^0p3n)S><(5XLv<{4TmD=+VR z%A5Tas)Albks0u-Fi9Mmo|(}QEYI=1#+bk^`Xn$UV~I3o6tqtGv8iYZfQ7a6bGTg}axd z!5tov-s~WUp#ymGss*0w*R!xhg}cmr{u%!?0(DMFL%MlOn)<43`xuw~Or%{B&1%XW`YRx6i>^UdGzS;XCQ$D5%wh&!0MV z%6kVN!`FzYrz^LSJBTP+n22}f;p00Euvs%R(U)pewv}J|trF8_Az#jtP|@>1Z*Xh) z;`%2bj_eT=({h~buK-NMvK6LQe!REN8r1ph7s%amRGJs}?RHCJKyw6wIWgx%sL47nS+?Y{qNLZ0gLlipf^6im*dn#UK_=Q?MB z`hWJ~iGwUH#k@r#2es~oUrg)6!SjeO&W*jbX!8I{+(Nb(^jR->s9vE+l%?g)jmqKi zN7g2zG$#9*96oT^cF=Thu}x3qNo>WvlVdB~^(i9PfY^`IwtqfkI~W?Staiaq4xe(_jEMB40yB;P&Ino*+!(bu5kJQ+#+6?PoffFOuBRxHR1YTfPtFBUE@;ewsvcCHen7>YgxoVg*9w0}s(eBf) z=VgV-Dzl#bh)wbE$jE+AkXen-gUJs*9avj^)*Bik-f$8kVdd~Wyy0(K}I z3(*j&RgkKuti2&9a!W|)D0XVyK#u%cy=&VxDLfiFiD+RLTgWMSCLRz>V&H723F;b9 zMVx+JBktzRt???cA@<$oI6MuujUzD^NvMRcG7@>$%I9~e;P?nPPx%QRu5YDvWM6aFsluG#G}Ckz)`CF)_yG7UKgJa z(&Z)A4DX>sYCysQ%&z<|K}u3-X!1ib8*ZLO##sSgmF}0m;PPfZ>j2|%oU{dNd{N@C zG;S&tmGc%%fO%4w@ErHjl@*WlH99%@9k;iB&6=n1A2s~!g2$zW>*A_q-abC6@VJvv zR`?H^C_}!1)i>UIYy}q`-O?>d`O*FLS9{jJS(jx4#ud-p(M3_>ept?=GK%W$YTNhy4D2Ny zhpl_bA$EO8t}}pa9o^k(d|fLCS04Y3MNW1O4n7%@KpVob)UQPe#lD?_7r(rO=nIwY zKH2ElA1K|1jPXM`tgbt#?!0$3 z9)C7WvcBzvFF$kj>Qx$16$#>oF+pS%1zR~<-&gzSOtSTJ?v_*NPc9m3PWf-N% zaU`@TqyhM4o=KJQ>*i)h6`vFxe7?=>O)I3lN1Wic&cQa_xuaTly|hnezc}mA_vGXV zuOzKZ%afpLZ3pOrQf1^+{Ax5n4J+q5*#k_2~$*xeAS`b9x;ByGgaximCUjFOY5L zrp;YkUX29&IDH-+Npsn-gDWatGEak{Pol*JE(4zZ$FF>=tmJIU32Z|0wFZOf6`lq9 zt8<_~Z!8Z;=E@D9X;ds7S?1V8oyQFHHrOWKMz#-8RlW+nPq~b#k>U4_aMLD);e_kr zlii?|+?yO#-&-0B0YhL^?9WA)A- zakts1Hx%U&DkH}l5!A)sVUj3W+$wOv?PMvF8vUOwVABdXvU$J$?OmsXB2WmQLZOgB zHZSVnKA3+GJB+_XC)O*ZZJ96ZGe=MWA7ff(V{Kgx@nfP=ZeI_6(E-Gn5M*7Pu4F)q z!t1zptz zm%xSHyLS`x2qEA*qEmIIl@US3sm9fcPEOh6YG8R4pXrYD^!Co}8o^HF2E}GQP_GP+ zoC0{8_MbnW0n+uLSFHk6*LG@%&=p#yyrvSK(>Kh1drn?G0lVPo+WY!qQc{G%M6FJ> zwziVYaJ40upXz`oO}^MGlr4dU-i<6EN*1Zg2kMfgp%!2<^l0CC6abu$52-|liDviM zcwfyq1n#Q5g^7LI*cd2?%Mx3fnzTVQQ~@-ztXsEJRkh|=_y-!&;Byke0bM}7*7$Be zZrN{ZW%UhY1x6Ik=m)PH2IV7kX>g(R2ZW1$9ALN$m*{Xx#4 z%&Z6Jp=h2n5l}EHr|97D7L)jZ$ex(9!7QFL@mB{ zzQx4nzQ)6hfkc1|881TK2!}SdWg)b!~XHaz*e~W&QJJY)dYH-nq50UwWCQWZsCkLBifMQ+)L|R%~ zKZ=i!C+D%5)g~Op9VA>{*l+OZGG<$CmT=0z=G_gN6kqunB?B&1;lhOo=<+Xt33~9r za{A~C++u_7(**bM$UYkDDvJSNpb2L6)Vx&q@aZ0&&+s<%Cn+GLfJ!MI3t7fh0Z^?J zrzJ|l8Vu<-upbTkWSCPhsILcn=ELyt@TG$HkcgG!&<&>rj}y#NOV63_DmV=Q3>>tE zB`I=tMMO&Ef6C`3Gv!)zJ3j7;h|$VSdTTp$HCBz{yL-d3rAj%%7qImu+9O^<5FoWB zdcwx>uy_Rp6(LlST+KyAKXYbi=&x?3{4~3bOJ!4pFLKrT4I82p5)#CwjvP9KenW{$ z5UHjeh5^=qlFy@~H_X@rI4M4tedM`ap%s)Xf(q5$MF++`Wk@mlGJmOY{ElN3dV<3x z7k)*OVu6Dr9!t3h4;2YP6%Q3Qo~L*ewTCbHk*tG&&oe#-!Fm1KwffUiQ)Ggkr>7_H z_!xQ>USijNUAf%U?RVKoaj!}?ItK>f+xPvl#Ki7tA+<*Scj-p4KkPh{61$V74k^!U z%d!MfmEXk;h0@RZ=fWuY-hT?c|LGqN=Yxgv=LE9Vy;kh=h+p%& z$GzcgmPyu2i-nsTs9qV|8zgyOJyskLs-&;~44-1}B1K8qSl2cGRT2HL9+b%P-vz&2 zgq?#!C--vMU=A8mmuV;?#ZxuYvI9}gq>t5(@}G+t86Q`Yl5#^0^vKsYqN=J&E60$M zzZLDMUl2z?ZD32KTC%UK9+bpRhtf@+?na>3%3`Ti{p>0si$FyKTpbf8;L)Ew~8)*Z985*Jw3hP#zn}Y(dZ^s0IKo4f4>Tx0l5u` z*+)_Rs-h!is80;h8_If)$zvd4%mX+@BjrHn?;+F_J1GSU-%)ay;}%I5E!sG41Ut&2 z(;3AMS;qt@qK130Ra^u8x4H3SfBp4WV{!>*cDD8Nt zoBmoYv#rkK3QU^>PtaJv2kE-br>Z64^N&l;4d3GY+|r^8v4t<-&xb_(!o_R{ zV6|@h#6+z`(|pGd-znD+#FF52tSN0ib?~cCacgqVrk2j_qvYT6UVh=99D55s~6zS>Pa z42reDD)g}R7cDiA*O8F+`CXPKAyX)7%wBH(tCv8!y`5cAL&I#G8F-&pEiSWX(4(Lz zE9;5zpnL%a70*rzG9Fi~@L0T*^JESV_R?%;;K*5trSlsoMY@w^jj%3>oL;W^jG}f7 zhV1h+CG>JRrcDeF3z+ZTwM!#~l3_QZPo+J6{20oxV@u&HGz#88l5UvHEnBg46ZFI< z@LD@7E!H6aZ`S2n#|-)3wqkN;Sy~-~3z|2R)I5+J!wy*-qY=)--riolzxL1VGdWN; z&^{aTj_%*ES1ZUsSeHr?2PfiR^Vl?~sZKL!W;?n*(^7xk@y}dqpqjq!6UrB8aQZ`!5Mc#eOTnmJ6^k-}p#5sOQ zKXSAvK*48bnVuxu>2A7pUO|v!X|CHJ}S7=nA0Yf;Ubm-79r{tvf(lmVvr2#KP2FrmjpiRWMN(S1O zD;E{1TNz^7s3}?xFFkoN^z_cwq}pkI(YcMkA0^jHvM8{Ih0LLXbTnLriqN?|A9mD) zqRf5TZwVnqrB%UC(YZvuQ78#JStxRM0apm7pp64Ge-^ke*5D4pD}8pn$H31Df2IVz z&hNyKwf~3>I&{zNh9wbkUuegWe zslBr^8n{!lJ(QqhXh@vYbK(>JaC#Z~wJoJmN`zqqMeQa+tTDh$Z$;!C7&wGrHz)}^ za!{HLM(la8unLGj{FbHtUan|g&fih^`w2*W))O3kY`RM1AzlxhMjA+*O1KUfM3m9d zMB*VB3|e@2i*~rVU@#{YvL2G&_1po7CGA~ZF91uvHZ*8jSXh9qPi&X_&7QkG_`f7Q zFSp-Mj$(8FPfn#=DEz|&gJ<~P|KXcK#P-ccY~CQ$bE8FS`q!_Qi2pPF$z}ebp{9Ly z2{(cK`dStj(g}xxj!uFNz%42x*^m>!;q0k9yhHZxzw34t843owVk)v}pG?9H;1@z( zfSyr7CDYN{`x;VKFC3%a2zLi*Q&al(dN}Kip_unVV}7j?YrN1<5pi?{;j>YRrg{Iv z@C&W>si5Z0gdS1=k&wcx^zy<2o51Drh+`Z>K$J@163)d_Y{JG>^7yB5w1;X)i1C%K zt%3NIv~%GWaI@Mv=dlVK`^OqQ3ks8PuD{%vF(STbLj-GZq8_nr&;A*l zY_Wyug-+ru;QW5WBHyg7=3q8YIXF0sQeHrJ+OE_y0s$xz+^v0iXhcjLImNY z0x+MyeEB8v>Qj)(1S8Om88H`1Q;R* z_qHI#X8SP?x!Z#%NVfkFPxQM2s*kRd6IvzyWX$Qw;(^pFtzR zVd1V09uiWXqYJ}C4*7;k7b}z(MrJSeb=-cY+BZ)0_*IX6#o*|>-rl9Dizq3CUE9{O zvzyMlC3u4FCaw3UG&D4x%p%Ge0!swNyJPUpp4}&?0#H|v(<;!T)rO$>9Ak9hf$)oz zzGiG}oNLy@Z`y(a_e$;t@|e0>e6YRPurl{^ZURXAgd)V$mnA^J;^zritBYaAIl%Zq z01Sn4_f`q;atEfinWI$lL)Bv7ipO37lrZMZ6M0TH?wJlV0L4M%Src|xQ(8oj)b|5f zZvXsJGktH)T{!`{Y{@xgyuJ$Pl*ef!Sj7#)Xqa2Mt$%}V@CRP<>akyCX&#K5o`D(Y zGCp9_uO>cTew;?mmhHwHR4MtVPYuL2zGRQIbfMUk>Pz{fnAYt+-GWjzNkzziHA@kNUWtREA zzH(JxD^^`KSc_^HkDVIEd#mMOiF(7?gUEx2-r^j_P*ApZh}l)x;1eUM_MtkIV*oy% zjaCHbGs-_PFz7w7cSL1lpN0u=!%O20wba zkAW=MeCQfMm}FDjz=}qXqRbbqe?DRTe7bi*j<8i`u%0?KRL5!_EbEH6j;aJ4<0DQe zTk=@XUz1Te4pbrdQ7ci}US##C7JSy;+MIkH!pnWswcPjig()1Cr`%RHjP31hZGDPS z3>^)7XYb4d-fnc<>wlV{oq;(A`u7sCYKh0txA)M;$A|PxYZl&k01vAo4bizc=PhRs z0HWu;dljhf%+@oB_$o*FMG*&W!`{B*S)UCRkq5oRH_Ap6j=zg7wuJJ3nI$VL+Bf1R3R*^6&6q_8hSt^yjRO)7C_r*$kQ`t5dF=hepT?J|2VOjWJ6L zcQD4zQGX#EwCir9V`}ga6Y9Cl`u>LWbHGG_oPY*$os- zPAo5&8$uAR7hrmGW$xDq1k&?pP8TvPT1_4$yl&m-)1{pVu#3a-R(jMEN0kJ-{@Q{l z4`MzZ+HL1$Th!p7e$HYulNuh%sEXp#V|4sF*=lIqCK?k?SBQkHb_LzvFh7)ISG>d? z=MffG20_+-ON+-AW0l4xa|VxK6JKfW%uOCAa`ZIAZD3DD3LhOL+h{mgty4=abE!eON)ZfelqU(Cp}5{US6X@ z!wRJ$RR<7WnX|SCSMBfTx1BlzVL%v~9}$Vco07%yC^Ym9BHUQE4w5q=&&)P$QW;Lq zn?FO8Iq>g9a$sUVBu@%1d>06U9bknBzk~+o)b93SyJ42j+yi}2&i~9!P5~c5{-{LU z3)F)F+wPmf=xR9&4Z|3N@KFt3>O16CE&u0@*=KEA7xOkFLxu)!NZB-4Uhw3}lRaU! z&3>GOrfP8-d+=c0yUjpe6BwEGS%8t0hi5{B5^+!Rh?4dQ@(n~?1=~j65@N$w5K)8) zC5A`VclGufH}5Ta?W8^?3UIn!bC%_4SnN#{#!1$EYVW{WK7*(OETlodYE_ywd#&1J-X@yR_3mnecHO#uZYHblf=%O1DFth8G(ux1#6@U@;*8;o`aPC6A;HdxjogO(0jz)h zNB%Ol7ap->1v&Z{;Q0rNT6K0q)qd!o$y$sBld+PAP6A5_;N>)tV6dkylQqlenmuKBxeD`4h7w39nr!J%{A>D8G(EhSuE8DSw}v@}3t^ zort5nA|fhCnVnnNTpy*5MslkB;H5(Ve1j)UoEb=^ZcHd@KE1kdM*TxP1}41%Io*^z zzeJ3Hh*Gel`?5$3;d$;KYgpWj44U zp(#o~S_SspsLv$w0bmY508zUart%hFL5oTArrVQMBmP2<6)gTke^D{9zQi{Q{}QTc z+nZ;o+&Yin(;#8U@;G2*)I4MHMf8NArf?gY(O{ zsfZ5IP>3}kp`%(f`%=AY@7}Y7HRt>chY6Bkn(-;MJ3!kvZ{Fm*Si0%4b4OP{o^n8& zR}|`S4K5B@lD@toc+ke-_9gK_bg~81I~rZ0ql%kJ7)*0GrB()a6!F90fCZ3ML~#&} zW~8`P-W*wz5CTy85&}Bm-bplq!{r=_7{r8v&3Xlck!XUur?;ZIC8ZA0Za^W@$SyI# zKp4>t_P>a02g0`6o~>?c^+9XaCmP&|)6I*+H5%ib8bSy(Ha4c9nFiHFo!GlYYzx-q zriDam0JXy_Dyoi7?D}LQAJE_2a11r&I&0$9RH34<#yJ<%;F#krW)qkFO~HAc2RSg^ zqW)oyiwBz$8(*{@BT2Yx_M>(p?kGTylUcSyDgL^(D0uVNfa0 zWnSN&pZhJqKK`3LSQ={pD;ycm;T#3ad(%X7GD!c|7#LXb;loe9kp+tXrCmEs`WY|K zfiQAF7an+3->wC?JY#k9bg*cWEd|x&^zPJR21*aEx?tVk%(qb!Wg3D@s1nAydGsiR#k@GfDEOBVhVNH_sVlp%Mf%-W!z5+>#n}4Po3coun z_#9Ga=;-Kx4pfHtAq0;#fc~AIx;zy)%7chvh&Rfhg&{4zgp^eTG(aK)WU^GKW{{=i z=Waz6x3-b(0q6;Aun*4XFL7G_ib$(aeSJ8opeX2ZjL@(Z zDtOc8W^AzMk@&^@n}2?hFnrJ}xTp}&`KqiiML@4v!s6%clUC%nl79czF(Si2{pO+q z6cA~cVPX!Br=awGI@y;iOmUS7d`oZUxje!&xExMueC~{~pSN{a=@z*YU;D?EWrzxTk>44F_ z_;bDjU6kp07pV>zgp}m#Azy&NAcPpq^-Ftt<3^)bfw)_@v2x zBF?G94GK_73P8qCypb$UVs10&DU&DX(;1WE-A6S zsHEfCPy_|;)4jTPfl5nGP6o=1^!I1+$|hO}%Y9S|AHHwLKVA?=2qVT%e;REEq8sed zguf+`9K--tcpU07;v0e>f`Cu*0yK>1d+T=df^LFnMG65xAOE&mEaBI%PGT}!49s|3 zHJK*Wues2A3{F6;uY=rx!}!wkxCqni1j3^%Ei=C+%}nzFqz0@%%W) z7Kk?#48th&JkAm1POYo|p^Ql%-0u zjNOa_P~CnO?7&A;o>yUI3Zoivno>pGg+;ZeYCm2byaZp6qjH^TzeVvq;FR#|e3vj& z0FfDDkguB2?cp~kZKN_-Z;{jUUfbr81g_tMR^V>bXVFXqeo}D{) zmQ63BsC@;A(PI&huh#9q-WtAw!os|?D#$` class, see section :ref:`blockstructure`. We first initialise the SumK class:: @@ -90,9 +93,9 @@ We can look at the diagonalisation matrix, it is:: [ 0. +0.j 0. +0.j 0. +0.j 0. +0.j 1. +0.j ]] >>> -This transformation is already stored in the SK.block_structure class. The next step is actually not needed for a DMFT calculation, but lets see what the transformation does to the local Hamiltonian. We can calculate it before rotation, rotate it, and look at the 2x2 block with off-diagonals:: +This transformation is already stored in the SK.block_structure class. The next step is actually not needed for a DMFT calculation, but it is good to do this check to see what the transformation does to the local Hamiltonian. We can calculate it before rotation, rotate it, and look at the 2x2 block with and without off-diagonals:: - eal = SK.eff_atomnic_levels() + eal = SK.eff_atomic_levels() eal2 = SK.block_structure.convert_matrix(eal[0],space_from='sumk', space_to='solver') print eal[0]['up'][1:3,1:3] # prints the 2x2 block with offiagonals @@ -104,6 +107,106 @@ This transformation is already stored in the SK.block_structure class. The next [0. -0.j 5.028+0.j]] So the local Hamiltonian has been diagonalised. From the other entries we can see that the *up_0* block and the [1,1] entry of the *up_1* block correspond to :math:`e_g`-like orbitals, and the others are the -:math:`t_{2g}` orbitals that we want to keep. +:math:`t_{2g}` orbitals that we want to keep. So we pick the according orbitals in the block structure:: -[TO BE CONTINUED...] + SK.block_structure.pick_gf_struct_solver([{'up_1': [0],'up_2': [0],'up_3': [0],'down_1': [0],'down_2': [0],'down_3': [0]}]) + +We can now look at the final result:: + + print SK.block_structure.convert_matrix(eal[0],space_from='sumk',space_to='solver') + {'up_2': array([[0.156-0.j]]), 'up_3': array([[0.156-0.j]]), 'up_1': array([[0.247-0.j]]), 'down_3': array([[0.156-0.j]]), 'down_2': array([[0.156-0.j]]), 'down_1': array([[0.247-0.j]])} + +We see that we arrived at a structure with 3 orbitals per spin only, and blocks of size 1x1. + +The interaction Hamiltonian +--------------------------- + +We now set up the interaction Hamiltonian. Since we want to rotate the interaction matrix into the local basis, we are using the Slater convention for it:: + + from pytriqs.operators.util import * + from pytriqs.operators.util.U_matrix import * + + U = 2.0 + J = 0.2 + U_mat = U_matrix(l=2,U_int=U,J_hund=J,basis='other', T=SK.T[0].conjugate()) + +In the last line we use the Wien2k convention to write the U matrix in the cubic harmonics. Next, we want to rotate it, and pick out only the relevant :math:`t_{2g}` orbitals. At this step we need to give the indices of the wanted orbitals in the *sumk* basis:: + + U_trans = transform_U_matrix(U_mat, SK.block_structure.transformation[0]['up'].conjugate()) + indices_to_pick_sumk = [1,3,4] # these are the relevant indices in the sumk basis + U_red = subarray(U_trans,len(U_trans.shape)*[indices_to_pick_sumk]) + h_int = h_int_slater(['up','down'], indices_to_pick_sumk, U_red, map_operator_structure=SK.block_structure.sumk_to_solver[0], off_diag=False, complex=False) + +Now we have the interaction Hamiltonian for the solver, which we set up next:: + + from triqs_cthyb import * + import pytriqs.utility.mpi as mpi + + beta = 40.0 + S = Solver(beta=beta, gf_struct=SK.block_structure.gf_struct_solver_list[0]) + + # Solver parameters: + p = {} + # solver + p["random_seed"] = 123 * mpi.rank + 567 + p["length_cycle"] = 200 + p["n_warmup_cycles"] = 100000 + p["n_cycles"] = 1000000 + # tail fit + p["perform_tail_fit"] = True + p["fit_max_moment"] = 4 + p["fit_min_n"] = 40 + p["fit_max_n"] = 100 + + +The DMFT loop with automatic basis rotations +-------------------------------------------- + +After these initialisation steps, the formal DMFT cycle is very similar to a calculation without these basis rotations, since these rotations are done automatically, once the :class:`BlockStructure` property *transformation* is set, see :ref:`basisrotation`. + +The DMFT loop itself looks very much the same as in :ref:`SrVO3`:: + + # double counting correction: + dc_type = 0 # FLL + # DMFT loops: + n_loops = 1 + + for iteration_number in range(1,n_loops+1): + + mpi.report("Iteration = %s"%iteration_number) + + SK.symm_deg_gf(S.Sigma_iw) # symmetrizing Sigma + SK.set_Sigma([ S.Sigma_iw ]) # put Sigma into the SumK class + chemical_potential = SK.calc_mu( precision = 0.01 ) # find the chemical potential for given density + S.G_iw << SK.extract_G_loc()[0] + + if (iteration_number==1): + # Put Hartree energy on Re Sigma + dm = S.G_iw.density() + SK.calc_dc(dm, U_interact = U, J_hund = J, orb = 0, use_dc_formula = dc_type) + S.Sigma_iw << SK.block_structure.convert_matrix(SK.dc_imp[0],space_from='sumk',space_to='solver')['up_1'][0,0] + + # Calculate new G0_iw to input into the solver: + S.G0_iw << S.Sigma_iw + inverse(S.G_iw) + S.G0_iw << inverse(S.G0_iw) + + # Solve the impurity problem: + S.solve(h_int=h_int, **p) + + # Solved. Now do post-solution stuff: + # Set the new double counting: + dm = S.G_iw.density() # compute the density matrix of the impurity problem + SK.calc_dc(dm, U_interact = U, J_hund = J, orb = 0, use_dc_formula = dc_type) + + # Save stuff into the user_data group of hdf5 archive in case of rerun: + SK.save(['chemical_potential','dc_imp','dc_energ']) + +The only difference to the other example is in the initialisation of the real part of the self energy. We cannot just take an element of the *dc_imp* array, since this array is stored in the *sumk* structure. Therefore, we first need to transform this matrix into *solver* space, and then take the appropriate matrix element. After the first iteration (here done with 24e6 MC sweeps), you should get self energies like this: + +.. image:: images_scripts/Sr2MgOsO6_noSOC_Sigmas.png + :width: 600 + :align: center + +The two :math:`d_{xz}` and :math:`d_{yz}` orbitals are degenerate (blocks *up_2* and *up_3*), whereas the :math:`d_{xy}`-like orbital is different. + +A complete python script for this tutorial, including some more input/output, is available (:download:`Sr2MgOsO6_noSOC.py `). When running the script, you will encounter warnings during the transformation from the *sumk* to the *solver* basis. These warnings just reflect that the off-diagonal elements of the full Greens function are not zero at all frequencies, although the local Hamiltonian is. In that sense, we still do an approximation when restricting ourselves to the :math:`t_{2g}`-like orbitals. \ No newline at end of file