From 74785679cb69f5a2823fa84bca0cb999f190a0b0 Mon Sep 17 00:00:00 2001 From: jeremias Date: Mon, 29 Mar 2021 10:45:02 +0200 Subject: [PATCH] press new --- kalman/graphics/figure2.py | 10 +- kalman/input_files/aorta.yaml | 83 +++++---- kalman/input_files/aorta_C.yaml | 172 ------------------ .../press_course1/images/catheter_curve.png | Bin 0 -> 62889 bytes presentations/press_course1/press.tex | 107 ++++++++++- 5 files changed, 148 insertions(+), 224 deletions(-) delete mode 100755 kalman/input_files/aorta_C.yaml create mode 100755 presentations/press_course1/images/catheter_curve.png diff --git a/kalman/graphics/figure2.py b/kalman/graphics/figure2.py index 897c0cd..8d84cae 100644 --- a/kalman/graphics/figure2.py +++ b/kalman/graphics/figure2.py @@ -72,7 +72,7 @@ def plot_parameters(dat, input_file, deparameterize=False, ref=None): if bnd_c['id'] == bnd_set['id']: ids.append(bnd_c['id']) current_val.append(bnd_set['value'][0]) - labels.append('R_' + str(bnd_c['id'])) + labels.append('$R_' + str(bnd_c['id'])) elif 'dirichlet' in bnd_c['type']: current_val.append(inputfile['boundary_conditions'][1]['parameters']['U']) @@ -105,13 +105,17 @@ def plot_parameters(dat, input_file, deparameterize=False, ref=None): for i in range(dim): - axes.plot(t, theta[:, i] + 1.5*i, '-', color=col_,label=legends_) + + true_level = np.log(true_values[ids[i]]/current_val[i])/np.log(2) + rec_value = np.round(2**theta[-1, i]*current_val[i],1) + cur_key = ids[i] + + axes.plot(t, theta[:, i] + 1.5*i, '-', color=col_,label=legends_ + '= ' + str(rec_value) + '/' + str(true_values[cur_key]) + '$') axes.fill_between(t, theta[:, i] + 1.5*i - np.sqrt(P[:, i, i]), theta[:, i] + 1.5*i + np.sqrt(P[:, i, i]), alpha=0.3, color=col_) - true_level = np.log(true_values[ids[i]]/current_val[i])/np.log(2) axes.plot(t,1.5*i + t*0 + true_level , color=col_,ls='--') col_ = next(col) legends_=next(legends) diff --git a/kalman/input_files/aorta.yaml b/kalman/input_files/aorta.yaml index 8181195..8badcf1 100755 --- a/kalman/input_files/aorta.yaml +++ b/kalman/input_files/aorta.yaml @@ -4,9 +4,10 @@ mesh: './meshes/coaortaH1.h5' fluid: density: 1.2 dynamic_viscosity: 0.035 + stokes: False io: - write_path: 'results/aorta_master' + write_path: 'results/aorta' restart: path: '' # './projects/nse_coa3d/results/test_restart2/' time: 0 @@ -25,38 +26,40 @@ boundary_conditions: type: 'dirichlet' value: ['0','0','-U*sin(DOLFIN_PI*t/Th)*(t<=Th) + (ThTh)*(-3.67949466208*U*sin(9*DOLFIN_PI*t/Th)*exp(-t*10)) - - U*sin(DOLFIN_PI*(t-Tc)/Th)*(t>Tc)*(t<= Tc + Th) - + (t<=2*Tc)*(t>Tc+Th)*(-3.67949466208*U*sin(9*DOLFIN_PI*(t-Tc)/Th)*exp(-(t-Tc)*10)) - - U*sin(DOLFIN_PI*(t-2*Tc)/Th)*(t>2*Tc)*(t<= 2*Tc + Th) - + (t<=3*Tc)*(t>2*Tc+Th)*(-3.67949466208*U*sin(9*DOLFIN_PI*(t-2*Tc)/Th)*exp(-(t-2*Tc)*10)) - - U*sin(DOLFIN_PI*(t-3*Tc)/Th)*(t>3*Tc)*(t<= 3*Tc + Th) - + (t<=4*Tc)*(t>3*Tc+Th)*(-3.67949466208*U*sin(9*DOLFIN_PI*(t-3*Tc)/Th)*exp(-(t-3*Tc)*10)) - - U*sin(DOLFIN_PI*(t-4*Tc)/Th)*(t>4*Tc)*(t<= 4*Tc + Th) - + (t<=5*Tc)*(t>4*Tc+Th)*(-3.67949466208*U*sin(9*DOLFIN_PI*(t-4*Tc)/Th)*exp(-(t-4*Tc)*10)) '] - parameters: - U: 100 #100 - Th: 0.35 - Tc: 0.8 - t: 0 - - - id: 3 - type: 'windkessel' - value: [10,0,0] - #value: [10,0.0008,2400] # [R_p,C,R_d] SEPARAR - #p0: [80,1333.223874] - p0: [0,1333.223874] - - - id: 4 - type: 'windkessel' - value: [60,0,0] - #value: [60,0.00034,4200] - #p0: [80,1333.223874] - p0: [0,1333.223874] - - - id: 5 - type: 'windkessel' - value: [220,0,0] - #value: [220,0.00034,11000] - #p0: [80,1333.223874] - p0: [0,1333.223874] - - - id: 6 - type: 'windkessel' - value: [160,0,0] - #value: [160,0.00034,7800] - #p0: [80,1333.223874] - p0: [0,1333.223874] - -timemarching: - velocity_pressure_coupling: 'fractionalstep' # monolithic, fractionalstep - - monolithic: - timescheme: 'gmp' # generalized midpoint, steady FIXME TODO - theta: 1 # 1: Euler, 0.5: implicit midpoint rule (one-legged) - nonlinear: - method: 'constant_extrapolation' # constant_extrapolation, linear_extrapolation, newton, picard, snes - maxit: 20 - init_steps: 30 - use_aitken: 1 # 0: False, 1: Picard only, 2: all - report: 1 # 0: None, 1: residuals, 2: residuals and energy (inflow/driving/forcing via ESSENTIAL Dbcs!) - atol: 1.e-6 # note: dot required!! - rtol: 1.e-16 - stol: 0.0 - - fractionalstep: - scheme: 'CT' # CT, IPCS - coupled_velocity: False # False faster, True needed if robin_bc implicit - robin_bc_velocity_scheme: 'implicit' # explicit, semi-implicit, implicit - transpiration_bc_projection: 'robin' # robin, dirichlet - flux_report_normalize_boundary: 1 - - T: 0.8 # end time - dt: 0.002 - write_dt: 0.04 - checkpoint_dt: 0.04 # <= 0: only last; else value + last - report: 1 # 0: print nothing, 1: print time step and writeout, 2: 1 + flux - -# solver setup -fem: - velocity_space: p1 # p1 p1b/p1+ p2 - pressure_space: p1 # p1 p0/dg0 dg1 - - strain_symmetric: False - convection_skew_symmetric: True # aka Temam term - stabilization: - forced_normal: - enabled: False - boundaries: [6] - gamma: 10 - backflow_boundaries: [3,4,5,6] - streamline_diffusion: - enabled: False - parameter: 'standard' # standard, shakib, codina, klr - length_scale: 'metric' # average, max, metric - parameter_element_constant: True - Cinv: ~ - monolithic: - infsup: 'pspg' # pspg, pressure-stabilization - graddiv: False - consistent: False - pressure_stab_constant: 1. - - fix_pressure: False - fix_pressure_point: [0., 0. , 0.] - -linear_solver: - method: 'lu' - -estimation: - boundary_conditions: - #- - # id: 3 - # type: 'windkessel' - # initial_stddev: 1 - - - id: 4 - type: 'windkessel' - initial_stddev: 1 - #- - # id: 5 - # type: 'windkessel' - # initial_stddev: 1 - #- - # id: 6 - # type: 'windkessel' - # initial_stddev: 1 - #- - # id: 2 - # type: 'dirichlet' - # parameters: 'U' - # initial_stddev: 1 - - measurements: - - - mesh: './meshes/coaortaH3_leo2.0.h5' - #mesh: './meshes/coaortaH1.h5' - fe_degree: 1 - #xdmf_file: 'measurements/aorta_C/Perturbation/Ks12V50/u_all.xdmf' - #file_root: 'measurements/aorta_C/Perturbation/Ks12V50/u{i}.h5' - xdmf_file: 'measurements/aorta_s100/u_all.xdmf' - file_root: 'measurements/aorta_s100/u{i}.h5' - indices: 0 # indices of checkpoints to be processed. 0 == all - velocity_direction: ~ - noise_stddev: 15 # standard deviation of Gaussian noise - - roukf: - particles: 'simplex' # unique or simplex - observation_operator: 'postprocessing' #state or postprocessing - reparameterize: True - ODV_functional: - enable: False - VENC: 172 # 241,172 diff --git a/presentations/press_course1/images/catheter_curve.png b/presentations/press_course1/images/catheter_curve.png new file mode 100755 index 0000000000000000000000000000000000000000..0cacd79523a04fae2abda4a8b31e16548916c949 GIT binary patch literal 62889 zcmeFZWl&sEw=J3k3&AD9J-7vjU;#Ri;O_43t^oox!GgO7cY?dSySuwK+Hg1d&N=tJ zpYPwRy0@#ky1J`bz1P}v%{j)HV-YMXEr#+A|J|!suTZ{!1IfR71&j0Q)f<$z@W3ZV zWbqK-<+X#nnDDF85yE}o4XlZfl+deJ6%j~}dT_vdMB8s_4zFIJ-@iOx5BA~%LtG1d z2MH; zva6Gl|K~CA=Z*Py(0?8=w8%{882=t}ai8kG{MW~79cyp?^C^ZAqCe7qf0^?Ce)<20 ztHFuSSS<~@BBWe9$?2SHe*Rc0DmZh3<0p~hZaNNm?;2OGg4+LEl`tIwkMDl>8tRGO z@&{pQN(I^D*W3@EPu`F4AkIfYNx~U;YiG^_iy2}6E-==d_ugmXPPl5`8oxX&gANd= z9_FF`ieJQYijr{I1D(1|bTgQbi#bzovd6W57-MFJ04p*PWa|sws~3Me`MA3*rl7RM zpSMDHGv?WP;}TiiYnpukdP;Y7LzD=j7{ zl4#52R+V=BJMz(@IS(qw4c?At*?@rR1tqB5NLTQcR3XHyvgO}62QH@v%FoAR*ay| zTKd@zNABK|2B&ve9J_$TWWt_$Oa1iZhAXf)kJT&)=_osV?kphk(K9bI?;H9t9%={5 zSz)syzxW_ri!7#KqEWp(K}^a`DqE5)m8mGpE{o9q;Q-_PM)Pnv3l*` zSMRkflCs+LhFYEAZ?J{ScEWd};bjZ%Oe;ZeK6H%p z`?9@{*NbRto*>O%5jZk38x*B$8`rY1sjcCSwAxdoES$)%A8HMt3)gd1 z58S28QnES{jK>x6*&DHCIl;Zl;K8@tM7hguzr`dF%ddKQ z&O%Cpo*Yy*)4$uafAxAYQAE6h!?cc;C5wd5dUVEKvnWc((K6vu`#5q6&aMli zH)r1q(*wZva8B;97~_|buF7nA0e|a^$GEHSgf36Jc|~`W15G{emEwApAp&t5bw3SO zk_I+@QPbl1u4qoBDX`P;AWyK{N6k_V>`Z;Vgfc+DsKB=01m_ByntNR~W*Mzaga7c(%cbyB>`s&$L+R)oz z@Ddx*sWYUb5c~nR-Vgzmh&QzPSz4s^4ezHDe~g5Ek{s`dgNm|<_Ew?oF{dpX)oLe# zn$mj|3Z3Q%7@Q5c?Xc7JzY7FV!_PKcRm41K0s;aH$68ZNNh^zGj46B?>(wfB8MZZ} zy_DqSA%k4dC8CcRpXgdSIkCL7%+KG>-=pn>ye64mZA%UlryT~xH0BGzlV|NxH#01J zml1Z1kYc?9<&CQj*POX-esWA)e&m*%=3P)l#O&d^#k#QtWpBgfyO5;?Tz0z^} zdU{A{eqwn!VOo~VJ_5sKZXIlM+OIa8f@-~BkjfF(ax)1!6x{#rn(ub@?T&82MsZV1 z&GUF5$ykWn+Q%w3j51EA5GTU%;Mr=~V*0*2oNfFB06MS+yWLZlA|CaqP zeTXte+=k?>*==cxZ{#U^l#DC72YN`5_YX85!C-DTU8?E)p!nr?1-w7P1|>|XD`hF( zh>dFYYUS{RnY2@v-|2?dDB~&p0&P^K9MaV&)CErDmw8bdn~(~po0qk8Nb4_3pM!JA z+?&oQQ2IGp!4 zh3#K34)w~W#pkq4XmWx~)KGcgdm-%_mSED=^A6x+V92pT9C~Pd1J1n!+UxkgB%k8w z-_5D4577F8>OTO5-w^IduHzCV^nTv}JjetBZJvloXghXqYHz_*TZgx(y3%2g?Abvm6eT#$U|qP8B_ze(0X77}nO-X0`It=(CzN z+eO*_d=Nfq&0iQ;Qo&$w)!Zmv9*HVgXxg_`p9Ka9?=YuW*TA6wxaKhO-pt( zJC=e07LgU-LL)@a_E8d`m+D<@h|>?%B3g-h4R1Ad6qzEek9k9XsGe_u9(YdR^$!D= z>nL?t);hkY1j`}>m=lv5Ug{71z7$PfPt7GBD`83tbl`8g^BJy+6}dGU5I&mrLcQZ) z*1OS>sDLK1O3d_B6E{FgOG;9^{8m*qT|-aEFIdPU9&}f7CBP~g*fO@K>jtLm5f~4j zv$8W~++?gG4}WBm6x5q*mNYP%!XB$}Gj(-1t1#J6`8U_z-+kYB+__Lw zZBwqkdM8yzaME+8vD#hJO)5F~;i^)Lpf9WB`UpZijAk`Je=Lo{uJr(!ygFP|n*35q zdGe<~z&w5@{E_f}zNCBCv}7MTZ;~XlD%YxgU4g&K8@$k2YOdUpXZQvuy-Xq=9)A~= zjU$ivAa7ZkmPNNhlO&`tY3B4Q)(T2>Yfe6G5De~1+UoR{CN?+L zxw)+}mKM9VAec`PSj?|4@qUZ3O$=boaR-QGwC%$m~qBo}-zW$9E%@Bk49O!KUn z^s=Xuf~ zBp1<>ODrGW7gk5(XxbhGlEAZz&|WJ62xBP> zN{F$I9Ii>4<+($90&m$qP|{AL+KyG7#4fh^PFgh34-*4po9bv6Ok?u@q!DHq|L1Ur ze>fbLfPkP{6a>1{8qp@_`#O}FKK^&}L(@(z1_eHM!EuKH=qNSFM(hzSBtq2mrf75c z^clgXjB9P*u!t}U6}T0kKb;zSgByKt;v$ccYVHr`y(B0#4s14?cY}{T4)N7tyUdpw z)`BiPQdoaUT2Y}Skgh_M@0b|%E}!omp_mq9uMP5~yMJb;y~VwGt^I6|Qc9KZw1NU{ zN?g9>>Tn?COoK^-sW>OxDs}i0*ZYs=so_Y$LL{LWo8dkADDE55`u%%ma|$C!^fRNK zuFTA#cSeYw#?C@rStKla1~jhy-$0BaUdFDcN$X;m= z4JOIZw&T5qOAP`yc@x+cy@_Cjo90Z-qS0A~0i$$HpSivOcW@dGP=mnWzIR>gg%Mn$ zqh5g#JD&b7^t^?~?wB^42fn0j<8a+X3OXIwKOOd27cQtuK}#x0z@O~Zugk~8_^D>m z^LkTy-|tgDo^_%f5F25S>g=2^8KD3#+YpC%EvT;^zzU#YXap!-6)O_u7Bd?8x;N;( zdCk6o+COesoM;ESHWQSa$(q)}416iZNW0qtL7-yjn-hMQ%=A%WhJ870E1B<)n&ip# z_1a$7O4Y`5KF?5-i-e}|-5@AY9>(eq zgof)$Wgzfk|4{Ey8viyz+G3|^KT~Kr6F>bhVyXur$zY|VdR5FK?4|-id~y;gDS&91 z{VXA$2m4vlOd=`+@#ap_!xLus+v>M!Mn;Ke%}0WD9&}e!S5(+E|0O3>S9(jgln+9{ zrsknz(M8CG5+G})`KB>{?i+{W%0Vf@k5YC72lftCwMpvy0- zH-XM2#a`8{4XCNFP7?*!{>Dwql?Yi`E#jrrdCGWOMwwu6w{pF;EkYc%511AuA7Dgo zHrcKT^D-VuhZ#^Cs{cX@792RbjXk@(6giwJndTp$NoyS%q~B#~eBz{i?%O++9eEz< z&K<7LxIQ}2*;C6SdGhE**!=+c(byC>MS5~OYPxzR*5z4?7t~uV%;M `{s#`!wW0 zdPZmnB`it85mEp_IyFNCK02S9zhVUjbdy!7g3PJ++HEH)pm6~;!R?>@X_=@(Q3 zz^Z_hG2Jd=Q+AK2T+rSecXMdOu8r7_Q@Ww^Cp?EZD_@^B!P^Qgo^;+wa6FJfIMMQGrh*>by**~cD+_Lu$0SkCaAh* zfbY`oI3e1#Rq^UZU&Kk?2!pqFW4%ZlU(IZ~=QY_+4>=v>LwGYPG&B18fd}d>{U^^W z>bzwHO7*}4ao~~g5DtzqoGt50^92IQGzZPJAVou zcR<6r83ak$DH4J@{JJ=n(Fj=OP3^4XZJ>%{o4ufUsoe~x65()+*9&v z_|Wh_YO03|E%$5+Uw)r^Z;KR>#^0U2SHyjA$%wE0IXnXaHJ@MieMe_!bnTec?TE<3 z`YA&fEfjJFO;{P1x|^blIISUL?`G*=e0p82TSP1d)br%MRgfU{edBI#iZ^PxFNL{L zWtoAX7|#ohCl}X)`jTEfxuG4%Tzn0hN08eLI<|k7{Zi)+a{kqMXH``3x)~68OiWe& zvcFrz{-U@0%XH3wJudZ;!2ggt+ctQ@*1mpLsqK=zjx5lbERO0wX~p$q(@2(W3(61 z1JA%eSuye>Q+`+XOnI0uw*{xunm^HR){PW@!R^{qOV z<~gZ+zGdY$m>=2LIQB8X#;Vc?P5 zh7#G<*VagR|7f^Gi?Oj^Foac>ew}U(-NaN$j8_>+EohaB7*C#7zxShdQL->gyga}E zaOIA8e5mrS{8#dwcP`B6@l8bz%fZiZ9%3 znxLbo9<3w(>dOFjcpIv8jQme=iiA{!>V+Gr7Eu`~xx{j4Lb z2O6OgBLod@H|K6JPqH1|7PhC#NLwG%LTfW3bi8uemZ1^5h&Z0d8roE`(2F2I%gJ12 z$DB_r6!Wm3E9`hUBW2xjvfOC&9hAEq0_>F-Ib^0SS;m=#PLF*ShJ~hXn%7{Dcy2GXfVTE`#ychbb_8K9t=~YFO$h^^CRwbs72U z7G9ybSI_kN68C22;r1eZux|v{(wp>xZcmia@n}viPB~`^o;VTm<4h~KXvyS!&68Y&HIcuZ*i1F1>j4}^w z7D3&y${r+hfPqBT<7)OpN~~qdLDuM+ z%i4#(h!S39Y*ni8!>=l)8Gc&1A1;6NugoRL&s zYKg2zQ_I-{^WnN#KJg4Q#>@LIWFXqyaQ%_$#rcG1q`H6WWQ}=@oZ@PzS2^ ze%$-9+Eh0jlf-idD(pNo=AJg?4DQ7;)w?HAt<{YjlRODy*Gc;u!uO7hb7rWHzppZ3EA@R zsHlS@&bB!Ylcuc2Hx{olcPC03*P9b9HOy_5ERw}DvrYDPM;0UDD_ zR(Kryke~EtcXn@fHeiMl&py!MNb+LPvpmwFno5YxZhUjo$S~3%%%QJ;RV$2+(YYsV zya(Kc%p85xq85CkkOGcN_6y0hGHFAr2ClMDC7s&#fb`G(n#ft$5A&QLb~u8$*fY># zMDJ$00hud0p6?kLm#?9T$U)RGz{DL#(;4N{#xHn$kwfppy6ZN~FbFQiVk`#AZ?E_m z(7a5%jnEQkuSl^VilMo=+30&r`gFk&o&P?$j^Ey{wC5tK_o`D2RL2?q19sxkF08`8 z+N#I7pu-AZR;*3v_jLaz>^@b0e?``V zW!($e*KCDbxB}A2xmRDs@7ZCJ9hAIY`@X_)8vx$gX|B*p!+jYTRW-HD90 zuKX{KCk%Cl6-%eTa^+zJ1G8hpa-&Q9rK*BX`To&o48PaUDM~lU^onip8UPd@FbW&B zkO{+Kd?}wqyk7?Bv>j;~pblz_9-+5)uq=L=*2@i*OyDM=@bH1Uo)8@7+yvH0rT;E( zK-E?P>KJj}SJmdAw{s0|o_l&WE|$_<$6J0pij-H_ye%U=9u(3xH7#2%J12Uy)c%5E zLVA1Tv?L96xwwUAI-2(Wvi*)m;pg))9wpNNGG*Tcbrfv1ql5envd5#OKT%#CGA;&& zS{iO@M)YT{-%sBH`ue`T=m>^9O@meL5E`x4h!l!M{bW{bt<-8}UYG7M+LeU6Lorje z9DMvLI#i@PmalgNB@EnY8ZZh7;wmojUMiWQBf(-WaftL zX^Vu&vTmyVQJJa(Yw!k90TR$BhYApkIBde=$k}Jiii--?eV7FBHLf8K}L-oSFG(To&7#4 z1_8}(PS zMJ7H`RWiK}$g}0bXj6V~gwI0)302d#db~T90)wP5zz-)_D`j8ho^qH?fa1u= z2>Ag(4uq>PvBBm9zQOE}s84ZMo2ZeljShabM|kJB=w^ktG>Pj%o>oVc9Vcq>M`10St5phcFJYC8)pxF%bL#gJ+9GX!Rd)PfKunrwRJ(L`k31@^v1~TYELeO z+fH~!Mk}l%7D2|4pPzj1)A(*^129R$l()JeEuv>&oEnR;nm++W4Vv&Xoq(kK$gHd6tx1P%6ACXje?v;_}GP@SsEkML=GBsxRK=H=g+O9hc1*nMyhZY|{!=wvv)~ zeC&9FK8!cHIjD`&LFD)lgq2@NG7`i=Sm5YzZgiUu9)~92=$>M*B(*f*o>Y`3!Jr%A zdaO!$;b+8WfTHxEJ*SC{-)l{Om;?||Xv=ZoWfmpU#X$X_3X$%Hon+Apx^^Hu91n-^ z>A1ls>%GIXsqLu#y>kqU^GpvJ1|XLLjkQ!bCN{ND6Ryx*&9zda)cJ0(*BQXvi{il=6TCH32XOOYz z#T0)vzLm;4k0rlqmOqB87L5@pO_QhRSxXJzeoq{3>x5KaM?&Aq6k%ZSuyM^iyJjw2 zh~$S%e0sjj`SXCNto4{g7D>`yqCU8;0xKmY#qc}Gc%T!S($*55Mh(PQ=q)tiu84!M zg(!$15^buvs!Ma$rsavy&i7&Z!)BLaNa1$MG=CPM?IvdVl@ckPq?PLkp?F4BGQrYR zNt21g)sr3SPMAAhJPu>m0$ah;3ae|ARfRlFwoB8Ka_&!GvRshz&=hOT=l67qzm0Tt zGZ=y|s;MJ>aKQK~@JaKwwV^vqciLXthA0R$^KAcnV5&J?(PDP49$9b z0Ouwo-~#tWCRonx$^1eALl@l_zNpPR+@ASmee*XpWlIwCQi?EfKZ$DILN5lH>wH{M zAkP;Vg#o4Q0VE6LuJ$USo{hxb&#lI1H2TaDse2BhaiGyySNFdzzZ$X^Sfxw0K1qE< zwOI+4bzRHBkw_pON&Gi)pdzYqnooU;>rpIJ{2m7o$lsK5EwS#3X>ZxLrr5uD{nxbn zz#~dh2uO~gaM&Q>o$bZl+ zyAIy#hqREC)nFRIEPs;o(#Q749sX^QMka8;=V=nw(?)9|d!6m;S|((7%^@O5Lb5fD z4&mfv)o{F+fu9WTEfhoFMy<78j>_RlyqA(nCLO@(GUa2TW3DDlg*UzPme$n~8^wyx z`h(p?Anp7pzBNOcgi6T$!{r&(&uTo%JVa0bPy0e<$(np|`GP)i1>g2>G7&>9KD#6? z8{&hR68ak_X@sx(T^HGPgz*>8Utvk$T?%|FjYH%5MlS{~bg7;>eb}8}XWg?*Yf3^R z#+#|q7KTc5rleSQFZcgEJ;}G0j?P7C*otKPVx2wI37GDwe={4{J3qM4CZ z$T|u**Pny{38+F&qFx9SPSBenopHauj~3!9a~UKpmykG%i&53OHv$AMn{zpna}N)_w9FZfz&72K zuo`O?LrXN}aPwf4f&V~b!|QUA!T@k5MhSy~c;>R1A7f*(Kx9Xx%o%IzuDkW^{MS7L zgjdln0Poq|eywtM#-`mst|X zba{u3x3t{3eF(+)PP{vs9szy0$=0Z~uC&ss3Sp2!^iKGp^uSc)5n~ga8J}A-b)|#V z2^yAKP_@d~AB|6XvIx|s2DC9Dv-JQR?cX`|j=`kC_F+?A-wzn8r^<*_BRTValZA8q zSf8CD;x7y_6(DzA>rD46} zwNntdIH0JgNS#f3%NYG-=!*b?_^9|a_OGy{AW}at>x(&xg7tx&rbye1!vHZ7byZrB z8a7g`ko}DSzxmV7eb&8iy-L-p&aK-5yChpVv;m5rR@jT7#(U zobR)C2bT=P$Rhw;O%DQ%syW{{ym|OM@LCxM__@`dyXI%zoqU(^zUq&3Gv;9B)W&R+ z(X{?Wjnwf%Jw`>xf}f(zqPhz@jdlXReBmn@DxPY?zt!FH1sgGTn`7@B+{)wq>pXt9 zwJm2L;4uFv0d#It^!eKZx=b|}IAKU?3q#b)P&+3qn6p-whUYvE$b~zy&q&u3@r@+f z{QGQwYVrj)Yvg6>bbl%)Ttcm;tLf}zt685|hP}qZGz#;!C8I_@krye!qXyI`^Mr@p zRxAv;p8Z_TD(ln=O@3iX$)2*=qC4J=h{ktQmZ%fV)~r$!?i{Z%Sqw{kyd7G z?^09%=(OR+^JSM^cdBQa!120be@A7lVo-oRpOq%JhkmqFZC13RB+q)J`t(UQj#4}! z-DB5n7Xw)B$9DJJuXET!{ZZlH@NnLCFGL7*q{nYeL#x;;gyt@186=Y-2v0ojI01xK z(<5ZpeB`94!XDs4tW+ zQNB=-wi+Hj53p+J7Yo+?FfD$v{+Mee@RQ<@D&vx#M7Luld0+rtM#R?i+OCQHob-g> z!;cr?VXY^BehzcwPk4;ezLsrpB&>1pJD8noIhV@bwRV3a9>YF#KKB?eUu1yANdvM$ zrK7z)P&YL-HDADZWXjNHnFX^MZ zS(T;y`)vc# zm;Q9iuMWy&6 z4&;OTs0H(`V5AI+CcAq04XagYYKmDf1g$YL-NCSaCgs*7tqm9EB;)ktB^}6K#Uo{tlr6*%h0$W3C0H0$ylR&TL14%ucch^zWIknkdv#0jBg z`8G#Srz%326zStAnLegBJrY}n#3^b_yaQ_B%?{y402PKD@p~^pVj!0Yg^|=`cOLoM zTL^fS5w<%6X@#GT5iq!rlS@>b+`)2|nyz!8EjBN9uJ5g?I!?XT!YS?{5e)0KDe?)N zjkv6^?%4jXEHm)TdzYnrwH(mJn+JC7@8nqc-1jF(a>V)aQnXaUPO|Dj9d$Po%hAF4rC3z_LXCaGMtncQ^)#pa)sR~Q!&v4ohFZ^uf1Oagfv$Wm%^S!mAiiP zS0|bIFuUD)|07aNt0}sVmEg>2%g?G&6OM0n1V3{mWY%UVDqsPj{~-w!NAYGoV5ulm zbGM!KhviU=7ulQ>d3|eN-$X3P@Oa|TrLrg;Yi?Ff-Hie(E0x~7#to&c_6mduxBzJm zAE8hqwM=0+wOsV~@S6NNCwub(QYE-#_AKBOd7`o{7qA!)yBj5xyJwHy zmkt{i^4>~Ev9;%YLy`(uoOx`O;b40nh?61Fg{ZABO>$jhV-W{GugPl*_*_C?(p=u3 z#qutab2B>q!(Tw7LT|EcmSqn|C8Ry22Fl*m(*_d)HUAWVE`HWqe(CAy48b&920CIh z?Z5nxM$tX`(SU+Kd|R{5P0qkS5p=-;&7o017lMuC6=Gm6An*kQ32!#v50kzlA(PfQ z3$(rSjju~gXsvFd**r~RRms^X2$ALsILPVaaYMioCbo~ac7%c?fpZ|E@Job-ci;20hnDqw zGp$w>u66i3(2oP0Kh&o4hBT!et}+9KfU{5jV-x0l?bn*)?7T8M9xO*k3iLP>IBP2y zX@*zk&F|RV=(i!6HOm$7ztLZJ3Cb=EBt8G>k6F21o=P&Sksza|FivGj*U- zKH`nJ(^$El2z#*-vGuS(?~=&+0Ngiy?#Jgx_~k?QEOo8pLzKg8zFdCUT1P@ciin)ACbYDwhun)GUCtC+2wV$q87;Q z2yPwcw+AgOyqT_@*HjCsiZ|fCdsRS)7&>-gmCBpyV&<|IPWq{k6qR`H=~_3 zj>EjItxfsu)7EG6YJGXBU1d3ggC%~up@^p}mTpOo?LHu!?c9hrZ*phkhse)TFRYA( zncr$!Ri&&j78P<#X31E#cLdD;{YVvQWQtBw4hzGJ+#UB$d+t4YI0B-~8R;V-Co)6O z9VSM@_#AFBDULe6a%= z(JxyRfPbMod-35@#e{PcuP#}5Z%1gE*~sp9C!Y82utmKQ-8M? zJulV#=diHKIEvD0#q@w2$yv-Wv$HFvAyfq4hh%e_xkqScx_z!_5xo?tsM47G8n*L9 zbe;V-sNmA&YzgJg3x5$_diD2}!*YFe^J*DKzf`M7kKycU?zG#7+tPi{Ni?hcIH#9lJu;rIbgT~Y12FBW{3KAzu!NluulhBO)Y)< z=H@ACbmc-*_B5h)zdO)>lK{xaW}2Uo#m}_`caq({Lnrr71tiba%k~?mBc@C@Ha52E z(HFZz7Q7h%VeP?E`JY{AYx*XIond(taUa4v51eM$|F?4ZFv6e$Q|7uUeA$WWP==s| zrwIqqOvzwx!6A zQm*tD1W7!0SivKix(-Mx@(r<5CAhlAFa4HQ)Me4L`TR?h|xNn3h8-F_YiAD zm(oB`0hrcvoOd_l>%Oo#(x)Zw3UwH5rgBs-q>mAn31da2>M#HTn!=Ym(aRUC*cm2$ zu>xFt0#qXTaWp2d5kYr=n!@6?dB|)1!ma|MgYZNUW*>q3cM7)qb9hC@-Oykt{$|^m zGV2&Ar<~5U69JwlI!=_1=Xn61GvsK#@?#}uQjm)s8@#O}BkJ*XG-5Cqg+BSP%RX>ORo1eadW;Gtm@@BX<{e!o!<8Bl%g@2KD-n%eh>Y@)ce2B&Qt&ip8Wwlrno70EteK19pGI>Kz{DC$?JeE&3*HrHscTm2C*<|X;mGuL&nja;6 zjT_0-&jJYw2`^%DuYe<@wVQ~;>{^3`r^4c4A=su(8vmh!ebK*g2nd5w;7XDgf0P>@ zAYAO*Iy=KlKlCH9MR!faAgC1aQeoRq5tEd9N+vij7E!4VefvS~Zw(Uoas_pbp;ti# zuw1m&!_8?t6bX1O>{HY^wwxzkovO#otnsJiuLL7?yCJ0Nt*W31XB zLzqi6YU7M7EsGg~i&Hb7x{|p3=SJm-d)D5_(fRuNRys1v!B)di4>qVh2GUSqNRF2f zpBFUL+;-QCCq};d3Xh!T@#@7i5=>va?HxhPZ&+_J%cO>>7MoQ)UZBI)RjeLztP-6^ z^F8?FYm^>= zxA}{8@y1^EFbvmjG^>;89V(@D5Mz-mq5b~crvNakFO=}1o-ec=E=%!brI~*4-L#J1 zac!?J$tk7Hd-uWFt;(qT0Ua= zEf*;)1_qby@fGI~9_xp}lU#gLGZj?kW^J*D+8X%5TH+tikeq~S+P4c5wa9eg_vxYS zT@BBAx)JRORrb`6z?m+)=Ty~wD6P7ZC(`pGL4hU+51NE`%Lm$7VE(v}KoR5>*caq1PJ#?!IynQLHIrIPS7lTAFf`haGl8XCcz(jNTjF~Qn0 z11{UDWq`ZWVb*#t68gfK*b(0B>l`ck^^#x(mclFdt70hm(lWt%>T_xvlFnmz>?#B_ z5aX0H_Nc&&T)cr$Hmjboze|sfO&s8nSXWqKXB{3{^b|1q{66S(j@WdeoD826iar_M zV6fII!BQh-#vv4YUS{^ROU7y?-zeNqu_pQb{m*>aOg{!FfKR`JM*Uyr&g!MrbDwwh zZ|;Ml(I-DMb{R1+ELwkfj0#Rck@JuKS$Ny}%$W>{D=c3$-P^l6%HiWG=Y|}AkhI^< zaPnZoM^Wq-6}I*5SW8}&&JL>8;+I|MqVceMQD& z1yPlx2>P_RE3U*Od*6iim$ph-Ak|N_q;JGKD$7}dwa@bXrbByK!Y+}o*KlIN2diP z#^?$c1MUNbCl0P|Zh8k(1(l3!|GQq5rw{52Qma;$vDNypUUSTkZqv^{2&RzL*dB9Z z0B3*jRv~C#7;_0JksO?ylq?D#ch? zUla2WqBL9Jjdrx_#i#T-O&|?Eivxqr6h)D9l_|*7& zU46ZgC%wX(fZ_WP3Rb4i3w@r7EKzHV$^F5&Y!VFoZUdmKJ{p}ipVWqEr-D4usWKg% z6C@mFZAGBR8gsa;7==l3C0NJcMRNG-B&PB8?BwYxlhy=80;r8@F}%V6z?JIaNWtm< zf^0=yy3Qn4^&1ZKl`=y5?27%&F)u@8CjH(rJy%(Bj<;K1LSd%h=v@Rkq>@QF%+}WS zw*omV)A#S7n1?@&o?8@@v=u}9Slc1L;B5MqvxQs?3(Dkk#4E88y@6v+l>-C*NHQ26 zX8hiJ+Gryuv9b~$VOXlJFNDysKFeF%u79zY|EACtuzn&qomxz=Sp8YqQx{e*67rZB+XA~IsnI!1wAo`1yZm+Tw5s!F{1VcUZ*5o=PnZNrXeGSD!{@8C1px6t zVunZ)%rc0%u1a&`E zKY1VP;wP$H%^34M!d+F7`P0YxN}Jc@kp`3MC$=oPN3m7(@AH$W=wPkWoT<)V&{ut_R4)Q1VNna^>_>riuif~u zw!d|5c=ExO---pY94HSm)h9o{#g}RRciImE$-+3=$v46wo0_fzy#0Zl2sCmZ$-7l>M-k z(xoskp5CGQ?SS_^j)hH?*^gBO&R(GD6HEJ6f3V=EDx;theSsxmT0ss=w`o8YoRtyu z^inRtehc*%mm5NeseQ--8dj%QleRM2mwgYegj_sRTO!m!9K zPlvd|8}cPTBM0Kv8WXC|i(fn~GgM5{YIO6KY2q6z_u*WIv{8Zw|J0)wMYiH(O?b-~ zx!{uNY3U4W#d@aiOdLah(;UO_bwx!*rDXJyc<~r%Oj%CMXh^~ZX|6zm{toz;f{OLv zdil;RRJBg>5jYt&vwJXIR7v$fyI=U+FcpZ0s=bobzpwIRO_Iq2`JrMn7)*l@P$(48@ik`RRfYN( zz&4Gv;$1Rdsu57XeHnmy#2Qk&1uEd?;VIxK(GU_>v&xVcU=ASsArnDmWpeNw*YrH@ zo#BeL%`Nu}DoqjN1soxibU=&el~r#y>rEhf&m#KrJ1$>A0OL`1Q=fFVmxc{_Z~kBZoZqIA6}(=c3r zkQ$MQqr!KN{TPd#sgS8&K8IHn>h=~RxWN0ii(roTQlWX%eIw}F7i=^kk^T7iXy~?c z591oT5(>!WUqC-B{TCfZYGub>(;OyNeZ`@VZ+hO4_)?}eV0ZKAYe60BjAb6f_`mo% ztGFzquFZ?2fTVPnbhng*AYIa3(%mhcf;32XcQ?`{UD6=k-7y>A`MvYa9L%W)@^EuM zd#`n^>-w)32Xm!VZv^-Pqy6TCqC4x{w|?;V%w6e)I(}Yj9;G3+G4poVVMPzpckuBO z5jjeY;MZ#J#&zhJD^z6GE1-Pc>-lZAUDL1ljca?gB;Q1t*%XGF#QVIHODi>uUop*D z*3I~z6~Kw-`ZEp|D~Xpm%+x`-Lk! zfZr3pjsBr`?My(C$-!0+aiM`=-}@1XtfklB=;~5+?l-3U(={xI7CNUd?jf~C#JIqr zaHHONh1Z&h%l8|k_7cs&B@2H| z=2u8^MX3rpoa)ryQS2=#9STa0aIROUI1ZLuHI~Ro2dyG)T(8^bP!^!I{pP!tMno;0 zML^zg9O>7b<0|RVD$db;dDm!AX!d=5`a1iy(4Za%#fe?IEK>(PN?80;(5eTU!tuzX z*h!5%`jzrdX>4Hdt`jLlg(?K~LDqc1;XZH7eSJXz)k~=oYpd)j`E+qc^59v3`kq0< zKfV1RC@nCWJCweh(!s+GqZ!NclSLl{_Jun1T;cDE@^B8mWXQ@so7EQ8>bw7Gj~B-> z=8$n5H$LI93H#%p+ctat8$C4~?(JJ8^N$}(>?-#*bQvX!De`1ma$DsbN+4RUTIUPU zWOEVe#NG$~!bM>EVpG$%z#!jnzk&whH_BiFI%ri#wPc3v5$_#>qy%4VNz?3$LVp!- zSSz;FV=5pZl7YJHlIsx~hK-ldw6j=@v3tEDSX4U^)fpIyEGM!QUMlfMNlBFYBfKXC+h{PDvw#S%?50IN%#&r2V zzIxu zU1NV(dYUy2W^~=@l``w&h|R0n+>v4#)qc* zA_!w%k&VlRZyFdbv8a#fZU0T4+O(R!cKu;x1Y0@br~1eE&*XbmjV(?insnw=0Tt1| ze52EmAXBwcvA620HQnp-9H6V&oK(+F!we6=bKPF8*SN9BMxzQVUA2o-bJVZ(H2>bf^8L8!4d8)xvgmD+tYKn%d&hg+h zaA7r=adFAjrm(^|dcFa^d4RD9xS@AlbM29TQA0Dp==ZDvm`8wuv~xRY(%eXY zv*J7yHKxZD=){1W+`Z7?Mhs`;AZXh+MW7|fmB5fqZ@n(j?U_p*iB8L>Iq_7IoF3O5 zr%+qP8z4K2_2TIa{4Oe4s*%|jd~(?vgq}U_KE+d@L361HANL)BX0|vO?bggi7^Ea3 zMDNqg*O(3XX1vvipd)34FaO15p(0Fv@Olct_skk&xlDu?EI;wgh9mjP3_?ssM8&N`%&}Vm=%{?cB91*HO>d`XE#X7_f(Ei zF>46td7Qzy!$>O;F=V*fE*@%Wk)q#m*V><-{QQI&-qx6?-5;mbg1B|dqSsMXu>e!T z`zx-m*RFUvM^^qx_uFe>A1HKVT<|B0IQ_9&Z0T&Ea6uvs4Yl$Lq2YQlz>D9#B1g=% zmXG4p{_*B-;U?p_zSVT=Tn@HP-LFubm`dXV$ra)ebV+bhYh z$_xu{IiS{;efHs}=kGhJYu9@|8VM(1a^JEd(6W?2w5{H3J0PR03Ko zy)c6XnrnI}6WTg5JFLfYemp5cYuT(}i{4pv#nuHT7}hw9ZMecliM-Zl))#DWIn@vK z562;i*C1?Wh4aH}BHQbScX4+gBL730pd`u+H}Dd;J;8@-+{NE;K;dHqC6v&E>jQfM zxMBJXhXA12mwq))`pTUo&qWHLEsS;<~9QFtDU)-Jxq0 zdHPDCBazrhq3;-LeTHW^IDlU1cw?e5d@?XtCct9uosGk7m(xf zN_E?!2RqgOJCi-s7*?jLd1nr#eSD295&i_5Oj#Rbu7J3PS%nPyYt5<%v9JZ)sEaZ@WQI zJw2>Xw5vf99eqS;tVnOnBNq`U->21#h2ot_E=qE8QS*OK<#MA&aF#mqOd)a%?VMk zi~Z}ZNYiXARTPPiY9DNIM;J{WR=dP0H z%{u=E*>Ml9*Cz}utp>fJ(k+LgdJPHr&&z-_7A_d; zk(5AY%SDAWtnERQu_E4zYc>V}mpl%dG%bP1)?f4HWna^%jJcMC+TL4b0MDC?V0>eQ zs+-+k8^^v^&(?8K*zapfG9Ml^e2=qwCEps3K>WM?7g{J)9q7L1)SuXk~~Qy0Ubc?dJ^zKxh9j_kDXQa|}9^!iifQg50S zVXSJ2QG8UKYb>EV2Ey@+6J)qzZ+W)P4B^7sjlkIV$dnQ9tl7oCB>5Y1-vexYuklx1#>rk2*UIs;jJ+B>p$;eVRtN^{Ch z>7BVnl;`Z&uwMozC8KG48o=rD=M%d7Yx_XNzQv(Rn=88DeRJQOhb?-1wcz<_aFuZEefHtd5|E{#lrDbCNJvTE3X~4XtkJmFqIhCug6pfNh z5<~u_C=Nw~&1F1p1~5?!zTQVVBBt>Gm z$GLd_;(6SbSH068-onA&+@!E?+o?;hu@=nLI(U^lVZ!Q!glV>Ng^6t8J8>1Trru21 zOH}&(cve9<(@`1*fo58X!`_vKNTcciHWhc&G4=*;W1kb@;}C?kR4_2uklfypo!sp{2#E3*tW7 zrgB`{zOklvj=g%Eh?coR9%pU=rd_4zMyr4M601I`Z5>jwpdad! z*#4n8=;{+BH)6@}JZiHZ#J0(}^wD78;y0L{wH#~nf^KMN=wD&h15r=#=2r1yvEft; zzn_t?cCR>s{Gw2VvxUqHzAGs;xt3tM`GY|$(ltYN@|_k;b?ukSb9%6fSaF(ZBIVG+7n-;)zySt zRP>xwXHztde2z>Qcs$rRj}XAET6YQe`soqoa8aR*NHCO5Sc*(Ij%5&wN?D;$qdHXj z+l!lG8m17|q&6+LhW+7{iv@4PkDp%peH9g&C>A8ZS~r5hzHrsqfvTgYo6Z0RPKs?q zqf3* z%6VLM0@BO)bmcInu(t;i(;^GY>3ZG4K0A|*s5Cj^p`Si9k>nk&>`*XHL>!&O+}p_% zg>2%L`l^n1j?J70<+#>Hm+|=(>#$ib^HO(Ics;rSyPW(;@b|bb?A)&UxYe6~wWfs; zUnQc!_}V}aR}>9I&F7u_WgSZ9ue+%KNC7VfyRUrp8$bJkJrjFm=$MfzPpNstNfjBE zDgze8SW(q(-!``2x(RA}&Ab6};Noblh;GiQZm!pcbIc@3`0i|ZI$aPVPNxxdNU+^l zqy4s_T2T$p*Es?qo3(3st9#k1-Q=j$Z9l>^`^_ zUHD*aFqPiMMTC)_IZeWLipo`OKqBgB9a#6g5QNP@FZlH(g@r_&if0()n@hW9MPmG? zjcZwn|92AZ66v)-t!W0iVJxk9Fa`xnA-I~%-kzK*iaDX6%12Nuw2+6+Q3)Erik;{4 zg?d?bZhdqdZCHs6_-m*7zw$sz#n3-ok`vt=@uxq?Q_#aMj! zt%i4)e&BuqwOX?b>nsxZ)AX6W7&#fjRPLm@?+o&7YzANpTr9a^oJ@V@;}+KrbP@29 zvkE3E_>1P|I`Y0jzVzu%L!8VGo6=LbO;E6YNcxW({4PDenO7(rvcffLkE_(MRP)<`a%{vc0w(qt{JQoxnYAcQGft3(17XPf%HD2_PF zGdU^}huKp==hn6+(2Na+WCEA;n395nhYgHOQ7m;;vRHmzVf#`)Z@Y zTVd0SjD2%-kG;c(j-YamI<$*f{$x@_UYv{c-BL<4liYTaM2 z)s|5!?P`}%|H@R zV@z;TFvmF;3JEXu{jxZ&T?*7_9OE9~_a3B<;Yj49(K{F+PHG2#-L&$$>9At<(yBDG zo*J`F);fIKAza1+AfA!_1q2C}U)i6Z+|0;C>9HjplIG*Hf4zv_Rf?g8P5Q2S%uA(T zcxr#~GW5V32n7x9X;*n8l!>Az$CBM^pHund)hDVKJie-S*u?}}x7Yo2^n_}k{b2iu@tXG!wvs{wLt zJZx}AXSwCK-167RS<}5G+_r~-FF$7c6CG*_hF?>g|4Z(Aq^xeovR@{shT z;^sWwE*5BLQH^%uHfuQR24#OAPKD=MP5C$uuljZ+QG=zY3d+%qq^`H(Hf>~6ah!ry zE$opFEwOU2F`IzpglgaH2~9`2T$8;v3&Tker>Ek{+sv)ojz(kui2w~B`6D!Y%crPRquw$6uVbA+(SxhDw#0vkmaL(3bK%QJ_>fW$t zh}!F3r*{+%ya?|QA0KIQr+>N<25Ef+r4A}){7}LHxc%&GZMU|&n8UDt)s?n4Wf!O) zz!0X-Wgmg}Cuy|F58ql})_hHoOTBj8-Jh3pb=N3YRvJNZ2%=fHhdnzBcdQR? zbh!9kA8{zs2In*v{X@pJslqohAC=$_zOD2T7AXEvvK%m=nT8H+dd0~A@w7r?i5Enh zja}2tH23T6Z-gQ?y3KDAwpd1wVT^5`lPto6gX#8OgyM$+YoZS8!Ye)xtqFJMyql*? zWd$1Vi~AB!ZE)Qt0B`Kc!I{kPPkW-)?XvE{^q@7rw572^BlMbsDS}pkTcj{&>z>83 z?bFFo9@N2|*G59SW2Girn1x5XYMkn&8+z_oMHrJF|MMuGUsZ1?-seYZLF9Ta%w?O%-9 zf^5DiAF{0@=WVLw@h)EwL~ZHq$lKpxZQkz{&Py1k0XmBdJ8i*HF4|@>FX~`%H(PL+ z{WL*b!)K)!-cWf0{rO-o4qHzkv*`C^WmxNL~IgTQx0O%x=J=^)u;OrCQDgqnZ4L+OjbW{uJM?U=}u~JRB(r^ozcIc1fHK&Meb6 z3?$JHF2SZL7TeH838DX_e0p_Bo1@D27ces%lA7?kun0t}8IH4iRKixDywAtXcU?W0 z&0!gWdJ4oqm7yHq@84hbqq*u#762@&hEY&;Fukm!j5Jk(LycwmYP@2>?l z*be^%FVU{rxo?wabmOYO)!FXe$omo!ObWuQU&bo$5sRi^1T0md$~iT%KvhazCwd5L zaJtEtcxOf1dZnFcw-a&qCK&Mt=bYjc0qa3ZuJYL{=5BTdxh3%+Y#Mj7pe6&&C7N9M z;gOq#{yLSz7*{xx;>J3)2w;u!~On;Hpz)Le1UIFGhpu)V` z@gc;SZ^|q+X1+v6!I7JSrrHmg7CcRKTV+PeC^JxZuCd^3*o4l6pbqNTZeQZHG#|!JX1Frs-X~;A1yxr^cdZRi zQ%PyrD{M`ZfRuOw5X~=4*0WJPo}eDUv2f)nW_zaqrH?+HER-d6vd|}8+S-Q zc+ugDF)k@<{3{`!R(-el%dFCUA~mK(d(_{K)hwCsEEX~EDW8ih9T`;aU8FiFM-xk| zWQz?5KZOZN&#B6|8>H&g;zF|}L9q>Y{@nfgBH1Mrg2o!^(Mu97dzh|7Bm}LB?oRa$ zjb>;ABY2+R$Z@qH0f<2}8G4HP z(32>E%d`|?b7<<7MD3*2=2X1L=Ww%~2iL$69|IDe+IVIi&r6KX&a*LiP$`WHhI$n@ zT%G+Ed;_JbaETSipV67Et7g@RLj9i~4LZ6b^7yA37^{`?p0>-$(@kJ@bI*4&xpqUo zc-r8hTW-a$-#_jh(|NE-W+kVBTwjk)pmoHyhq|!=_xBdfzWTfJSf?-7Qu~UnP=KjG zzK!^NF)3?C3$)Q#zu?38!mFf?N(X6GML$%{xbD1KO3=Zzq(A*DApcii<(eG2 z$5cDU;(^7~j3I{%WA^d(ahK^ivu3z~gU@jOkbmh=IDHJSLSXB)1a$3J_2(;gX3&+W z=~U`YqO3*F~^mzD%a!K(xtt7W4<_qhV?ikyz@esKPv|clJuO@rs3d2x<&- z_AuqHRFe#4ONPq;V8-%tYP}rkn`6~*4Qagoe7PbbV*Vua*{>J1ET%a~CtpZALtuVC z6#a#?f~l3XqEx`71(p$m$tiGixT_2O~D(WBNgg$qqFse)(Z1EDg`h zCt^s9C8@~o^-@5KGHKov#^7fXzkaKaE(qOz{2}X0Wv8n5NB1-(BqSv2Jh6khO5vIF zEYJHZic-mT?`JQy&_`fG82o^`$rwpUImrLnf|G@NX*DuCN|neW*>Xs3?*hyVVFT+s zfd^n3^Q90i;mOSQWWTuEHJ?&Ff3wHoMgZBK9&XC|Rey`kS_=P`UwvxN*S6C0bU#BS zqkA?)YQgaXBJ5)<-u(Qm7tE2^z=uL@B51?~$SE*FPVCG4k@}xcT&(y9Bj1m7+~ZCg zaRHl(G0(C40uHE}YD7Fv7HOW*?QWN5YN3V10L;+`1mTURiCjsb%-#CI(G36bz-0HD zi$JPI9+Yn3tRrkE`*bpmz2v~3sh@PvxJUhhY=F$N3heK01&@0RNV%^q*1l{j(p^1c zEz=uRdOUER;BveZZWYxpxi*c(X~7QS2>(P`FSK_ml4r?hCs+y*(Npx=h;N!N8UQ8` z_)yUyV3tDB!kl&Db4P$J#E3d=gHf|VZsr^bKqSsYl4kUb`dI-cfy$| z{7ILP9@!042sNzkPbHC&ktN7DEaan07CH?-@VkqYrMY1NYiwtxG`R$*fq4eszOL`6hiPwSt^%YYo@i$sOS+JcD6+ygyol1lqD-oQG8pzz&Q9S zQU*|ZqN906Qn(f4U3&nO#E3f2JD3;!18HXVkM#I;^1*VEkB*Wyb{`k8nSBMHADL^# zWMgNy=AYq%5dL7LOhJDfUGjqQhv^40OBDA@)q1f<=y5;ba-%7cX8bu8Ltnhq7&$I^ zv(9p%KYD`17oJv>x>!sE3S#I!snu7+>~TRkZgXwHcFyadsRUs??><4tE6UN@+Pcxg zK)dqA*vs((VENGGDf69p@cv3mPT6H_76dw;=nFOs#1E%1C~hqK*MtY}%VShi7p8@0 zs{OF_;x@i!sJ~o|zd^K*14+ilj{yFUeUf+;Vr9-u43*I`Z9*nn<_}1eyL1H|R{Y@R zS2r@FQ}~O31HoKm1g-Y$=!`O={XJUvThG>)UUVYQ52|ffe(0n@iSnz^w za_kOhwyNp9J>Leezw2acL(L>f;QV3F5cEu@wa!Nh?|(4@fvBh zL3_TSk#emZAW=7wNi)W#oEl4z>viac!f94AwCfmz)v`q1sKJQx#U0-kJEsU8r(f&N zhSaE3knzN7?p0mTl5odUh<^@6JtA)qNwVTtwruFEvHVKumEL{Vrzbm%0?Eb1W zTbxB}#74lB7h~|bQ&Mj0^?wV?9PINKVByW58?UAwyts^?B(uApi2dh7Ol&Sm2a=th zFvUvHjamrRF^vk#=8CsBj}eX`FyZ3^B}?J@-7Y;dwQ~A+p+;_v`Vb)q11D@Qxomds ziJxF_=N{U1^EJu>&4Cb9c_LPw8Ss~5*R#LpSpPncHJXyS%BOwZ-}KAK-|IucIX#JH zF%z^Y={7}uA3#y=@`IcKg{wkni8)d=^su!(7AmOs9gPMgEW%VTm?=)jFH-Yr;vA*u zA62FOWH&EjU!BMurWhnA-Jflf;jV>w>Dg(Zqx}$sD4S~wOtMQ}ugQ07Gpv|@ei3RV zpFMll@pWQ9{Uwd60!7OMA3zx4=r8?(hLycgUnff_W^IzThIHP467QD~NtI0#P%i_$ zzG%oIJYEp!S#G?OGpEU(_|Qsbr0ff(UV)+A$aFTAkWlPeS78kI9)@#rZ|L7&N*v9VXME3At3O|ajtT;T zc+I{Is;i$47V2d!2Rt1n2^_F?19#3CpkP?=c;vM?>YP2L?3(HF1+hu(m>lSBV9}ou zdI1@=mOZzjfvYypk!EuTR?^;4NmLe}&lQ(y&)2T((`FkA=IzH9Eho1}!){Q3P=xHH zsA@{DXOD3J929mm7W?H#Ib?=$eqh zJ)cSoYXpSNXQ7g8sv5Fe=YgLsHXeMy!R)?h+A$M0w*_IS(udasqL!arI;+2w_mFwa zLOHk_!{teN;ic8><&&`Jj%-!mtn6;+I6qRd#!4_BHr`$o>>j2{OQ`kSNY|VQ4uyhO ziJ7fo<^H(;&KDZNZpjS(#b0Of5>)<<3ZdXbD~hrypYZ5o206K+y;DM1;3)8lt3pCK zzl=LX4-YaXFSfK|h3=dMZfI02>-M6Riy)_D{YUCL$0*Der_f()og2?tR6rrkZxed? z{7T!JV~fjgkIUCj23G5As;lKS9N~YJE-^`&wGik$y5hPOMg&hJ0(2~^GX-Nn!Sht2 z0`=bO*h?vIdwOL2O9ZqWIXU`?!z2zZ$FT&0w-r{Chf0j@!AnYQM0rPp zSo6K-!=P36d=~cIwYbPiq{=m&IiQGCwcO1K#s|0M#exb+RMpXD&D4b0$&6=_<_LN^ zek!!1THaGL^@=e2)vwVhEGPZAYFBwX*6FG9j;OJa80?N0Bqv+c-O`nC}*AX{aG7Z6x&nmb6Bb z#|+Lm%%K0cH10BG{=sXBAg3*3$2z|RNj|p&$-}Cc-059+K_I=@Z=EepvXs>6?&!EZ z>LGbAUfDmu!OGGJ?%jz}ga))HktK0yzdqYa7c@8YsOzNhsCPlHJ2!c!0jMvoOU!*y zm3HAUwEfj=eJ}BH0PQEkKglbQsaUkK%fsX>5wgIl0x+V1zuv}0ArjEOQ`CX>+MR_n zx{T*F|5Ih8u>4v;KtOE^W{Xa^P!Tw5Yhg57T;xkp<=e@-*yPO5+x^@$%sPWi8hjaV1d)?z@ z((h%tnex=8p0P5+UbeQ&yw^$z$>#;SsBYjANh9Nf`|PI(vkW>EMUtN_xkovq@s(^< zzF;?m^zZ_PPoS%x#$X{z=k1T!Y;+Kc8!OilP)q0iGwg4o+&u?yQR~K1xLzd;>|!Ff z2~P@<7?x=n?T_huX`O$zG((>(Bu}(zzp6D+Ymkr}&;l5P#|)ocrYbAqtpU_v=`EBa zUTcKUifnHU4lA?$Je!O6GtLJNv+sh%p?a$-sbH}~jyT>(A$!0PwJf-+O!Q+?v^2ua z`pX)NZqrnIC&Dqv9o%M}V>ai>GOMO$Tiu=2LN{^kyQ*UELU4Xd&AERlz9p(|013l?EEbt$?#ygJi)7HJMKCMq1n5aSLYOtE9re zG-R~iw!2BF_1Ya98bSHt?779pN-ABMr0{tC$bkD4aTorT5p- zS6S^Ja62oO%Et#2cee0%6wQ9C6^^ZI_}FocoNMK)D!zX0^BVwp)}}G;u8$cdZKKQ& zLkLDw6@As3;uwZc0obUBotTbJ3cQ?N5Xg&eA z-@jcR>yEeg$pM_VDwrPKkpeZQN)ttMjQU$Ah4x3+(fIdb*(9-vT7)qDpgCy`dtHo_t^mQkJ> zR7d@dMVOPeRJsPli}k=mJiz1C3|7wOG?erCeJperv+$tAO1`Fo^P}qh+_^r+~X7 zN_AIJg!{UaJvpM`Ov-`)Wo-%)dGT&=cGHR3h#}KLQ3gkA^MQ){^ysBT0cSvtu8q$b zBD2cnH9RWy3>2q`R06i*0sfw{yf13(r*-|gXAQ@`Z{NGVGf108EJn}C*Q@)-#Kh43 zRsysUwb1n6^XC#M@*Fh#Z(2A%SvfuCWloLVs_N+=EA^HC5(wDuvg{^(_jpjBV=(rJ zcL_-^n*er^MN?-(NgZ0JOQ~nB-yLpXBm#8}3b6Z*$Cf*sWEKmKuxFo6*|>PuQsl$t zUOXLSM?!(|c~Y`gpAQNUJX~B{`X@B+gnu-a@uKby9RwGYNJR1EIef|uKkHeT4jtGF zDRu@JJdvCj1YZPP^4`GJn=Oh&ZcE&}M@M{d)Yc=Tf@m(~>y=O2cb90{5vU8FBa zhIg5j198QI?w`B-0!2Kb;E>Y;xUB>37^hrh@cFQDfAV5+lL2o~jl~~@cz}S8 zp+Joy7Bzpv*mB}=4z(*P(I17buaq$oc3RN)HN{2qbwnKpAk`xgW^-m=-Vn^WJTZd= z$e;9cO=6QY5ePMH08H! z7iT?X8?)IviCQH^Ur+8d91dDwwHM4?{}~k=>7e_Izi0kDcZkA1+s@tL4+`vY7bog9 z+ieyt$Pao5p^ZPzwsO5gDXRlx0xv@#AR4P&>HYC8U%R(<%ZjhgG=dDvbP9{7s(pT4 z*2moRBRz!jy-zYncI)Bm+{{Zli?>ou@$nHASQVLv(4+0%knS&R>e#LQ` zh@AjTSa=yNCaxArsUMgps%ozN<9VV0+&TL9`nucHt|$s z3o&%>?pPBsRUkAh0@n)f7hJv^J4&VTsv(&*d@uYr6e_ZXY;dykPuU{qnSAi(Y!gFs=Jd&QMPb%s^H#e7I`?oid%@j>0 zRpMy@AhOUgFmh*Do|pCC&|X8~zNOm4`)B33^8S$Si?D~hOU9!ie?3LM4Zh5(qzVSGR z3-t-P?d|P!)stXxfn1fLILL;Qo4yOaNyY&Nfs6q7il1K+Ej{mdmIko|p6d`uPDiJdQ}cu0s`?PMErcd4TZb0anXzt)uM2&G z!dSlFC$A$!7L{#om&fKDdt7!ZX)alORZFe$4h*$^oKIhh(Fap@8bXVFTl0!to2CF| z-;(#BS@oWgVIx3h{ese)`8hX^Hq;Z6>?WAdLEjvFK2xYjsfOkAlR`3nb-gESuFcbZ zq6Ynx?`6dTK&aWJ2dA~aM0vhCH7P0YMBaqCq}<^SH!uN-=VkW!PuFbWyF#CPd_bCS z<3L$?Zi9dz398&I=T2lRf!&)qie%qx>)1fv|0Vy3G6=xB4+S-xtib3JCoAmF>~s$n z4w#rV<9i}d9lBfb{*c6B&}{gTjaf1sZ)?^A%zm5Ct;XEOwGa*VD>nMV=28dP=msvh zVxm(8^1`JQl=9UzwdS*fj264v9upt_&lewjvK+o|sDj?j;T{g>Jq4nG)a zCbpx#IQXMD({r25ZJ+ncgPgRV*cYW~my%w+roWa9JFSe+$NxJJ%RzealJZtR$(q>l zN%uH=CORaSQfz^S7CEIk;Tu}H7nSu`lfy3QLO*;kDm7`jFSY()_q3wjA30P*re1Y) zo(~_sTH-4?ZEB%R20R@QWI^O(Lu?#J%bqn<;naDdw^R2I`G_%+i2NT4&>nEt>xq!H9=eza(O6KOBR*==q@>>tT!sWkAHFvJz z?)#n|x1&ahBbN^zAz_btcyV-7kB3eLr5~-E9Scm(p$wd1^oAN2s-&g73lxZNu{|FT zpJV$Ft{66%9QVGIq3B5?m3-EAn|R0UC?tJk63$NHG}XR4VrE0qe8#XNi-pFGHGF%I z0ucz#HVWw|6bKm9oPN{F?TG%5Q`v7Pr{WFup{77WZ(xxw^&nb|J>)vr`${Lu{4{J_ zPQ>tMb8|DD)^VyRwAeY})&-zr13^aK%V64hN4EUvSABYT&?tq`2uV_w;9|%o_Twaf zEFK5)K9-FCZ1C~P(4bM`?N>z1w|CHZU8R&Q#0)3-2a0Kf!BMbxK9I)A_Sj$DHD^Y9 zsbq1378v7mNl0P+j&2CQn{DhwYZ-oga|mr_eE6)hS2|C8`KlD@9GJt2UPZ>?sWNMc zKBVn%c&=}$Cd|q6>Wzxq+S-~)H* zP+ns%x1}EtW45h7PFCXmsJNs3vxU+GH`zDb(Wf}0KN(NJ<#ZJ)bw?B5m26ngO3fKY zyho}N=5iP#5NQVmxe4Dy{1<_-g&4dEoWIR?iQ4rgM(3-nhC`9}a(Pl;$Z`BH{UR6F zueVNAhNHQ}!}hRG&m|os`A0WSpJ5vo{`4#3{rD90kGpZx;g!B zAMC;(>YUALEPcLAT}M9W&~i_RY2GW)&7__@Q(MHc_N-lMR(jZ2I*dIAj7*Vvs!2_eOKn zJPwI@>61q(SoQFC`t(`f>IZ@7Uz5;l-^jzKUq^{THy+o*4#~9|b&m8lV(;Chr8mq` zk5NfGD$G`lYUBFvLWp%_#)v|c0F`|>h%ApUN4aYN=A{kKtBs8UwD?#7FhBgi#cSoC zOdUUYTVAze^!a?XJP8#esgqE}AnJ1rb8t>PN?FLQye=^kJSN?e0%~m5G^M2rC6I_jYR}C8}sX+;kt4$?*dZB`hfzKM`3=FYp$MeR%ZX84?^3hL7r%j(zb z)LG!cq`i2-j;}y0z?wr8|M=`1=i5e!^XpiJ4Z{XEcy!MZaBedLcY3GsGltqMg%dvO zSE&d}rzs=l(V~APqkim1W$*X4F$}@-2)~n!Ex(Ln>3BrF`-n>W845a)yGn=-y77?H zkDcrPZF-Fs?fTqHO@u?t?0lnMn7t>$N`ctbUT zhZDbph3*|*RFmI@rl$WH_Wd)2$RT@rFNO25J|U+ja|i*a@}yO(2a4e*&S7o}u%Dy5 zTI6Lv$9Z^H9hZcP9fVZ?F=?-xOgy%!ZfLPIb}xTYAoFRAQPMKU$|f@J(Q$d}33aX) z`6j&nh0yi%JkRVae77&^Fi68B(Vn5n!clYb+n@~v*w+6xr%6gL`;5#x0ua8w?_~QJ z&M-zlO&y&ANCrki37^k9Yn*Mo%?ES!-p)6poIESPVaX9*G_fX{Ayhgp>pqx;A1Wuw z(|>5SdwaxREACw6NBkOz^yJ6#yW*+vP3ubX(+#x;qF`bl1{-3PlX0u{XH^S*Q@P3z zJUCQ2)Z|*aatUr)mWc-p+^QJtQ}P{AlV3_|f%_kh>IBt zXkU*Yen5)?7IBk8$9;77H_gr|eghhB7qG@JzQfp5AIXW=N%Hc)oyId`eBBaBVj}{X zquI~+&RoF4FHLfxi_(L`tpP`SOIGvZ>3@Whh8oae%UAa)yQca6V)EY;hF~mJ+=UX9 z@zQ72(MZfg^3h<3oQG=0v(+zvX2?Tmxc%c}m6PBUqj99;{>(tK_Nk#)r~ezUFZCL1 z3SmVJX}~UpTGu47kgKjCph;j!d45?h(lHoEyEzC?uTl3!qxvh@FmfiOp!wx$cT&vq zcVeYM-yp)2^NJ@{y608W{7=XAJmW4#gUXV>(%k~JeZTU;PLT_KVYl1 zZ%<%Q>L|x|O}ReYFatME#Wth$%_XC6=X`x0EUq)Yr%};E3TzpL?o$utoD#6tqrU0P z%=esl`!T1EcN_uTBYtR?kg@dwt>{Z(R$Dap_3=u=5X^FeEhND0_5f_lv`RR7d$s;f zV@!8x8MJsZE7DJ^;){zK zG+c=|q*U*%?0RabnBBy2RWM!h$QNX4%Kep;D6_=sBi<*Vj118|x}WUkrw{UfynbdB z5Zjo+@=tnIq7M1Y-KeJ>*{=WICAF*}3`}ZC~AJ6XS z_n8D;71S#-WC&=2jmbTp>zHkVtsvL5+uUm0>I=B|bk1kfrE=>pVlO^VMeC(WFJ?wn zZ#Y}g(Iu9ZkK8n!KfpR8-Ol^B6D4b9&2aB%nRm`845vXCGWJr^j}EWoz|jO3B6RkI zVX?OQO1Nzw`&=E@7u;RDZ;eY9-2EM4!)g|sC6&{2HJ+oS8PEtQmVWlt*q_~y4tA+} zo|F|Bu&@7H9ls;DRpShHA5#5E!k)vNO@@;(9DkQwj5 z)sycVQJU+dZQ#6RVQv4Mp#tEvNxH1;r}yszgTp8r^`{0wGc-Y9No)nxeT(h`Sy_@^ zIBF?J39M_mw`ChYT#%1Qa&|{ww9yytX2^@)Pbw0p;`U4aImcp_^RmzF^XO zetKkqD_r?H(AFD$HdGtR%nRGNkI&Sg>w5d3=Dn=&l8zc%Hik`ajaVprmLXi}2xy(h z%NxXy%L~SMs26aB>LXk+XuKElQ^ddEWCjev6^uPtq2;-dZ@yM!Xq+QW7??BVI*dcZ zD70U8U@^Oy+?O6a1f~99*R0C@Tx$ol&KUi|*Ejlu2{K=mJICJUV{k%=`#`tnN|+9Y z%Fj37U)b}dtkeroR}i-x*kCHQ<0Ozmq;e%}V?M=$#YbF3DZY%xN8RDK-b0jEw2^`O zC$sIjaCC52_kDVu)*AA?s9&WWQABk0tU8@ff}a<5W{Nt@G#x@$%4=0@j~ zC#53uRazI))^y2P?l!;No0j}Jl0l@3LM9zgrv?mevV%^jE-0La>mlV0I4GQN)}OoP zArWHFyFX(C9+_P6(?oUUH!JlGYC) zBytj~yn9=DkktH(zDY>>;_<-p-C04V2jZqFyM&R2*$@Inv^S-LtZh64`#H>qxJf}^ znm|E0zv+4K(5m|Qf}WvUHsV6|rC@CeLl9r$)i~XtA-k}r-{BAPnBU>-Hd%%gZ}GBa zSiKLL{~=YfL8828v%oH4DL#`whf(vmf`rfSCIo}|v_PcGq*koDw07wZdvg8>Hr8Yu zwu}H~y_nQ}N+kwE6l%Gp7XE3!9<10UP)T!SGvEI*b0g9OA8}>Z`ZPaHlcm${#e(g) zJwo6b(G$=N*PAGrvJiq!IhY@;?-1;A$L#&&aqHT1mKvWn@$pvjg|t1SMo36l#p`xK zh4y@}A>L-*|K6yFXt3#JOPNAT!jAjPh8OB5#yYDdnti~f=x)mYBvR57bi}83q5P{5 zqRDF~Al>*PFrR&qZM&ea;coWT)UGq)@r^}aCn1I|*~0PKO^x&T7%k??p;n|VeVm?S zU~s2l8ep$hx;08%GFmy@shDg@K_L8eQ0Zv6 zki*?}CoI~dXkcqkndp9$SH_)G!&^j*2KS^a=gwMIOJWNb_?{!qF%?0;Q37|qT%(Ns zGB&(qLi&4uKfcB(o8*QNDptC9S7JJ^?^Xgwg{uR^&uB|M*6x$fub4OcO}&ucd1TPB z;#}4djeJ^)*_I8T*?C^zZ^7|=<#`Y*k6P$TDb-N|n+vWH>?isL^|IyYG-sk`pUF}D z)5;}Bmp_q88k1<*2_hNe`1xc_p+}C@I3R9W}*q1AAh0P=d*2G zwJ#j(A=~}T^)5TO^Tu4s&v^EVNAGqIg-dZNw+eD>m4=cN;G-eJ-w=zx;~Pu07UZN5 zO1EkXXW!E9XPl+!{TIyqp>p4tVh9K)ivSXTC7}6`pTC55x_aQzGF!sg8}6O#16_zh z<8&A~S{$S|zH^%Be1M6DPKH`>W$f};^HyewUufn14EH+~DeIS&YQW2^YLzw8A@rnJ zYeXAPqra?>WzeiF8^udUD|;6(#hIn^hXVQj*sxxdA-RZA${mu{SD&>WIz(oCl3@}ds%qlXIBgZf%sSDn6Caf9QD~|TTbIvTsi9U& zc_N-FO7{h{FkU{JoYN%HnbiQQd{Sn!YwG-Rczdd2_B6avJT`Q0b~|n~k{BaO6#GyQ zHx_V{-SdTe?o$OmH4!l*eh}wiC}no~ZE?{jfoSms=ISa}zImGq-;#}RfBV8Gzx9mV zzn=!fp(EoMW`J_IH&`U%a6g}KE7f^9%IIsp=#Gb#JA2O{lzy<+o{^%;0fG?d$WBM& zR6XRF;L+hDJ|;DoaMudegeBjVpmA0<>D}!z$WrAo9^7b0$e?6%c57;`ZR51-P2}yX zHttgWo)Jc$)L=8#C!mMWF6FolU6A= zay1aY^D)lu1-cvorQ#-be>-@N#Y9^~yof7t_G^cHutRl(f*tKE8L1kBmt?^uaW0)L zAQmFd8NDvvkO)I^)$=w_m828pk@68o+LwyI{Y9N}`^jJxcp1a<!7Os zu-hvw-JORPIHYvSp}Rp+T9EEm^3WjNjUe6KAkrY+-O}B6^LuCBxifd>{>=>g?D*F6 zthGJ?!$XvNh?K!LsJ-PlJCWUNl1cW;G$goP0k{y~K()98zVViq-lZm|?J;oh(27$r zD~_^a<&6r4A{`7INSP9bN*tf`M}O!eb}$U0alSt6DFV3S!BUaw55$go-Bf0&a=aXL=N zw*OFHTJ^ja6of(m$S-3*{Mqh;MFYrUO`@m6{Kyj@u!H^SW;Cy z99ibT)uKel)j4D;ofP7;PiR^|q5NcDNP*`P;+wlXkiV&hNJ;ku#U=qT7TzL&F#qCFOSt*-UxQA1keic?WEgG67Se z^f4OPas$i77jb9!o+M|)D#uY~MG$A350nLb0<=OYJW{~WXr1mA-_gLl!~2^@l3PKY zx%YEaESusR5)ir$WFwV$5#`@(*%4iOn3A6bsRv>n ztVSwd84B+Naj>?zfZ8l{8uayVIXk!k2Hn1Fn#iDJ$oTTM;mlEN2b*$Gt^Nr?z5;1U z^(XPIpS35EMA@)nBqlYe&^R)Aa7_4uYc|^{(%OhK-YXK#3JK#ew z2&GA@N6IfFSYb*0&L^Kvk`cPrpUc^x7EvS}&HJq6wdyRzvnK};3lQ_aFc&e{y>%jIrPmm^KT-H}`6rH#xjxl^j-f`%CZt&iEjZa685T9Y#MBCIHge zIpuB2Zro!h{ddXC0qYK^vHfAhz(B_3PSvhzSmtjv#<6zRK=diI84*wf*e)D)+JuQN8CHeo!%?AgAyJ8V{@fUUM9`8`93v&nmqqn?t_yOM6$h*dlyb? zI=Kz?s?LuPaFrjhnRwMBq=I0^nxAdEU?%w-0Har1AUtL*Kp-R$l49Tv*q#eGQZ!{{=A~n)9H<9cdMLec3<9&BrQu zN}1_^O_y#uhY+U~G881EgCyB3Lzhrkr)`b6iU1WCvIv1*y6q`d9gKURrGibsB#EjW zS-LzT!7Db+phB}RSc!yQXh2?plv>lxQoo^9eko|woL4v*`Xo+@_+S{u+7YNVEm`cwT%^Z?j6nJzD^AT4o893mVn zoz&M96`TP^dSL5CmDycuyZmW6FJjPbS4psmKn;ZIu9r`g`@sA?&Fj8rtM7E&YWmke zvv&8{2+wM~OWx&AJyJb%xJP7MJD5w=#K*sUeplQu4JV=)0+V^9D?X2izAuFWdSS(f zJr#rnaMGl81W^I?FcOpKUTvPYZ0~K3ijEomS4G}ojLH7zv_=MV^F>mT+wQ&GO4St7 znFeYXQ?Ge)0n#9v3zA@0rP*mWYGD0rZ0O&7&GxC61E&B` z1h3Pe8vX8c_i2yYu!*nzRPCoX&jLe92cofzjFgm=PKWy$NKi;2`px7W;IVLUJp{#n zx#b6rS{3!ees{;8`riC^q&~M2Ns22K%-K|15~;z_1uu)A`#qRI{^@?3fEDQuf$cXs zj$(ptpTMcC)$j`h z3ni_#-IpKza-=Um=tK+sY}>>C4o?Nf#Na716sNVVL>v=(v&&*w!ZLuJ6xltMg-@ve zB$oAuB&UP9kX@XZX@Ek32rz>V%TXozEb&*Dk*b^?l1y*Y9ptjO0GlSjkk=1YXwBrW zq=cqno-F|4|7|)=c&vwXWU;&0XO8$8Y;8x7ukZ<5jgr{VxfGK> zyaY?yqEDj;W2%(XNn-BVc2AxekJ}CIAy7LzY4+hzI!UZXOD&x(A7pAOz#c|CRESso z9C$|($r%3BxnY|u=mL+edB{Ritg$tSI(+(!lGI~tYV`Rwx5tvazTXI25=*(#;~c5y zh@hCIFw~5_O_n6I z>WxP`ix>KUmKQ^v|!I@!7{JD?60d*+T#O$d#V6X}{0ne5|Fnvxt6AHFzyr+`WU0X%6efg5RKy zEiSUx)+WfplLqd0&pF3Lg}NU=)H3CoM=cSPy4*|x3UovUia@Wtsr+H{dEtishEK-2 z<u(VsI;ly6*f@6oqJBs)Q78L<)CF+TU%$%aPf@E~9`)W}#jt#w1H`Ft_pHQ!G zVtl4$&eQw(>cP>pi>V+Xr=H)~a0zY?1_8uNDM=1B0f4Cq72&FIcy!#1>7Ytc$UM_4 ze+Z9>IFmS0?J}$IC;9?6KRa2nI)Z{+FEMEbq;OwyeqenomttvV)e`dwa2Q@bdtjI| zFNVL*S$LK$QN;3VUDrV?E4&?G{Bq|Xiim~~`Y1(I_3T%O!yf0YjqSbvp1k)W*3Q@b zwP#Nv`CrezNU3)N-08OcW$O3=f1#P+FE_4akS5vL?`Jin){G)jaR5Q?u8Ol=a-3{$ zRXtI_XgyO}JTnX&EF>`ivxWZ+cJqdGT85}NEvy@H!Qj9?R*;}owq6+lyUl>}@Io6b z1-h@Ut%>g~lFjNtglAS{XM@>u2-?{D?6{2e&%exU&<6{s3DP4vX<9OXCm6QYi`pX2 zpx)V9=g@vsSlUY|iC=o(9O<3P;7A8;8aM6$Qc&)iqX<1HffW+F_b?Y%EU23H>ohfV z&!f>*Aw|~AR<&bSf(B!%Og@`LKgT2~JL(5AnzHO*C5{ES^Qw&9^2jo43Ik#676!nu zB;*}jhW=05B~U()uTl)6SJAZ&-zCerT zX`l-|jPUC(763%nvP~<>aZG6-+tR)ps`{xHJCE4g4$O``5dG9W-_RD)c*_pE)5n4SLu1uZ7VSq$=msd7mA>FK*v65 z#P@vbD@-G&1i!C0?6Bal9AAA_x~vWreGTZByNASd>tD?6;Ub7d$V#_$`@5$Wy*y#X zy8Ml-2TG_vn9F2d`Re;8fZQ#^YxF8%@UBc&h6v;jFK}ciR%V4Zjz2I-C5X+l`TMQk zWpb64qTl%4eT4~9s2p2zh!k8I{Q{%fP4Kgi>dGS#P3G-x_s_Epy{i`+3WP@|{12>0 zq7#8QmY&vF!1;($3MPzF+!NGoN9}OqPHBo~d$b4u3Ma6AZV$+Qo*D)57fje8X_hW# z*}f644L1sDtioUIe|g#oxsUVSR;bqzi`J8VGz5&cy28m=y-0G))zd&LW$;ZBlA!W4 zAkYsAQy6c_pg{2BK*ee8=<5950V>^W*)!;9^n?h0#5nGXGIv?m`Jprmb!;OZ`A%m? zv^+X}5S0_EP*tvs!{^_bdr#Jo?8oj)F^yo1CG1P0M!%z|HyQ8ncIv@Wf>xnZgkE42={?!TcVOI?2V zJ@u9{0bP6R*Wud-f&ddpb!9m8XP0rW|M3U$Ufh;2Zl$gM&;)nq-q}xnl+r7wn;dLBm4VTxqgzRw|?=NZ0kER9|!ir)#= zTv`fer0H_<4P1_+oY!2su?Et{(jMp4aK5J>cN|M}wON~9#{5DcSKyN->LR9bg%H^FX}EW5G`6ncf#4yD|sZBp#nHT8ea($qG&{k?V)ArEbO+ z@_xCz&_GP7j6*+u9)h3jIRaEL^l|_SkA34^-${Rb)XNn@p%nApoIwZ@^oUv0XBAD1}dM83`!;du0a9>4S$$pw- zW?O$3KnR%F2nI6V*Kf>T))!V@O!aGt2Dho9Z~6#Zo~nIse_->@3Yn4Dvtxp(cDHA- z(~82m*@;C-FIq#>Y|ngvu%qqha?#1XP~I&si0~EucZ{c#{hKm^Ehrk53uR&4GBH<` zN)#SXZSw?7Gs(*@6)n|t!?P+M@`s@;eGx(`e_~L~<}(BTNWA2kx`7}U=37+#Z(<07 zwf%>C{^ZO5%xQf@{CVlFkPsY>sD#JPJ*1mzUdNuxg)3e+YnB)0=yNtLn!Z%}e7 zp;=|gnlX__g%;iBB*N12WE5_Z`(dEYkmLbPF6DyAzMZ44dLNTS?f?M(<6HiIp7Q@q zSZ(4-Nv!U<@I}R{RB&n_Z|TVCI1(gt|H4Z6xK7Swp{U#d+;yjvUK2f!&gia zXZAj|hs)XqeM&o941-2+7RODB7sEbdjqjSS76DJ#2BBPp5hl;2+H897fT3USLetRJ zNJ~!f+YMfaNtSljv$n7ko$(BmvnfMwr44i4CDB7k1VZHAd;@=v*?qPqGS?HD6bDn~&a^ zGd%c%^;+GU13nXAgIW4#!`@mHECwv<_sGI6QGk|~+VzGGfT3z-n&iE>fwb2e)R$$l zZb<>2Y)(^+>2Wl@KszTdndw+lniY#pwiIR;g1)`P<4fVQ``Lw%-iMF#XD|5lJ};ca z>N!h2uQBruFT;von=ao3rIY69@ry|Iw!;$&ol<-A+B4&m&bFHmjl-y9(*Ix^DGvgL zm=An@*^$@BBE(Uxqh zInY;`UUL1zVv;J$3OO?{AiNmDz8(JO50wxx?*6^hB3Nh3rJlj>V)`YrwRTrm`r=tj z0-#E^U1qgb;r5IVa5<>4?i7)rNQRT}IKa{iG8oLPjigh7<~Ik%k=!eIV{o-GWdg4{ z5`p0*Mq1;Zq&kbm-QmrP1h!OfEGA7Lrs@sR z0r*#xmCCPAQ!2F#1uv=xMOG@1k^&YD3-vEd2SxZUtJA#-*fcaiYC8a((Ysi6bETW7 zNps3u-n@D8{w?ruOjlJ?b1^D-(LUzIDy;Cl)CceLCnekxzXdLXI&p~LWFuGmAHe~9 z<#A9)c3QIPV^%LWR`1JxxD5d&af*YR@K4c*WPhDYW=wSWm7cMd=h)erb_@N3)u#hd z%+V7O8b9(Tv@f`B+kKHp>(PF;g*TqK)d)h)FCzX(QP7(T1P!<|KO7GZxaF~rf)xma zy<^b*UwvcYD6S6Wx!4|QkV=`an@u7XsOx_uY4QcRsAh9;LEor7V@vO^;8G7k954_@DGhpn7>V6#o#llh)HBsyYa6!EFSmL+Lr zMd#RmW}H-@1Os9`ZLy64c8I-JZ?;Y&a`?>fyJ%NJuQL>-DOf?svGsIBR<9J@%MpD= ztJW>DZ~GTH(~Ly8GuGgaVV1oIO-XC$v^ z7@)n`Gl&R(`q=gcpJ3;`!`sZO4MVpyD9U`0?CgCj47+hXs@S6lx<1m_nNUgZ3AUdh z9J?aimz&wvV6L#1^NpMP2ebAsIjnoJA6O-o-I;|KejSlf`Gz2&@2LjB`@3h%l^gFx z?DbKZ9-&czH>RnOj0McvcnO3MAY-t6mDyMg$7F5dti5g;TQ!)AVvlGOIFJGancm%o zuNFJr_)S-hSP#X-dW<$1bT7c)XgS_%p2Wn7DS`ipCFnVE#Y}J>tbJUFRJ*;C_G+DI zYCfm%dDCIYW7N{9Th)U|;(LdK$n9l((XD^VbMzT-Dig4K_1q4nN=#BFIH>fpK;{5fp?SlDr?d^06L+*SC0C8V=%F1;e8$> z@a}M*HnuJFIWu~-`M4#bXH(xiD6Rd=cmeb2f0mGv8frp?%vTtDblRoVVV`v30lOs*_Tra8#ry85U5$(Mt_+Cyv(L!O+dY-loadwI_A9vE1 zd;_Guy84%lniUN0lgMuD#1Vr_@T{sZkPn*;MG^aOmb-U`S%4}tNc?8}2vlM4<^hPl z62ZTuEu2D=MPd~`jbm-TK<+J46AXvz_mjKk1IFdPUk7AB&x^06)h`AbWdnU z*$!Yizkjy#SbPU!m7F1GXzqeJ<7bac*wdhbqPRHZWF^@8hM-B+sm&uB7X|J)P-3h6 zg80jp1n5>@VJlzgY+dq9Wk({hL9=icc|;av(WcZ_`O#Ndo-}-w_^w)2V3~|NHZk8# znk{t&HCn>)X$I8faU}9~Q9STuN%c@wWQTh<`PSjiq#UnTkDTVczD0@?@>7E<>KPnw9R0@(ZlQ+ z`BXo3O(S5x_D#=xx?IwusTKe14UhuaA41{Fz`E*v5&5VsgdhGbjL?88ap`;`xyIAg z?UL?fx?v8tA;=`7xroOKyPaW0C6kY_pz| z6D1&sJBWc7lXoM+Nn&st(lSGP0|nZ-;Pc6PgYHnwqm0`nY_$Uf+Cg2RL9J|2OoC#C z5B-Q|HJ}L-gQ+*71Lx@Tp4eNErGyH@9rmT{BeNLYaq*Fq9J$BUDN-DreP4VhS3MaF z$t^AcX!fInkbq_0IVC6rX zRe^Y=ZXhf38|Ptz!M3Lr1rnbDI8CY%;1KZj@h<}5Y-W|H8cL59mNK^1z9A7pr8}CT z*?>s7=sSGEAKHmJ4;Q{TAUiDwOj2{TBv54Fx z6#9Zt79h~yslzI><2mBnj)$601R$2>x|;jamvs^k%hcp5DP|lx`521uZ>7oGc^`w2 z-FH$fw)&>0HPA(I%`oS1;kYRLcX1HRLngH zeYFjnz}9Y_K2x*3!h?;N_ym1AEtTOYU#^UG^loi=W7wk%c`$S28p6GyI~E5b?bI(( z|MSLv)qRO3klR{|WsOuABv7+eiXV}#*2q38 zPRM=!Sn|t#LV?76a)do6`d$Z4&>fW21pgUs`4PQiERzPq4C@k}X6hxH5?g!4ujpws zbnScUix@?nMHj`hisjv{cRP$w_nQ_Q-c{pYt(mhiI3Fr=7i)s#!xQpzN=dS|n6x>w z3Fw6PqKa7)`+_JEH^)-L9jta2Eh&|xQlR5m^f~*fn!H!Lx~FUh%AsE{BcCScXnT48 zDrWFAMi*r_J|Qc8ZHXk|F$Uykr)Otk9llTT7F>|)E&9}gH$BYifIVI=>)>xPMqcnn zyOw)l`%$K*mr)F^2y^*iDC4QWU8ia7A6z#%vEG?8`~40KaUXPJ9}%qrVa~I5;kzxI zZIzWwuVX&5qy0~X zeq}Kf9nbfN=Uj<^$f@Cx+-b5zi&Z|e)HJr?4FDJI(~_cAF7#k6knDJ#wmEh5q3wr8 zF}ZO)_~;RTBrgO=-grx7WhyPXE`@2f7E$+`e#MjVQcb!cI=sq-R8f#G{UmnP%LO-z zHH5BIo>EAWrP%f=)p4wXFU;)J{V2{84wYj#e8ePaKH#8Z3Nh7-Qq;qbJF*6L z%ih}(hDq6~k#B(oHKu| zFAq5q3d>nqnN3(&S_V%OX#BKt0)@hbfI0$>*-#VOZD`FKlv^Hm^P4lumT(}KNGcZQ z${_@Mc@Lido;AJ5V=>%pvifdMueEg?l-iH)!;QpxdJ@j8TO9XNhhmWQ)5s&cw?crv zqBi@O(Rbp1DGTigdWtrqE^`s(#n&qUDml8foGu+#?bojX!61iBAceM0*Ax4o=VwLk zSvl6-@loZ3*$0Gv!*m#kamNw@GdZXNOVbnTRPFA>*67=u`WkyzL zf=Fw;PD>#BifK@3E7}v3PF~5TwseC$Hr0^fkmG4hN z!r)_C3SQ;oxe}q8pdZXW^G_imd2m5@!;Kt#9=W5KUo+@NCK zg_qm+gA(~l0tH}p>5m}d+)77paXSglGwJe5c+x^{p(qj~KDi0WmP_%vOoumN19k5DVN&crC5kcp4=U61SZSd2(}`7#%8J`&J)G+dd9nze zA71NBX>pCYU(|js`sNG9)esg_adY3oY^wRe!qf@ro|~;vEac}P$I`D2UrG0h%v{SC z=@gW%yoqSM#p6`p!ZpGam*$RAfq@TOf*~~4-fYags;wYcyuJx?^!bz;IJbugZKI&B zK+)9vA}TBRmmL00@q$h{SgIh$Zut?zrB!pAY^v#^?VQ1`nD)TPylHGgV57nN3i(Vn zGdy^#PI3%KKOwbMBGC&kVa2daa8Pm8jp}~*9mGyIS8y;b@}b6_8n5>GP!NNvCZq>-qaLKF1 zsZZ_`)8vK`R$AFh2q5;9J+#*hFZry@ju!kvC@_js?&~>6N7FE+QHFppu(12c{}q+uyKBe zIuxVlt_RU!Nd#FYLBM1xjfL{3xJ*+tO22K_T@#(Xaeok~Yt7$b<9cW#&9}4d>0D+m zJLG6v2hTB8;q&fsz$t%5;CLV|pUikOc9y8ZDGl+oxl---w+`4BeNc-FG9%@d>Qv*K zW&gGRS@9#T@|jrNVeff}!g^r6?Z3E)H2yj22!!ikotVDB=^b7O>JT}CR-c#DaQ&bB26dP$91j;~81|(ED@t(4tv`pD zt!c~FaE3D`p|EqR`i6#|s=PeFb!1H&NS7z`Bipt9pzW9XYa*yUPPYAblnIY7E`&KE z9r3XawZV1?Vr)E0F#K`*SAw$M)_ zd&;|kY;5Axo{evFIthaDK}RRzVd=G?^-RVMM524`W|MQryC4v!%sqAsa!7O3x2M<0sw3;H>X` zd?0p|GFraRk}99=bM}aQMuWCU1lX?1-f25Bz_ajAL%v>N*gN}(@P=DZ`}H?AXXPy= zuAX`%$qqLu!x8GSDK7cyke<&ljtr}#6fT-DBwMV>~MjKpeD)}*3e?w=@n%EMus#YFV~#akUs;I=CPOn+&BOT)f+hc zUC#FK{kT5PJPH22Fh<0z2&|L}_S@i_6)z=4XJ#n5I}qD3Kp5;OJE?d)LuW8dVD`s< z4xvG{WUuWFSdAIYg(AH#@x+|zgOb+F+PBkMe=HCto$v!tFsEi`$7=Js8NG3?2K`y1 zz^BGI8J`sZm6Y~~^^S2>=rih48s9cZac+zLktz0?t+EQF2h#VFtl#8OHiDk*>WMntH0E;~bfQYqqWaQEl^GThU--s#jNhz0Bw9fXgMJF+LNOpHiaHF7Tr- zYTBSqN`KJ<_=ZMN{2pb%ap-Tg!G9$KbKED`orxoM+CfC!A*m(AgLrwYPGBZWfG9gR z(8kis?S4IP;jsQgFDSvK`KR~T#s+k zdD6;p+#j%LX%_AmOq1j307n2Op|fAgRmOc^J;Wlm8(F~{tJ=QUIq&#xuqzdyV##BA zg6}ZmVdYOcUfNrqdYdMj4r9;lF>}Pw`bk4F58|?H?2j~2!6B9x%18-ANgpQ{5}vP$ z8;JB6JSEBWq+6O1gmeu%1HbjT7ymG~2%k%(ZKPZuC>~R}siSd!7?TY78v!N^>w@3x zJ^I<7J0t2Z!d)oe6%lXF`z+S~b;5VTz3nv8b#!Nwwtc|i+57`zLX57Xv_9ZGWMiw8 zRxFIfyJ_p@yt1*lq}OLuX>{D8j7o1puemfLzcy+dRlk$ zX+fMDIT;88J=f(+3CAlvkh?13$kGBY~z~qORG_ zfa7-7cd_9DQVO;;61(Kb>QIAQ>+4HP){uE*Ih#+lHtL##zKJ2^cyHSd6+|%J+r6GB zvp2uLaNLkedavj>M#}8ZXQdkH$Orkdn9>UIbRv4YGu1$b?^$K!HfSA0&wrBb}Y40 zoGRc!p8T2>1}uij#crrt<$8=s=6PS?DK9OK?#o1*(p<#Q$oSR4eHkuJ1-{Q{Y3*=a ze)mrIyAbCig2ji+{YghQQN6EWnFv;uRID2=8&?2$rDwjzVw(zKBUj^5u|D=i2+I`n z9kzrbqeFST3qAtYC6Q4+0EvV-3=9&?)AKSui;AI ze9bc?Ps~QS|5Gs?Xo&Lo)G~3Y@ICg}2%`rnt%GjyKuYqESGn#rI}raVA#>@C?G~i5 zaa`V*W~~{P^miq6mbdE8g6}4BFIOMaDp%Yg$0W*4*=05ErG{d|p$$Bk_mYT&01a1| zm=x913lsg0#Jg1s-R(Q-2R&m2QcgleBx}~C(69b+fD&oMNljv0h}%`mpf71EV+dGuzp}jFo^jPY@8cq ze6*D`(MjcJk{|j<&}PvEIm7%3JcgVD&P6#s*C&_^KsAv&{!O9-Vca|m!8g4(fmdbzTDkc88VpnsAds?+4K$qBk$z6|LD72r&1804w{!1E1wG*mCIKL`) zC#s2Iss!(ageKo<#hgo*XjMHBe%%a`nemm4oJpDEN^dPC`<OwFny^T^33#(ZzY$+3~5JBM%Jlsyk6<=nai}b%% zv>ISHuutqjTgyFTVoFQRE)Ag3tV{Mx>{3ENQh*|=K2DRJBXuT0#ZY9jPR~+denj*H zLhCA_AZmfE{clbX2gPrIYkW^;Qo_>{Q{lS4HfxtBm6mvVc{$wPx)YE+Eam{J#su|u zi>9r{OMS8_$~{F}-j2SW+rBTD4}BAf(K; zB)J$TR-~B|UyYz%tSSMEeXs2XHi?oV?93(MM;m)0};iGM@el`1col|H@ zugtIaR8(R*Lml7V*qm3LZx0OtSZ}kLzis~2g&$^IOR5VKt{i4&52-zWAp<4Y1i3g| zceopD4u2F>Rh-x+b|7jVH&~wny}vW)dy&AGcdZRydyc8;9ZYZR=2{N``IBAG!WEao_&*fVk{ly8jlgE*xA z+iqSHB(nHHzcX>0TqMPvn!gC|t>k>*cRmUKPHpjd`MXF}?kwK$-8pL@d$Xu1-(xgv zjtCRR`RA9L!LY=D0N3Qe)T$OrV`JlZfY$Q&=Pxy)hzsXfiwmrt$;*Aai@3+fwR^>s z;PUN!pr$>TTKMjCc4w0P9cZjORV%husfDb2v$iYfb(8~PPkg{5`l!d{|If=vx+A!< zLG`_h_B_w;7&pYO&$|wKjr~mwq0EwPuATH4=XIn9wPbq_vF8Dwtb&in$Cq}EW~G)s z-~eiB+Sli@9GFD19*6w){FKWF@`RdZ;+|7gqZ-%Xaj`S7z1|bD{daCTDy(NjVHM`x zvvgbK+n4X&JKz8IZ6mvHxW_w*dxel5fO)fV#o;ZL;k^Q^YFE60#&{7*+h9nM)2sgFcwCAD+*{uXNz|8U1cE?N+7DIKnxP*_{g@j* zuge%o1>|#iS+r`&biX;#!kn;J(+n);`M2FNh^cP`kkylb((J=Ul+Y&9rzH^GH8#dN zxzXZ>XD*M@NwS=Q=^$Aq62OE7KlVlmlGp>Zk93$Go#UHp6fv=4qZU&)t=5pY#g45A z7E#6?k}7QtS3Gcg9kdkwc>SAWZ{mXf38}n*UAv_#=(czA6Anm={;|E!-RCuA)qG^z z^tfRUmrFr62-qBVx`8`a=`i_PHkw@H!%}#8Oet0Bl9Z_5|9QDCK>N~DD&u~(@^s5* zn*bBb?A1hUR%SG*A}Jv;U^bM<^50tJUc{$)t^L!_<1?+1mf-`xu0`NXGQw=N=H9aG zVjeRqECOnb+-o#7*Prhn9EIW1Q_Inygwk**kZAaISltn?={$$3)vZo(%a0t-x4VU^ zau8GU2!TfBLD-72oA1#|Vu`G5rJwm;f5N_2@;<{W3NQ%}yZivS1AoP7mTjKGo(q;m z*UAjvws?v@=*&?Dl7hjUb?yck&50=yu-2aL8p265J^}BTgY8u)eP#$f< zo&3tLt^Z^8ipaV}1N>ws5)b@Kr_W1W@Lm|s`?DBN$hDVZ1ku%ltdyX?6>O6#zx(@{ ztbs}160hWKKtd+syJrR>ynOvxOg@zBjCTn*pk{@|AjPbsgjpqn8S)aB5Xr~0ehNp% zOf|y}haJ5;3w#^C5*}^JA5!R#=ajNfDSY#$cT!$ees@;Pt1c?VjE%Sb1Lv2i#;GGq z>UI3zzs%5aWzb=%2ot&4e~0cdG7dY(CxO45Kt6n(QwE1Zi(-CxzJYxi2iK(zYHW1- z0c8mGgQW2N9u&wkOx5%6;2Oe)jYd*${*iU-i} zJKq72RJC-lyj1ZsGZ5A+ikgKP9*jf`MS1C<<~(1!eD6w8R>o0S{J3e(EUf*&^GSG7 z()ZuE8Sae!-<|otN9?EVBZ^J9bVsN$ULTQ4uJ4Ds!c@4^w)=HLvpZ?5MVJMN!!E=Y z)`AWAhQSka&+H4@cIxEYijQTpK;0JL01v|kpg>GxkvYi->(AIb6`pO;+mz(*S2;m@ffvZIP~d)351b8XkUcqu38Jnnz*{J+%@S(vQgwbKm5 zAxel{;*Zn%!-9CeBjG~l@tke$1#t1KxHZYYc=0;d859WqCAPdsTS^nd6bnDOxg8P< zIObk34Qs-Bi-rXqw;nS-;!5k{AE4RHo%q)+yCwtk>EY`FDW$ zKj@!UwnkA1z6=e#+O`G%pl$=lvGM%;JR|T`niiFY5=bbEZqX>4Wz49cW;EikGDdS? z1Sp5DyfcOdSUvSg%h&S#sPIPGe;b+KFB{qY!i#xdR5!LdS&KGg2x;2Wcfahh z_pZ9&n2GutWWlt(r05MuPcChC%LRdmS_)u?65IFw4=@?Mg#xqv>xI>)Aj+aVNR)u< zPA5xWW6Pl(>u7^!Yo&YKqQ)fkD?x$n=8I+(G^@{ zJTVrnnB2WLF(EcfjSrY2FBP2JZS$y`g8FQY?Kolm%M#twakse{Xt2Qht_Of*pm2?t zeFnRgR>mYCB;$3+z9+^TOjMnrM1+N`1e8LmC0&vOl)-ue=`26DB0kDM{HC$eVjT84-~DL+mr(2zee3PH1Fg>K69=oTB-m-EfEV?JT138+iW?o4mr zbED~&*r~h z5BMvYCgQ|sY;b?;tcr$6(O<}qGC)_5b7wzZEF`cM%RPIcmZCJW(p$)vix3HP59MqP zyfBWLR!WkN=C=jEjW##d`UKBc$yTpLy*-Z#WWWs)uN@rxdfk9iAP_+R7~`4p4?bSn z0_?l)l$D=3q;Qsc|9A@M1Y&DdXzaWTcBpH==pzbM57WfPG2&ehyhUXvU6*J$m$ez) z82$HgPODr<&(kSi*EzUv4F*f6 z)L&joXMF!XFs$E<9!c!edmwraK^}`T5dPL+u-#J#68R$gO)WxGWWqUFILXKQWf!B* zO!wIDT#7W69#H=Nh;@CaO@CRND1deR(Z{2wR^{W&q?i7}v7?PUUYKd;6UY>viU3 z9}s?W_M61cXs+_JqjFBgg?PFSKSEZ=qV6$Y4z@T!B3M6+ww&E#$5}bFm@~k3??k4` zYWjUssYXen9yV4A)9YK$eI57!H_Qx3*2)a)T3H}RxYIGv$wrgjW62{TYay*;{E700MTi$QcYeeA`Gx|9r9fjAc}0=R`6)*}QGtlrZx~lS(T- zYxEWOV#eZR#|ONF-q&|f`zSBvak3(07HMK03JudcQw7@5Ff1Duj+8L4Pxn`GUinA! zHHd13^07a$@ZaG9Lo&FtQxuhqpKEL@j7d*2&-(AFH)}mwriL$vxU@91I+zre2$)9w zU!SezT}l+k`}eZI2^6G8<&U`B?9!7f5pI8f`Q__kJ=$h}^q5;rW307~8H=A@#vp?= zk-xwHmMR=CF#nSUpwL0WFgj1drE)chn6!1iBrf={yNt*{@=XVn(;Z2q5~F53^!9OTyY9q?t4HmwKPgtl+#oDA6 z!g`xwIhBqF5!RR2m$Vct5A4@Gu)F*f?6`s6{Xbyr?H#M=dBmp>HQp%TorLK_U6=uBYa_2ahk@9 z2Oh~B9{7o0pRwe;`i+j}Xy=GVPhn!qftxHwPg&zKhCmVv6krKFQU^Ku{>G;cBh0#W z1A|e*frP1mEMO@0EDD+K)& zOS9@r{l$jd>-j#0@>!=adhx-dgYAq6(IjtnlULWJm& z)rhZ`D?tRUJfM-YWs8ydXZ^2l6$a(un-tEytvO{|Deo^QMz(uQJ0Yg;fNm1Js+sen zGaC1(_uD#7ALBkFH9hAZrtd%MXF93>cAM9Qq@-~G zmDU8J(wk6`kYeQYo5g8)rTBnv2zycsd1YnAF(nCTxel1c(J5>J_^2iJMx$zMtY#A; z^D?hM0%4*ca&W!$JB;5BBb=^^+eSZqR1gjwi`0j(8cD*SgpveM=jk$CEwJ5EBUMs! zui*as<^LX`4!_F*O;Rq;;@Rd#1e4oA3_p4(B}ncw?H7haQw~)uKVSi|eF398-{JRj zw+GN&&q$x?Ml@%tyu}S$4fG4}9lh^I3@c_5`M+dOk$;@^ zJABCi{&V6i5O2?r^y?J`4NzMoD-7D$fYh0S#6MhQ!rt0CLur;M!)Lf+wN4Uu?Y;2% zo)ce+nZJ-g>tB^w$vk%Cg-TiFQ|7A(pcIoJWW22_U>*oZVmoc=w4qG?%gu0b{f!G^ z|NU8Fp#4WuOS&5SSl{xHmaL{Pvv_Je>GMynx`Kc}ev^U`vFLbD8Cf(v7MID#^u?t( z)m+siD)n*OU-KT-)ix zgoIUum(w?_`_UlTS46!pxU0h1|D{=$T8*&O;-;iC85S9tAm_3N2k`g6N=iy*XuM84 z5!H$M$}v0Xa^ebv^1uLaPxvC6@5Ot!9YWXuMrca{|Be_P{lisK9U$<@*%>}6y}zz; z_=5dskX~18+BYX%4D}k#_n{0KsUx~0$7Sc|=Rbyk(6HBfhz>Z`OA0jgrQUbW}rZB=@b znB`Jq6)>UczyK4Pzw}xo0#2|`q;(N$boFf2)gf!8``Jc>*$CrOoNTsv#kzaLF-vL2 zHIe51iiE5zF4*n~TW#lGuum2i;osEx1rAh#>5xrKK%ZdG^AU&}RE}f4yJgbU?H? zAkYc|&bb4Jt3EgXr@60;i}LH*RgvzJX6O!4I{)Gz4HBYsseph(cStDRjUYL6cXvn) z9a2Mg*GQjx^t|UhAKowLywCePXMdfz_r3SpYhU|X*Ip}CZ)eO}tjjVekP8)fU|S3B z$CU7??`LX_^KB5wl$xZH<%|${Fi7Pl?DqTHwgF2KubWeTp2Oc5!(s`MhA+UswpEBX zfLnhF$m&swjz0m5h#D!@6{;t3B`Yy#<|r*h<{nC*5p_T!AVokQp6k?T2klSPP#qcG^qW>85d*eG~t@mXBUL3N^ z9ROg6lCJ$;C>lntN)Qw%2Z}1P?TcnUxdPI}Az{x?d7rN@{LqkO6SFb0SG{$%tFZ#s zwNrn@o+cOO0@ce1-3I6Z@_?8n7Bw17hC@`)` zX=_uB0j+Dq9PHVV{e~bx9c~@&Rd!LicNG-{0g{!q%bz?7PMW?d)`u}F z%4WTxD7WTU+`mg~=6-N|ncT`udi>Ve_&CssjS({UEbn%Vzet2MoMfGJqnHg*tl5=v zCY|sswS0!xApf~N+P=HC+%d+fSYRrZmeKJFw=gIzN8Gps=NAGHmZo$i=X<-zfY7(>y!^|Z}7>{ry7Ue)X(UPsOO zVy&zO@%vp?DO_gHsKsRvoHpFcao!W&dwVnL(ISLrvAE0W@_E~Jh_OP>ghL@SGl5Hn z_6MtDO~8A&T^_&vD1*g8*J8BT9#oW3Q1OMlMgK0BrjEaZi<4>}ElhnN10)~PBpJtY zQBY{{Az7YL&@e7KHeow~p<7>v6}Qxsx4pTR-}ZpPVzg*f-OPE&^x9yoiC8V#6|#`I zP&<~1FAu797_qumiB5o|s9W~n%px3^;87mlH*!177h<*-+1h*8Me&b@t{SapS`wak zY9>$rdQ(}vqtMLa#o-EtT#JrTE&5!-FV4~wL&VvJo}@uNbuG3AhnY+3iL2mW|L6Yv zHc8NL(Nw3%!Q$LJKYXKcdC}6Be|p#6VoGd*A5J>toO{Z!4R#1Vkv8qA$=RxTkx8=M zsFutDI;`49-9wizL{VOD>{#_$ zj%_R_6Ayhmg~mj1A@a<1OqdNzX_HM~QCW<^2-YhnQuNZUwRVa{(8NNIqCS7Wk+t~g zt~r-A*h?`)CWoabA$Lk+6LaNe7#G@-WTLCVKS#%{g>{z1W0;Z!`^RGHY9Sh8u6I|v z;2qX`>k&9k%_}N@!3~ez7^Y?td0uBUO0JilOZUdD?4EI`o%tFQJbSw^#$=|9GTy6; zB^a=3kpY&!ma?L6_X=F2$>!IrAlbeb56K75TBrok6;WYl_=v2hA}@^jKd$Vv56&LtPR;$nLPoR|3hh`R@qzQLmGY4XFd zowiJ9sbnn40@sr*R zWv_qcU@lX|{_5=SI?LBjU54m!uJQHSXlyB?!yf0W%X;io;ygFZOJcSN#tlxWLbD3# z)Y3@I&77@e%Jx5K+rxyf9JF=+HvZ@iT^yvSI>gKml>E)Tn6w(fqv+jMK^SHXfr@{) zz*#^jWPHqeV)GB(ov@2cAB4VVg|@S_H3fZENM=Dn@PxB3QQpPj4nt5zIBS1^dpX`h z)<$}OR-ZrZ*n(Qv{q3oO#)mEUs#nnPn5D{NiKf?$+3LAv?09OXx?@_W9LOcaR^vnB zNp9G5Dss17LD0%W<%@Qa<75Xcp}U!MlUtTN3%fL43ET!aaVEnFh8lg;7klfvH|!-7 z`TcFxB~%2UFH1*YD6B&~R9Y_a<0-%6%eW$Vo823=S|5~^2sB1L2n#jmjAl^;?Di_l zukfHdBc~ioENsCY(vbnt<+pn}jp9MLugfk}zIuB3&KB~Z0(V;w z@3s#0A+BL(`bPOoqZ7pwUVXpf%F%P!&+gZL_Jb%~VLNR3} z;m?UsV$NrcxCW!fou5Vjv;i*pNYXx@5y?!^pVBW9rsjAv8tZH$v64~K-;wQ)zFL=4 zWzT^wX8q=0j%kUHQj+^wpN5q2)bmOv&5-hJxCzzET@*ileC=X98_Xa2?vBvFdH8}o zk^3b=VmqAI>5lS&o{sISJW09`tXy}_kdfLIMkLe!6tA{ZP3g3Vn0M; zuyB(Nd4m&;Pw8hw+l~g+vD)|GvP2a+7H-1jcAw-5oRt0#^_-=~`>!#p`PUbA2?Iy9 zlG1aJX}f7|q>e6&!sNV+saY;;1P*Rft;3e*{t0f&Prc^9{>jJX){On*P>P5+w<#d% zRxInq(A%f(_uQ*#3&FjfnU4mplRd41v|wY-gDk{N zl1rhj_By6w+*}JYy~vi(Aa&zI^5?{9XQjBh>lEs1b%*kDSeR|}u4vUC$*<1S>Zs!a ztLmtP8}3kC?w(_N1{8N)KDIw>kNn;gZb_sLl6ob*!kDcJcHU>JfgrkRRAel)RzR#t zw^ePZt`X`t6kKUX54FH&T^^-EqcbpSp8T1qwFcbij%DwyLNLA#c#yI;6N zUn56S9#%hBRLR(%02`h@#bmhf@)?K3S7y@vzw0L#UGFs9*EcbroNp#1h_yZgSb*DE3c3} zzj(<9IR;k+!M~wOU3TPkVv8RTlLD$y&9h*U^iRzSZ9}^lL4$6FtAx@wyq5i+^0d3{ z8nSjdw8}YlM9c+;fkW&^Yt0|?-|nO6+L}cR)+qEMBV8P^b5wt0-B;a!w7B*nbZS0& zm0K;f0GkHylf#jeeZ6bSTgzDAxVX5F>?p(R0sOQ$R93CVJ)H`J-?%xnit*6rJGvtixq2}%Riz)qc}h8EL$36zg+DgxVXLa{>LgE?ku8cc}s`Nksxg8C%VXEg8sI?7XPGhdV#aL zCn)zZUdf8Zn@8$Y38B{SHBe{_y;SURH$qn54pvPEXqSKUlO=ZOQl=Im4i$!K?P5xo zq>PJeU|UNRm>XCW1n)@XFl?;Zg;F3S7~F=t@t=!!_Pxi$v$DjI+z@HV^rCS$r^H!H zlW;F4%Y}N6s zYFH<<48O8U_&Zp^m2Wi5qf6i7W1900iQ1{=k6TsOz25dP(%iW8jqzhn;2xnOXHPY&xYeX7ig@R20y*Dwmp4p% z4!3@n?#auof+gXt9IT*aNY1=}|JuXsT7H1Hm(#Jsl-=h1fwCv|cIT~{0u4p=L}9eE zp%5mXDRbJKs4QGIm$_6Bk7&n>fyOmmjIjliW5j*JIYl`p`lUxUpCf;7$bWE6<3Lr0 z5mSRu2fGrNR!{zY-vf%Yt5ZDnxzb?;)mRqZn$PR8_G zVzttm3N5Y%H6LDZpJzYdmP%0RS$C*_$hV`S0DDc;hQE2_N4R~QsIm?#4Q^1(SoER7 zrl5EBrK@{Bj2y(yber9ev+fSp@Ds*UAE;;LcTbFDSb$4u6cmb36hA++v_P!Ov5V~c zkTiMc^K0sTnM`HhZ%t!pjSfwPnTdEn&s5)M#!@>PoFq+ZoF+BOduka9`h9`i!LwRc zt#ndKtEy^+z07&~-IQaWeEgKvwZfem`s_pzzeWI!_H-;eyKPb(G+N_u^^cDC**B7e zcs7|Ib6m=#Qif)kYs4q-yr62UwypmFhjBErt$g_qG2UX@MJ364Im1171<)b(>Oasx#B!KckvZ1zT3B1f885rJ8Oe=#8CQ zARbSEdvbE0r1Ljce&1O8EzACNJk?1*o>Y(RrX^PnFXm!+mp4uP{Y{gw?Uj~j#`Wn( z4)1G&3KCbB9}~X}?4Rv;YAZ;2s)WSx6w`7~j^1jQPf&0dw~$!2JbYoFji_{p6ZLL- zizSrsKvxoOda=6axxF?|=HSJ0Z2|`l8qa)SG**@;jktEg8x}luoe17ow??zikk@AZ zRcBrMN^>KL$8~AmI6RSYK!U&+VhNi6L>1ng6yC?fhRKe&v3-B9mVz0xG+%})S_KK; zxyTVyM;LZ{4m<$hxGI6;9R}VNDjk0x;js8PjL zrCIo1r0Qm{v?5rxFc0)ad~Ri7^f+;&W!hp@kT;!!ANR`r5sB!)+B3Qis>zfi5;5aY;l^Us(`Tv7b`Cs(JzniQ1*FVn8jJ2WZ4G&;YQ>WdO zaSL2BHM8YpZj}P~qcf0+m@HTX$+;;-`o{)p%+y&4IO%()yT0clB3OG4YL+eTma?@_ z&_Bq|V?5fqT{K3=9M9PT_%r*hq~; zb~J&l6{>-@H>RzE^0#OkZ!fkSytFUBN(H|9sb9kt{$S>&Yob@pYY%mSd?M>3MR2;a z<~rQD+srsoAm;Cl?AIUN!7C1mQD+F}uC{h(OFPxOGADKx=X^9O=!F#XLdPh#yPBN? z#Gy#D>|U)rdH?uofbiYNW?B{UYcJPKzVsx1)fBLZQt0ke!XfjTi;@!oQ1%!l_rK~pnAJsSlZ^6ck zj9^E}cQ>P?tuharjl#t+jt=n2I62=1eVEY#=an`zVeRYBuKG1^3k1%Jq3{ef+x(t2 z@HjEScHuwwHyh6}7jr`4>1Cx``9gLjR2+czBeO;DlhX}I&~TDw`|7cpmz1#%Jl@7N z2UA?2)N)4+8}}O4Te1|8tWlF9hXl$QzrUF;s~O}IgQZJu6XAxN2#Hmn98a^Q+*5?3 zmm#JmWFBj%1Ga4U(-ayG(A(qQzhIbOEp}&MgldaOx{^@?XPv4wSOm_~1%5l!X<#t3 zzl&{Hu4Em9s2eYZhR$nZgs}(MjrH#|1Yk~s;R&p$X z6D3E?G5Sz{2j`V{Xc+Q6@ftS zTI;7x|Bk=I5BW3v5&TfZQW~o@HyDS_pFypu`Fwygez&0 z>YOuv;zl(%BU?kGNzhXC6>3rY&Px&Htn^VrCX}k~d{#hq81W)U% z8!>tr%HTVY3zP7=&T#Uo3pccNa?|ofM8{$Q)~|Z6CdNwAP;3gUYgbcXT@QxhDga`{ zU0)i#Pm^P@>#wUbn_>E~(RrH1$vHS+8y5zG)d9vJ0Un0VC3nl5H#4HV6cA8kT*HBK=VOJD`zf!t zYQp!a3D?7r9u0R}jq9Lmrt@kGKame8OS;ViD)vu1#okh{kA=y9rzJdgDROV@ZJ_<; z(MlataM0ORg1aX z*f8SUEj|u16>#^Rhxv71bxuXe0OiFwx=EkmCZU5Jt5!r_P@yHf`;_~e(Lt10Qe$cl zc_@Gqh3#83g~s(TqzNlRNxB^`3dxkyz#Lgar6go*S|Gg z{i0KB*rEG8-t=(Pkgs}?-Qzfh7bv-mZNV2+M|c?gO5j9kT?UXA34i_kzb&6gLFiQ) z!i~=4HsbgBi62*A(O9GlFA;tPY%aBge6Miw6{B8T{Vsl-sArLp29zgxk|S4K3(9cO zLwd71@+j<6ygJ~D0)yXW%!$_6gJXpN>i(zjFU^-(3CRz{YK>OYdTl>N2*uyh`oCY_ z&U}kEBnxY6$_e1hfeR19Z^%?GFWV0%3}^NPC}}KPe1{{zA4CWTUJn+FxK%d?(Rib; z2C1{k9Q54Jwgmr@0`*Qlje5$7(b8Jfg=kFSyNcV$%WJ|5Y_CTNh@8)2-Vu2^M@kqa z^Lu(uB;)AsjQn-j9GMGS?Bp8P1UXGJwAPCY{^>{VE|KGto`hv~<^lNLix1Ok2aN9r zz06m)F^2^^%YK8%LVqoXC8ct%r>qsYSc?mliNH$lvCW99?&UdqH^>F+aT!`cK4t}t zZQW2oQEyu^-p@LIs*Uy8_dM*g_+moyveQ3?)VQZZqkz!8<#9eQICc*}S}_pn)cJrw?^P&r5&25|OB4tM1kRA{|8W%PZzl-< zOW~#Z$Vj4NGD%lU3Xj$C*alxW-;v^RzRyFbU5W~%zym%UKPJ2;TPp7MW;DoWFRB_T zPJ9ZZN?6s1qh!d;jR0IrzYhl2UZ>eAUWfg`+))+p-g4&Sg4PrIMb|UjMcsi1fY?jT ziJ`2FeAP*yYhJ$Ch_}EbyrVFJwf~)Vy8>5d4TH3o^1J#ZkyVu5UeP2e+LB*;1o>|J*}J+}CIoH##( z+HVnVY_*$yk4|Izd!e*k=$TAD6A-Gl2_h^1e8r)e_#AJ6;}=^T7N@%=Wbv_MYyhe$ zI^qtp8rNxmbh&7Qb`i9dk96@BQ5}etkZt^yA zAoLf{NwMr*)OMEV<5C-7Ir{k2@s~h%dS6)CXBbP0ypfAf;13?azp#bK9>L2AEhVm3 zx_dsvH87aGAx%Yd>KwYHvi;HwSmVzzW^|L@;N~OhS(>o@VDYMSW9&^8&!>AEqD7)@ zi`d{Z2R;B2CNb-Gw&32Ud$u~c^wi|3IutF>oVYf9lbCS}WfcJyZ66bpzDOXhw*a^` zyxOR6;)wufZ_d_nvVhm(l@PBSMP0zr9 z&HS|!f}lvZH=W=~P5a>+LFm;_1&W{ulLvW+yMpkfH$ShIx$KI7p=(n=X4BQfW^F=5 zd{z@cY>t4?7HOaUfqFb=#l8R#1sZ}=s;GmJvq;U z=L5V@y3E|K$YmQz8Kz5i_~=dlCXpW7!-O1~kzXx%I@9@qD+=8FmORzF8P8y=h>+UI zQ$2!iqqi_Zr#in~!r2j|`~WH+*7P+bd8$m|(gTgG-*R^2tj)I{%L#WbBS9DE?$VCb ztMwS1KipA5owE0(ahfsl5Re>FBCXml);BxgxZf{jrY(QffjUskuNorh7h9$%>}W(* zb7j@m55lL{aecEGeJmV=PzCJjnLN9^!l4vOTwbsn&4SPG&La`p_}#LA-@WqnkMDy5 zaXJC(BuzOZ5W*{IpGZi~+@8S|CetN)Lv`Mz`%<)B%@q@HWVhdi_k4h(5&$ICzWpJ2 zi^PBWoIOaGR(XZs!6rC4i4d-}r|&s)ryX|F1go&qofrC3szBo$0T?!*UKGkWA9tR) zjHWt`Rjl{SI!@g^KIM=H~lU*@H&@8 zZ_h#1I2yQuQWRrPk7{d7M(`xzJXT%#rq>Ih*PxZuzjc5ILS>!?wG8+eSoU( zkhc>-tea7^_!-+Gu+@=2{BF8`p44b0dEC8X z*0|hD?|zJUuWDn26kO$yfUQy=L>MQtAak|k*t44fSX zKCaKxzFg$-!4qcb2d3nI5g+cNclkk}#`%tbSnTbX)SD-O^N=a`E$%H$ zJe6mNFB10s$*|u{M}OQe=+WHF+{EJK&q!s!;MR|{Hi;sm>GzLYE|2J2-q~TOp{4Da ztJplMHClw=9^aFJHKlevM`6l1Ajs2IEhf-i(GhJwO z;MO)$mSbmxA3+ZS#K?|`4Hu7Ausv%SbMof{Q0;j#5EY1#QT70<8RYI<@ZJ+go1d0B z9@aTKfRBU!G<&||4%M)u7Ll7cRgUG3O*Go-3$Rvd2ifR#SDSFTy&2N5eQe?g58tw^ zYum<)Q}g}7AXWfv+A(yaFAaCHZ_O0a<0PrFbMKLRHR@w~>NI|#_hh)L{TJO3a!UWL k3ZZ`oHvFqoLVIwB@40jz>% literal 0 HcmV?d00001 diff --git a/presentations/press_course1/press.tex b/presentations/press_course1/press.tex index 1c14d84..785e997 100755 --- a/presentations/press_course1/press.tex +++ b/presentations/press_course1/press.tex @@ -25,10 +25,9 @@ \usepackage{amssymb,graphicx,enumerate} \usepackage{subcaption} \usepackage{hyperref} - \usepackage[normalem]{ulem} % for strike out command \sout - - +\usepackage{tikz} +\usetikzlibrary{arrows,shapes,positioning} @@ -98,7 +97,7 @@ -\title[]{4D flow MRI in cardiovascular diseases} +\title[]{4D flow MRI in the diagnosis of cardiovascular diseases} %\author[Jeremías Garay Labra] %{Jeremías Garay Labra} \institute[University of Groningen] @@ -118,13 +117,40 @@ University of Groningen\\[0.5cm] % \onslide<1-> +\begin{frame} +\begin{figure}[!hbtp] + \includegraphics[height=0.5\textwidth]{images/aorta_healthy.jpg} + \end{figure} +\end{frame} + + \begin{frame} \begin{figure}[!hbtp] - \includegraphics[height=0.6\textwidth]{images/aortic_stenosis.png} + \includegraphics[height=0.6\textwidth]{images/aortic_stenosis_1.png} \end{figure} \end{frame} +\begin{frame} +\begin{figure}[!hbtp] + \includegraphics[height=0.6\textwidth]{images/aortic_stenosis_2.png} + \end{figure} +\end{frame} + + +\begin{frame} +\begin{figure}[!hbtp] + \includegraphics[height=0.6\textwidth]{images/aortic_stenosis_3.png} + \end{figure} +\end{frame} + +\begin{frame} +\begin{figure}[!hbtp] + \includegraphics[height=0.6\textwidth]{images/aortic_stenosis_4.png} + \end{figure} +\end{frame} + + \begin{frame} @@ -145,10 +171,6 @@ University of Groningen\\[0.5cm] \end{frame} - - - - \begin{frame} \begin{figure}[!hbtp] \begin{center} @@ -161,6 +183,31 @@ University of Groningen\\[0.5cm] +\begin{frame} +\begin{figure} +\centering +\tikzstyle{block} = [rectangle, draw, fill=blue!25,text width=15em, text centered, rounded corners, minimum height=15em, line width=1pt ] +\tikzstyle{block2} = [rectangle, draw, fill=blue!25,text width=4em, text centered, rounded corners, minimum height=2em, line width=1pt ] + +\tikzstyle{line} = [draw, -latex', line width=1pt] +\tikzstyle{alert} = [text=red, fill=red!25, draw=red] + +\begin{tikzpicture}[node distance = 2cm, auto] + % Place nodes + \node [block2] (mri) {velocity data \\ $\vec{u}$ }; + \node [block, right of=mri, node distance=4.5cm] (NS) {Navier-Stokes equations \\[0.3cm] + $$\rho \frac{\partial \vec{u}}{\partial t} + \rho \vec{u} \cdot \nabla \vec{u} + \mu \Delta \vec{u} + \nabla p =0 $$ $$ \nabla \cdot \vec{u} = 0$$}; + \node [block2,alert, right of=NS, node distance=4.5cm] (press) {pressure \\ $p$ }; + % Draw edges + \path [line] (mri) -> (NS); + \path [line] (NS) -- (press); +\end{tikzpicture} +\end{figure} +\end{frame} + + + + \begin{frame} \begin{figure}[!hbtp] @@ -172,6 +219,23 @@ University of Groningen\\[0.5cm] \end{frame} +\begin{frame} +\begin{figure}[!hbtp] + \begin{center} + \includegraphics[height=0.6\textwidth]{images/p2.png} + \caption{\footnotesize reconstructed pressure using NS equations} + \end{center} + \end{figure} +\end{frame} + + + + +\begin{frame} +\begin{center} +\huge{Conclusions} +\end{center} +\end{frame} @@ -180,6 +244,31 @@ University of Groningen\\[0.5cm] + +\begin{frame} +\begin{figure}[!hbtp] + \begin{center} + \includegraphics[height=0.5\textwidth]{images/catheter_curve.png} + \caption{\footnotesize Two different methods for pressure reconstruction (STE, PPE) against actual catheter data in two patients} + \end{center} + \end{figure} +\end{frame} + +\begin{frame} +\title{Conclusion} +\begin{itemize} +\item \onslide<1-> New approach for addressing cardiovascular diagnosis is presented \\[0.3cm] +\item \onslide<2-> Consists in imagining the velocity of blood flows (4D flow)\\[0.3cm] +\item \onslide<3-> Main limitations so far: +\begin{itemize} +\item \onslide<4-> Difficulties in obtaining the special type of data needed +\item \onslide<5-> Practical resolutions are far from the desired ones +\end{itemize} +\end{itemize} +\end{frame} + + + \begin{frame} \begin{center} \huge{Thank you for your time!}