From 22ea568b8f08ee0b059d5dc5d8c596b882ea8b0b Mon Sep 17 00:00:00 2001 From: Florian Atteneder <florian.atteneder@uni-jena.de> Date: Thu, 25 Jul 2024 09:11:32 +0000 Subject: [PATCH] Evolve: align stepper logic between all ERK implementations (https://git.tpi.uni-jena.de/dg/dg1d.jl/-/merge_requests/189) This assumes that the first row in the `a` matrix of the Butcher tableau for ERK methods has only zeros, which seems to be the case for the all the methods have atm. --- Problem explanation: So far the LSERK4 and the remaining ERK steppers worked slightly different. The former used for the first substep the initial data array `u0`, and `u!` for the remaining steps. That shortcut avoids a copy and appears natural for the LSERK4 algorithm. OTOH, in the ERK stepper we first copy `u0` into `u!`, then accumulate some substep `k`s into `u!`, and then call the RHS with `u!`. (the reason we always copied `u0` over was that we did not account for the above assumption on `a` before) The problem was now that in the GRHD evolutions we sometimes alter the state vector in the RHS computation. This influenced the results to the extent that we could not run the (newly added) reftest only with the LSERK4 method. What was happening was that when using the ERK method we never altered `u0` directly, only `u!`, whereas `u0` was altered in case the GRHD RHS decided to do so on the first step. Subsequent substeps of the ERK method depend upon `u0`, and so this broke things. --- src/evolve.jl | 39 +++++++----------- .../grhd_bondi_infall_avmda.toml | 33 +++++++++++++++ .../refs/grhd_bondi_infall_avmda/output1d.h5 | Bin 0 -> 60280 bytes test/IntegrationTests/refs/testconfig.toml | 3 ++ 4 files changed, 51 insertions(+), 24 deletions(-) create mode 100644 test/IntegrationTests/refs/grhd_bondi_infall_avmda/grhd_bondi_infall_avmda.toml create mode 100644 test/IntegrationTests/refs/grhd_bondi_infall_avmda/output1d.h5 diff --git a/src/evolve.jl b/src/evolve.jl index 854b01c2..c0675e5c 100644 --- a/src/evolve.jl +++ b/src/evolve.jl @@ -225,33 +225,24 @@ Stepper for Low Storage five-stage fourth-order Explicit Runge Kutta scheme. Taken from Hesthaven, Warburton (2007), section 3.4. """ function ode_step!(u!, F, u0, t, dt, tmp_k, coeffs::COEFFS_LSERK4, callback_substep) - @unpack nstages, a, b, c = coeffs - kim1, ki = tmp_k - N = length(u!) for i = 1:nstages - ti = t + c[i] * dt F(ki, i == 1 ? u0 : u!, ti) - - @. ki *= dt + success = callback_substep(u!, ti, i+1) # +1 because u! is value for next stage + !success && return false if i == 1 - @. u! = u0 + b[i] * ki + @. u! = u0 + b[i] * dt * ki else ai = a[i] @. ki += ai * kim1 - @. u! += b[i] * ki + @. u! += b[i] * dt * ki end - - success = callback_substep(u!, ti, i+1) # +1 because u! is value for next stage - !success && return false - if i < nstages @. kim1 = ki end - end # for return true @@ -266,23 +257,23 @@ Stepper for Explicit Runge Kutta scheme. function ode_step!(u!, F, u0, t, dt, tmp_k, coeffs::TimeStepAlgorithm, callback_substep) @unpack nstages, a, b, c = coeffs - for n = 1:nstages - tn = t + c[n] * dt + for i = 1:nstages + ti = t + c[i] * dt + # assumes that a[1,:] are all zeros + F(tmp_k[i], i == 1 ? u0 : u!, ti) + i == nstages && break + success = callback_substep(u!, ti, i+1) # +1 because u! is value for next stage + !success && return false u! .= u0 - for i = 1:n-1 - u! .+= a[n, i] * tmp_k[i] + for j = 1:i + @. u! += a[i+1, j] * dt * tmp_k[j] end - F(tmp_k[n], u!, tn) - tmp_k[n] .*= dt - success = callback_substep(u!, tn, n) - !success && return false end u! .= u0 - for n = 1:nstages - u! .+= b[n] * tmp_k[n] + for i = 1:nstages + @. u! += b[i] * dt * tmp_k[i] end - success = callback_substep(u!, t+dt, nstages) return true diff --git a/test/IntegrationTests/refs/grhd_bondi_infall_avmda/grhd_bondi_infall_avmda.toml b/test/IntegrationTests/refs/grhd_bondi_infall_avmda/grhd_bondi_infall_avmda.toml new file mode 100644 index 00000000..a5265ad9 --- /dev/null +++ b/test/IntegrationTests/refs/grhd_bondi_infall_avmda/grhd_bondi_infall_avmda.toml @@ -0,0 +1,33 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 0.02143872868864732 +polytrope_gamma = "$(4/3)" + +[GRHD] +bc = "from_id" +id = "bondi_accretion_infall" +variables0d_analyze_error = ["Ï", "D"] +variables1d_analyze_error = ["Ï", "D"] +analyze_error_reference_solution = "bondi_accretion" +c2p_enforce_causal_atm = false + +[Mesh] +range = [ 1.8, 20.0 ] +n = 7 +k = 16 +basis = "lgl" +periodic = false + +[Output] +variables1d = [ "D", "Sr", "Ï„" ] +aligned_ts = "$(collect(range(1.0,10.0,step=1.0)))" + +[Evolution] +cfl = 0.2 +tend = 10.0 +method = "rk4" + +[HRSC] +method = "av" +av_method = "mda" + diff --git a/test/IntegrationTests/refs/grhd_bondi_infall_avmda/output1d.h5 b/test/IntegrationTests/refs/grhd_bondi_infall_avmda/output1d.h5 new file mode 100644 index 0000000000000000000000000000000000000000..6f9227a415288f4c49e200f88648973114f66f4f GIT binary patch literal 60280 zcmeFac|2F${{M}TIb=+Vq6`^R2=DzdMk1k*gfh=08Irju^Bfh@Kq>Rg+Q!UE8YE>1 zjWRW<<lbGlhI6j#{ypyRIp1@i>+x&<;r-cbuf6vAyzKG$dA-+;&{Wf4q+3Hr_G6-< zA)_Q?{^jHHpM|B#`^{p>^8E60i=_|bOEG*Y=1?qtASYV<v5ae}oO!AJA*#g!&82cf zeRcI8AD7!-j!TPwhX2bI&{S8}`yqkxk`6RWv2|%bTbli4U;FX-m&(h9|3dt${a0?~ zVxvFz-+w{e|BF`rtNk~8={&eB#qxRd-?jfz{CegUieB8tXW!**{fKmAG`}7*i{*bj z|1Z}iC!=1hx70!MrFO~57YhQHDlW$r#PT0YSoknv3qK|jnG8Ld%)*LaX2__Pwo5@l z@k;|VOYN;7=c8Y&L0mX~H?c0Y@}sgQ<0Tt1GVk9l{m*~=P1ofs)pA_y%RlN$y>#6p zXC^227NY*O!{lTuf6?h5ua9nvS7fr-rJ6slkH1Rx^TVGg@K@)Tm$ZMZym++!`>v0_ z->=o|F8;c)-<4C`2@TI&UhuCRH=cZ}mYX5YMvLdEAIFoL<ZtJ!-<4Asc#v3@SNMJT zolrg2<#G)J19h^hrBfsME^@N6A4^IZZf7pm==^2+hIu&2u1Wngk9|6>@Zg&=8EHJi zBh+qzoD?`PJ3K$Wg5>&6B+mXa6-lgMmr2@dYSJ%j2d?}z_U<XYYRyhVieq`;l59## z3f#!Bj)I$>w8ksSFNKVeL>t~DvH~%YMg=_PxG9)PT|c(F`?px`xj8c+g_)GRPnI%u z2Ma0c*}W$^+^nQ0rxVU_F|H=vW0pAF{bn_(a2jrpN3)UA&Jok&-(raJd#Bw3?4-jR zD0{40*hwZgGIG{j<se-%v@jYw$4R>Q^s9KfHy6nt(I*u1aFc9)nSP6VnCj)Z3%N<+ zlXovj8LlA-dOIC^A<aWN$}c8Rv}rBL;d<*=TPj|X*4aMuz#LwZ^j2b;{Vno7iJ~4? zT1Wa;eZk@T&^pqrn!HL0zMjPQ>4Z0K!bdt^y~}D%8y|@wPe1wU8vaGyNWaAo1~=_Z ztocb+YNhF8W&ETenvm0^O&dt3%$`3J%-let*ihIWxN0NGowf2CgXPA>{-63SnpUm8 zi{IEt>bTBBv3q_aX{YO}iK|Wmq=`l?^~iAnl6^R*)Qf$Bq{rUtj~u-wNb3IO|NdKS z;Pc(SX-1GF#S&~+q`HZ8^T8fN#^OyR1CfFoGg~&3?!5C)(K)}F1cPC1k0&;heti9{ ztA&1zn@<m&kCYQ4r5u0R&2>hIv|6#O>gb#hsWO&>)#Ru!Y1_Wh`Qi~_(iQc$5+oUs z#s8nlZ!u;<W5w)A5z_aky4wQoijaISt>y2>L>Kj|6kXI$cngXChkm|W7WK>gE&ifk z^Oi;ZIK&q9J0-TLADQ@~ey-w+`qhgs>NoydJn%z5dI?h75B<y}Na{cI>y#jU9KO6k zW1}Q#Ytz}(*#nX!$C{YJ$`HxL^^<>#7R*XKt8*nusWi$mk9#CZ4_|M-eu`F#B*Nae z&qqUwWcLYmY&t1L%A9N%eS1}E@%vZ)7L!!nbTvAqNL>n^whSMnNObE=pH^;>CLPkG zkI1u^CZ&tad^gFGCiQOk$jR3(z4-qd_${`0We538N|P2`uBBmITS<z0qnRi6Y$e6| zLCs*~R#KOlpmbEzRua{IVQqtnt&94R$&mi}zGv}nX8An&FWl#FlKnC*Kglf4{r)~j z>zDhS@WloHC4c`}kHtHefB$`;(2_oXeIK}}&;P&of&W20c$Rj7zq~L0zf6zye-AzW z{Qh2#W3l%?@9+OhPHg|>-T15f`&RZpc6|A|r1!_-|MA>1uO$;%Y(jFWrRDM=>pzzN zOZWFe8-MNR@9x*h)#L>Jtem*N*A!e{@VolN^M>)J<?`j{YoVor<$g;peV@Rk@_&3@ zUtWJX?ue8eiS{XlpYf*-Klg*bJ`Qs0hE`a{#ACv#n@zlQ<1u6F@p>hx-(t7D`=xzs z@t8T!MuA-q<1orJ-PwVVI4qvUqdjs<99C`?&bx0g7F#o#dr37X7W+17QR)%%TYOd* zUOem-iwSYAl;l#0#W=Mm>+?UwU~a|oN6sFO!D_Y6lUY8C#>!saW%?WtjWI{u&bHk5 zTU0(-1Ja_=Sgra~f3Ley*pU;*)}Jzp!dN&K5?(!w#Hc9h%xbJ6u@cmn$Fw>U8_x>c zEdK5+wtOD=8GkD9Q-Qx$AWvXtjRBntepZ6}T;oY+oUzMVdp+s5cxafqUzNfc*RK5X z>Y2Y2KBz9gvaQ<@FSw*qX2{`)U+tFeX7_Z!6SoL+jgmRw3gW4Q>$>cJji+3TNFUDG z<5#k?g!XQ+$44Y4?P${N@a=L>c`lIK;S|@>)|vX+;t9cXGXiUD@nc8NNj-UI^J`=) zXc(ezvcZLN3-o2|Z1CbXnao}9tnqj5XX-Sat?`ez)Y`pYt?;oZx+}cbt#I8OOV6`m zRydjQp$6wNCioeD>hN<v`0L}KE_<V&W`r6xM0cZ}j9(4oYeH<wGpfJFK&3uKz9v=7 z!o{>U#8MTj%f1<8I;esTtxM>P+OL9Hl*c^`c&3aEaf^yOMJQty&UrQpM}CWekEq;2 zm6Wl-qkMI3Z%J5eN`uFWU=mg*N46$}g@~iqo^hFB7+J}o)3J8}v$)do**5OCSl2H| zn%4)cj#^}wLYET8r#`iDfkz3en`e9S*g_Gj>+;eas#d^A0s&tlG!(EQxuAKRMFIO~ z&hzi$fI^~+XVae?@W-?3KmQ&qxunle4)|j~SwSu-vADnd)MH7Ht$znSe*S*H^k4cJ z{#1d#;(Iy5f2{uR=6k>6X5^mCmww&Y@5%|z)>^W>;Qz$;ezzXEr=;xiYJdE_o+S5Y z<-~UM<d+xxt{sBsS?pLYCnp?BKQ#D9zl@yrVhzouivOSbX?$V1`QP`G=#O{CuM2+p zFoj8GiJH!8*8<<h<4oE2YG740LzK|FYT$Sk{_4y7DxjN{n@QxWf{`}v_TW>M;KO#k znO3M0o}8z3J5OE-gFJIvP7v`B3soBoRKRPM(F3B16;PQ*soV+`kmdcYyK$-<K8#(o zO)4%2n>A?+ZBga$F`k`b?qoTLlcq9m`<25CL#$@0eK{C~v}6_QmP4Z)RfEOGa%h?> z`Lvdz9ITm-mM9ID!E5~VTzX#_v`N}Gth-hQeJba37c$Erfsb5Nuci#X-=eGvEG~l^ zDz7Tu#+Jcuv72sJPnSWH)h@Fk$1=D*$@xemuMDCnw{;Y<QdKh176<DQqclc2Q9* z1F^lEC!R!?!DC--MXU8?fY*F<%G^~3$&s~UmlevOp;DfMC#4LIQ}XX<3Mc_*IjQ&T z_e-FuZdhkvTPYl+w!AhZR1W4_*{-E{m4cn@^G}bx%Ao99G?!a;8K?(bsC1er0gkEB zu)=Glz<VTpuJTO@IJ~}E-)CD2oGwM1JjBZ(`pqn_Wp61o^5xYO2v)$1#;(9Yp(<bt zB4c>?)rH7p#}@EC;E3?_2#QOW{eC-U`}g;(I@%-F()o96B-?Y`;>M!`9gHy(KAi zAQj|bCCJ_kqYZEG*cCUx*O+IwPi|_0wB`>I)YOgemEJh+$gW0E(rK)_?$QW3?4qqc zE1E&fbyev~oqM1XfWLlX+XpIw;cs?FK7#5L9l<{K#~}PrLb8{m58^gU3Z9qfgZ*V$ zPCL!|;ERPIXlOkFYI4m}`<q_E=V!aRST9b&Tlan1F*#FUw~Oa=>)~lgyBWo=I5P#$ zJJ}Oml%^rQ@=4$h*=g94^Lj)&Xcq20?Kol}MS+b=-@Yf2M1xiK9&QVcr^VKLQ%#vX zro%|mXKuBI(P7-96Mj)Q=rPZd8~Xy}>98Hz1@@<2uEH8#W(MtB%Zq98X`r4DB3PdF zN2T89+pveLghjWN@4)U4uNwODb_W(`x?<0o7-cNQggLY`&Ipq|wn5I+#~0gSba=zY zfhepaw^rUoB_5+M`pl-I5sxYQtaFkkIN-i5Z13;N#A6%=XSUuTI3QPq{f79#IE<VG zk&i;-Fb~$5ln#kFEG%kjUE1?l%rtjb)5)u`*!?@^k0nxKF`1pcJuL(ee7Seq3p$Tj z?CZU23uYRznB@l<Dz=#zY<Bl9jzhLF*zFMm!FB!7n7w7WHdR<OMt8F_HB2)aE6Jwq z(BB)4jnzr@KG_nDMUm%IoV*i-h3&s62*y#^6n%e*I>7<A#kn}yT1R3Ml$S?dutZ`$ zuX3pG5gagVGRN-uhhWqcW_>W{O%T#MxMqFm%K&uiY>W5pBtLZ7o^0X7ljEqWa<|9o zb;r=fFXNhX4o4Axq;)`XnFoqnsi9Rt#G$93%`@mdkZ7eroPD_)I>lkH#IE6jViP1K zwi!60NYkl)9&&rso~jf^5on9jj$SxTGiHOjFG{Rp9JE1(Y3neNZ8qrs=d|<e->p!L zR<rbb+ast^`*ZjFq$LW`G}-BN-2!bmaCpl7qdEE*)2et;z#Pf1K}MpDhfu`Ji&ZHU zhtS31cg6L)%+Sq2@iNQd1IYfY&IC)X2{L%;bTZ%02-P2-Xe8;FpuR8KXTzHf5lz-0 zWt6QU@)w_|)Hq^*B5t0tKfc)rQNN#AJGI9ExsDg4SE?JJK-+}0PX_wvsda*X>!y83 z(TPqfDQ+(^lPwOSNYzFCF_&LuegH&iIO~fN_Mre#J<b9rU9>j7hRoSh8+po-wKp?c zpsRAuDa9;%5&xqTYfbu0kz2iYz$RyX#F+gw?7&TJB$x^1QmOmV1kHxx0I~x}IrD&q zVXGbbeo@xDV9FoOJGZoyzw}3&k6l2@dxMao3~P{^T>xs|;dpo`BoNh8dL<tB3PK}Z z$DhUd1R_?Sy>9D2Mx$VB*DFGIQqU<M+L4JB$tb>NcumLW6jXHW@WI1T$>`0RbE7XP zQ;^>Jl~(%UDd<s)k#yCmbErW(JHh2y4sxVzXf6mSM5okV&+;-9BI8Sv!R>-Y=<}|v zr`22v5p&yl5%sV_L?_TNrkhrXMyCo3&#{ysDeWSWgA)yCx~!6tsiqBG!>hA&$J)@b z2;0)>FYSo`84Z2-zIOB~`VB8-bvtri;da?ByB%rKO=;(JcO#}7xw)CI2hgf;sn>HR z&(Z1$`P{RQo+CMp4SI5KM~P?_EqiSgo&4f{S?}Z+lIEXX;k0`ctv6wDUg7)-HBwYQ zmEZpfO-8o|FiTM5qC3y`KJaD550>NdY+pHW;hQ<#TO>Gfs@HoxZ#-Lv2iKNu+jVd| z9us#g=unt3?xj*$RFP?q`!$&E+gRz0^E+{-%=tOv#*q|PABj2RVwvs<+f|+M*xg5O zucC3r<({U#Xo_*d@1!4ma^a36ZZ?WS%-1{O7gP`R8u>Wj*EjCeTupGm_g%_Y=<nF$ zq}qA2l6HH1b=gKT?Ie5r?w3an#!~jUFxN)u^{IBacHA~SW=cE!D7o~UueUAE$T2mg zvDOxM96v-G`N{?_>X^ew6E^tu=lSbO2oCrz&0x<xTN~UkOkqU=!2v(C-m#@2IA9$m zD_J7J0r?}s9CCB5@N3+GVVwjAl(kYaE47z`O4+dcWbTqM9GdA)31Sd=%A!wpuP{)m zr|UA(Y=ZTY&sW?n+6Z2TlM=1k8(=#1%?sfwen2`!VM0V)AMsJkRf->shczprH}e5| z&CrF7Z+T%_n?jwuod-%Rc7&`FUjt4w)%zpFxS=%mO|WAd7dUl|;FgVCaOcFuw`PW1 z(BDPhT$sQKW`j0CXU=f|yTM?eNIg3Yzb?F2yOJGV-(lK)yoe20%lCZtJHQ6TXEe#r zDzSlg{Ps;rLTsQ<9qdIpzZ%Y2d-`lYxEfw?*ByVnb~TU%ok*ySWrfRU2jAZoUk$5N z#HlseSYTrB4W&DlED-9hx|OAf1s(w(ML-@aG&w8R#HX=9mOw{WY6A->RF6%6iD7}( z!gwn?3s%@W6_y$EU=@sIVQWK6Rsog#NTE~{6J&GhUMK&;3S_?ar{W!%VY5iY3yV!G zFdsNPB)^&+#KvE^)dDj}N?&uY5?Ku^PxyY#Jjo0VMJ!tMlFV>fB*nK=c{T9N^V4$n zv%*ELSKN1;)&eGiXX@~6g`k2i-SjaTI7KacPTEf%R$gZBWQ^Yq?E$ueYsuxnYx0JS z*jss6<?xLud{7psBCp-}oTv?#ulR2KHf#(D`!(<1YBT~$KfkS?#14SPx;RnhFk|p? zr0HBKZ32XukUx+A05nHjdEXv&5bWC&;ze0pfzIrX<}QxoP-b{SR;kbnj>t?b9NFp( zHWuH`<jEh0cVE1G<&GQ&-ay^eTm6qiuvXUDp1pn`o>=$7lI<*fKGyiErz;+)+h2SW zr$_+F#(f?S?<N3)Wd30)lLYAP*uL*jW&*Ib8R=b#N`O-=RHrT1rGV9oHbr)wt8n5% z3LPKyEl6a@;8Z<(3p8(9E+|SDz?RR>aSx~R!Cynzmn^FQ6ht*On6&euCd4~;ZhHyT z4mG^B7i)yU)wlOXi}yle=fRZu@*&tU;id6OdkiK<G1dGNW3YGHG1oF^0@^r=#h#C> z#5^cH2g+K-FfwhF(Xl|nZg9;>jHsz$%H7A#Hu9=rgFYs@SH@JasczxOPw!Q+#9GGQ zEv>3pUae`QfvYNJRG?$TF{FYmTxT}N%v7-bTJF*<Bg&ZJE35waaAj<(k!eVuHF3Z3 zeB*GAr83rk|CW2SyfP-;ED$_8MZzkNBA(;`64sP%y+39(32Tqxk9}r>VR{|oDQfot zyXx2aUN9W6&cXEXglND{MGxN<Gz2VS<ME5pol2N?TT9Z(wMv*tHh+;F!2vDBi(l@m zQNZk~KJ7;83K#?V-bN|b|4j}kx%3<3PY(F!98hlgH`t&4H{w}z1-YE+(zEGLJr?!Y z`FGIcCkI^mFZ~REs=!}yzyhT|R{wW%!2iA9x8w{Kd;a&juKiK%f0`o&UD5dV)uRDM zFR^~B`C1R;BF{fL^3=o0ZDn%XH`f6cW=UFg69KR29dfyVTA1LkD&5go18-j6Zl&|B z0o&0~UV0+7Z4PA_mac(^*<oYJ8P&jUHYB{_TNTXQYbn}qRt4)`Osu&3rV?B^?^5u^ zRYF^j6>pDoC2arJPG;g&2{7<nG}fUK@bh%##VnP;JpZ};=z|InCM&$F<XZu^Jie~= zFs^{ni0xh4+!dglzEe%=L^&8w((XA^TMk!yV}`0YD!@A_V&~{|Ib7MBGE9|H4lC!k zw|p%xhb{SQ?`>x+2XCVrU0KBX*eB)6)C*-GFK|>f|6>_+-B>kcK2!z{yLYT(uq+2J zrgw*G39fcr-}G8+MmZ#xpV9C*Qwlz;A>6HZ${|A`lAMvT3J#pXKzgzQY$L357?aBZ zVH@eCjY^@Qz@ue4uN>MYFKx&@Rt^SzRn=_-FPy9}r`o%N=%0kkwhOiuaKJzeJqV}( z6qK6NzpnzWr7E_0NR>hMtg5)Wa5Zc*-aeN-QUkvE$F}zqyi{@5p2mZ<6%b41khnLl z9Q-pZs^#@6;8E4Bf$DwbVAn7gY}Hx`E*V|gtBV`Kv_fs>@=P=Eu_sMXjWvSIh(%=^ z-UJ+n_iT1M+6Z0;`VS~DHo-&B@g2$xO`u$ornh%xJ8Yb{^OLH007qs8j4T5lL&zzQ zlVbyY#Ph&W(JiHqAv?u7aaDaEJm^<?&IXSmqgvqgm8Qq=RG+VEpX@LkE|!VX>Uay1 z%6^>(PfP>&vhXA%O@s9vg<Ku&8K69qrW{Ai`%B;8jcl2QiPfi{_-f8T>xg^~PyT0+ zt2C}>Os2%v)7(+GSV@bS7k>{?Y^28$z6tADZeqYfHjXPeT&2ZYl8w|~ve99A!*@$> zrZQkhT<!Ko>#$*^jT}28Ci$>;jTb4H-$-HUH~dmLU#VhGnwly(oeyG;H*+r0U9!cV z`=9mAI&O;@y}B>*cB==*L49HK;qYLLQuuX^8FdU6k~G*z85)n_CsH~0#l>S{vk3_^ zhVj_CU@nI{!tq%B(LDxQ;_+B=`(x@Fo_MU{)Ap{NeQ{VS?TG@x({UL4Wi3ll@i^?c z#TTd7&tkD@2c4@{*J7~}-ixX#@v&G%(O}{8xL8anCUGmtCl+%)%(pd)6pKY%$KTM; z#bDv);+hVGAF$Vf&G^gHXe?7idt+^2G)7IPUfoRi0rg9A=~?$jV-=pQw?B(SV=m9W zou4E)VCIXxJH-eNXd-QXM3dlvoxX4PjSzmo777znJ?2PEWg#==8^HnjF9@j!lO-aZ zl}Y>G_Qs>uTcVNIvtrSk1Kw<2L6In{ocu<BS2((9|30-mJPciT8?EQ43PHHyeBt=F zV6^)^=l$*T!N{^kK1W<T7{S%gxW;H8YB<7^>=k+vZ9V54+wAFwBGsmKcCPV3ci;M~ zl@2|Q))>svFm!n#>Ah~|pZdH|!fCSP&^=yAYck{6$q$~WET-Rs(eWty(s|KWc+efC z#cwAax#@=JllMFJKXyed*%9Y=K6gfDah<{~hAyaOp?5Vkw+njH_-baIf-`dIz8A}X z%K;fu<hvD**du|oqWBF9Hb`(X!v4e)TV#ffxdn&YqK84-W1g+GMYenFigHV>(P!kI zv2~v{x>?owV6WQ|bmdw_SM&iZbX`#T8=aa3`j`{H<6HhA#4|6fMmBDXG;H#_c-XB` zp6T7O2Zq*YR#Q2}joJ(yxil8fre}`kOci}6LTym4?Pr|n@gc-&N_O;%zB%f;+<oF| zgC)9mCs5-3k;91YN_Rnosx@Nrv+#(~K8#q~pQ}EKb4O#I*RPp&1R>>#)$Bdu!6;WS zb1R2+2%_D-$BAYp82J~Mvz*uyf>g2$&U9D@qh@{kVAgG+$c|aHtY|U;1%3_gC|Ed$ z(oAWtoIa6?_UjZpd;anqnsM!E-5{TeQZ3`RFHoFAP2^dj{SQ*n7UoJ?*$0{E!#Csb z-2=Cft<(`Q&(0zgY<w$F#;OQu2;O8FStvqPY7zm;Dn-a3P19!g)gp9u?(5UtY{jV0 zdcSVDbrq5*pCl_ee;c7<H@)u64rKDgH+i_S12GD%F1y#=iChkO-Anq`fxhTj?zGz6 ziLTPV5rHqAXn2q0QT=uIk)Cu;)TR6(w65v+J2C%J<k+3VI{0`DJ^SXpJ!j`Qdg#Ym zsIzhuk?#xI#4It2l5#7&^i{{vpk`+9d)jH_wm(PmLd{piYu6Dp`-KS~I=4Y7Bz*&( zQ7IfYhLrF;!!w^<HfrGx;c^BwuQYMB%_32&HX7hE?=DN`9y*MF60Y9uo8yQptvkM^ z$khekA<x%s!R3OdM<(U)o_5BM)vWkfyvrH?ASBw&rtFMgI;Vc7fYurJ;Ox!R4|c-4 zx^C2%JaWXZ#CS@{a601tTJE}aM;-7vZ|qYsxdWa#eVoUO;DFzw3OIv^I1O4wVd?gG zPQ!kMJ0kYD$ZnICdBP8Ps9Q_hj_?EC|9HA}C*cRw)ho2OApC%3__kxBlQy`s)o39f z5$mE$KE7|X!Q+Yr_!bBbXvOfzB$D8O-<Sd;e{evpJKt{{BsieYu6sTV1P6R0r=~JM zaKO1Yl6P-9X+kwf4d1+{0l=Cnb>^Nbyt3g)%}OT0xOw;6tU3h{8yMNbJT41)Ce{w= zaoga)3ZFG6GG)MYbhp|95!<#WY^kf;3LMndG<qwfA?!`k3Qbx`Xf1ENVClaFDy80k zsl6`@>52)D+c|{be8$b7{ym#P$5_Goi`izNzdX8YU-u?Z`kpjHe?$<f8Tr>`J=h2o z-5*QZ_HF>tww~kOviy+qbj@4EFg_4ncU#_a<$B2ExUg<0Vm(YMKS;Y8wjScgyQNy% z*1_Dl37*6=YoWG@BME0&3$Jd4oo5f>hKJM{C4K#CpsB|$x~+~Ia-S5`?$6|g!-9{; z%c8j8JSF!`;x#T%)rsVBqu_+k+2&%>DqL`zyXV|jc{Zrh8XUe8#txdU8@-?Vtp&z? z0{#@u+;BaGq9l=u8?+NgG%jyk4Y_Y$HVb63gLZ+De(Wqa<XhV-VkQLNTkTooYRC@m zN8XvQbKnH^-D%k@&p4r%XGgk81vdn+zMD4K%mL0TT|6RsH$p=GU75fVSy<4a47(Vi z0D^w8oB11+fZ@d(ajjS-c-+`qQ#Y&x2dUP`neJDF-hpemb=xqAjvGmvpVf!wVYxyi zse>?nUxsql&J?gqp%NTN%|P=aV=Z>n6!xtj5PfTV5S%|9uFGIK2*an_OSQ+2K<<~= z?P9c^pgq6VCfCFVta1_)PrvpC<HnPunl>LuH@Mh0QR58(dp=!wt>**g6Zf4DfAIm+ zgiOCP*MmVPH$M=?#e&y*7CGslM2KnaR4qtO1hI1}f@_?Uz%k<<$JEzEka3-o&t91X zElkP2mW4@hZ^}+$lrICgqV41TrLIF{F%QLangXcqDdke7DFpe(%1Bo4Lii}(8e07Q z7OY;opL>H$K6H}!tTC(*e6(CR)op5Ed?4_OwsH$FT$-SuVj6(*18Pz|Q8U0Yl=ez> zmKoc=(!NsSGcWcsW@iLh4lgGEdhpT)zs;Dny1roAd0A{7xw)f4hzj<|s36>|T@CX- z%SdOKpoY;Vx}JWxSq;mRzf)`USrv1o8E|wXe1V+;&Y~Xesu<;yqNXt`Rcw{#f!O*n z70ja(8*ngE!S<$iv&s;Dz_Y$46ZZ%|;O?(=K3W6^v}M1MO(LR8DlpA~GIpno#x&#& z3A>ywYt~Ho0gIJ$IL!z@-~~Nivx$Qkmh|Fv8hH<3A<@&`&j~-E4wKWl5h6a_ul4CZ z;Rh@oc8z{Ya6tZ1ogoi`1Bwp>nA#H@@aDEEp^t<g(Cl+|O0Jp$cC1tDo({nQ|4;sc zf6OOZ&O3_Tz4$D#oJaM~&s~eohJSxPEotdF?7!PLwye`#oj>ZdT+Y1bkLCYC{fw5L z0{^Rd+wXic<eoIJ)Mc`N@U_$}?S*87m*zh;t`Eud@1~zx;F76`==kq+CHzQay34}; zl^=;@zF1?~kMxuKEi3jj{+R-Q#r<;j{;~Gohx?Ii7kmEq;(mX1JZm0Z{5H$S^U`-C z-AgOCTq}4Rn7?~_U-fE+xOx7=E5(}N_5J)&-^K>$_GoB*6ITy4N9laH3+h1M<=w!f zXdPr#c<oJmfdloW?XultIJ}OMu2FTvA*bY+nO8P~)4{$WOgC%6Sz6}7&ZZi$r1LSK zTw4RT3|MQv1XaUc@op}<*lIA&KIk|XTn)V`6d$)aRm0_;+XqgzRY6p!Pn2jl!4FUF z`a<=N;DcAh?#VV*LV(nIebbCesO!)dav}I&U+uS>2Zt+w{fcFB<fjU-WR+t|w5kNQ zwIUOh2^GNac&|aGssdcDO<eHUTLG?B<cz`g6_7w_NOrcp0=U_{%E?U$zW63+Zh+v5 zcd?VTdNLKz(LvsRTZt$)>?lr5ECT@xXJ>`>N-*wn<k_~i45H;(NSb&lc)w}85_Pr| zO0a;$zzt<U?)5N~PO%auLvP-n4kmcxc_Te$!Y8Q3dv5JXzDkHvO-^N8UjZtw&z@_H zFM|vBr)=|D%0T~t<JSnnFS)}!w|!m<2Q8X@@n@2CAk^Wq{p6uaFtxQ233^olF%?_; z9}ZT4kLhi_33CJ{)Ei{{U21``cK6!DR1J{I!(z*m-2%7RE+nhHzYX*BHNBmp&2U=X zMDhk}Gc*jYGiTe|4uzjB3lq}r!Un^%c7C0^(9L8vk||95&`k4EW!(EQ_@{qQe-P9U zLLxJ4TIBt(K#Qb{PCS8chqq3qT^t6Nut(R_na7~bFu0^dZw$Qky<cmIy$0ngiftQ` zr+`bE#->|+24oFxJw^v+VE?sq4ji4caKGNcDn9))*camz4`aWB;!UNkA1M}~L6Uxd zI8K3WFFq)3Tu+VRTe995?xDw$T#IS!r5Lb*XCEXa1DG*3VbZ<XL^h0bU3lg7Yz}N5 z&)rh9z=6GV-&5(O#*3xr(|cL262@H2eV9xLpWmay){GKi2Qhe9VUfsl4BK8)|1p{= z7)x4RBjR^043o}y%XjB$7<Ovl#j7i#5!lNxdyjK1F<8>H?tpGsJjR(;XEZh)kMT|I zm+goo@__kLtDl&}W7h;4FFluz$KKCn+fd8JV{MX0>9V*8E-1sm==vZITV=jE<MgRG z%$s}PlY63ZSf4`%!w8WVtVe4`8;~1|y`UMVY#@AqtUV)o-=bqNt7K<G(POb#ugH$d zbRsX<H-}8{4Urdol!d-bh~R+A-naB$5_!Rc!Ml#85qZHgUqvbn4WhA;LBeif7>$Kn z@8i`c@`CSAWO9xXdBH}k2?lCJUT|pHHNBez2h<XJCH?JiB-V2(M=*>v67%g#ED<9( zp!2l8DIZ-1%4FIvI<ol!YK<ko_#*opQjXcA{G=}l1@r26N9V?)qf#2H?j49lys{P! z8t<ZzfbNju-quK@C)fPI=5{118q#NF$T^Ens?2SV4T(Tv8!nDYESy5AoB0}jE{CE- zjZnosPlM4{SB0|lgdn7Y&Gqh<2}1AY&OTum4?<MB^BSQof#{`Qinp6#0CFnzG*VdW zkL>oWNhmGxMZ?j7pEO>2BOR)*T9+xk(Rtt3dAE(d5Z%@CE$K~O=;Jmz&P+Kk)OqO3 zrK5evP)Vy)lNzrlVt#)szGTb;nbWtN`TW%#eWC9Rzr5NVjn*Y8+zoI?!zCu?7{0rq zgUknt_fff_!W%ZpPj|SYkq@0@3+b+C=lh^0BNux#7I1az3wIY(6oQL|_Bo(W+$(}* zZdxM;*;k13o(+0X?-jP<iZvpMt7o6t?SLq1^{XddxS&V6shO(JZP11LwOs<@4rr@x zqv-7lClpaS+*wVuJLv9P!J^}U?hp4`QF1sS#*!VfjY&QTQSM4<3=Bg%vb0G33Sr1b zXkmL_Pbhjem3CD3VJNE7Ic~E)C<0krt_~~`i9mwuMA!uH#v)ysSdD1KG}OOm=55-c z^Jqgs>!9QMbd)f-Q+XpzI*R@(@QC4H7CLmjW!sgGt4Li^W}C;}tB7YTC^n-g4}Cfu z|AuS42(i2wSQkfCf`kSPCFd>|BV%1gWwG88G^WmYBmG)6+H{C>oHefwZJKzec5tK) zWmOxUU?Xov>%6|d+;yWJu`(yTxoF*mgd~hLj#+der+Gh}Z!~vNVngiOqk;F)XM+dl zO9dYzb@44~(Rq(i#^sBhil&2zM^AF@-OUknZuYV03YKxCOeT5aoXa?3<vNkh!uST| z?V%6|IyH@=-pKIT7|bB%i%Ii`KF=V!{=U(K&`*eqL%~9$mKv86c9J>sf(z%WG3$_d zhv5}oYp(|d8sSH<!ivV1hjC`*Y~vq(xjrM^#QE3OIQQF^X5ZuO@b(MywimBE;SpKx zAGJQa;NkCQKZXpu;Ci`5ykfL2xUiAjhwIVK_{RL)=S3>cxa!5X*9x_r@pBRdM-DSM z<E8EC+qwgs@Xh)&UQY-=ph&Tg(hQLoEKe3b%thn{M@v5`d`fV@Kyp7f4<av^iEY@I zjEJ?hJM?B#?Q!~MLAR$wUT|Q6@+~TY16Dl_PYotGU^cs2XE?zDnU1bfUP*Akv?^?2 zj_?6SQ8!o{5&6M8P37vi2tS~v`0g=tA}^SFb>Dt*;`_ZSobTB}aKO3TcG>JNR`^!m zlRmEq4tT^elKn?saA>>gY|Kg{I4NgvdyC&b2(9oj>`c&ylbU^MjL&p|gFje$?3@-% zH`bh~?^TBr{M$@q5>z2pD{_)QM+Giw=UJ3pRe=?s@(c}|l;L&4xqX{glfdkrZtosd zMd&G(O35gfgDTCPiUVu6!#W|yGxl0CphB6Me4=_QcrWakk1p8?vF+as*9sEz+L_}C zm!-gxSzG3@fFwANop*2ECI-x3H`w2p5`pg>8O*Q7gkgR*{PNR*%`lwT?O<#n1eC*U zb^c|WVV}Qe0;|Dhc%dRO(!F00hOo)e5Ly9H)H`h#!MYKySDJKD3T=R$3VY`lZf*d9 ze6o)(*f)R;V_KTy93Qlb`mdg^-~*SadHxk&_`u-U1EZVKJaBpIg!r2We4t-g^ZiuY zS~#Nmq``ro3#b{U+rG|nf#>*!fL3EJ;2D3<&-9uXgeH|*CEoEt&w-fWc?xcLiqBkc zo8*P}pDv}SG_Qwgw_MF5bZdcDHz!PKb}gtiyxYeu#0#so%E&*j5`)lBvnAYLm7qL) zm8afU5<C+5P#=;<0)@`&tO5-rsN>2?mTFJ~+J<g=lS(yU-uT=;q*Dj@U#;)Wczysj zg<W;7YBh&fL8BdMCgxDJu!14(kU4N}xWSdxY6T1jY=@)moWTCjkScGDBM6i!b+|@& zfbyICIcHxV5MlG`SRd#IA6TNg+!_2p`P0n%k^U0^4%JUyjf4Q_2+x}ojc}NG<UJ^r z9uAuYbI!K&#lSH4N;$gs39#}QndK(lWT5D}|2}I~GVH&9Z$R5A6<8xU23oi>VL0qU zW@kY*SZ)|@xh8!X(ltM?vS7*uTe2{Nd)ID4;s;jy4_<{(Ss)Q!l2!=Sl|fktM@k{- zhR63}k6PGtw<x)1uolj0EPNh!#KHGvVZQ>Y5#+)*ZjyL-ANrFRbFvv;LSok$S+VD= z*x^!_AhCC1Shqwq`M0ZbSZ>meZB6$TF_r`YF1%Y2dtg*>K~$ZDUCmPsqE%DH=p8A> zuWVAsDq`En(?ry<Z5bxTtl?@{iotIG2V!cN#(j1>UBU;rlSis`|5sJ4chCNp9UZFJ z+dyaOlEbRlFjg~k`iCFzg$`#wkr#YPp)u|p!2!<&Vqdk+C}ZuqQJaH^ykO4cj4@Us zwq%BQMk*>}eom$)FhRlszxY>D6C6-vV)l(4!2#t{w$3{ddBKV1Auo)Hyx`)m4=YDv z04s0JNezevOi!4vXSD%f{M9$u!il_K^Sqrd-9%n+&bR0J*9i{j<S9OXvq}New$L~) zM&t#b$awGkBQN-W$^la?7M~#g^Bhof>ACFh!2$a$|ESY42jn{P$MXLi$I^6M-12|7 ze!t^Y<ep^uOI`jW$0G9jES4Ai2laFPyXmJEz6A4#j{i<qf&<bVUe@m?-&<DVXZ$k- z{)+E~FO8qKyw1Oy@BMCnC0DB`Uf#?<?8kkj%M1SbJg~^`mizTTAMaMvZ*d!cFCMt8 z%l~_1CgrF{1U_hPgQtf~D6@xJ;Y7oLBCALX7{Jph@_o&4DE5%r0^SG~srw$eq&0x& z3I`S0>3SGtbGv7mR1fccW$R1R>f!W|!KK{Y^>8)j68~gR9kjQLN1ha_gSyYMev*PX z+`jzwP=XGED+)?iZ7gbGv`J|pFt!%riofmeCpcEV>BB3#oNB@M?Bqk?lp3gUIkbI- zxdu4LjZa&oRYPCT%<;ySDxmCWpQAZg4P{COiCHdH5I6ni9(_X<JbKVFV6IjL@h^Lf z{CTUuxr{G&fw~Ij@^y7tZdHO<B!gc6%1W@LOvuUZCio!fW2pP#DmbL1mV|$(g!{~! z7|c|vKux&BCtSJ$%D%*Ino_HPe5T$&%14#pkUZ$7Lge?Zrfb^6vZ)fr<|RZzi1t|& zYm==li2UiUI|4OptKriD4Nd!Z6~I+Tbyj+T;Dnj1^kxN>P~rF8f<Ctbpj4`<FS!<4 zr=C>J2sOYaGU4aK4fVj~vtp!Ls}VY!+iqWSuLnUr)9|+9YM|C>)UUlz0nfJHd0lzG z9@g{;>f~Q<2Mym=8A_&hV6jnT)nn)YznSQ!ewKC^?^H@zkGDXHc)Hu!Z_OY-)=Hn- z+ya?4j=Z0A9zay%`p6YrPaw_z^}}QA18~Rw$jN|R1JHNn^qBwT6R_V%WH;G9fvi2{ zpS;+ffXGgrTY)M=u;GrY(1*d-ke+Mv1Ya`^do8LwpS_rY;i(r_J$h&0xJ^W4e&-BO zyj1q)DVYJDOYGcTF|$xl!6~tl{1fO;R4*7dQecM`ZeQTqPmK*nl0U16p~0xeZXF+v zqQP7ZGw2WOV!#BR+}L}Pj2@$u(%WKjZ53uz{yv~NgcX~(U^(C@!-viDzSY4GiegUw z^u0H@Brx|hw!%DfNsOhR)SMA3hc(jHS?}oBi@ko#cICLHBQ`z9WZtBA5*su2J)bgp z8uR+Tscy0(3d`wC)?K+a7TdSto%1T5SS-?|Z+)v!9Htxl+$Fj-9(x+TH6&6e0ZaGI z3O)2H9{avmqTDz=9y2KJbM4WJ$GYXLcJpkH$4*2((|REhk2QCm4JTuX$NJ~xJofd+ zVY{^pdruPM1Ma82rSDIS56I-QDco%+7ISdCNkJm=f)$9}(=|30^Tx5Q-UJWy&&-p> zj>lq@x$^Ee2@a@Y_oDXJ#~5r)a(0>Bkr=GZoOin?;Ro!GboUA){D5|vR6M1GA8^gU zMR^*+4|t&Q;jxnh2ee*!@|h690RwMviun^9a7>Yb-|j&q=DG5*i2&gT6zfWR>P(Ce zXcVkk%}a2=y4UApf+VjXp;_USw}aWpBjzTzlYJI)aeq<#HSq$v{g7KdK=d5?IQB6( zo-7$<z4)lTus#u`cE?J~s3jm)gIjs`wGxm*+TGn#A@RuP#6sFm?l`ped&~6Y>S$Cu z?z~AS>MXLG>l(-rk3h<jw4Gymr_q-+Nk`RZPa%QNTPyhbPa(03QVJ1APN5rTCe71E zLeY+%5e9kn5L76uEV_Fr2)UjxEQ!+(L`DzTYrnAsqPq;LJ7o(`qUr@cv*eWlX!Gh_ zt&^|((F7m$b%i4*P|g0W=}8}akuXDZbi0}lYRXJ%U7>Rvi3Ozhyyo;mRK{qR5q2C6 z98FQjUU;FVw;kbUwj4(r-j{nP*LkA6GK0L44<2Yzw~SiS!xP>0k&D+k>V_Bw<9S+l zxg*gdt{-)tdLXNykf&9TJ<!Vlwsao>PegYpj#<3qD5`6M*4q89s4o8P<;)s)^!!uo z^A(@n(Z)Sd?}T0-M{GCZw|V;op%*l5*e9|eq<)fpP*3F)DhSTx$o~+5Ec@mAa_@wr z%a*RgjU{1dv$;Nx8(kb4F3c`QUoRlz+p}7DSOz+9`D;axTqbhdUb?VVFcWRjuDg4n z>LOxN8+uj%7tp4MRh}BD8EEkRixb|g`DnDprdjh+33{q?i{qGkDav=P@6#zLMKgCk z<QfN-qLCYYWV>WbQPwd*`)!X)5S6TZ%#&-C1pm8zDxCW^VtjBQqlmE+y(*H*(EQkq z;<4P+jJ|H9%6^sNB@qWR6+#XtcOxT9-tvLRcTf!11&L6dUPRfL(7&2)2$3IZ%N)oX zK|893{mj@#(YGp&(w=prs9KSR$0hUy;>_!FVqh9aNqV6#d&*y-IE$8ZyX_{Cj?^i) z6|24><+J@e^K58w`A#*npa%@Nm*t(_4PH$6ZkK2&Dq~Ju=T>va#SJoe-Sw(jQFlGu z)o`3Ime&&ha{jU0M`t^{KFj~cYC$Kw=pcT}(cT$veY=~2&%qg&y1Ie1;NXJS-8r}S zaD^+re!lBwSg|V}5#KoZu*U^2lB~RQfyV{s(W)~_C-Q>Jm?Rx|)t&Li`?hyn!5Kfg zrvsPx?u4Iw`RZ_TkQ3fokBr3#4wzj(zWzSDBVK$kR50?G1OAbF<)L*8_PA`W(?wFZ zJx=>9bdHtafU&b`ib84jxJg)u^KF6yvR;*py_#l+Go62$J40~5>NC@Z-+gTHhDJl( zg9HZ*Fb=&uPw>AjXSOA;CH#PgwmjY@L~uasQ{E9Q1P7E~-FsPu;DBGx?>=1QWQ`Ah zO{d}{{D2YPbjFm4@c}7fX7s}d4jAz`ci-J<ODO5fR}8MPfIvg(H<w-=1haa_(h@Rb zpzW}~GB98OY*sUTTgUVOuZt5C)zXDLyQ5C5RlDI8i-u5U!){=!7Ho1b(FPysRm}&h zHNh!Umd-6x9ZEOdbZ@?+0^1Vn8TL1lz#%Wj+t3|=WNFU%wp<D9dj~v*l9b>y?x}FT zSrHa;POWdZ+X3}j@7P>AWMO*!rk3eT+rXtFjh0S-E5J#Qpp>W55LYM_+xl7xGFB=M z&90Y*BZ}(1t4yUJE|b}X7n6kKub=PPhlvA6Z9AT$xdk>)+la5J6$ZZ^Tsm4t2;Q9X z-b9-t3{IPw_P^O73^xv49n&~33^@gvYaV57g7me8sV|QUfKZ7M|2NxBu=^hN+Bu0I z6cu*gF}c41H2qeym=_6vi+Y~m?Q#M5*zvtfM0*p+ru2VWQ$hF{6ZgHZyUz#v#`jBJ zOV|iB0h|w<+BZUo@V0Vk2T_RfSFq?lv<+k!Bj-(vWZ*bk@fH_(C2$%0aB0Cp6*|;a z<)|L1!Hd*cp(!#IxJH)y?NPEWn3L0T@=%*YmMDcLw%Zbh+OX|Mi!32&my{$M`(Zee zw^45Dp(XHqX})mU#uAv;4X+*;ISlI))wIZ5JVCX)<}JzW1bADNrYxlRLrdR;+5x(g zP!~QuAT;g|>J-(-O^W;>P(nybkk=o^$_|cKIRpc*Q&D01^;lqj&y?A_GYNDJdd7Rt zq`;F^7n#?4rNA?5|Fnpp6bSK_>3F*{1uoJavpJ`f3J!jz?DcyuLJ#w~`%Yxn;kE2T zqg``1!Q|je;pXmq2;wYLm3fm7Jr<5neZCX{t`~9jTvQ?C-JlJ49bE!<>4M%3t*?Sr z+P93>Yq!G8`DriPn~$LXsmPr%nkSGNuNu@`JOCk9QwOt}#-a54*t<=;R$@XEpWYl^ zBZ#dJAb0CX*p7WbZu=gSNSO2OF|9|$xP1(TE;I=}s@UzWSl12rR59+n<coP}YFHx& zS=`%Vbxcdwgllt)I>w(a(QL3z9s3gYVL&QH4T~)ZHDO(+hE=Y3=;HrT6>~Po85*2Y z#XvpE>~@_hHZ<3hSw;8(6`z<T#}eDwWt369mEeGhEQ#}ZgUT4k0mTO+r<F0<xWFu7 zD`l+rzR?~Tf&*^KrOGNJ#s`cR6rYeHIAFP|$lQ$}5|;ROR^kr90f$}nvONhupgF}Q zx0xT~14bQSxlH%~jU~RGd_qJwQscY?!2#cz8V@ZHe!$sB`fFzi4(R+X_3K@8MGW6_ zF2;xO10IjR!PQQ1z{x%tuF`*j17><IK3V+dIiTLs^Vokk$6DrnL&yH8({efA@jsUT z^YiQfmY-5DVCk9g@8bt#CSIGjH2sbX5aV3kUD^x(lYT(GkiVOLYOzZuAfn^H*Y8d! zv)8h&%kLZb$MMNXOK)9_U8?w#7cMLIGya(Zf5i(ceEwK_nak<@vG{*H_gDV6)Z2?q z{?Y$NjMw_>NPNrhzg&)s_5aZa|GVRc+*9xLul@Xextc}zpOq8mDYuB_1;49LaQMJ8 z%jM)`w0}B2AVbVz4boCY6+=xm;w7XMWMuS2iItFX`r-lxGKQrZA6uNI$%(li3jgEL z%$L1bNp7j`a$qA(?tXkWd|wx^=R`?)$#g<b>*rfq9E4v?tKO~eSqnU<Z_1b-YzC`V zv#I`|CXm@EaP{h~M$pZlv^yo$2sR9pT_y^R@L1RLV^LoN6hA?Yq1p{Fy*DNDyI(z= zo!iqTSXu}5PuAi`b?QJ<y{~4E0uECeaVBEYIQSh`s2yU)Vc*Gbd+7Oa&>$1Ma4#Rh zt5_S#NSj)i+r52nn@%n4*xk@xv7-joQ+Iltcv}r`hbiExa1C&5iD6qyQ3K~WWn;42 zsv$$d_2Q+o)j-jwqkW3L8m8Zq%f2SY)hnO57H4_58Vq+@@>&j6f|dNq&|)IbSb!#R zOqLkWFs5ii={n&TEaK(ArSC}8_aWZBr&tAKH~SA8O%UT5zKiv{MYM0pU%s(Ym1wtj zbBxNWDzGW}oELki61rO6$)8CthYRl*51khy#s`#mWFe_s1q-tBWOE|b(EP>iWZ|(! zkk$2jIyuk;jy<Ixw7nZZM>|ZF)YAY1oqSt-HR~ZX)O5aabuDPDh^-9es)DPVnJ>;U zwt`dI%Lp@pF4)=mF-yd}4N9u*d(W}8!zu%UDJr(XXL^UXJFHsafc?&`Nl%+$y4C%q z%e8Jeu5jc_WkWw~ciy0PY5!AL)g$`SZfXEJh$O$q*PlYfke<u^^-rNAnY6~F^$95J zGcq?2%TIUdky*Nqg7>->s}iRv&}N|!EwY`3#33@-P^DQgQrh)2@8}1x5exiEJO3V- ziZ>vC-S_ZBrFy+n+<TxI8Wak2pNC<U<8l=GR9O1$IcXmYdhB`4SsT$qbXaTec#Xy| zBla|E)jJf#hzYxTu_T^n!R#}XQoie~!s@t=n7-5F#%Q~8uY~w-#zZB%`&1M)u!Q3e zn0lqmFwxUgdV5w`VYEk{vFAuwVRXm7RhHyCU}{p;A!HwpVf4;P!mMP$Sj#1bUCv)3 zFy-R8u{qONto&**<K9E@*rj7DGTpz$V@kKw)}LO8$D+|u<Ina9m;lrCY~a-d%s;04 zQMF+L#@0!CGd~rN<#Wx-#GQ@D>b`7ePTUiZNsh2npCSB!Q`@bXb`awO>bq_{`iAfW ze&}3rjGK5*2ye*gr0>CTSYi0?T7Dugc#<8RHhUV2)fH%EuDed;1%Le3bS5em6Kptu z{0I)H`tsnbY=Q$8UNfONtQ?C;eq|~N`y7K6TelXBS;k<oBI16J2|wT)*1hzv2|wWH zeS16&2@aTZSozduBI*>dZlxeN;7IG#{854f9yk=^X16~IQ#^fn?h(NO&!;}6lqALn z980NWJxBNf*Xx}yFeNx3h3Vr*v(wj+c;CIev>mxfo<SVDNq-f+G>|lVUUCUhoSMBe zR+xbbALc8(EJ{ZQ0%)(~=$u0mjoSr}3#6j&EmiB9Bva9wHty*`$rLnFUo+Y5mV`{H zG_rw|fFjgS>7EdeLmLhr${tOPLF5k>vK%O*kkG0NOk;;5k<2bVuQ9bq#BNVHT8+-4 z3&V$Z%<eve`Z;1}>F=LLR4%GAEeB2^kCo5Ar38f{F{8o*pNYI!dNPf9ORW%eKhtJs zgg^)qP``NTZgDW$c#HOh-|HZh`)M^+QWuD9VglS815Tn;e=Y`J{3KfN(i1#tdIBA6 z<J?Iz;e&Xt-QZk<`J%5Ebq+Vl`XI@-6Fb|TeNk#(qvOy@UvzTK>+=sLkD+RvSG{!F zKIl$Lv~|WiAC&57{Vcx08;uF?kZCRQM)&M_gHG6aAvrz{+SWCmNUX)7D8JbYDOu3l z;`H8V-}@K4##>LJNj-N`m~IT(9&|%w_jnZA?aiO8R1%Fctl8;isG|{C<D>nIbdjh< zQsQz27LHu>-m6Ayr=b^<j^dOamk>wonnwl0nW&8JRzpZx7Ao|%=y}YNg+dnMRQ1a; z5t8J!4nr5v!GM?c;e}U_EEDltMp_BFontXszqTCduQ{T+@S+SAbbr<9JYSBqIiFrU z$ybiZz2Ek=%$B0=#EzE@wPk4H;N1l3r*+6Q`pR{VD;;QRp>5ZBsk^A>iuuBI^1JBP z>lNNUw|me7ggc74^dL@(vy<|(cad<K<8%j04|*oBqj*Q=WAxrC%gf7f1idjEUU66Y z1xls8wQk&J47CEhNjHCiKC5dtcUZqbrW5XOZd$)XYDH-;WSm|hZj<NaCRabAg$&IV zjGk1ur`5Ws(J>)B$7w_^l~)C4^6xsva8L_>`iwl)GEfWOEBb|6Ygi8#KA|xrK6?N+ zr}S}~HnYN4`ihY|S~%crbQ`^dBwg@R{442ilDXlJvWtAFW!-R=2T5{1@@{yUPLFY$ zp&QP|N@te~Zn&DsD>rX~2ad-0-KXz!!H;Gcl^>vT!PAa4j#39X<J6d*5kHX+%*y!W zRW~s{pq9mJ8rGH0ILXV|D42Lp$Z(X1CnNEmkd!I@vRNW8_{wqZo#5esS3U7jj-YVB zH4-=)_jM3?!TNdeVgv`|4t%%yJi!59WTm{Q7PZIiof`Yi2tS~+W7q~$f&&V@(sqv^ z{D6<m&v36HIG|%nwB0>|16sdVBEyM@akssiCGvuUG~VY+5FF5tv1H=$gf*@!+;RIE z!2#Wv7*6;T959zLvVD^911d&{@mwG{V2NV4nGBCT(6ua#dfC~6m<E47i<}kE^IHxH z_E^BBzT0CZ?WSPVIVU5OX9RukVpdbh?t?7>mg!8)hOlP;3Gz#83}I#6wXag#1~9u# zETVhYUU=$yoL50!7mn6ab%k%*1uF3c6{?OqfkelgAIhc<?|5anDxB3I_a-^p1|u~n z?deRkXjO$ImVs>pt}0+iX085Godh%TUfiKgN}z3Uad2n#4#;>-Dtp-_2a0LUZ!bp6 zftf@M$zi`780xcR-?}6VOg!#BzNXv3G3+Eo{wW!V@#9y#1Ja=4rgdmxmoz+`EsNpG zmjp=-yYT`VaX3e{ovh`!IOK_BQ+)cg1v=l7t)?;)hi9q!UK+DvFxULBTt7$%I$y8k z+tjxOvYbPUQf9=UB{`>^*=GyrZqFRkcG?2xUK%!k3=#%9f6w+<gH6C3&bFziT^J-U z2=y`~h{Db6!~_n)-<Mi0-xGOhCxiuaDc(M_6J#hizOOgb1eSsDTN!nlaPpOPAopWU z*igU9IBQr9w5^}FUN$`d#;460q`zB%>RZ+BcmGd&Zyrun-}e6_LzGmKNQMwfnP<*r zo+F|V6`3-RDf5_6W(_KtLLwnUhR-sm3?ZoyN+c>2sg%B-?yH{rxu5&FkLUP(k7v5S zzia<-oqO%I_S$QmR{OZtYro&`lfo8C^p-iYZ0%r!$LzqaXj`1Ck@}%+zb)i+zcT6n zVhz3LWPQ7xoxwn=r|h|FAbhtLHr#kH6b6e2DQI*<L0m6O)O;coN(L`Y9u5zMqP?_- zJ0wG3#7uYo{q|6B#90|*k}=?9svA~dkOnlqSwqR6F2Y^vxgnjri*QqMLh{+M%i!QQ zW%_YC9k?mQ4cQ<a<Z=sk+&Y>L`s;2l>-OaWLmkgs)#+jolrcJRY`P4}bnO(RDN3PU ziK97_stk;FKB0amRtBC;JC1Al-vv>5NvBKlci{S#+zUVYnxR(SGR-;Z5rhd9?$aAy zhEv{aR-3kLCfb|03^zG&6ALP|nIl2m#8L%qP>SRZA}yrJgpN|27=6VvxTLB;7@qf6 z*D6#Y0_z%frrD_zWwa{Q{7*HA<;C#UB^ga3fxI^KR6&#Yeqx~KoTnzSHP&gr*kMh= z?uJTmeTN2dz$oe4Fs}xY5x{UjE=rw9dRbybZdWJVS53DIepDmQN=>tAFRBr%4)<5@ z!TErM5Gk#Q=Y-sS5`VFBM3s=EO}a3HeZgiunV)8cRS3-^+U-xV9+1&<Q~n*C3;5C3 zU_cXdz_mL3tt{9VoaNvYritf-^gS7$Gr)Slahu^`ZOj3G-1b*<#(F>%j>~l+m;*W+ zxUg7Y{->&3mn@6xoQ|CQW9$pwdHY(EB-R5yk!l~8#T?MZwD1up)&nvu_Eqk|`G7Aq zxsJ(U4k(gZwmkA*Iwxe`uREB3F8#!RACEKp)t~>}JoCS?A1S}?O8-~$$^XWFT=-|` z$DdsAH%s}m{(BwxM_jPy+~02fpUnj|jkWht(BaRZqFDdePx{UJbbozZP5iHQn5+He z@w@(7{LR*WE&juTzq0=2oqz0_8*{j$CBNUGn3=brxgQfeir=67`+R@b-&;~s?4vkI z!HpMxt@(-j@Pqx2@-H>Yhc5NsvdW*oG5vhk(_gcKe>2wKzdb24{IL_!U;FxN*`i!L z)BC{Iu?PCU`1NOcJ%ndWO^#fPo!F;q8Ru{F0Bk9OJZKi~!69R&s4qv`fH7@m+fnyc z_#pF)^GbgUh;ES=(RtAV2L;yiY8-8Ws=bZ#rX|e~vhTz0zR@PA-dGZQ3+vqo%C^ES zZjIm->eW;6paE_vYLtOf1Kj4gV;b$&0JOGmU%Ej9OvJ@@7>PGPiqHjV?&Bnor)*ak zln{8?w+if8sE4ln<*FdtdPrNoW3Am(50l~Ej`@-GpcuAzHP*Bq9661}?H6i+t*p;) zR<jOTR(;^4bE<<_cEb7^&L3>QZtuyTS__KZ!<t+6)qxpD+1Ts2I`B>HQt&UpoUz-b zZ-;UXWNw-}Torm3sA|^Py?3aGtxgYwMAy|otlug+X%#%z#L##PO;t5`e|$yC#$x?n zN5X;o3YhyDYuB}P)PP>A%a6zc1hv)`Dcp>0FnxmYxZt%GaADPBQPXGu<-F*)AwvXm zY|Mh%mG!`$5&yV)Z7mFSGmC9uZ3NA-Huf(io$%wczS1ObH*7eWbHw80eK70esrfN` z4@CCt+VLF!UAEfo(9t$)hgL1RSf;)X;PmR5V*36Vo}3@FqwanV9}9NdMNW-C?lF<4 z7iNaxTta&G-m+0p6YpI=x9TN?J^E~U%jyMSPgi?g!Ap2Cx`R`}U<UZs&<GEid;ldE z>RG!B3y^R?=|xGx0whV>FS8AN0;_eE=cpQ%;eFV{1cA`65dJ{n^rxI<VBE4X>>Cpm z(VY=WvGkFSc&29}#9h3WNM*iWXn1@*p<ZS_S2VYQ$hV-B4gAPLtoLwuB$UiTr1vDe zot9=N#s__!iUN6w`K9x@e2c2Yx$ryFB~?d>Z1(=Fm=zB~IdE&m^JE;dwV?dr?Xf_D zY3M$A%QBGo5<0-dw=0w|k}KPjToOrWymIs0x*SVL$Lh7e?MftW<|s045=bFZ%dWR4 zE~XG!YPGF*=u-(v6(yat)v3fPQ?B5d<rLzYbLi;H$`qpO$)&Q_m;*M)KB~Hp^?+aU z=A(BfCKGR;$T;R1Clke|&D(c~CKI5oE)gr1Oq^K>E3{`$CZ18<r4_}#VBu%!=fb0s zh?l}ix{TNt?0RS-BN^)fvm50iU*deg+(DBIK3ETEl-~2e2Xnx)Y}!RUcwX@1<<W*i zSP!WF=zwxI_66S)mB~GdeZezSTLh{w2bAWcI&%thKw8sP>MB?VIJdEQHUo3OE8AlR zzvFyBmljv3#`%Dw3YM#=F$Xl(eti|5oFx{Lsgtd+9#F_f-d+gn0j=-qc`{-S=(hjE zgDk7t=tSlo(UTW%p?a0W1ES@H$daeNNU|UgA;nw~s;xQ5%{uo%N_i&M|4mfwjLAR* zpYP234OdZ(vUqJD?^Sd%-%@eQ@MWZ^e{w>z{1RHao!-hWa{;-n?%lu>l8V?SzRyc6 zC!;eAJxd`MlF;m?S@$C^5|RIpu*X7^iReopt%STyBKq#6vM6Gnfc)%_Y`*mFEVBIW zwC1j2EK(KWSr)o`22JmN@ygLF8qv*OO=&leMtSE1Pka%KM*W(6C_Vc$(%>$B<<}mG z#19UxV_1$rkJkpsH|>o;4mx6voZ6vic0MJk2J7!uhSW?(4u+wU6yFusz%Z2Wuw`k* zFAPaqrGBB74MC3WB0pqK2BI8Aaq3m}p=h&YQ38*95GwTzeOLSVBzmP5Yjm#a6iVGJ zsHhTo3caK8`0l3}fK-*m_$ayp5c7pt#bV!3w6wkO)@tzt^twZNx2$(EDz(hnRkSk^ zU8|S=I@1)7aOmMu2Xj0kXEU~TE5{+6IndvGI|)S<_~c4RXQ5N#5=F^sx#-%Kl-|{E zvXHCA!bZb)naDoZVj*@Q6Ma8d-YK&)3zeJeZrwYWgE%VC4xX;t$RPBBM^t+iGND|i zd_jE|<!Gx}tUFMJqWV~zEYj~HO$&Fb$NBZh@cLnc=HNOM@A~>_9u0|%2FAM#eeR*` zX1fLvw}+_y<b)x6U@uY)SlOs<jk(|UgXwpzAEEVgQ9%zho+2UtEcsB@C#dpBrMHyX zQ>2?^C7Szb6n(n>EzNZD75eeel9SGC3h_|P%pYHxLE|)~T0NBS(J<SBIz!Y4G&TKL zOv>*AYCdRiRdr<vooCTj`RGMO+K%2i)*`rrjDF*oVpF3*hIqT0pPV-(H3eBX@1H(Q z3NvpVZy~HndG}U2D>-ZOws5WaG5%wu4d}g<qj4bLq};y!!N-|o)HATI<Z~y>gqha| zCwh>Q{IhYo$b($CP~Dqc?m>n;8Jlr0^&n{;Wi{<O??KK6c{lCS@*uzF(|c)@xs!W8 zjlIn6b|W?S8M%L5<3>sbO{Ujj4mf{bp~8)DC3_%Qn;y>%W*KDLf0)vh+)V8-q=7jg zWBuBXQa#S(np-{k-#0mvLVhp!#IPRF{EL=<6rK~ZEkoLC4b}rz*^oTjFbCv$z4yHd zo)_$0wtjOA)&pK@eD{G0b3knM%#6kyke@TP>LliXUtO6;-?KT8N0Uk)ieVj~=EI?k zIm`hw(Ydal`G75kX|z)~A5c?yD5e~9z(+R{y^mvGu-%RKjKY`$@>d(K2x2`T{Y;u} zFy?>~UY!FTb{^opc%jEX!3}H^mpS{LonVW);d}&xJt&BsZ;aVv12d7rR_E3n0bTov zNsAYUfR;!U4luR=<%9aVdyiQ_b9a`RD7879xL>DQmvI1+4e2!wyxtGG*Y~u)Wi$fq zt;sAw^#%|z+kQuUNe|w3@MYPjVm`RRaYyMDU0AbSUdsEd4ydsdPJjNS1y_PQ<{n<s zfU*{a9$2LYldGP6;0#xRu-2~R&j*!3+4<{4`7UKx9~fF*M5hc{=Qmz1W+i}PZbG~6 zhXSa}@LI4Pl82;A?U~`6a+qgMjI#wx!oBcoPfp57fe1yWWmJ(Apd?R6mJunSk9|Ua z>Zufzc6HNE$Vq_6^#!nF5rc$hQL?g(lJG%t+`aIN7!;qaeQ9=VH`MGL^w&PP2abI{ zGG)lT2RyAVYvf6ag0cO{DWMCZFuX4KOtGpYjQl878osRwVh)jyjrk2>?-Qrbk{Wwq z#feipS3@7xQ&H2uO40+x^(Q5ApJ+pKSANY@{yyN)6cdO`w}H~tr(6t#9f9K>d!_jv zdpOzk*j1R-4p>*kPWVsQ0`G(D3tz45;H&Y8dHRjVf!;HsCbI1W-1Lj1eg7&HE}v$2 z*tk9%HhTrEecm4mH`QJ;xF&?d)0w9?#otB2z>cx@<dsOcZ&U4Cu`3$3&NsfHluQQl zHaqu^ix<J#fFmzx-(@(PzWHL#w#&fyZTF^)0U5ADS<b_ED+j_Q>pZUHW<zJc%Dv*@ z9N@8hOa0Qh5Z0b|^~}~R1@=N7j<C^k*dVg?DCN2;NZ3XfARbu_7ouK{KCeblATzX@ zc^<(I{e<l_ds`t!gTB<AxfdkI2Nx){zku&**SOot9E9bep)G4tb`viCLUZ49C5fK# zp)lAlL-?v)``D%}Lj+S<;9oN2h{e9OV>_3C_!_3LmeE(0=#bIyD(282F6WrCoV=?^ zT(1(h{=%w-=X1H=GG){vnDxe!c^R~bv$>qEn&X;;-k>{cxr-*@DaKLZiuHhdwTC!& z;rB<#P+P^vsuMfn?lRZ$suR!F)i9;cs}Xf?bhH5H11_#f)jrg$Mu=~fKmGQo8u372 zCNUlRf{$phnc14E5;2z<H~L~dU?rGLZpHb4REM@D(PBMd$Q{L0dCURl-d}mUgnhx~ zCo+1zVh%XVelTYR>j7^LJ>2SvIUsvjG>0MPfXK$;tvKd@3U*wPl*vHYr|45AU=C>R ze?RLd2VCPC>`{UBfT!FQyBFA%h~6I$ikYnyiASQf0@tt}@JN_K@KwwK`A&<hyMZ~N zs)@cP1<n2I&o@z+<NE?e3K0Q&4vGT&bREUy4;DU(p9yS?6ieY!dniKx@^q2=t|SEu ze!8B*`--v}#m^f$MvB0!<a7%DzdRK*3ay|J|I1TKE5o}K#9y9TR$S<!NXJjtQS9kv zd_nP3;>AdD@XCWFil1sFMv52rvi|K&5q|y8$A4`H{+?_7!EJy4DgN($bY{P_Cx3F_ zpS}G1dC^CI_2*9x{PPN5MR_#;*OmS!2mYfUH~$&>@sBuPb=I$=|5bi}pZsQ1f4lbo z`_g~i_tw#>KaS>~>wD`TF~QLCKhEZ_3c$ZTuRqpuC<{brmT`0SflBvNi?{oFA*=F3 zmqO!17*n8je*Ll&I0AO+@=!m3W?nPxJpOj*EcJ}0R=5WOZ=b2s#J9n<Ybm)?*V^Dm zAm^r`AFXg*f~F?_12)A4bh%tsY6184OOKv?YJwU5tDP>AO|b5?V6pC|CXi5Knw`De z2qRiapDAk^VJ&OjNXd~#I3lc}H~gppDi@`nyLmOh!cD8!RICl4L~j(RXhedi_n-@N zOFeXw3;Cx5>!IyV!8;rOdN`h&^7!~rJqXV2_157fq0oTTUt3rQc|8w)Xw~D~!NA9w zKBfqa>U`-b3Tt7XyzcgL*;=Tr-845{RRi}ahjKhwYrshG+IN*7cze_S?^hQQJYrgM zLO7Rj^Yt%5>I-+_#Pq!X(~LTpZdH!ynXZ9FqiKaRxwQ~0^jfKqt`7Ec64yMXn_!K| zuCx)_`{20!Mv10zE4b#XCR5kcK}q}zaS6LdxLoO{(e2*`L5g3Tzd5(T;na_DU-xu^ zSeNMjuCG17{l)4!PuxRr-Et&g%lmc^no{VBQF;h(#bTT(W}kw??Yt{ZA%kFOpBH!% zhJZ6IzW$-f2ry~+J0)a|ffT=Ic&_>=Xt*oi@+cXFi?M?rCsijP<MGWOb4Kss#`T*$ z+l=Pmc;)x+N3h<H3~P729WoC`qApTieY5~Nkt_0NH!MQz<<v&;{Y!9gQz3T<-wH_6 za#OuBT}4zl!`Hz(bc8p>0}t6jT4G&&gLB0^4big1;!Vag5DM4M?+V<=M7&}UUNkzk zkr1lPc=a`S8*xU1WU4IMO=utGdbBD;fiT<TlVPz-fjGf^IK5|;7Lmbh-D^!8A}FRs z^j{rxA#M*QoLju-N3bTkopc)xB~<xjD?WEd5_j&`5duZ0iRe4*1NMcd38kk2Sxd2J z2+shm8-o>R2^qr7C51hS*d((hz&#>`s7qC`j&VsP>atrtY*R}mZtuF3wZfN5JbX31 zu*j22Y<H2pdgNOQ@mb>aR1@ZaIX>U3bc|DoAM`VvvD3*!%Vb<jNOCgK>=3+tPA{2I zrWTQp{>lH?mhMSn9%w6~7}<b1pepq`jb~4i2-7-R+t*kRsG}-iOtm|S=r_>oamV?9 zIu~6}rDF~_$SYl#iTNLUMQy(ut{=E<)(B!f;LbyDw`yVzm~7N>kOOnTXdzoq7VHaN zlFV2c$G+fx7j|C-><hM)+OqpK&Ic3~`Zkk|^8x20eKtE_4*2k!hIDK9Ib#3aQyW6C zFPOfqW92c<2jr1n>sgFBVBFX9P0pADc3Ql*;>H}X{J7B374>p7GJMCVbh#9XU!Bah zI#7ZR2M@^R<rE^9Wns1OBY9|COg<qsI0tP?afgZbS?Jrx!VHOxS?H5o>KiF;oHy5G zkapWN6OpF5Bt<9o6F-j04HUYJj_weOUs%3?R_G=6i#MmCk7?g;dUvIw3(B8sa;H;} z-6@XntxTyXs%Gv;oM{STNHck3LYa)dQ!>xT-od==uF{$&`ULcrhhfXPhO_8R{pjv@ zrm?7pwr*~!Fb16yl(8Akk3ll_at%DmG3cyyrK9Rl4BGK>s$`|_G}>bPB$2r_5;1Mo zxaGL&G>Y0sdGYJ%2&9lAj7O)0BiguTL?0WD__?OM4BNs`3eSe+%>m(v=|{-fG3Q8R z;31dKP#ubdZwud6Y7If(K4|Eq2Ztdau#A}E4@K};J$EiI6j{?8JjYTRilVcAFgC_U zqpHxBl}N=@q`*9-WOpzP`DtmK5gLw1Synv>u4Bn4>F)8galcE*bCDx?!R8XGOk2=a z{+@{j`r>Zn5_u?Xr@CUBLmt|w5OsK;U^cpGdT)eJJs-tL3z&H~mZMvlO}nNvD^TCz zD7N$C6)2adQs~p2YILdB{6(R6Eoz@r^lzWKiwY=rNF)l>prNh>Z%Oec^sZLQ<!t8z zWIx~gi5hhx)yb*iCo7#O#&`AP)#`4PvO#hsW_ur!KAIX@B-V!_qY9Y~Mf%Z*Lv5@o z;}B8|j*V@1c!7i`wd~fNoj_Lx&(SnJd4-azE}IY|<0$>du8kg4uaRCs(?&bN85G@T zRT`W+i%9Xxa9OWqbU-K~uJX$|Qez;xbq5U>Iov5gX3cV9eH_c&8Gm7t^HjjWw_2*? z-HL+!D+2q;$n|PU+8KvQF7};QNG)5EyX?lv%~u@A?L5WFO{=kPPH6hVJ6b$%cMYFc zoTU@#u|j{M1<&iKlIx*n`Qc7(mhuhSx#B_IbQ;#`7xW}AhZ*?2dg(#FotkZ^t@9v{ z>})7)sr4Ym##;7Sq<D~?)BTPGIOsNq@$d%zQg?D*XZs0}hi;_!+M(H%wQl5F=>q~3 zL9XOPa<6g>=6^b>i(i5;544QPS9*>)U<0dd>1Lb{80qUD^R>&F6tkf^QHS#Zt2Be> z_;EfU*Ne_HOUwZiI@C@RnEzFExLW_@fE7{X7E!DRRBIGD>4P~SzoMnUJ?smfJaBb& zFXn(-(*wA7U=G;cEp{5P9`MjQMMf3O|BjSY)Stt(SDj_@F6MxNh1V5gFb5RW-|wb{ zeZe>1-cVM-9FXt0p`t$KfG>!6=AXV`KHdeF)0hLYr1M%ed-%Y!&A0vX=)GZHI?6QE z&jU)%Uy-l7;0(tW?Qgl%+Qa7?8=|6YjsgF$Cm(aHH5_K(8yiYI3Utydlw#M8!k3I0 z_Tm#qAbalp1dFgG_�J+Cys&bgmcpokL7P@p!{J@z2JP6hdtg^3Vu4XXQWbo7{(e z!jzp0bVi`0MLRsKxexBG8s0Tswio6~FUBqR>%#T@j*rsxwIRqwH22=T21FKpioRN+ z2K$%wLyKwEK$EfIeIw?Cn(~JN6?;`7OkQtAyIU2s-$qoPl~#sxyXrZ@lN2Gz>3C=l z8^HAD2j_Ja<$<ro?gg8?9Bgv5D|(bK3qcd*TuJsaz-=pgbtA7VeEF>GX-cC2d_*Wi zPJ<N8HGj_8B_IjcDrx+7ct}G`)AJ+yy(FO|Lws?&uoScuQZ*Uxl!80$@tbwW2vE5w zVRh{r=6-1h2fhdFhl!<kjYs_Vf~LZrl%3|rF!fH!Q*6{6eEi(`;#<sN&BcP`V|q3) z^~6Kb^4W3FwAQ2doOFh@8&03v&+GtVtxK%m?A&0DKbq}J2!PmrNnKzIhO83S3&B5v zVX*#0luL6M7*yA8Qdt`btV?et2WcWek?ZAMS?LI%m&|(b>`@#vFFfLi)JOvXi=qP~ z-RGfhLwQo_x(g7`I&ax^<RWa@`NUhJ^$M8o_GcezxB^D+Ub)Y$&xAvZ*;`wV<$=JS z+4U!Ri=nyV)FHaPWf+0k!kM~KSX+SVj&|Gz$((++JM(35divI?@Zw5H%Dy*Rg@0F- zfA_7o<xV>&@z&C)DUQMuTk*l;pXQ+5TllyQ?L6dCAN8I6K|##8MQsTG%t|y&wCQXV z-AN4F^%f6MN)eweHcIklC=e4VF>|4x2;!j2NFv9yGI9NbcTMvY&d-Z6&TxC6LZF*3 zSdUy*BSL5owtpct2<NGHq0EQ1h^wA81q$_A1l8u#O!f!02)c`<7CxL>gvx5MT^Sr& zgn*+Y)!|W1qPFjxf4zeyagLtSM7~Xf(2F`IRfsvDajB+vcAPpvMIRq=luMn^_1?yQ z9`nBlkK-eGm<OhP<}+%<95AD7%Z|%9A20>RpIpa!z|FFL3$L*+ShVaiB;b6&L47~o z1w1eKww;iOIp%+(Z>*3UuIv{2->a}MIQyQ0&Kt}Dc}G-sbzlxibz3666m!6|fGtCH z*cUv;c9$Dr4%pvWAZ(EcM9%>?gE3rVSWbt!Vh$KZ&;ErGbHEPKi927g9#CxDJhc<+ z0ZC3Sy(gFhmh0r_w_^?%xLM>!CFX#d+Ulwlg1==E_^zNzv~UOIdP-_nJ8wtsRe#(j z;AKj&(4XoG{P?*bew6vI3ld|S|GFS~`<EEu=VnsMhJRfny?*~+7i9MR^2Pq#Ot$Bj z`S9n09Op0JFJ7P$6J7t;Ra9bPnGgQh#h?HE&+5RRdanPhPX52SpZ~}|`CDtm{p&m6 zx7O(IIn5s)#sB=dL5{z;&7VB==M}h$(y{i}mG~!5{i7e`KSMwM5l>zH`|AF6693se z^^e_hD~Dwp9_@Jg2pp8XYDC`j!D=eQj`Onk_gGs`!PO=YLGVT1c@2?H@S5L!aH0M_ zggV!-xs0}hZQOAKs=e)a&PIH;l1)2kIC#bI^xT7&S;LkZST{I+obUV;v_kv6`H@V| z7Fh4v>sZX!3~sTJTI<T2;Prt}BR-ENNK&Z1pBLH$RA0)ex9c^*2kDUI`iVw3m({W$ z;oS)KOjRfJbsGRCHKWD$H-P2|=@)y48^DBa$v#NA0U`p~KfSbXfVr;?ta1({?3sGd zx`7^hWrM60^?xAPIjI+6`vyV#sxb$_;96MI;&7?<R6UqA8I25PAuyP{CRM;d!kw@l zMRmglFzV7@a}w{Tpw-CXmIcnGp>$UZ*j)!SiF?gEsIXu8aW2c9YfTVP*G}iB)dsxN z^Fgk(Z6H#WSastC&Q~n3qMM-VgxHg!>3zzrAZ#Aeec7Y|Mz+f~(B*Xl_eVA#KmR@` z<rv&5eXS41k65Qgo$dqSC>m``wLb72RX(3_ryD39?W<zDjs4lOlk?INPr!GTP1x$9 zF?hG$(SF6_1w4N`pc;B@4A_6D9GJL122+)#X{_sBz%BOR4_Z89z+iUXqgd@FtWl`% zprL;UtMgYn&DuYJildL&_rOn(=HfXU;I{y6Vs;mF-Y-BC%XFn+>LQFi5HHxnvINg_ zK9zpdT?VP%b-eTxs|cq@x$PhI>4?BjIt`1QjD#*<#Av27BjM~VbbPHLBQc}q*r(jP zhM;o~+@7Vlkx*6be;7H*O04xv*u=d}j_96?keTf}K-hq}r=X`TasFg!#qy#fF?5rQ z<M@Ij5u|K7YufBeaNO_PVCU~cSh2CS(X<5;nf;63QtTs$uJ<(3(udCwS@XL(X5`Ki zO?wJ1meIu%8xz)sRIH9CUi-#YvNR<S&uDXGZkHqxpDBiI#_y*P8#Wdvw)LbE#xhP} zmhq{?dZjxfDe9@j!t0{PdpPm(T?OmAm;;W_+O*VAr4rtGrcZ~;Qwa8TH@F0hQi!Bo zhTF+^$%M?Iujd%fB@<Syk)$W~1v_?HZPLOVa5HyvxIXp=3)L#!nr2BRJmbS8e6b#o zZd1VbA)F6rd6r%H6V3;$e*XGI6`mKo!GuwEt{{=HIOCebkNIECQ=5hJxOPzTUtGp| zKv6qle<St<FFXlL;KIJ((%$9=##j&NE$Wk~i1mPtbEA7IgX0Nr1Ez|#nEzeoa94ef zd0=u`9&M~pJh9M`L&t@E!HJo|R<|(+OuhZ({m(fev>~<Z+psT~p~HG(E%pU7m_9#d zfH|N@4)gjJ=}NTSU48St(mQDU?Sr9PA4^f3f88C*^jj$4K%=#H!wu9>Dt?AmAP>bC z_E0c|=Azsnr$wur97IhX>+z|~K|Q<Ro5Msl5*_9Uv|-9ZUV~jMCYc#%K<T4p+oQ`U z-0En{7IYEusyQ#-BQBtv4Li@pMw~}Vb~k9=rJqN(2Qt4!d`m;@Yv^q)+fz~Im-aaa zy%Z#|O{{qJh9qQFRV2v0KM~#76=2?RFcFo^Z(lmICIPL^bMoF`mVmBUG|3wzo<q+b zDxWUFobb8X3RB=9))77v*p+ub2IWpXlQ=aTh5T4ULKAyWqspS7fDGSgbg4{br}IEG z`mrysOLP7V3MfgoHFy(^f^K{;`f3}2XktAZxL!se%iZ1D`j{tX)+h?PFvg*o)#zE9 zdJ?jZ@LkaUo`}5lotFgmrXo7o+HM`Wbadiu)R&+g7tp7=IqulPWc2v_2z7{i4pKjP zm^fl`1Ff-)rM^W`i0l%rVoK-=QBKPYo&JGBq-E$YDYcS^vY(AU&sm#`wtw<|>x1X< zGBe2y2hh}^9V4nGTA}smVLsEf!+dqfAUj;$hqn&Bx>iuUWxfst)SgPvv8qM;Hi%5> z@sQ~3rX7M(oDa}3W2wUL{XJ;C`Hhko#(p#z?@P@j(}zf^BS%(8^r2Vc4!$*0edrAs z7jM|#k2p66JNle>iWXMy)ptJg0<F<|EXheffv#!u99MOiMmH_}8NQvEK|FUJD@_(m zBX^NKo%GqSkwbLwWh6g`WImkQU37L4@vc<YUKQjd>tdJ`hD``E)qCgJsX~2HW&>4d zuDcPbB<EAc<z+;^cR4hBG|z+-9M3!vnsJEy%q#lh%t0ISO8UqH`VaQxGk<kA1>{6- zoUFc4g6DFG6!}#vS-O!mr*dW#kGPRiEREfDChp|%2eH~&SoddEd+}S5izj*g2Ax{a zRZlWvpER4Vq$fGfoR`@+=0RR5$?0#a_8`wx7=6pD@gS>n*9WCwJ>X=pBkeZK0XLxP zJx^}Alf2LMPkHvak(|@ZR(n^wk#6_Bs?UVDlB>9%f1|)0uuiq9?Gvth*7rV1U*$@^ zR&?nV$9llti4vtY><f-R7H40CIiT~Z(zFyjFW5CqZjT3^6Y`jSf6!p(aq`fGPe)}j z2W%f-%|e4YV3+yY=hT=3x}P#`5yyH!HHq(a5tsvh*FJMw7jr-p*$rjcSP%G$^>XKa z%>Njd-?R|8di3s+IEFc(PO5dz&wN0uP3(DxFbAx<bd5?JbHF?+Md}ln1M<`<F?`3q zV5g-x+KZS2+DjbFd{BH6^k>^GXFi<(F_|>YupVz<y|Nr~@RkQ$oZNi(W4bdO*zs%| z=sG~@CKg#?IXj#S$aceC*cPNB%hWcj+JcR7?5b3)V_>VZ82Xvh8k~sL2faR9fwHfa z)6^LYn9BJeQ@U~h7P#d58$?V&X#2RFhp`Dz8og0>cQOHUIjzj@()}<wDP;Io-x!80 zx&8Mi?t_Oeaie*z`d}_E%a)d)3tTM9jB|3jaLy+=o#wVSlm|t$Sj=lfMg;W}Qv)sd zIK^O^uB8sz8|M`R)zyKvoytI*R}EaY-v8ipf`C%Ge4S_=Wl#^e9=5hZ8ANZsPh6x{ zfs=Ey58DP+K=93j9rhe5aQkBIfY=%Z*r{jNl4Ym>{(;mhkNg#ZDdwqlm75woc{<-> z$gd0c_>Sy;lBoj^6kik3KMWz{6pNbM$N{j@ZkGDudH|jeg$s|@8pB6J|L@#)?SNm| zUtY@14J^Oy3k_~{2L+Yt39irXkX&ifP!i_>TNXu!ExKLdyx*pUkE~8`U8vnQe7ip= z2&|6{$&3WqC6AtqZ=zu1Mzb7=w~?^dCTpo-JQ92chS=CHMS-^6<yvOXNDw&BS2P+J z4M8tY_t3pbg=gvN7D^nKA)>dV+($GMPT!+3Nj1I-ljqv1-)3Be)V)bgO}rVffmcEE zf<Y!s6qKxU+?Wd&>4^tL9mO!M<v4t)x)dU^H>a)*uLM&+xpg!fszGaPfbZ5oC7iTR zx%6>cIp~WTik{1=g=SVdLgjW781UUcG{ip-1cUy&oa{}+!$yVELxsG=$z$czFW3YL zo5jkFuek+@;E?=vKT3p&9D>@{P+Ef6{Ep>Tps^gGzm4*3t_eZZZgX_Jyj_)GVk~HW z`BaUtYnyvsnyOAbHkZof$GLqYC#X-pQ`aDPdpiBx(=`csoIy}PsZG4-A28=})h2cv z_O*G`r$x9tovBVY(IQ^?+?DIXx<HA+Zc)muTEx)6YMIGbnuPwDrHCkdO+qy_z1F=2 z>j;f%%00O>2%2OY?yzWeLaeL(QYnu*aeij9a5v_F6ZI+ADX~9z&ce~wpjC}9`DjE= zV?Cf?<{nRToDax%wS{5`b3m4x)w>xn2i(8?@$1)E4;ZdNs8C@3Cpzz3V~i{QKZ2kX z%mMp$rreFj9B}8(AS)s43*OaaW9*7K;JZ_;<TuO#_sX5#Q1lRp4D*L!r!fDkqu!4W z;QHnkqg}rd5Ygd^>MNK7uD-$ZjswpN-ki*Ft`u`X2X~$f9_$Mqk6u^24|BkQ(Ot{~ zm;?UFfqva!{aOCq4*WgO`One|I(Gc=fq$z{*+}u%{p&DhW7X)HfA_flr(XPzxY)m^ z%lUb%bNBu@;{Se5=JBh~|GPQaZ=NW6|8M*B`*n|BzX1FX^zY#xNA$nDfB!0Xj?(8$ z>#sxp+dRnqUfk`!SN#3I>;EcF>+$HHt$*rczd{A^f&Z(%;yk~jzpl=ou3d(|cyRyd z+D~!(`oElJ=HMS^tNrVw{}JDt`ZvF5{nzllUkC9YmhH+MS}!O+WIhJwCp{OxruD;? zUDs(oCG^4Zmk*0Iu<L1b!whTrnQj<(yVOY^&;h-)b{}rw{J@h5O6k*?_hC+8oojN@ zeVBEU=6$*MJ{-^AagtlI9dxvXDDJqn!TJ1W@olQDFtIssGlP2zIPTG{lY7$)E$K3i z!Y$2kq4C{nshVaua$R8aarb7p+9oowV%`K?uWoi7X=ntoJ;x?$HJZSVhJHvYr~&+T zk3t~_3HOB~EU$hdA?C;1W$7NQ6O6myF^~Pl`JZmdoaStV_LS3yweWWO0yY~nChNiZ zemupCqmAI(<n!T9NdvgrbKCD5LlCx>GO^?tg5lCCc`B^e%Vf(G(4(yb=WW*rnudBv z5kGKuvb7O<_{2BJ#ytS>YPy8^qHZwM>yV<F>jJ}2{u8W=U9i_^OANhO2b?g`mlN=6 zh9W<7Dx3WtSf+ei7?OSjFNJS#-LiQEUgi2)=Ouf9r~kvcwufC%RCGgX7k@V-ZSYFB zmFt4ci|w*=hx_2g<;}@gN=M-Inrc>g%=6rOsB4Zhyo469+)`JCG0>(>lVK1TgLU7v z!VPN2K=*8fR-xuN1hnqVsed;GGG!a9r5)#i{p%U@u3`}$h{QXulUxE?img&dct69T zO>Y@i&wU2%msjkh-M>IPEsKa#!#8+1pAdd^CnYg#Z0CBnh=!nqA}W5?HAK8Y=>Sb2 z6CpPImEB#5g)p-(l(&^&CDyaW3@@9r60NSM8}9mY5$!BY9z9{Q#OjAW8JjswiSBdf zyDWtri0sq|J~vh`!Y8swe?ZcY;D}o+a=(6x@ZmFZD!p-vIHWYRW2a;w(Py+G=}r+w z#O2zZZ&!~Z5=Pd<C>@9;3N+SISji_4J)w#2mu!-VyAev|A$-Zi%ivR21^AN*@8C0A zMaGf|ss8Z{FV$3n{tK68tauu6qnLhHW+jz~YhLTnmy}9a6u2>L(Mly0&RFVNai<b) zy0fE|xMom0Q%zuB@T$&*0|{j*MB76i#Y9s)CxrKTPvBHC(Uo1Q)PQ}#g)AbK?+ua( zA06tKVwnG}8dN;-lLwj^Q#_)@dcYmqn8F2dK4AH)T`QiENyMbal`GVk1NI!)UC@p7 zfN#Wv=%_IVysm9#?T7haL^I=nBCbojE8J2s2Q0a^E>aEa0q+Gye=Ws&Kv5~-Asx&C zlUGb!_Fz5W0Gm^q8Rme)Pl~m=F$e6u*TZ)a*F@JvzMsC}#taJEI?Mrm8U!m8F$Z*C z6D%i*eZg&^yPuxMdcX!UeQFif0}h6rE#}1>kX_WLxPZ8ezP^3Bt3RR&rE^RLd%zu3 zJ83dJB3O!yBLb6{_ZOof)5mJ&?+ehC8z^3;Di7K1)(eSzaUG@ZU~YQw{yK7WwmvZ; zdmT;vP|NZ;lY<hKwx~YdorS7ZsKXXf1{(PxtD?<+1yLCFuekfCqh;f>O`Ojzp~LBM zY9gJN(3&l#ixuB4BIhk8JF;8Oqc7&ik~pc;@SNV{RRUL1kvIFk)oIO1=#A~!(&HkD z=<60MS))UVNKZ8~)=M`LxiW<_^<GRwg$}pc*m@I@Wn!%OqUkxb#w1puX*?E{%kNoq zkBCR#RQfxGXU`&@iUm9SE3s&O^6J63xLDN0bkF?Eu`}o`WK3M=ibkJT9el%2Mx%%e zN5uCBo<-JEYL@!b7m!Z<<gTv43?!$sEH2b~4NYrzQE=YBhQ9kftV?CfKywIBqan{D z8omj`9okn>Eyem5js1n_N7<uq`K-lAE!`nM{7nJ+RyI{ttD28k=kk_%-Sd$U&1$9g zt~``J)0|=8QH-wjX>{0d)}ZDw`up1oNVI=tPqL>Qi30P4+h|zo5aThYg@DLfq(1Xt zJAX|ba`aoE)ah+Nr-B@g#z(fHMyHvh8#eZ!-Av`xuiGA>Jwxq$jF0=#Xt!7$W$|NF z_Eb0gx!nL#No{=BJ2Qahoc-K{cMPI;f!n3j*N>sH=%N6(J+F{SGTWo_g()N-uDe_G z#SAK$Wqn>*_z?*$`^qSkd_wa4S$t38J|X{)(xDa4SCB{U^5z;VPV(k=ZC;CECGxR+ z<N?WoeI)0#qoMtV=49L|&%qS0BP6luj%n^jYx1#E^^3U8)@0j4d84qq4XJr_AU{>z zjui1_iD{X3B*o}Y+MLaBAss&y8ZlpUCs(FN;!UxR?@NEc;yYGP(tX9Bi*}1Axy9l~ zl#`PuDeU{+&uiY3Y>GReJj~`rX71fGJ#-7}0I8;^Kn8EeWkmV*g$GIX=ytg*=6^en zMZA7}*MltEe@N(gvIn{BuewH0)q|AlTBNnWb3#VvCIwBq+{ook)1^t+7ffv}M=gLk z;7fUasr5J?kW%2~k`U&APeU(z^I|>VsYR*?5_7=dwNokDcuvTtr<2BGcuvT-aJpe% ztOu;S5TJ4c>j8D1bWL+%{?~MG^7%W=10Uwu3YTGD@WDm@ZC;oIZVvBxLyI|J6wUT+ zp_l_+TQ?d$a>{|cYvyhlf;r$`5zTj*nEypExu1x`HSm+9^%CZQ8Cj&uKFk5H6Y;b^ z=LN4-RNQNaIbfD>znCKCfY$mx-`KDo@X=9M;b_bOzieu);^Pm3*V9MC+js(SUah?U zI!!-#?9>u>?zA^BgddYS@9qIRI%LzF{hZ-xGFJ^A5&{OD7gg=j9N?mL=E2W-4lw(v z+ft{-9zJl<T++U53x*mb$FoXnD9d?N)jM<;?tSFueRuf~q;~`#p*~{{UDx8%oGZ=1 zf}1(6>V_E@Gt9>X<NU$<HL+>WcbI}Jzuwf!n*HEkk{hyzFoqqu2T;)k1Bf|iBFm_! z2h{6z1MgYt0q3Ia?FXuQaJ%X0kjoW42)r`(EbO`-yi>cD^B_VC=+}OX=vY>VRZ&km zi(0f`qG5yFfRiSK)CZe=3Q>oPpC&K$xTph#@<K>%5}wB~PBnaHzY6f^k9r?VQh{%; zb*(l8Xu#A5)*mADCcypu!fb(@B`i<XUbQ6-18>nLT|?Ew(6(aX-xh5J{-)}#ynF}Y z!p8H9vg?lmxolR_tn2~wBA@N8A9}&w9gclwZXOUvnh9NJb_eNJk8T<{%ma@#u9XUP z1Af^x9w+sEK=H*lCQ*oh-~)waE-#{?)gb(7+C(%A9Qsh6q7VfM73w(x@sTioGp36x zBMLfl%S>2WW5CeqjGOYCB;ZL>zgB3J4)o6izljTE!un}OKh9T~a6Ih9AZvCOT$Y>( zxT2B`32(Lq#-ME2ByocxEc-f8t%gtMkKF{D5}EoBx}`7?ZM59yPyv~|+z(eQtKjLG zAFW5!8eww6<iLR^Eg)ZU|A={X3&e3h%4Z`xf$x?g#i`VJP*GrJJo|x}py3fTzN;%h znDWulC_NM-Jbs+9VO%Rs+<umv89OXPY__CZVjY(u65M(c18>O@gD)7JpW7-CuRhP8 zI?1R)oWAYt_v)k?@rf<1#voFI__%>mWqeSRD8He~sco-CFvv2bIpeA!bGb}sN{bkr zG?jGj&?b6bUSJIt&>;><y!ldaR-3r*TSyu{(IWO5Bhj1#TEw&MmLSTGIbZ^ZbO)Oj zk;HKERQiY}k$CW$^G-WWV)tCoQ=N7^LQ1U6aqU(O!d*h#$O`*{zm;2#n_&)EmV5g3 zI?VrGxgT4li+NyVO0Ezu)&r_Wo-T64`G9OAQ+gek1I}#C{~&-l;BB$Oy^k>m{5bZp zQVw&#k+-@NikSa}YmX^yz*QtJa;GllfRz4BYOgT|tksCkp~4(+yW25dC(HpAh!sP? z9MFwRm%<fuKv@?ic}naDem@>~d<b(uRvu#VDb@oX9(tVMg7X1)n&#c;#CkxE%zMYz zVqY-J&g0RwSPz(S*j7IRbHJB|PvWcp<beO-9n$~bH~4#A@K;Xvvp4_uzG08iKYIJe z{(krQUvog+|A#o>=b=B&@V}p9o&443f6TG|=<~_ne}ey!`kSN3KL6W3{eJ!A$lt#H z_t)ir6~{U`{;MJXP+#}Ip@091S5f*TKKqxy;lFwv&#Mp|{=MSAv45}s+4`qG^$X_V z1OM0kdl>2X>+1a7#s60-pYp*SclcxuJpupq&T$(}24KVG_U}2{`r+t%>X{wYy<n;C zB()hw^b~RC_k8c~0>7%&@jq%i;M{kiO=nhh0K2qZq4DwqsNS^sskg-gpmpha7^`s~ zEF0@Xg>a7D_YM=ob+v7fJbr^VkZ6TX#U1(f6gXFKtV%r{^FZaRk1l&ZX$B3Q*3s8* zo8iQng*D@LEfD(c{<G4hCfJZ`F~YS5`)G4L1vY47UvaeMkfv!9eCZQu=YQJ>Y+`4| zB%d_meCfAQ&D$DbyTf78Qaut!51wmzCxLZ^`~8To?oD97x`ZW290BjEjkj6{>)}Xt z(5N=9u3M{5YKkF{v3h$YuLVH>Rc!qikp|da)+fWC)Cyd=GVkltdtmpgg&5h-U0@rO z)2cCz{m!#Df>wWSf$jcwF%ElLz}n{%b&*gT9JBcR_A7laHg2507<6C&D(6#8221)N zhM!3&+W#TkKW6s2%U}Tf%yr}RRE9yKT;qFD+YkiNUQ)4g8wFQgk)a))FM+@RM?w1F zICSjUSGv(3>-<cNTevzWz%jeOdolej*aufu7;KmWew&(8i57Fvk<sVoJvk3@qK6&B zJQpDEdG*VNdrKg*IA$+I`voTNpZd-rML}!{bDj&zq$I3Uy4M+Yt|F$_?C6beT}4ET zL`dtsrz0e^g0|a-(h>Pwt;a6OG7;nABpMB3CJ6J?^b#Klg3om4>$VJi;)ZR1=O#uI zLLydVwT^@dF`c@Ng6oL|u`;8`VspWcxc*f_u(HL2IJUqlxw^oQNIP<?*UK-M==x+O zXE7Z{7|i*sPm+!x6wNe;a&TS!_DU5+Y$Or&`m4s-_-JC&DHMEgGLA@&sg7>&OeB1} zoGW?@k_pxnV?p1|sYGW>R|ju?D&g^|J?PW*R3hZ-P~8|?8u99BonK0J8nIovK=E)u z8X=?Yoq6bcDzU%I<osv~*8APc5)Fh@;?`w~D_gdw5=y*0;<{Y8>Ro!1@;!y1`MxrB zu_}e|7Qe}|e_smmaLp~vLo;}OaC?Y(bwV=niH~|kcyBUc6tHyI8~cGNo^fVG;mYpq zUUz;|GEtG?<-mY_!8hx@Qukp!V8+25)(Ol3RT=grZN(hWdDC1GKjwfJZn#YuVg6T< zmA=gz*O7O~gaiA6yV%ogxiJT1O-McviaB8V*qRS#@SG6JXFe{o!|?>s!+doo_66r{ zvNt}Aa{&`l{VCFMjSdQLaL02(w0aI*6UH1cnKPTc6mvj{=BsWmu^y1_N?q?&%mLL^ zjit0O2lV~AA)g0x!0VHb)HW&Cp%937TjN=S_%E`rjl5Th<WBf0zIa`R?rNUZJ*0XY zy%%%KHq<Rb2hQFR9n`yl<QGz{vqcL~V6K?;<DCU4g{9_3{GNQ|CD@;I>ezMEcb0m9 zPCEx3@_*}OVV;HJLe~{14`rZBnx89q*Iz-rLN0Swr!S**;pdf^TrQ)%+al6;a$iQe zr>tH-ewvPWrBl~64PQVmmp*KNbT}0;+Z&SwpHtCgKhp{vEQT)0=C+BIrJ%IDCWXnA z6qMrn{aIvaGU8hDab|K#L~dl|kjn8SJg=kGn%t9&6yh#+EIvMq3<v7d*vI0~e8<*Z z%*k<xZ8!;PhvE?b<B5!xM{($5@AvZcU(ceiYnU!KHKZV#gye67$Fk5pbx!eBH?z>0 z`Sp**Hl?G=H_^u?=h6{tH(3|ZatS4h7=?RJUO}~&x6iRn<)ibdO!Ngfr*D%{tejd| zF}hN(H{|Ag1F1MZA)D&UkXX!BwL-^g)Tx`VGCf;`CZ28&;l^C=(-7SY3)2R4_&`Xd z5ZQ>BmuFlRml1luT&W#$yBX0*y^vQac!2cF^@c>~yU~RUK}-o{-AMm(6P?(%9(0#< z=9u602)W;URdkSQ01X<cd5nw>pia{(CU~|rQonSp;6(ZuvhnC$b>j64q_pWhA+YZy zqT)L3zRvnJihCzg&p$DNW>urwF4(<8&&ng0*-p%(KBh8?$ZA$nyN3G<gQz5F=h{4V zT2_f%yOLab*H(!fdg$BJ(yK-$KA=7xPqmK>WL>nO(>5b73biuo9XU$IQ?+OF(Ats- z8uP+6N%rJok8h?=wF4=j>=P1N<v`w`nJAcj>_}!XzbX%qawfOWs^{v<x{*T$b|Y`F zU-!BuwR%U9Cu#AO%cfJsi!2rl2;lJcB3-X-jP~{NB0Xm+POeGzB0aVqOzgVkMJk^r zr+VbQ$SStr^NuB+<bgXDeteh%YJIavvl;OqqjIlG#bOTFnWHv+4s*a_J9UR8{5fw8 zm8XN%J;)h>E1gDI52&ZR8P@l>k)G`ntJ4_V$Zv8oGDX-Ie6?XX%o_WFr-x5;hT&St z-M227%9Yf(`OG0A*o7>KlDo1O>j9Zp?&n{^dcdx86j3j*9&qx?u6Oq_2Q+xE8vFzE zzeu6nEkAw1Vig5@c(ESv%FZKZrkDdpOI)S?fPKN^HG3z;F$XNHo65h3IUv)UQ$OUe zFPN2zcIFx8e?l(FG#_xCau71&!}Ef-)T^#x#T+o@KIPh<dO)FLHy1Q82h3FaR;__K z;9#!#0SU|jja6$Ze)@u&-fR@ke-a87d$&nQoe6<27F7I|>;bUII3#Fccmk%VH`DgL z^MaHYeEXDSJs{)UgO1c}7pR*tI6xup4AH7+x>ML04sIL<)$!wy<#xPcG0Fid7QcDi zYPW^9pyqeg$~N$Qqfk|P;SrEcR()@)Yl(fU3^&ux9Rf<)(7MZ>hhSjU%}cHe7SOl7 zR(q$i1$^7tzdxn%0Ceu5mCPG322OPY>+;q6;iLi2%~kTo;Cx1ccjqZ1Xzu4~h%_^T zh!ZC_cPQ)w8&b;HS6LsNM;?599<vv&2kYz>%`|{a<;+i`Wm?d^cVe%>bxrJN2y|T3 z(1iS*5jlQgnn2y^nxrqM3D3i_LLK?EAmPG6RvF8E5ZA75@1S)Qx+5;_QsTCO+o_Vl z>K%tbhpT}rA;=1PkL8`7S#1gO*2nqzHynkuF@BfwTP`rh!}7U2(-+o04vlPD>jN$u za}P(ldV<NpLi34Xe=Mzh&v%I~9G)1vJMYU01D?gH*>_oyK&#DGQ#El0#^X<2aC43W zEAKrrmR-?sh+F>FHHrjCPp$XdRhtG|+P=NnnsX5%Lr<*f+j|K<X*G=<H@pJ$?_a2% zczg{iqdq=%70m{rn!dEE?reAz{`uy`l^Y;8O!qwIXc74ONGx2>EQW!IS2xWTZb8-c z@w19OcfeS~obPsB8E}r?5_PG$3;6@%9JAT=p!~2#D!h9dE-~=2mO3#IAvsbnr~8-* z`*{7KFEh-9;g`r(X`iixuV=l=-m)FU!WA<Q7Gp8Or+D04L|lg0DsjKqj7gEW_dYJ* zp#?#l<{a3$%UPKiU0h-ic2Xt=oGg<@c~l5etm&Bc4OL>`Wx%1-Y;|HQ)dGD=o+d$| z$ja)frA-{AZ3!%r)gknE<-0V->kw+E-bma!r$elsOaAV0Ux!FeH%@x;Mu%8ARUl|> ztwS9Dp%Sw1oHlXCSvxfVbHH+!xI>2fw1^1p9_iPZ18Npr9AU;B&}y{Rz<ErQ2=!VF zTzg!TD6>$YpKsP6jCuoT9k*!^JGrQvBVyEvJM0Um8MyJh;FI0PSFs;B{7&!rN?h~J z3Z{)))QImFKk#qDzTlyzwy(3;7hD(0;n;yWpr#UkvoYp?+AVoLrkDfrs~O(r#{7?? zrBCB$K46%`JA-1(0rNzM7P~M9EE`U9oWs6gkGX;kJ1_@aZ)&!d1#>`$YcXOKJwRxu zx8AYC9Pm=YmTy1v0XLdm*AT?>f>+Zzw9;V?IMlDD@ly{t)VTD%AL{|BHbv5PVh$L6 Sg;c+dIbg`5&BZFr0sjwfP>1gT literal 0 HcmV?d00001 diff --git a/test/IntegrationTests/refs/testconfig.toml b/test/IntegrationTests/refs/testconfig.toml index cef00ea4..8ae07561 100644 --- a/test/IntegrationTests/refs/testconfig.toml +++ b/test/IntegrationTests/refs/testconfig.toml @@ -126,6 +126,9 @@ reltol1d = 1e-6 variables0d = [ "D_err_norm" ] variables1d = [ "D_err", "D", "Sr", "Ï„", "p", "Ï", "ϵ", "vr" ] +[grhd_bondi_infall_avmda] +variables1d = [ "D", "Sr", "Ï„", ] + [fv_grhd_bondi_accretion] variables1d = [ "D", "Sr", "Ï„", "p", "Ï", "ϵ", "vr" ] -- GitLab