From 8cba417317a54bcf8e264494bb790d992eaa9d5b Mon Sep 17 00:00:00 2001 From: Thomas Dehaeze Date: Tue, 14 Sep 2021 13:46:28 +0200 Subject: [PATCH] First LIGO review --- journal/dehaeze21_desig_compl_filte.org | 5 ++--- journal/dehaeze21_desig_compl_filte.pdf | Bin 1124696 -> 1124654 bytes journal/dehaeze21_desig_compl_filte.tex | 6 +++--- journal/ref.bib | 12 ++++++------ 4 files changed, 11 insertions(+), 12 deletions(-) diff --git a/journal/dehaeze21_desig_compl_filte.org b/journal/dehaeze21_desig_compl_filte.org index 5824965..ba1fabb 100644 --- a/journal/dehaeze21_desig_compl_filte.org +++ b/journal/dehaeze21_desig_compl_filte.org @@ -113,7 +113,7 @@ Practical applications of sensor fusion are numerous. It is widely used for the attitude estimation of several autonomous vehicles such as unmanned aerial vehicle\nbsp{}cite:baerveldt97_low_cost_low_weigh_attit,corke04_inert_visual_sensin_system_small_auton_helic,jensen13_basic_uas and underwater vehicles\nbsp{}cite:pascoal99_navig_system_desig_using_time,batista10_optim_posit_veloc_navig_filter_auton_vehic. Naturally, it is of great benefits for high performance positioning control as shown in\nbsp{}cite:shaw90_bandw_enhan_posit_measur_using_measur_accel,zimmermann92_high_bandw_orien_measur_contr,min15_compl_filter_desig_angle_estim,yong16_high_speed_vertic_posit_stage. Sensor fusion was also shown to be a key technology to improve the performance of active vibration isolation systems\nbsp{}cite:tjepkema12_sensor_fusion_activ_vibrat_isolat_precis_equip. -Emblematic examples are the isolation stages of gravitational wave detectors\nbsp{}cite:collette15_sensor_fusion_method_high_perfor,heijningen18_low such as the ones used at the LIGO\nbsp{}cite:hua05_low_ligo,hua04_polyp_fir_compl_filter_contr_system and at the VIRGO\nbsp{}cite:lucia18_low_frequen_optim_perfor_advan. \par +Emblematic examples are the isolation stages of gravitational wave detectors\nbsp{}cite:collette15_sensor_fusion_method_high_perfor,heijningen18_low such as the ones used at the LIGO\nbsp{}cite:hua05_low_ligo,hua04_polyp_fir_compl_filter_contr_system and at the Virgo\nbsp{}cite:lucia18_low_frequen_optim_perfor_advan. \par ** Kalman Filtering / Complementary filters :ignore: @@ -535,7 +535,7 @@ A more complex real life example is taken up in the next section. <> ** Introduction :ignore: -Sensor fusion using complementary filters are widely used in the active vibration isolation systems at gravitational wave detectors, such as at the LIGO\nbsp{}cite:matichard15_seism_isolat_advan_ligo,hua05_low_ligo, the VIRGO\nbsp{}cite:lucia18_low_frequen_optim_perfor_advan,heijningen18_low and the KAGRA [[cite:sekiguchi16_study_low_frequen_vibrat_isolat_system][Chap. 5]]. +Sensor fusion using complementary filters are widely used in the active vibration isolation systems at gravitational wave detectors, such as at the LIGO\nbsp{}cite:matichard15_seism_isolat_advan_ligo,hua05_low_ligo, the Virgo\nbsp{}cite:lucia18_low_frequen_optim_perfor_advan,heijningen18_low and the KAGRA [[cite:sekiguchi16_study_low_frequen_vibrat_isolat_system][Chap. 5]]. In the first isolation stage at the LIGO, two sets of complementary filters are used to form a super sensor that is incorporated in a feedback loop\nbsp{}cite:hua04_low_ligo. A set of complementary filters ($L_2,H_2$) is first used to fuse a seismometer and a geophone. @@ -843,7 +843,6 @@ The source code is available under a MIT License and archived in Zenodo\nbsp{}ci \bibliographystyle{elsarticle-num} \bibliography{ref} - * Local Variables :noexport: # Local Variables: # org-latex-packages-alist: nil diff --git a/journal/dehaeze21_desig_compl_filte.pdf b/journal/dehaeze21_desig_compl_filte.pdf index 1fbd3cd967c818862043d52c87448cf5752d0687..7e092132387353d9880408d42beb366715a53d03 100644 GIT binary patch delta 23089 zcmZsiV{o7ithQ_0wr$(?)?3@=wokFOvGsOq+qP}nw)wsDo%8cd=E_X+JNHa7Nmk}0 zs^%o(slZscQ*$3sXaU_x&nzCa?)T~w(|5s4NfML^>m&uL%vCc>kBh(K6wnA~4jxuF zBsKe==gnOb!HQJ={_jsl%`~!wYNbN{lrY@yiY?xrSpNIZkB_T{{xXkJrzA_(4Xs=D zcML2uY^q=sovl?|-S*whRa?#O?ZVyBe^ygAZ<}hR*R!QhJph_FXet2CbQ7V;^P}b5 z$AM>^Z+k1I>Ppu^9G%C%>kCza#T9i=8``))R805NaXw98WS(Y_lnyRgjyv@xiH4}i zY{i6|4pyu&J;nu@;jwj%$Ve-eT5>h|**qE7Yej#^JN`8K_$P6jX&`@;smgF}&3CQ9 z;y@ zT+m3}5fGeySb7>Zk<8#gy5spLuWq10e_p~ObMIxJ%N^F`A@1p0)3KIWWqrX?VK1m9 zn#I`UzmLs6N)_lV)ljY&e~8L}4Mm5|I4aT!I*HH23gDE0|DA`mq&re|9R3oJZQ5k~ zYJRB}Q3r$67i>5L<$G)j4u)JTVMPA(WoYH}16%x3RyfP@wu$Dvh=_#e_=qd}=YS=}#cF>F`hL!an+cWyt_?^hf=y9@i*=D?j;SKLF*D-)?NIkO|B2?*tBXC-2VaP|n%*SJns4`J}>yaIWI@t$gp)mMkOc#a2h&){wiB zd0*xD+z*nx6%gGh7%`x-9G^t)Nk=u8VIz58WT3<|OrHV5Y`w1aZlOQf{^l|c; z*G;ZSsOS7aS*QZR`qrJ80z<7IjS4K5DIXCh4R~|4Mm%-I3z@V=%1Di~If+ttGiLMk zdEm#;Ol$y)_R;1N5Y5)bK%zL}b(82M29hofBJ!j6)f~?vUWf6y#vl8BNKaD%nC3VN zcAaK!RE3dnI>k>rl)#x*;8@d@pbUO+vEtycF~8^Zy88Jurr?k#-5Qt3;XQUyKjC*g z0<;ss)&7rMR&l9xX(^Fxa&j+%zm}p$9F{D*`8b_&bytFf{_{8Ll13AnIYO1L2;3=% zAt4WV_%`hV#-n3)aR%x}xD{r6fGt8a`#YNYwYiacbN7Ci%*!aVMOG|iz<_ZDR*LYs z;k;7DWMh6)O*w?e34vF$74OV|bt9P8@_l@=WfypGz|S(6XFkE;}b^rBh))U>Z>@UE$#_%W3p8`H||Gt@{%;H z)L^R(^1u+J^S~Mu!?GIjN`*5DG+5x1#~_B}i{!3X)i{r@7Re2s)wGRmqM2jR2v<_c z4Xa|3ly)u5Rbn`rR=ukNga7#s4TtRE6J@%93 z`6_F_IjMABMXP^x@i199A+(NZrKZ%?_lgMK8i~YXKW0v@n1HyFL%y9})4UtwM5I6m z+{o=f07h3FKl?me7M_;_r}BCwMlkyPwbvilPM$ZkjZdt6!NOK-O9*WboEUu<{;a+~ z(f`8d{Mgi`7%Nf=pk8L(wGPl9B~Uj)aVk%uj0V%rlr-@HnvtKw2^O1NwkvUL(V$iF zeDk~SPNC$KrUHzJt!;9)^hb-|C`Jd@76T%{&ASejZ<{}5A*p^CK2j ($cl|h`^ zE^|K6Y!bJt6Bi(4Yk^`w0OK>?;BJ?VGN3Mesk;`2a_IHb~huXpiH{xVrK#sp$5+Hq8onh^;U z%Je~=?#ou-A3Tn5_DusD1Lp@8S~U%7DTY(=g?73$8(8k8jiqO+5#)$ieM*alsG~;%H9D9cegV0qu7L@`$&mQmav6s`-^vMEI?2 zZxp5j3IE3OqLYy-?bGD@>W(J&yE;9rvxQ&c3!m)o-|wdsf^<3SiB~k~Xy8FC={ls( zi&p~(fF)&=NN`9(PxL4+QDB4^R=vTn4Ll__47-orZ&`gh=B%2Mj+|=^Z&6{9nQa<4 zQ59GK?d|+lX+`+s?rLIq@U#gB=AV|kjHmTa;b_L>&Vc3ICF%1Bn6fjJsG-Ov+ooR& ze)fKEQe?+;?Av@+2k^L44LcVHd=bw;7hLF5k`0y>Px=LndIv#N9K(zteVFB- z1p_7DVNyCRs+#24d>^*=<%{PHcI&tr_Y?l^{G%v}+{1?9ks7UO58vont(&uoM{OZ_ z+IWAd)qrtV>>yWfB|~x@Hmx|dp#7@2T}>BOokF%ZyAWDS=Cg1F^~yVfZ7h)2Vs_M^5oA9b0<$fR$_k7jCy{*?3!TIV}*K)LClf}|8@IEck%>WHa7aVGR!~R zUsCh98WUD<)iaIvv(RCP8_ZQ<^01~OB?E1>f-6HtA>XH6W_zT4!w?UY%+b8OVu@z7S&=)?6(C2aOrQSb*~6v!(=jnu z?=FuFI`mZM#CSxO`AhyI;5P{tZ}IkP+vW$`Fy;I26djzqY;8Y9?LpJtDx z1uCwi9H&T^==AwZG>S2b4BK9g;gv)N3m`VEX1%h=YRMlG2!<~NO4H;#{JZ!UB?rcr zB{%e{4HU0F21}6H(@G<~gYJNcYHJ>y;7sAaq_RlcSH__kaPr$FaTuA=-E6{#gE-7K z4mW>tq;Gq!v>)EvK6S}9b@Vn_cwRqc)1gV7Df1&YKcioM+h;Tg!RsmPruq8-RhHM& z)==0_Dw;b!`>1=@G9uK_{5VQRU;t_W>Kf_>^M-(4#6*uP&H)$ZKMG($kv^M3*C4v0 zsLX_*t%*NK0SBs)#EJ$Lj~;`2Gi^tmbfH!di*cW*-|?Jk<)^w;h}8XWO7v`S>&0=C z{K$Vr)+sfW9Xo$6TO)U#Mawh&W5eGdok_WFH7RYKNkMns<4gkTn$u0Vx?=N{VH6ln zm4U;nU&tpQ&{Cst4G;jzfCFS_)eJ=``rt>R_p zpYxZ3h&skAf64wy*qjxpguY&$)`oyNDS>8F>Oj4C;qKhd3bp{fM1;tGv7^Q1a=$^| zQ=+nlOi1X`+q=CFsr^Mc+e`gFbgwD!78%8Bsb^T|8J{{k#1Gg0I!SSLqus!X+dz^u z6O<_VS?t~>3dB|X@#au9p2<~I$VLAL$6Xzt6F5>{XbZYpvS7gD;iaMI=ktC+VcO;`%6 zHgaN3L?I{ql_6goRn2G8_))g6xD4JS(L8=CPuW;`5UFAMB(M{+SjdI?BBIKz0Jswc zL(cj3UbndRFwdWWkp!q$pDdidH+1VRU&l9{W0&1@$92i`wQ9(s4Dtd`DaAwuGpXdw zAhU<}WZYW>wjU4(CwGvjA*rv=6c4bor%74mv^W#p{B32=mBhW{ZXcv*7ib>T+fY-; zH{D@S$@r6|OOn{S(VSjxSlx(80mcpY2wVxIWXH1+doBy%sP7b+0!7`X`9t3}st^)R* zhJ^&Ab0xw~Mh_znSIj9CFf7K(_ReW2Q8>6>rE}mkTA-Op(4`IhqJ==k8UtvHFh*{Y zQ9;4w;1&?A$}u2%RfDZ3STkrvK}6Mz z5<%fl8S9NIlhd?oSX!)WrIJh<6m(`9iA>AV_5!jn8JK=(i-JV zVsSG!jHoyA3>(Wc|OIQtuO|u3^=sL=eT%PmOHB*R^bx(T2x@!$%EG zLC9q>9yJ+-fHc-<01QOglgWs%mj48Tgg_#jpAs}tVjS>z-~CoJtUcvkRL0B{vNYr{ zUc&m39?ra~WvJNoBFP)|?>)rzi;z%`ufeV4RAts!DGZU&qJPt*MbFPi%?Hx_tb))lF<*fC!P=zXuRk{WLm#Ki7O-?ljdx$mt`N{;Z+I>fc|sL^+GIRb~zw z*9yB|H@3m40yv95hfqcRV--JbRxn%Lg=mrHbzHts-uSTjaj<&P5 z|C}9E33T@OVf;PW_IKd+9R|DZ|M2)OXw}%Tr0^Ud_~wnLg+hoDFWi7`UY9r=HS%#B zosI+P61t9AsO|i^M|I>DKD{K;{3f(^V6H}14|}|E1diX106FP*z5?!B zS1_l5A-3Pm$oPHM|17Vnzh8IrEW-SAU&&#sc`ZJ_36Z9udWJu~4Z|ii38A985yGTN zhCc_#67E{xpOUuW!*{dWeMZhaIFHEdh2P^nH+TaMNB8!KPM4?Hrwr;&qw@U32mT49 zD*$ajp}9U;yOEMk>&xuKSGg2mWZ9|bQ*%oJu1AI$%2~rV(ut^#^E2+JjQyZNh2z6a zSeXT{CR=Wg`|K2RNVRI6G5)I36X;dsGGaI3H2<2Yx&A={)bl8lC?EZpx3*<@Ja;Cat+2*l z(&!*{+kZHPC6JjijOkoltQ%~lb|!^_hz&Xi!~)Y*cP?L~@P}l~_3n!pmsp2Iv}5O& zvO<|5hv7g@jB>*^xV?SQw)aeV9tKtBkMc(s`ZaH>$a4nf>QqvAL^B~Gg?v)^a_F=n znvpd@hgG;!bY=l`11cxB1yDsbWkFw@8}Kyj%&bg19lw#3!N+hYx@2Xw<}he6#Yt}g zcX!M)O___oO;txah{l@b2xCn5<|)<7<75-{&#!>XpVa-_KPyL(c!#vg4Hyb|sWO(% z*;KUL#k*5whinb@0p47>c4d6N@QhlN4I4}#fT})ooqNf)=3}^3rGnS9KUJPP-K{~E zw6ETu<-Vt!4EF#MN0JK@tjPvzQ~{~qr^a9pYb1Rd1XLu%` zY*$Tq%Bg^X<(|-D7=4$i=fj%3Q!Ir#KHg5Xh)Xn!FYa#sO;ypq-d8c&J`_R^o<~L2 zG0M2}3AEmRJLK`bW(-vknb116?|`bEw(6{7GO|e5M0?p-lSAjlUA4DZ^D;EnWbi?r<=!xj^OVBuj z=ENn><$}mXW9jiwDMrtv_Rg;0j<0hzx%suUTGD8opS(y<>|lC?uju})*qs9V>0lwb zbv6C z5*7k_i-beA*V-$?+6mLTY+YsoolpN=^vo0s`YnqOst(T%%nr>Cvc{NRWQ4D9&fa## zKD&1~?)j1L;k>4gOiTDv-LWV9i+E5dup|@2KV6lXmh$16t53*|d9W+6B>$h3qz}`a zYeII^gMPs=nIQHlFE?P58lGQLh}==2gdqVVoJJEqK?@s3uH24OV3O-lHb=A2j3_`1 zx7IdhxxRoMHsd>@$Zx^OgeGlR3d-R(T8NUkt3uH#`7E-er$sEc{9|wTjYabM{Q6`- z8X^*9xOGJ|#P=pur1I#IQQ#T0!tTdU%<+n=9WdRTPSWj?@UK~U|O$jubMU=){Inxn!|a<8E0gl zkm*ZGl!>bi*$KF|a)Otu{3|@Fl?3-`6eaf3OHPRG!d5GF{e|3!R4wx2VmALS0VRLF}UKNEn^aEg$I9t*5^PozS%Vm4k>GzMjk7Q-hj1d@)7#dToMxUI2+{jS^A6t~*@77#g+jyYm z*mkaZO;NWwK5=Z2Q-)~;wlQEl-a2I#MALf6q!0%!;;?R6h5>h`Z~lP|Q)7fwC2Ksc zMDv#A1t*ZHE1*T#`WuEM1N)P&-=H;3PP3^ATksgc@wCnHw4h6jpxmh{Pq4JWCEkCH zf6Jp^mzZA}6rrnkx7}>m>L%4Bmxw2K@TGVLbWG(fOl@T;#qr&%F2#1~php)HfbM|+a$xQ+7IL;3V+9Ye+a8#guho{8cb4V2Xxy%8Xd> zgwIqB{g9)hx{h0`+u+m5^Klb^cVL&!I^t{h9leaKo`T~6$_|T^MN07D+Sf}=MuP9016yp@*DDaH zOr!nJsVD9v;1WSmqAy?$AdIHezZ9`r(~>STBcUFmRU|ve&_3=$lH%%1i>J1V)~3uC zd0e2-UZXf*2Mf_ki3vEN^(JBuCY_Oi#2Um@B%(w>)y=Y3bJ83u>OpVJUR$7YPB!^n zDjNLZLxu6%`&MBkuo#rMo22()bu5gyRP(g{&>2BTO=b_U?~X?RiL8Qq{n#yMG4TPFJX@A z@@=#bcrW!|VhPDHST~*jrOMc>h^sIu8T_lU^@!;^zycNu7khcoMcsW}l3>Yd0(rO7 zH#&*Is2f7KmIGNu%M^{B*%S~{zi@}tTO8b1aiO)!yyOW(lcSr ziD?)L`d&yhFwoy&N%U_r&kpT2$vM6Qr=uor@JZXI1$NjyP7<5thrKX5Yt#B8PUA-9 zuLqL{E>oY(z*Es$fFMCmO*W98=FST$1%3%FX>}WIkf^|ToXY!FX3}=Z`ZhgL%Ru99 zs-^bTu;T~=1(H@5{XRq{=dA|3=O4@aoqjPBM6cQupxE`M;jKqC0FM9UO#q83N$-A|?t=Y9&A=|tpb*_6`v(Uprze#{sC0vwb>PJU zoRgj1!lyo+&jLsb!O!}vlV9cLBJx=C;)%EW$4S(p937yv61#cqkDQN7ejp2L9Y2F( z>S)FR3{paKm3C=lemqowg&*&oYK_#s9_Bo8Q@Jq_&q`G=9Vf}K|9w(yeH;Fa@$ z1w~$tnLjh8;9Q{225vW}ce3yI=?5S8~ZMpvR zegxo}O~q}lSkjnZ?ox&y>}awfg*!+0@-Ju{;3@dA`{Xd^GWXsGPDN`aVpbdz@z;4R z|JEvF{EK?-g`%*Po%cJB!GC98ZI)bC3uebTPuSlf*SX0Bsc#F`j#V`}=WVQ)5Z)KN z*y3EH&3z1$gq|OIhmc?WLW(bM9$$hw?(gqEh7X@Z|H^t+3$na|ti7h~9Vc z-ZzxBuN191&PL-wu;6hDZp=Kyw|Dh=BjOL`hO6bEcl>;0oR+=-m9e;M5+m5{{Lpv+ z|EpELh_y~_zl`y7HB?raCP)b5o~a`QXzDY&EiDx{EqL-%%hFBVG`&2=dMq>WN8+#X zJxi6&S{tENisPfmUK-vTZ6*=$v(KAG+xOLJ_-7iGk2KvQXJ=#NINzNdAd-B~Eui)<73~3wv>i~Nyb{AEay36+mVw;@f1xOw`K{5kh&CfL))}#TKWv+VEQZyA z%~sJ-;)B{x;?{Iq83_$h**DJUa9nf2b;<|tr6Q1|geuR4@;di$3{y9?n~uZ+r1 zg#D#0?jRCa`n%<7gD>*kg?P$R0CJ0VAjqoQaZE9)37pgDtHxeXYU;=K=HTgXS<-=# z0+ccCO{Xu@MjEYu{bnd&Q<|A9`?3QOWH?cy`!D@UAVmJuAR9|QpbsT&AoX@|6sO{g z-ZC{s`JZCje(*D>)ET3B7dLn z&wqa6Z?VAau+K_fCp-%U0x}t1UJu~deKbO>M)ki4cgZsBFM_xQ*G@94u6d!eV>7)Q z)0Y(|aMJ7q)UdOrR|AF*kJfNtrhE);?XL~RV; zc2m}j0z!peNp&Y2kMSF<$N7{W?Ny(ab3DXCq^N6VM5MKkQgOR_lRpjvjYctbE%0`wOqZE*r3}gbqEW;X_l zS_9=tRege`0WNjy)z@XvzN!pax3O%HC*Ri>mo!mfV1JAY2ng6$g>E)-=zR9+Zk~YdN87Tn4g-C*Oi&yZ;!&FbS znoK5LW5En>QjpL5y+aaWB?Hevl%^^nZNYFtlh#Wu0gP6V>9U@OFHK10A-ZN*jb{aS zoPdN`O$+o!1xUY`n1Q~K#ef=Opru7QJt;jMrY5#BbOX zj&JY?!4bVJj8X@9c?|mTr+V+9cK3{6_Ax%902%hJS+@rW?M>1Xg_;5IqowD-Jqgnh zS`HNR9*2;-2ylOT8T?-rRL3fDqB#-}b&n|Mq|<8Q3Fm(UM9KG-WZ9&wz!3zTI@u54 zf?{~&$xN1-zoch;Wc`B)I|^_>{h7%J6U__~U;`KM*#sJHQ?kH=z0PHX6doe}(a!ZM>xe1AI}bvPh3I%v66b;sXya6lrr%>{N#T{+ zIQi%8R%PLw(sk*X@9fCNF-?@2+iUK{LxxmremuSP0q;UlT%fx3gzDEaKy_So-k;<( zDkjl3Fkr!#`KT({xi*m_Nkv*A~yK5nEUwW)S4B|?3ZYL~Iy8g6|v7I0j3b9V2pdYhDv78k)@=91$IQfpuZ?lwAxJi|16J2TZspSVoN1>P-DHWKlm$ZE zN6@XhZ?C#M^X|v$z%l~%EnS?yerAp8Epk1F)FrM*LOLT(|hOOt8Ft%az>~o z<89}sx6H!JeuywHn8hv7SDl#;s4?4bzWG!#9OsUsTy&1BkbHGAvEQZGoDG?K_$_KK zud8rtxaiHQwU(D?!k(WT^^0&T_*Xc*r8roWW+KQ@&rzt+Ck;RVKmN%jgIJr8g_S$` za?xk@D%}>P`3+d8&)o?Z(V$R&aBd&-57$vzRpW48XnX(BU7HV56kJ>oEe3PXh~aWq z4Rplkhu;s2?xu^A8}eM0`Z^}FjdV$|AvnxLV9|PN9TsEdDk^!Z^#e=n(L6V&Q`O7! zS73c2!T8q7At?$2G*!uBXj=cf_5*<^{Dr4FB+sa}fALpEN;0W9HA-P+UiS6XO-(;p z8Ab*?Jc&CXT2|6`^A043c7nZaqO7)|*xg`Zem;+hTE{gbM%cW*SPECX4CY6inmmQC zZ8{_KSDLZ+0Xy!9k218Zd{SB9Zi5TqKp2T4RzBJJOjfoQ5MKYwLw73E#D{OX{J!oC z;l)h^(a!uM<`N+rDq<*D-WF)5 z!CwE;qN?L-yOl~l-3u>=dna<%wE^1Xw?6$@a*h457Vb8jKNE`Pwh)gALd);mk|);% zCB{h8etIZ81Kum)xZgL|-y%d@5`4HM(~?|z|BW%vcWtrQ;Tkozx4QXiuOQT1pkb9d zVr7rW_i$}q1)l<;jU+fCFz41%vB*a$D_oV7_}8>;LP(nxf`vG-wb?o6+eRr6O6xhP zb{Ov_vLqP?E=ZRobI-r6QwI{m5vX?)duDU1R|Uc;0Wb_6-Iha`&SKO2kLUhsOF-oX z!U!bLZ% z5jN<+>4%S2tr@=piss3UX6`_KY8yz3_I)W6@3}>iYxac~-VGi84Hy0fC-Vyeg%!jE zH~^9H>TnGDFA{{vgBID_6iHcT1dm#&4mhDiS=i!n27}iv-ozANz9Zsiv2l4@;sH#* z6$$7~89kSY4EzNiRKan-${mC2L>YZg@?;?T3f2sY1N$Kp_>SVoq2GF`?d5+J5xV$f z84b)%dR!_;R1sglueEl%U0Ol??UDyx2bQsN&)yzWxKz;m@~yVdsB~F5MAt3j#n%5c zVJ&p+o*=F~mf^cTRMnJk5QfE~m7T1NofC-VQ{R}*pv_0e4%kJIriEuK=@GyO5JMzl zX~$O_o!pT$rPMGISA_1vHc7;yhxbE)1-YOKibARt*Yinj%4xW5uf$0{o9hBmdJ}jTsIg>YbAyQ~2sPqFyurD-fa5`YuQbV%;Z zlN8Ym?3|2m-H{A>CH#R+&D6QXKyRsk=)&Xy9;hc=*w7TP2dM7D;Qs;f zACUh6^&imx0n@A|340d!zl~a|`7kUdYwX{Z(q^Hic0v0jXzkyW92ljClm^es*9aCg zs&xlLel@TYASB_h7DxCp1e0+Cv%{(7!?2`$0K~tuh_N8)4@+uw++WX?N0(-0jYw^f zFpZHx9>$tg8K%;ZkKlqjQ6YV!RoYM zo$Ue{a)d`Wt?IP>04^vI8{>T(7Eu58 z)}|)coY)le1bF_Lsm^%%_0zv3Rb)DDxKa|qRhw)i;u8wO%Mbrta%(Z*?m#_SM8|zm zADk!r#}JF7l;U#JOfM0t5aOR)IgM_Z{|yp`{4|W5Qem>6 zyE|l^LOdwKO!}jWI5XT2ya`(j3k%PF>pz+fxk@xGFOMvixJXe<8-fFLz#2$lR#vVED`hx_F^)pOcC$xtE;uC{9Gb>9Op$rO%zHfVD zd#IjGB&Xu)%(Z{=B2=9GU7f`$y-iOZmdQt*R!Ft?oFwE-%WpH4&JWSr1z#Y9{>fYr zB4l6N^Ll=w<)n<-bslnPr^l!|ZS*x{{tf2#R<1(I-)^YP)M{F(Yee#IdwF4xfD;?5 zAXbq!u__h?u04Nshrng=c514H%a{ur^&uCDXY~hWDNt0r`XyR9#m#NBqbY%1ygUSZ zbg&Ul1WoqP0XqcY09o{>VQ`?+$;Pi=;Qi#irZiGP2lSJV{!zo>f4c_jdzZjXewIQY zjn>wMoIL(+f~dDj2-absKkG`@!4opoGF@a@-qGsx4=Z9WbxcJw;(;xxKud(hZ9Ci+x4j6=##UfU~aMtt(Wc320eW3B)FSl+=~e;N|f z{@rTdytGnFsW8!#2qKq0(;gqfR`|ZPK$G%z%htj6sH6!RM2b)ZsQmHt>o-VO14QRpmqFDg7_Vjm-$?E_z& z=wGqk75ThM0d4ecHJ#0~QlyGp0D?gOP2iq7O_NULj9dwDU z#sSVSC5+DeysUFc5c*uo>8ct28ZXYtNtGdZIv58nf*i0`92|QN_>|UigU2&14*_k6 zA;Vlp{YM}7(dO=TLE+=zl`M95OKZyxl24vCHSvqCsn8yQsV}poZa!GbtVx9x0D<2E zgc|?;bU7kV1dY*y>FW=-bu*1^1Ics$S|G2_pB$zc8I;P)V-_1!eUeiskqGflzBCaV zInXxK;rbi!CrG9c%p;7sQ0oZgSjY2`OevA!_irZ{YCNN4Di?g}S%WNnQj$cPmrcgh zCY_A2byCdUYJ=RM^{$Wcc})y|mwU(x7mj^b~Rl@&LqY^I1>?ipZ(R{mJGDL%}hT!ep9{tf~4 z<%@;=Q_@}^ zX^?RBKHO3vQRgI>>r3@Uft&SJe#{t(wd$vexf{k{w|6Ht9_ zZ+C=>3$3y}cR3jf^20qSf@qIJmqc<$m$JY+t_|ozCs0w@Dy8|jhm&yww$<)^Q2lpe z)D{n)%L0ir;57LVg1b!eq_AXyRc*Ay2EKb{*+DafWhL6jH| zd!v~WRj5XrZWDut)cH*IOCjzWQ>{e!a3ah-`)TGe;OjRMzRX~jf!dGXhs1=*Gbtg< zL`8@5V_L#!)yz&4uY-sHXEwAs^Oi}Zk5MGQGZEx(gqkyUel}zGMgBI*R)^35d-!k+ zMso)`r`*2QPNDr@-FRnv;W(tpp9-jCBpClG^^+5ELry5a?6s|yeYwqexait6f_yP`b2sm4@o{N8;nb5O&$B@k?)hPh4 zUf5ildO5hoskK3e-s||HaKD40I((^Ifq#-LnVGX=s;)jexr|BUt6V2N;XEcDfG8sA zD@@L^GgEcs-mCBR9c;gM$>?5!5;g=IwCr+RPa9Y;w=@NmE?lqV(D=BG$T*X18xn|} zu72Br$2veSFJ+Ggnf}9cD3KERD5t-SAEpJ98B|=}(z=enh$vMWUH1Z_08ZlC(^%jq#BOUdMeh9_461i`_1+KvSqL zzJ`3oy-XzF&M}K~76$|3Aq@f29wyF9e1Ev6x6FzaWV5MWSj0`a4oFFWw zhaP%66!5Tgk?NzAi~~uY)dcW$Z}C))5=I|&rLc1w0<#RLl}j3D61-#>_qJfDzH^Vp zstAm!wkK<8iUi$h92m1vL*{w*AJj7qUFMa7kRxnbfZ7Oj(NLh|%lUiH$JY$jwD zQ}kpC5onOLGDy`-mv=OOFikLZtRNUBM;E!>hbZt^ zF#naxC^wU1Xi9Zla9~;QU@ezjA3~%OGYh^R&J_XKh%-S{(%SkMXm?5(e2D}}$=UX% zzRn#u^2{9}{hldec-J|1hZ2(avHDqqa@JCM&Dv7(aBVB&UK6T8r%*>&N1Y%12kHyK zH$vF!B@-sy>nN~KeGlOg?nPw0Z5x2v3FS$=sz|Fn(V&q@4}xJ?fzh=sS;fhTz!W2RR!JLK zx`_MpW_&~y492#hc2^bD69+JoH}USAgPkuDUY>7CTtQ_Ae_o&QdvL<@XiI2rPiSi_ zOI-q*>>;tS;LO2Z;#p#aFZY+_$@j0uU||dEM2X?sT?Hz@`{6bISao*lJnt7IMz$n& z_IKI4%+baup!hXq#x^?%i?GRUpY5B=A80_|i-|T7|S?yEpue2Ep2?@@VQ1qe zJ({Vm?Mvv2qdQvh)$&eu*wi8~jRGR(5N90axad%&o9g;9Lji#~4D*0E>ZIs_>b9 zqd^P2ib1gO)O!yN!^miRJa95?L>k%XYecgT!N#F$Pr8PW+^3EIfbkN{wSpO6Z~rrt$*B^cuEu9?wpZr*UR`l#A(+<%3>`DZjL{dg zw`cF$wq{sU*YCg9Jy2IL`^vX^_hGqSK=8@jSK{9amq52v5-*2`mY?uRz$U??&=Fg9 zrh~}LouQbF#}u5PXb4QjE(d?BcVK}JiXZ>2HKb?RfI_r^BmS%%^7|HQw-rG)5+nv> zW0AV8MHb8ZRy7LN7TcWf+*E?+!Xx;4Ep6@E>s-$;eoBjFf3d+Pt%74A32E`^`xwdK9G@dSx~^2XgC z?Xzfx7+&$AzFRkD${+IQT=?`nR2BumhvrrsqPRh{V?prGYrpvx@{xQdp3iRcU5}ia zHVYo7|Fw+afgI}F5OK(eiF1lIO9eXsPhAG@mr%!D()}-ego1I-6Cq(!8w}xZnf6z< zHV%yle|JU0a%H>yy&J;R=?Lx0Z{AC+h@C;>^5-2Ce_Uj3VcS)0wYM6obVdWx!;4@R zg)G75v+3mZCya07Lmt90miihyirjiT`|X2xN|LP7V~I4RIWCt|X+t4$JS3rVnNyIq zJR8TM6;~T+kRuHbk?#g5p_8P2gNwMF6y&k*zZ_cwqJI{t@6xF2O{PoT4|#uy-Bb1c zk*#HZ@&+DX&*rTsK%!ZEC;k8j5&l~pf!IGx#yWb@OmkOd+yz$4+vo+b6Qz&8yv2%6 z`Oj3!AT$o$%zEJ5ka|eik(A<_8l`4(_emumc5)1>;FvMRaAxMt`*l(RX;}MsZT#{* zBc>wGP=@-jE+@t~ChPGVAmI9Oe%x9<(yAR-Neq3X`cc6`I}Yj$S?mC~i7?oM?pW-V zUK^WK3Y0n@`#OgmZ5FQnx zYT4;tIO(#byO;E^)f0he?%{wT+Csh$XUG{njB|UZK}<7i_txIs+f1*jJb89$GsY;m zpsB>Iv)YT$oZ215k&~c-wq_nib5EJy{R%O;I+Su5yyXc;%!rY{(VR2NK1FBu%j}mA zFfBx#mbH7y_9O;m>9T|M%l6#Ziq+_rd@i7Q{yl;7#W;N2|Lp|2n;pt2g9)<2tvaHh zvVKbXSr!X+a`3k@C(@FI49%n%0Y>s05DaZR{QURuP6k~iy3v)QtNv8qN30>Z?(Iza z+2@w&V)Ht%M4bkAXoJ`-jlACG^kM7D^)227h9+-U-4&J? zV~Q_QGNI*&3;SvQlca{ru4w_IrQeA1G&py}q4(7mQo zLE)L|+fnz1n@@YR7*!qc?H@1PLOu~&4?TKXC>*&P^d^6=dJ*HHHDP>mNrEr0oAlJm z^|*#qTv6_fpRqt@m9TZHprjByJd6e?l5!1p5OIUl3#tGM_HndBfikc-IB)In;@iTm zyCCNGz^u&R-8g7mp{L>b=WsV$=P7Kw|JTTOM>WxP?@I5zSBXe(p@d$fqX-BpU1`#! z1&|I7#Q;hdBy?1oC?y(tPeQNKA_xKj0qG!BO1Ql5SAO4j@4D;WKlV9$_Ic*atT{9L zdFHG$(;!_{W3Rkn#QWUjcc)&lx?jzf*P->N@a~QSQR#si7vyFMf+bS0k*=1ltjFS8 znR*f?=La&2MWn*is7ce{EIV|O^qB3k>X9er!vxug;m3vv8eKL(=^>c_j;+v#z+|rC ziRkH%`x+;!JMmpju7#}jdY-E>-i64sYit2)z8ul49Y-L>C?p?Fo|!zdUV zHYep8-@6k?@HL;%`Qd{@_A*<{T)Jac zb*SiGZzXGiWD1(p#m9b2^(`6Wk|m|PD~ryx*pCk0v=)&*${G_Mp-0Z;IH=7h()@+D z`dE(#9jsuxDfSKPi0|ZAGk{th!F3xoA+M?P>h9W_M;`uUv)0JpN%+>ltpCzDAk#4u zUaaxLpfPj~a7AV~$#|TEiX&v-K;Bf;b$_{=hvsmSF+OB+YHBw>Wlbs?8x?pgrg+FQ zgw8(`SKez$A=z3t(Kji_NA65jiGBV~bXujlqpiU^nUQ)lH?BTa7Z6*0Kt0Y-BcfisuEPa;;Xt!q63#6_rF@^vM>yvt$Q}8`6=~aLy+N~zot6Ae&%LAHgi?~ zf{HN9S=FZ)x%p7NrdN7ch4DOoSzm4hNoukz4d?`<-y+6*l#wQNRL?H%eW2?z_yt_G zF6<9ANzQTjVbzhW{cM|~?D2Vxq+wz_x7DA{#Wk4c$5;zrMaO1w?9bzIdl83+`aZId zyUXLM*ctRy|4~KoF6)=vx_$=ZyG%~y7ZLLNW-2OOApTV{qVl6OEL(raiD`z9n*oC~ zFoS}@%y`*5vZXkS5!qVFulhhD9FNGH3cYi0PV2$y7@palKDk@kR6j!^Fevfz?4H39 zJlTPp_-EJJ+OqVGF`+vd*WK7235|yGqzfl*-*pqbmc$W-1@`9+)!*r?M2Hl;%B zc$ef9uBVuspIp>3S-EmV-C26#AychrIb=g;dUHUJ&a6qn{#{?2X?f25Wb+=l-%>|VR;-1jN?UZAb` z7SY9iPv{ZXh`O(ds+o{c@L@#oXigGIOOh}wRVYnzz4m|{pl7$=rc8kqWG8M5x2 z3VG=KU7ouEXB$yYV@F&6-ff|bIf`o4Oz9lAUqojL?u>zU481%UaPB zlxLpxPuM12nOUSfBKy#}eeKOufI+7o@Fyo0F)TrgZnD7Z2KA(ZU%fluGNb#d)+Z=% zff4g*DcjnK*6xl1T0s5MUOK>G=yJ1^PGdO<7O=LmjT-(b8 zMSrsja~W%!!1?~@AO;VZ9ivIoj7ngvexH8Hsqq;KcE2FtG+AE{0Pr%Xpiikd#s>kzrEgG+0s`Md5^uU& z&+mnQl6n8Jsj2z2E$p;yXM61|YrA0WtZm089xcWp@tUq*UeI2EmP0_;GB82{RIE$L z?3m`uo}gR)D9ytG8&G9_SMyQmnZ8dE};~Z!7(I(?JXKqrQ_dsX&$n5f<%+wH=x7`oZVhapXH`0$YP6tb!@N!# zHc*;xgww*&U($x%0tOn{XyKSwxW>TQ##LGiwCJi}Uf2p0QsL_jB>S82RKc!R{H|(1 zw_joc5Oa;lVYR>p0{P@;zZbLoQ@FheQ@)06%xcHk*o%3_Y$1pO(UMa%VfIp+F&EHl zv}jQ9by^IAB4y8lVqRoAl+qaE+(2l7|Pk@`U5#97vT2)~8mJhYCMtquMIJILqJ4*~ZOE=n{ zJ24W)y>P)-+?Jud&(`yBB{V!pBaCr4n0VzN1w*l`bqO%%kvf%L-*IvTL=4^mk^zqg zcihrKp+W)4&lY8il^e2^%bvR`E{?P&`PZAR&SBlwUZbKXruVjdh4)d9jLDW8m0FQ;Ts)Q%*b@$rd?ep&CmZzG zhT1Sv6sdsAz&OvHTZ%D9yl3fFpQMZ`bcc*tmmhysKB5v(>-dze0Va_MlpDr2SOlk$D-cIG&xL`~J@{)udSRWJ&Wt_9{k z9P06C|9WBmtWR``?3063#nl(^tbv#?1Qix#Y)PkLPvrT+2*FX;VLnAVD5-987+SLzg`ddm*M{KCx7Tn1Pjy-${MzM?) zFbk?`*2i0d8^I46CwUP^p0)T-0j`xn9{GwW4+;`{IcGD-A2O=ec6%?pE!abOI0r2D4xNF7>*$D60OS$Dm*Ki1VPV= za@>9B=OQ8=vy31#M~hn0c-Tv{kZstjc}6B7DC9A$Q)+e24Xml(iY|cRVzvAxe>`F z;&s~Q_T#R*zlL?a+4s{&qfPl^KaitML_s*-**LOYdS290M%{(2P{~CMh}M&Lu~$f& zyH&cgwrICHuEBJ_$G25zm$?AtA50*HV`(v7Yd7aHJ|+ErdCwy_eZGlfbyz#^?4Ji7 z|Me~m^BH_vdk-%$Tlx!~!RwVmWgF!I^RV zNC6I~)GN%Lm5f1=?E5lhPmKHWWV)69JMj7+h;N|Nb-~%_UzX7s#3B@CjKwz%Zl}KH zTXv!k7dbK8Em<%`!K0-}?Al1vZQT9D=1UH8Q*^LA>-#nutd$vKK9|-tdS<6Td=jw2 z&>0F6_BoQUMHm)@;*XqkR9gRf?wicu(o1Ui&=72rS#VX-B))*9yY#2>F&bg6UQ; z%;J$S1Uc(fI;LEo-1fUy5@-?qy{4J-LGh2;UEYeALj-|5T5yW`|CEn3uLx!n7!f>6j*KY85O@;XFd|IW&rL9uL$ zzXDyHuz$II!2jO(f1tPS@Xw!G?0Ghi_goR>&CvK zX5kQ(dGO>z@4d8l$~k)K7yHJXJZEIPNuZSy9&bZByh}+(pM86qg&^ zPb9Ml=Yy4?hzHtyZue^Ea=C89rxqzU{>R@E= z_rdkVFmL_>6-?NlAHfjEpZ5gdio6WG8YV@+v2TypV4O%{n4aP3)~sA8&atIGN&Z%9 zCZAUej88(-I=Y6~onnl!SipeG(1O_A3H6H4OZ+~=)wUr#jhBH#M-YZB{+h#Vov!6! z<`(kx_uFO7ff)JI3^ED2tC-lld0S<;&7Zp8UduE6^$?Mo|HAY8|GgW0o0owTU~v;0 zFi)%`w!ziTfnvD3?YRZb z>YF9-j>hXQ?I#MO(<5fbVQn0htL0~{pNbA|1S9}HIvWBl(cgAS=Lr3C*pCL? z^~{dR?5LT}N8@%Aw9h-Ujr&*C#!k`n8e#O$9pkIHQ#6XKw2`qcG((&@5IGOQ4#t4X z&~^M6;sF=_Y@*ZupM-Y*IS#77L7ex(q2BfB=G~ZkAVtK4Ps#Xw?|ILG!8lBT!&0}eE8OxuD=SuzpkM!=8vC9*nx3MYV6Qezzoiw7T<+g zyp>%dmv!>ihW7uH;q!m9z+>yG+llNr95*P0;}j3~=*vcE_eOz)4Op@V=((R2*`D?Q zDe(6|aYsu)>Fdx0Q>14ODw{Dp2&O=k_Y_fBQW2L3z-;0zEwaD=VCDG3;1@F&D-mQx zhLw5$x$v7w8#Da6N~{8Y`jQ)1!BomDcNtHSsAFE61md;&tDmdQt;Rk*b{rB$+Gb4d zwOImf8$)b34%_H5C`UmDOE9$s4q&j2iG%X+cF;6X)Gx@-Mx7*vL9Z8_?b82*nVS6P#WX86|m(b4Ar6wAPM+Op3+kGG0HUtaLxH?Q>P zxI?Sd^5ucb5{w_poE#RVj$M)Q@XxMo9QEIM=luIk4_Ic*8*LzhLlZ*L(_FB|nXJtIR8+WAsJd0*vheSPRsB1KEi%H#pZ$zMdFNOv8r4<-0d0Y~JGmS~^}T9BNLIorL) zTl%@sjyFy1j+0#oM`e%uCR-bOrI7okhHVOVT<&Pkv*zEN9;>h`pLDTEN3B0%aS3ox zdo+EnZjY&^-8<8uETrDsV_OU@7I6fpmGT*B=0wYpxW7_OGQn%Mz?5LARHw``Qi>1m zz2!4Q#YY>FXj0x`2^R@|xchP8A?j1K1xYpKW0paYCy&3Dy}gBccSF$q1`Dys)PvoU zw+}SSqs?^BGfl24*aNC^Rx^oWYXhJ*^(kml>*GHoRmFRwE$Q9xotJT$KA?Fm@2)<% zqt(I`w$4ZK&GbQ_ZLL{Y1C=k^vFFWTqUr@mi^0;y`A->aLK~3UjX#?2c_*HL6Hl!l zJx-~)@v6r+KEKtvVQDSql-f6Y#6j8Gh5-$f#s-0}d5nNQ6HZ zU|^#T;iA+zWe9)5>&tZp_Z47Z22u{S%8B{Gl#hD(;-nF+rLG!ACIMr)#O-58=jY?d z+hqzQR2|;Y1!}}LLqx2-I>3LXQXzHp$ke3BMvi?_KX3 z`20gEI?fMWX+?RWgMVf9lcULfBSif_68ll!!KTvMG~hpm^e75q_MG}PjTqg}1x8|% z3TU|0+MgY6PX4l*TAwbxiaLy7;RV=#!)yOVRWXE@rF zBdBhtvtgZcg3fk(WWr>3cBL$@X>1*zb50`*OM}R+{yy22;j3X3zSpecz9@7V8o!RG zxATqzPS%RF&wneQ;6ZJ5kgC<7;W_PK;v?(B^Tp{f%F?(wf4i){uR_{|)DbcV|KX&V zy&;R>J#e02vjR>UF9izfO%l{4PHTSkU>QfStM`ECIVXkf96S3TBK{mL&%k}`+O8bk z88$DwnuyUW2maLVV}wm6a2g0c3Is|5uT+9Zn#Hnw?XqKQ=my@b{c&&)X^{ zDYQ)D%$+<&o<299A6Y-MiqOa{d)a=v5jC|17x;-S&UJ z;9Mg(o>yqcQrYM}?4ufwZIM&LZ;j>@3Mv3%*7x0ipC7tO7^LsuO=?(6`$jGduDGR7 z3V>KAjQ^ZZP-?zuRHwmQS>pcN2alJ{ENK*z4;q`L0l=Mmt=i!S*W*s6ujr)=7`?p@ zNM4kI`U^Q14ehuWUGX93)5&r>7VQj``R@3SGo#<%tXvk2=y(L>Tq?A4SE{+MJ9AXF zy06Dw%2=O`O{DlFWAQcio2E6i#Uz&}LPY}m(^K?&d=h!*em&b=)PC$SzFCn<2E1-DSfi6<#c>U*EZb3rbY1KS zrWtKN@aw+oc`(>^aO3?n)LqMC7y`ymAN~p?q&^8+n;?t~P)Mi$WnPm9c@ibKL`8v8 z$>tg!vEpS)mMyqZaS*DFw`#a&sxQ^{2wV-kDat(D*HM{ z5XZ5w|F*m})te;9Q=Q?~6bOG$;*AaxAddDGSGy}Y!*~mU{w&QMhabaGuMnbxI5R-T z8QT>T@UA=7BM!~>m&?p-ZGS3Os7Z8?>)fZNysB2tQDKnedhMWB$G#rM?)GT;v{zlz z*l}EAJy(Iw7Od~wVQIXHVO=WHHEEFL5gZo=Di)zM<{5?RJg5nm0s@wP;JGV1+x-fF#>os|4f_i5ym$6ZM^$`)6B!2HOE@__%J* z&mZfLIu7x?gQMi#qAD<}6|?_Fml&zJr5#eyhcrm=%ln0~gHvL0p= zkefdiqZmtF4uYH{x+8;P>5g6R%;e?k7w>M?^iSqh7?W;jBxV?FGpy*v#p@3jzh5Zu zfZ;}M;)Fa(&mnc-#m52>ex5g1HAG@tgu*C7PKp|N4>3G#33}j!aAdh3QMfRvnQOe4d z0XhA_Ry@9|N3xrgI0C?E8FyR%KFRvBfG{EJ!Ain81ISzk08kZNUMFD&k=61H%XhJF zx@E!i3}b&hxcBnfqzV%eJbDMa1|h~F6N$DKvyxF}*)^uL%4c5gz6Q(vTvroGS;wA~ z;pMjPSj-hs2*Gd!!{IqT6$V~rl)20Zzk`WWkq2koy5l4uHG=~qmvY91`C<)8XbS5r(dG#y)10%_BZ5L?J$w06HQS9msk)7zg-d+CW^?J(2k?SJpifq7+$WByU?4dhXKuJ zK}dg9^>!d-Sr8xuw)Tz(Cp#WHl2*0Gy(HqHXb#ITTwF)$t^>e}G(IF^5In0^H?Z#M z6inC6o=AI`jC5P^R?u@naz}boACBz$xR3}z+uIT|GC@KL1>ixhzO45`6S^El-Uwa@ zS;173E&>ddgEu{;^N_KK*5txEpHRtni zFpO|{MPIE{GFrhniXId@i((-Q3`D$Z1!+ONsCAZEd=!{SZ@g7h>yWKLp@ZJm81}YE zlIJQ!CMeKjSC=9qUuY=PYi+_Hjq{~wy_-$tW+R>WH;b>eJga&rF(p8M6~#`ab;5s= z!hxmO1$#eVvi}}r>%4TR;qK~kj-LoaS{GoQXtb-`gPU&x4xPg+BBOvHUlB7^?uEbMwrsQg^HIV8+rC+lr&Guyez@Ai)CTioY zTHL8n_zR@ruBrs;>!EAtRi@u{=Q#)&w7u_|Co^Ek_CPSne6-w$<3c@sc5$4xoM%_zqO zs2CRRJSSGM1QUP<0)l1K@u!9QHm0KKF)u!?d!cHr-8(027; zdEL-WeaQt{gnY91<+(y|+XfK#B3FRV+fyc7rNS4?fpS$fazA(Jc(*I+h%~Rrp#B>M z@)gH+V-XvWW|By}4yLBPg#^Dv(HU@#ebF8)4>60y=iRZMU7&Nuh8%zJ0|6Rw!h{Ft z{Cr8VU}ctR?k4ew`YhnZ2|I`mcXc%$nN;$9{yZyOu;$f%Y9%0g+9!dplR!v$Ki#Cu z3aDGsog`ACLl7jM7V=y@oRQ(vD_JBpSSShXn2qB>GEvB3*2!3Hq2vN$CYKvw5SPje zH!LbIm^)WqFb^r{ujGG)z^Mc={q@)rb}#SWTsXftcZ&h=X;p^QOY8b^YN70fd>v|i z)Sie^RTE;-SP{GFman?6tqb1>)7;M+0Cy+ptw=m|>05@9v#Y7K5G>u&f}?-a(22iyRWOVnRkHxY z&7?;d;^kVLXMt(29QKM$EhRd|6Hk>WQ$^_G!WVJv`8m2q?y}U?juF0H=LDfd&z08~ ziBix`(t0di1}M8&(=nOL!?|e0Wv-6JF{N8sC}l4~#!hD1vIpN_T`go8>{&CuQgY|d zKk_Fy`wVT49%_Fx2tEG#HnU;`q?Dreju%YwODn>th52tG>(kS*NOHhv;qeq-FoFwy zxJs?&>;pjqDGA>YzH{;!&38~SqS?K|)ZTd?TsHE&m>wP_h#{g9>c;k<8>1Ec=YGkI zLx|Uuc*_xUP+76AlS`gbtg@kuJ(0d)CKG$o$LWp}Xq|t8pS;ELVAf2Xg-c804RfOf zaS!~IzDz2vAJMv&U#j@G2AHE3AGh~##h#yRS!pPhFt_g*@~Gvv!YYe3it}yL&B?jz z$V74GVn42UhDycXrs>(briYNm_HuY{OEJ4#Gt_gd*Y!`qZ zHrEJ9%@y3@XZJMojBCmt%KPf8dkD1(*DfkC3AyQSJ^a%Z2cTjW0y0`OE>98P5jszf zYxWLQ%q2U&e?Hemgv-@MJX^1}Ciku}>f8)Qs=|K{P%iK_4F_v5^T#_R$l)sB;s-yW5D;1ZvXs6S(V!a7Y- z_uLjk*r%U=KR4d*=d^|i|qN|k?NYKhk{rODvBK=8Ulj!kxMZ6(<`=SU>L7b3r+k@SS# z9M_buG06SYvU|G6jeaslCxi>``GsmlD23Xt=?)Io@l(9uxw56nqYPZ2^R6$JUvz)q z={tLT_nbYUK0AAjvHyM;?P=h-=cs?ODomAaU!RbjU1QNHm5k7JVTH1ORXFV5qI?mD zH1*{V6v-p$cf<407Ze=O-*72UE?WB?f;rM@{8a_SJaX9@<^Qv=4<{)gDP8eG#L2n) z3`7(ILubA_`h^rh;{9>a$WQ5#gxh}~-@f_j*-XUW!UUjBo?w97=tHtrnp3~LeUFSH zr5Jyoqfgm-7h{I3HDZ$_P+|5*eOQ&}fiZg7hm9legO4N)5*2UaJc!-r1?MllOtSdm zL#yQ%AFsL|l20_loW5dG+e+G}1|#bcyII#Jjld>~x?%W>OH@bt&G0*(9PwtWHEIq7 zpgR>cis;+5Kii?-BW&LS+-T2>{(C0<&FyZ2#5OQL2T&J2e$LoR12(@dt!FrISe8J8TtnBo>cnbLxB#k9u@%ZyZv zJLZ&Sf-~pDMBsIrdL}elrI-L2?l=?ZFE}T{Tz^{0qy(ANF(RNrNrJitZRtP|92J;4 zVC1D#w&zLO9P9SV-=fPMH8-rX^T_j#f-^0tk{EoitboZ)SZQYS~B68Dr)DxHFO*OhEYQ9tX+=%}BEzDTWL4F`~r*u~IwH`JZRH@o4SLH(!I{==ckg+cVq1!}*H#orGfuLgq4 zVer2X&n~d^{c+2?RYE;La?c{Tt(=KudXOt$K9vG~hyc?5FYAFLQDp1W&5`JDwZ z+;La;S=}zjd}Q`#?(>d22e3SuT-+i3+Ip68RpKjNIO@GXt)D71DTHpQ!B{Q zr~#(d60)|BTee*I#Bt*yO3smeq%BSzM23pDQ%*NL>k&Q zu?d_2;KKoEaT@NdS}sG|1{i1)%BHHbx6pG{LeK32@E67PB3p<0ZB>tm)`%qSc4#76 zD-%!k@N;NytMZzg9#pdq5abPnUCg4PHZ7Kq5n+1R0BsyZm@uooWW>wR6KRyJ)nYV4xjt5h64mf)-g~UA5=V88s zBtyw1{8IykZ}Mi$N{fUxZq8v^=gq9jm)>V(egl8r+%|1kMd#kTOStwh z<`QI+)a?N^K1naqaa&+=@H8&){d9rWt678|*MFsa)51<~!lIassB@rJwCN>Oq7h9= zqY_0Mw^AhD^2fO%NdwSg545EiN}UJ&c~^UBT|)cZ159#G?gH6wpAGOTH|kl_YV09_%6%?@CN3_*h5oe?;9oGDSTg zQVy|RzYBi`(kF`T6k@vy?ylyQtq zusagdJ0))up58NSmYZ3T_j=6U)|d~W+<$uh&fKzmesdUp?q%iw`uBg3G{9&ivk1&M ztf$;8mv2sWHGNm95rNu~yBYw{{<5wb9?)FW6a2Z<7j@+tbN`Z_64{}HM?2{H0yfcU z=mq&YZerQP(4iR{Z|GbJckzafukdYkgwhU?Zr{*>6S2LaH!|KRy4g}w)NMdju#rHbC*F3h7pN{-=(le2siS$gQ zXCgfl>6u8+M0zIDGm)N&^lYSOBY!;`>DfrnMtU~Vvyq;)qgITN&q_M%h4?{2gqov; z$SmRq{GE0&h{;0TRHc^4bQ^d+q_%aH&$dA@E}{33n|92!9B~AlSOaJ}`G+V1`m}<| zT=&Jlb@idNT}~dZ4}(2sqaCNHS&D?po4^!^YNNxq(=>Qin`6X~azaeNOxOyL|q!rzthl=R% z8Z*oGxkSE=Ha-cWpT~;WloDAPDwYEIW0`X z3Gvu0Rv84rbO_f1HWp}wbbr!gu{I7Z=a0<{jD+5?LPs+&5|Odfk8L#&5F#D+HXPwG zogH_yjR2tv#K-dSlo8W}j(TrVN-Vi_zWE%=RaMo2`)ro4A;?~=hk0Z-fWA6X=4CZv z&Wl(TiAHzEJn7C#3ztYlz1nz$3I4H&mZ}^I0n-3esR;WF%vn{x4}ZwziW2Cuo1)oZ z?vrqfaxmFdzFfwF*di>)EQldig%jAJM%$0#oN&Ve3Gu3$hhXJ!0A;bN>IAHqYBU7Joa*g;vezexItF+U!hHI*!#;T>$tYOSDcW^sgfWXfG0jnENc~ z-;32M)L1`DFZP{iRTrV$g6$^Qh+sX{3G@I~rBnfRm$<42s3O966T4>*iu=TJ{df&n z(LPeTvjv4PyC1+Az^fdq(;~0uoqUUCwJj!hC9;t#k2=M&N^YK?;weZI5rLu>$)b|# zG{|{aodmDVq`MvFv_u&LPJg zN}x;A3s8Pa-_H@&8L`6zRcouyV|uVmN*bB)9bjqGB=r{v@pyUue~_JXm+`j*6N6_d zw`VB<*ewAzli{2xf9)FEa@)A|-Cw~zm7ZEy+%MbdgLk({z0FM9#It=!^FT|q&50zs zBo!xr-*W(vA|=R@CYem9(+3R+-~i`-!J%mOkkRb*1>Jf6`+~ufW-Mb=5v5F~tszF3 zY`(g9|AA%;So}4k#5j{Zd(c(Jq$bK~cwc7kF8*=x)9u9%f3LXB9C4cK?Djrmj1eVN zri~`ditKigz0d!C$!%UV#j5PeroH@d`z!3?qg|R4#v~l;!T%W`n9VeCMp%S!L{gW{ z7?Dl|C>EjI=p0F;cG(E%UsTZqK0&NeJ-Rk|kT(4~oWT-jFpw6t=h4i+H_aD+zD;S` zP@)Cj)wFc>e~6|h{BwABBGbwd&P~AkIw1x_tPm41&=mEYXsWY)s0mU@oaH@}KPHah z#B8O&E@xa=morLJhs`u4%phruIHe{WaQPw!N#mm+4>{(*;|(MvnX$yMGnh$m;0)p5 zpA!y#7XjIb15G$(6AtDHd_uwcgadv+2W)F-A`%9ae{Nwc;c5#))4n=VB?CDzPH7S$ z*Tg#aZ7{MC#)LdUW5g;M1*)ee!O|unUl1s~vTOuIkID0p84Q6&DzGT7X!0mwh|uC> zlY+BSHU*xj_0*=;iK7%FObr_3mouRt%!?Po4m=I`_!$YN&sIMCNDX3l1`HV$KnWbB zEcg%)e-zwaf}a(4U`YA$lI0~5q;mG-($e;T$R8( z>@S)xVH#q(T=ssc!A&FhL-h$u#}dV`4r~3ae@Hycd*<=kEXpPkKN$zfi@^)~SZoiahO2N~ z_gHiRQ`Pn}5o5bqjZw?(%Ci@vw;9vh3``G-{r!kkl1d6(L6OHZuw1Z~n9Tl?DE1%? zTQAzoluI3QWGSewX$H;UN^3T1PEebA?L=d#ZpF?G1HsP)G_0VR`=ar{ePRT3 z#CGg1_d6KP7$HweNP~pN7<(h5em&csgFp=js zc?6$|TJmiIPDmGOjcw(68c+PY^EA0BYeKCZU6oVa7-px(tJ`ZdfcNQfX2 z#lacXe#3^**s=Mkh6$SY&6%fV-NiP|Zi{BrAb`=R>EVW^huroJPHWgO4n>p#o#5@N zZ+Sjq4?r@Lw-w?D@8##BALzVVE;radX}hNA$~H`TLJ{=7n0MhVf3y`Qi%(4*8Y8vu zlj78>>>g`Bx2Tu(!o_q_!)7K{YEHS>9Mf)YTIO^~W zti+}U_BQ=X2?#_v?_#%zz{yhO*I{v0_YIjq-RIC6Y06?5rkB-y1Q-N2>8=zX{i*9n z!Bbpb?=)Z5<>xN8e_@cw(l=IOF9jIeYLY?@bHw&HONrem1LQtEEvtEfyTIO*wD~V( zTRi|8Nx3&MeExjN@m7c z^ZQNfF>4x*q)uLQ}36ntxI9sj5j3*H%ieM}zkK<2U0eJYW z0B~FD1<<4F0?b2GAW5I9&O?GV9qbqS?puJ4IMGZNWd}x5H|Iq)enmI2 zGa01or`{j$e}Ago_UF~nJ^pCE*EjuewbN4k?dC0VbOtf|&DfKufe38(sM(1CgLM9S z?TLDOeKJ*o?NjCV{tRI%h7$FG^k(ms%R#v;D)$DB#6=SD1x>SQyD%7A5)9j}c!*#{ zx3dowdll!932!xZ+*fvOc*8(|&qOD^N8bAdB5C5`f6*6Q2&MSMpw04Cgo|i4u@m0( zXrq6+zV94(I&A9IdR{j%Yhc3hNL_G1X9FIW-U7e#~&eEXp?LmFsR7W z{~+1lO5=R^jEG;03{jL%^K=ou znLL#B(__8H{Hwqf{l_GhuhE?hg(GS(^(KI}eyIxsUp{G((hp9LzprR{P(M(Q6YWbU2K^vwzQ)k#oxWsPf5dN3b$G7I27^WaIV#9%(>-mvAcU)Wj_XqS(}7CoZpSba~i}BuWkmwR^;sU>g9?mxNV}YQrA%K2E*& z;T#sE6;kA_4ffx8ABQ3lBi6?NwWtcnqpL9RehPpgA`a(<$)P+@_*XI?Ivy)bpuQw* zHo^yg1oSqKo1stzM9*=Lk3g^tl|4D0b#MpUOAOOO9qx&6KVjITv7lTw>g0@XC;tIg z1Uc=85Mcqg5Mcs~+Lu6A0V;o4ZByGw68`RAG2a$y56%1A)z(%K91^OMT`F*!+_98P z*aEb%Ew3ae;oslWk_<9|!(q9te1O98^h{6d>6e-59tZE3V8MH)ExrQtL`5bA$RL~ZDxwv$j%dwn zB-${iiMGr+$cqkei1y4Aya;%NFGd!CFVYzBMInKiXvLx=TC-?~HY^6BE%X7vC>$w9 zJsyofEP+X>s0TVrfFuD@0%erYS}E#LXgKQ8XhiBUOlnO%kgNeUN=K%m9!NHtcp%w$ z>VafSa1{hv!mBjdlX`zJJvxVmXNpW2RfK6_p9tw_K(dxhc}YEr8knVFewq`&uS zlLHLoXh)t9jrv-XWb|t7!DBo%U=b6j5eLykG~nBa321=XN@_$j6dG`cRZYXAfms&f zoij{3-*^aKhUpw$qhp{U&kXD`a1M^wf`;IE1HBv!VH^zL5NUtLJB12dg*8=q7kS^qd`6a60gXXK=*NTz+a?QMlDR)rN19%pVQHJHp`{J)Pv!1@#-HanP~7mQx?F1yI9y1ZgryOV$7EZZLNU&~2W=JN}-^{>67 zZ4$LuOtTpzGePhgRKu(qm-(VD$|m^`?A?G5--Cg_XD1*}aO@-bFq#D|Pfs^c>jl(0 zgId+q_o9}%OD#$oymuR^wAXp6{9ZIN_c_WU%~^dYJH4incPivc8hH-2>0{tj8grEY zHOaYa`|y8FyksGwn|Zq94XrF>i5ppN*jeuT2?lqc>C`G1)LQzsE}#|oOFGN$T*Nf^ z;zjWVB6r}V9c;9yI82i+o?>{KUR1CpbW0bnv;5++hDAZ{R%3z;E+ETVf+i36L7h(X z@$UR$ni1e5nBzBUyv?f$dcV&n^(A?xefNGo&5(b(BE@Z#BnxIW=+Mxln>L&4w!V^k zx}h%FQB;!{MXB*Unu!~|l<(_>n(UE!#dW=mSm~v$+gz`;<=l4}g1vJYlKa2cWk%jo zP2IT6Zm!z;dh6?kx@9L?DtA$bdLK5?&4-?DnEY7X^ex@YN;ec*dLYd}YLt0#fT~Xo zlmLHe;z0IS;qH81)V-bJO+NpK&8DE`e?n0pPWfwo!1qpML(pC^uGxw5DzUVu1)F%P zse?|0P;~K-xTwIli=Bn>BgS79#|7guUDV+opZXju!gogo5m*YHbv5d}3L|HZ;6^nm5 zpN$q-d0v#Q`JN7@LBP>1T!+$-yX6LktFn*z#d3U^OE;?O<>YD&@u&Q(gsVG`S52<) zbeLuIiBq11kSIi?h=~m>y2K*&`JW5KJjnuC-#o05c)ejxCE8&-vpTVw?ne`8?}l>@ z!q_4d#?g{k5s~2-ToE`o^fj9Yu|k0AMhkzp8lm^J zY?SJWj_%sG9{dRP`?I-E}!KOF8eew=-&*{ux!FukA}7nVRm zP_gK{!s=aC70Ys*QCC+&zRf0i+GHNo;qYqToVJH)Ny&uqU547MhYD|{8&twV7pU#G z@4iF#-I^jg+kJQAt{uOI3Ppb=%W+n+t;OX0IQy_86C<#KQ5@`*SxQB~7P?m?*2WN} zJKGxmQwb%7RZ>cW1n@_peh}gZ@VzZ253wAV>0*&h;DK~n!A6J2FFz#vh@Z1`J|`^4 zmwCnT-(va)`z;JfAc66Ak;`Q`U*_HmKc zmsOUsPn{5l^gwOjRgGQmGt?;&#d%rPdpIuj7yT;v`V!EubNAY5w@%nRd>u|nDo;Li?Wl?97vnhYuTsGkHR|rpj0I%^wxEHx) z^gN|s)OkjCTKtmd^lLhw@)@7=f-kUn8(#8?*L=AdxTjXuEB+b*#^2C6!mFpV&=Z^E z(PdPivY6S3C4+3xU?m!)>rNmX{ohfP+^Uce>%!5tj3uD31hJHap)JrOiP7JBjqqPE z56jXDWp0;LI|C4hM@0jNM@0mOM@0pPM@0sQM@0vRM@0y?M@0$Ofs^5!DVImb2^BG} zo!E)DKDN6(Y1#)^ge9~TspTbJ_t)fh3CsW>i1}vb8%P|@kHpc3myO2> zAv#}!kM%@l!6W{&lc@~8`Bw9GGSylK-2#MKMM1kTD_702qk{T@o~%{DtX{54yTtjn z!Ua0hJdx^mm*~d{Fn{gVZeXwd3Kqq$nEY{fc{G(e%;9sVQigdP(`UD~YU^e)jq(g$ z>RL-pQx(|AgW{<&OE)tm$!EBig=x$TKE<`2o3-I&nnZbUKFOk>@A`&U@F!U}e{s*1 zcqqJzaZ7G=^eikLcUjZndfo6wm{(6^;7=7#cYLxqVmMjqwSP~fUf9;P{=3EY7cXtM zsQtO%Cq@wr_1l2f(p7kl-hAz#wA2_KCiRDdi}d>7K>WCHB{yL}-_Ykodb3bd-}h}- z*_PjcAle%6(6_FxLfrfDIoNB1TlkVF3$uj37z4qvOlIpXr@_i>jDq79rrqou zyvciVbr}7oZQZ=G#g?sg@cH7y&m<}cDYJe+9f_fx27kTZ@%0lOY23Glsb0GeSF!p~92Ew6p^WhB) zj|7jQ7zvU;Xs-8LJVxyFEuIaw9uVtFC2&ZNl`TZCC$uG)MT;EFVp&o5&CJs98M_** zhPejc@_%g-C-|{wBm+RX_aTBO_ElkXO+s9)Ch%*;k-7U>P+_YamXo!>XT1o`>3Q zofT^mG`2K86kuoYJ_!5l2gR^x+nTv;d7BqEXC%M#zHh1)5(or23eYkRbslefW?hj( zow{b$yN>)45LWxB4~pMKEMy+1@RZ>{81nc$-Etq9h9XU2jDJ;8yLVEAQi}Yzn{}bg z0e?tw8cLBRk&a(7p43UGWjw4DZrcs8mVtFs__91yDO*?VweC4sH-dA&U3k}cgoBG> z8GIVVkVqW)OKU1tz6b0-x}saq7+Wn2pluYX*pB#IdF0ov!(%7WbHk$)l8;YWn*kUd z!E|QkQ?V^=V{j4P#Fd>MG$ONwiunUr@PGEH?d*~>ACU}b+(Hl@MqHyCbJPV8t@EcI znH_CDMnNt%5*7U}R2*fm#Dq#zm`B+zCWK6mMV4jL30NtJR4QP1_)>ah4so7!S1;lu z$wEj7hj5|+Mrk6qaT4vt$rS=3iGttUa%r2Tsk|Lh5&Y^ED!A>-zddb5;YT-L(0};* z5msPS-8k!$SZ5dzARGfC^#B4P83|UH)r6C#E=hB5e0D$`tq=|Tcf`;RL1TzhVBfc7 zxSN_nNX#KjQIPhANBqS@(UOvIPY{u&eVxXh_Hj$#e*3pGVYrF(N1- zA%LWZEFD2y$ULIh_Q0KIN$`_hnWo!BFAc$t_+LU;7_g_~m{5tDi>mFKemMp+uXhT8XMawatD>?0 z_G-gR{DJ`$Nn+CuoD*`>p2;lwgoEqGnteqDViQasZ{kp*B?h82jdn4xK9Law;=K?c zpN@l%z3~pMg3m)Sh~kIF{^FDcti3q~s#TZkt{(gtj;;iokCFt4S|62Z zoX3G5S5Urzl)^KS+`gqz*MAE%diYaFU8=VH9-%VoOa-ePHh*6+L8FsTgHcLXLc~Ad z;F$p4^(>*w^jNAejujY12HA@dS^$@p@?a(xNhp(K+i1x_qbVJ+(NqF?+DkMb91O&= zyJ3`RiiKJ3NIx6G^WQ8#VfhCKrL$$>wTEZ%#KV{UGzitrHLSRy zf6(DOCbeT^OmO@Xz6(-;Qi3!XA`%X=+y!HkWKUSunKybY(tlUR6Rkp_quqc4UU4k) zI8S1k!53cMo}MX{zg!>y9oR&sXVA;DFjaaRd)Z#>U2~ z`D98Ur?j1LU#LGy{@5KQf9xI(<$TPS zYj0i0SO}_hPZgTplxqJOkqqcYT8(q}h5+fn(;xNPM6M)wM6CBQ^*_Tu+i#%)$AkBT zz_Y>eU-Ug-5sB6Zsdx^YA&s_8F2!D(qr>Tl%5fd3;L9Y#9OX$jaHJC!U-GbXi2iP^ z#UgClvVR%k_i^1%Hnw#4m$d+C**zQnJQll$BwS5DKe1DU8w7&uLCEqJ84^b$}pWO=! zp^lsS9)pWY|G9=#2Wv>B!X(=vdons$Y2eSVm-s`6zk!g6XC5d#CC*Tp6{z$P&}$zuuB&N&|nx~vp1})=!YFQ=bM;A z#-!3wNhE3!0wE4AHr_ps6Dsr9KkFEQed4A`@JodkphP(#3%h34n1_y*yJZ0T&qo0Q z<-$jtUrufjEpJJVJ9s5d@D?xI`!@T*#VITvPUL1)y7Tit0*t;ti4&cw7wJosgo!wm zm1LHNkYed4z$~LNavPmzcW?d=l>$yRlfgU`mrK(L76LOhlVKewf2~?ubKABOe($gF zBRvz*#DxTzH`{R>CvlTw`!v&O`#?*ug_t5$l1k$K`m-01qG-`m_w+$SV0Qt6*bjSy zqJ#T`=-~XPX!Ja}`RVAbOb@a!OJs6zb9WFYVUno>BjZqMeQ>il_$~NN=synsxOv|X zJNRuDM^SJb9?o=>e+HiqWg6toyvmDLK#@ss;W;1jl`C!Ssh^%vCO()65oU4Be>-r#af`g1TM$)B(MnCf zQcujw1vgGBep6i)NYF^Qv$cisBgwNLOLyzqS1cy7EQ4a|OuriUyz5+@gTl;E zVKS4U7AeTkN*9%TEuv6}C_6G?n#m-AImM3pMk*JP4&zjIGlaXUo@rtHhGbz1ulA6L z_d=pmv5IwYe*(MHw5VMbXThnbUb})7k4kCU+Is$Qye-Et zY6Y{;b5^iUShc@Lail>H@Dy8$j6~E?RA3VGdn99HfAk9xAmGfVi5w#!63NlmFsBv4 zt8q~wzA=ns7{O|uLqwW}Qm0dhP)Q7cd4h;!FCyOHfmc~@^VO~&XcJd-;FBbH?^$qa z6Hp*@G|0QJ51#+d^RFw?w_W+e3ao2eHmE9O-+?ix{OG%mn9p0+;f&tlZtKpjfShZp zqU!?Je|hq{wsY6?Wc%mFtvkapt@ZcSrY;GapvLNspV!rfkGbBS_g}V)z^+nW^#qS6 z*pLHksF<*!vzJ;27-6*7kBzU*H$ViM2N6aY;6gqNidcjRxP&P)AqAjtj0v+B6QuY| z4HOTZch-F^!EjqNs4;x>#|roT`|Pd;J%IIme`Gm6{I&*FuX3`fNOh}qP4awDJEbub zb%M{vx@P6(+{xRPg*F`!P=fqT<0VjUw9Z~{Qx2jbJT|N?HX@+N(*9_f1Yc!`s=u$j zEozC_f#*9eKKGwN*hT8Xwqz~mZAAwFX`h!No-Ay(OK^;hcOt(WraGwISIgc1mi`pf ze_$-FYqxj;Hqiue?n_R zEm9+Gno{=y(>cl4v`^EZQ zM0&T2QySjc$%`G9ti%$D)&-d5Cd9zeK<=J5M@SpjvC+Sog%r>tP5mei+0{BGJel)*wRr1 zMyW5DZkB?_x^dcS>%o9^^xz> z`963htAUca$36xctLoMjc2-xn92B=q2JhzD&BAieHW)plf?uQJW$<&ge|0hi%Xo7P zXq-|in`Z;iTjc9)cqI+wX;rbo>G7LnYJy8UM`SJY62i`cNA*~yDR?dz2{~OSvFe~Zs|?*87u>PH$5=w$_Gxbn02}bH$Wa(A?S4_8Ek)d?UGyVvbL6+gV-Irw*ZJ{y3BGlFk-N? z10$w|?kWb}kRdo#e@s=ISrQrYs0KVLWQ}}hP!!R&F36CRBRLL1BnXlQ7(hTWNK|r= zAd<7>*aV4!fGBYQ$s!1Yx!3YI=M%un|WPdp_){7SmR zfAdVoiM@UTe@LyoL03c>py3|J7g?goamCh}ExUW~bWBI)bhRS7^OJ4PiiHxRADw(G zs%!?o{+*J*gQgF-h@fbNd0xZ%^~WIOR0_c>e()o!^_^526%_y-pkOoal3qb zdP4xzsCOv5u(JC~vEtO#nME*0r0f%89KQgItTLBc@?WwyukA3ujU#Twq94b_dr6IJ z;3@lRh!whm61snft52D*Y}F=l2kXTB$!g{;tHedL$Le-h_q5$7n33W{hvon^=ySj~ zX@Pb&pUs|?G^n9`YfSq6AA-o~hXeWA(fY(%-|`WmKN$=q$=5=^ z!cycWjR?UZ%l@=2;g9#cHjDiShFS$g;4Y41A4bI2l77>_M^H6aOwlTTmUaeZfA!2o z`F`1~wx%S7*S_JpkBmJa{%oGI^Lo5MU-=2Ot!kW+8Rd;OPYST&m%}HQbZ_MS0Ku0y zPm+W6Z#E?0EiaQA-Rj4>8|*->XwTPY-K#hsOFOmkC0{oWG%W85BE)UmI&8XxoPCu; zt7IaE_n1H&PFs)cx!moW9wU)vcRtE&4&jbz-QN zw==Ps7{An&yz)>(_N(fX?BjNkc)o||TJ`?d-*w9CI&~TV9hsq5{D^N)tlv<3yI-%T zhYSUK97{HREll28WP)iDKu8qRMfEOMYccB^yZ0UC8o-cSh>vvRo*^`*}1F}Wj z?vb*_W`dY$fv3qI1^0^jnj}Qq8Jha}^ys{ZcA`HL>Ss`G!4#IV-8VY5O7oBe(LG(9 zceE-jtc05bQFNKstu`_cpA@jZ*O90oYE()S)I3D{y=k ztXjnfwIi18j?ilLg#Y%ISOqm$xU#t$3u{&gh9aX_JI$j*8VxS#^gIG4S+=8smv1$SxW0wl&wJI_ zSf`N*;-za1I{X&1j>wnR7gl6%DDWB>Yew(CBi}+wv`Yiahc9-S5SK$n;+kxD; z_T)Wr__2K7<>jFMBo@9P?#GpaXeGf~Wr3fW4-)Ve`}u#$g>3$`B_6**vnbj&!CRxx z7l|aHS@Qw|F!K0nY%!juj42y5ED2KmKR01aopfs3oV70MG5n^?N^N+!vE9PJLX z{?TmvHtYQL$SOW}A)=eo^Y&s*n&pREH}0~ne%j-<_~{7EdstfCv2)07!& z<)z&_L+_rmXx5VN1fwe4ppHs4qS4wJ-c`OYxQ~tRaq*bW`0Q-yN)0c}2*1Q{rD(Vc z0!XT_23XZo$4`D7wg(v%4S&w1o1v#qFf7M3v`CvkM;LBe$?LD^AcMw$xxxl&L#c-^ zcc4eFd*@nDo(_#m&L!!4m`64f;ahig4omF;L2%o){VfXr#UzQ6rrh18+-#qnmvvvV z{dUq_A09qyHE!=coE~xEyY0bu_`cZ88~9So>pbF#?{Zq+S{)hQB>4oorxS~e=u_A) zc0&WqTQSsHbszsQn|<$7q@^wM4LMKatBc4mLRxli|N_^_ku4eL?q4>inL z3=$Wr%-}RvtfiUHuxnJ@^S7V&ZlqshtIt4&k9F~77aK1!Z!cFf^(Ue9*zD#Z360~< zsSRbuFZ+Z25nudW=@Ks=C`P9j5f5@PXOVF0rFGrWB!I__-nh#ZNXprY$y|}M2~Lds zS@=b{FA!pk2xG(;Bf%Ic#>g;6Ugy4W+nt+00468^6&Dgq-5>%p{NGU>DKWT%@QIk9 zxSWWnh^U~LI80VSKtw=E04gsotSBolC@v}jg-WwL{2x=E5%H^AyV!c#v4{x#-LcIv z=~u>NWHzxMFx#Z(T6J~V%)fo?Qs>h&Q#-X?+q8Xrjd28RW`}Z8+K&>(A@4pSOtSIEq=dpewc@cRBjFXC2{^&A^uu0`Z2!Z0Qc+fYiM-KE7Cx?9(OQq&n0j@Kj z4ll9~cXmhDJ(U_E@tzXwS))cijR3oZ@GzasBozlO-E)NY)JIf^0siZ1Kd-ytrc)Zu z`;~&UkZf=7m|EB&NZr9ie3&j17+;{|rL7&~cL?ZSkTWcBZmG2;E4Hkx<(tfwfQNES z?l=T2x^eR|tjE23T`lEZ&!`{gsZ?4z$t-S4SQ_cPU0)kTW`Lr1rkG@g0j9+KcA2Bz za|;r*W-i_vq|PUdu z83g`uhx7zF{=7AepaO~cCsKJ{4@OkcF6hhj2IW_eMPo_*pxwIumcar&TqA4H^Vb2W z8QZ;$h9&RNje7QALlvzEL<{GH`7Tt~3cC?%wu{pSqYVLW-f${p>NUL~F+Wm*MdgO2 z`Bv1tKv`J*)utPcO9x*)Df%SoU5|c&-~mGnIlJ6@H?Chqn{wTI?~phxf*&qflw$V< zqxUU!NuKM{5{J7B+`>ynUdXBhJ@BI5j;Z2xyC@I7riHr^_pd{Mh;bzG)%XI8R`|*n zLUGY2OMT9%vx?4r6FnNfS>oXv5ExZPR$ef;2_aFw>>W~=8XIl!6K8uSWC)j#ymn^% zb~YdBm^(C|qLM@Ag4yzy25#Yceryqg(#Q*e4pxBah;azB=GZ2&x|WehI+cUK>dq39 z!BCPx28ZymJQ#qDyR!7gR{&yub}P|RT9c|+j2bU;4Eddw)&z&$8ZmRRNTD=@)MlVn z$60c=XeN!T3*8r;SQ@w&5-KYuP>UDY#j*n*!pV2Cm?z0IkITnI1XXM@UMFIKgzp-Usi6#+c=)KVSW z9glq=3=piD{rFpnNq$Y_iDxhU1iqsse=}Bq1=Cg}jU`Q4R2Bl;8>ut_%CsbQif~3S znT1{L;W36#>%% z$0&5(xa|Dchh)ZxP;4$E*t{}^$fZt~b7I9RUg;s|3o29EjSn17o-`#)29562rG)rP zcm;2{E`gq!^jYxv)%+TJYtYmq;i5O|AZFN9XM8i(SoNc+X0ZHs5Yvn&@w7t=yH_i) zN4n1>NVsuvp3QYV(n%h7${Xd-Q0Awn<~&eiT;zB?cE{Y~%Zlv4NQS&fki*JoLCM+E zsL^3V;o_9B42_z*<0Pfa%{I>l!ZTt-q8!dg`%6@HF2<)I%x@k7c)Q4^U?Eb9y)hG} zwo&PqX%--*T*<)P^PZuPVmkibMfrR{IAfQ^bSbA|(hh#fA;ZOD5QPAt4g!mME!eCMeTvy4Ef&_ET#i(8E<7}=7yC#1YXr0CKS~Pxqozy<3h#E8 z37h~+$W{~%rnavP=5I4bvVY*rfkceY(>9ZJFJPL7YnYQ1Io)j7h#THuhd7L!W}dCT zqwt-mR1C;e)Bcc#t?}(}$_)j2>|JfpgFHv>nbWT5aujJjVcooFGpWzK%i}4aE_wL6 zg4Er@JDS~b!6-8-Er_H|;%mI3duxCjgc&LC3rbiZe5kQr2ehD6+a6uG#zfMbwFjWR ziSo9g@X#tZvcY0>cn;&62=s~kna@1EqV3a04-ri_Y@3BUrpG)Uk5}G3c_omV%K#Ql z?P36ffnK;!7t0@aKXo|{pcX0g8;FmobqpBxDSAw_+pWZ3BtD>6t zmW4z1TNBHE;n&nV=8prp{{m^;WE1{Kf-dv`IwKrcF2gDJ}O^`vvOos<(UIAo6t1cA-1!w8?s}msQ%3uGkoh*b=Oc zdY%@_fuKxjf2m#JZ4ky&(?ma7G*k0|ZM!_2#!|{DBpyL1ilUzQUo!zf=!-M^D%s+f zSA27M(zwV$qmtQp&t^8hRcQd21W}ItqlLs^p2sx&pY7ch#Z*sxGB5owBx;X<8Aahs zG=#L+ZFN-dhG_}h0+lo(<6{}9#6F#NcVoa9Wse`n>U+4sOS9g81$t#q{sw^mjQqd3 z|K-_!%nk_@+;R4yKkIIbvD2bIE8YfgYp~y3c7wOo^1G#}z;rx9x-UNVXeZ^e$ty-~ zsL*fot{0#yhdRtjEBB;MKeLKna|Ob_9zzlilZx0h6fRQ3SRux>#~UP1gE56C`0 zhojIdO#kU*64ZYuTlIWgdc(6Zd_jH2x#-bcU4BUuXM&iTQKnIr`+#%7>iM4WHLF=X z^ySC5SFtKobUdPZ(;vC?o)j+Xb;4s!sTz1*>izoEqbH@Fe{(6@ZOt2IKNlc(_hLH? z>c3*O&m%U8C2I@`tKO#$s|MK<4&N!@P*F=vVR7bt1}=!ly&qOxsjisn#tf!Y_!G{c zG)|Utd}^mb^jS4`S+*}c1)q<6a!rT1`r#knW$nS#H`n`Z{{a6?6ydJE7xAOuyV79l z05g~#_^+|PoJoUIVdeXF{R3k^$K}uO>-1+a0CShoTA`W8ftY(!V;;hEK3=0c+yej8 znTY9r3E(XWBsIbvIVp=+R-LS3AGKe4tj-?8r_59JQxDm-rGTl2`9p zsVnH8pJ}0F$SN<2m#Wk(>3E~40K^^rb_YWo#|2zNeEts1AJYt$6i$(R2&&q4Z(l-h zq~-)}lO_39L=XZ(hU{~@Z5H}Z0x=7>+LK_!w^})1z%b=2ojq!Jq_7+HVEey|8Pg~R zHwu?_XJ3XLoKyGMrE$89+T zx4G-pwmTI?cj+@nf`YWnK>R0E&pKSeR~@%y0!UjQejS>LNkEt78p&5h+K>*=XNXW( z#54ljUVyt!_nh**?0>u{blKL$TLB}qXuSi9x&}ulesZ@=^^>$+s_nNf5xCva3z>w5=7BEff2Np2nA4Etl_)MoF*G#@8(qCiR&yy@QhXu^^CqfQ1kdx1|KD?D7 z0(G?e@g4ILhvi2K-|oaAxjb2~aWR~|B!HbPfJNwj!kv%#soPi>6254Cwp>xwRWi=$|+XCa;4fem^^%@c%8NFT$4Ud=4kvrd_W zlZFqIHHO09cxkM?n)!d^xcVL%`C_fkOT(iA+^}~#6oGo2sy$KDuqijXa-?28-Z*RP z(!J@hsMd!dzq>@{%&+BgIP(zeC0xJ!@piFpHtPAV$Jp(+Y^A|brd5flW3Q^pQ;q5d zcn__qU4=`7)2XJ*Y~D~!_itF(42YFp#eIoQGg9m?W3$l*SNR)#Yug-QYH9ix0S#0j zbD_sH9$C<1gBu>{W6V4pv25XAo;&~PyCA?Ans4?Yp{AjK($vQJL6vB#QEOk>iZ%6# zkYkV$Q%>|@y4hx$*+3rknvkPM`oHF7_VLbH>mLfe4oI(Q?H@L`Sr)E3NH;p`FAKDx zejrSonjX>D-(hV7EC^M(r5n}vm+9K1NC>NCr56{bi}dx!+UOsCdIMK(WFJZT)l)HQ z?9o@&ul?YPAb?NGrNLYW7f>*FNVbw zJ!Rsr*3}jM%=57BJy3RsurFtP(AgmhefS<2YL_9|cf`FMV!hsLBiUC*@wjbi1+xw+ zXRzT>*b@y?7I3TRBo_6IoQ}?Qo22Qeyp+w68zd&9qU4I!709#DkxdMHifc|}M%B#Y ztGD-QUa!40%mCMeNSuF7NBhIPHaHoYL=;Bp#B=&o`_l(q*%G$wWgQG)EsZFn{sI{H zIQ&Kb4NO&*@=YY{1DBnqt1qb!OnOY*dE>7JQWtU89~3kt?#ML~{{@Rh@t~o3Osi40 z4@A`~x3#QYOxbVc1ap%uf1w?`<;M^NC{vcy-;eE|rclh=dA%sK_*1ihZLm8&j6o5e zz9Oaaxv>+VED=pkl>g%k(Uxbjf4P$z@fsynbN2u4k;V;Td(qu$V!? z>gP+IBoB_4@COaf-g#0M$5l@kn7wKuTA3_z=~ z^QK_<{jhHTT-`7*ZN!5Y7+Er5Hd2*SR4lQs5n5a=l>qpQt{)Xy3QlNzuZsF$yZUG- zI5eAcOo{>I?;U#NM7;JB{Jo0qgDnITNX+KLl44-<7mYajD<@fxXE1USun>N<%kj5y z;I*~-2NND2w+F+p_`{_*r2;hlkKspKlWR=a45I#`PDc+9cYh+#RhIdyss=-5@!1id zvz_UshDZWxTaO<8+BMqR_(?HR6peW{MFBv>g6axvS zMHfQ%gQ?KgdvL~=fjPSDl7>647kakt8gIv=5xoCCn!2F3 zBA6{N!rq}d+U`llWFv3hI@D`}fz23RBwc$qFCF?@>V`Xy!+zrsLpj6H-`5F>)w~)t zH-dS~_8X|18>aiiK$sq9`($QO8^DtnWPgsYrBE%y z#~i?)hhx9bd1IpRQX0a*=2Dzq65PGhN%E35+&wKsH7yJmc6;ismHDqf%b#ee#J`&4 n|0_u}g{{nY3(~FUrp^FeHNVkR5I2~WP*jA7gF{hEiRk|T0ScHp diff --git a/journal/dehaeze21_desig_compl_filte.tex b/journal/dehaeze21_desig_compl_filte.tex index 1b3dc46..6b26667 100644 --- a/journal/dehaeze21_desig_compl_filte.tex +++ b/journal/dehaeze21_desig_compl_filte.tex @@ -1,4 +1,4 @@ -% Created 2021-09-08 mer. 10:49 +% Created 2021-09-14 mar. 13:38 % Intended LaTeX compiler: pdflatex \documentclass[preprint, sort&compress]{elsarticle} \usepackage[utf8]{inputenc} @@ -74,7 +74,7 @@ Practical applications of sensor fusion are numerous. It is widely used for the attitude estimation of several autonomous vehicles such as unmanned aerial vehicle~\cite{baerveldt97_low_cost_low_weigh_attit,corke04_inert_visual_sensin_system_small_auton_helic,jensen13_basic_uas} and underwater vehicles~\cite{pascoal99_navig_system_desig_using_time,batista10_optim_posit_veloc_navig_filter_auton_vehic}. Naturally, it is of great benefits for high performance positioning control as shown in~\cite{shaw90_bandw_enhan_posit_measur_using_measur_accel,zimmermann92_high_bandw_orien_measur_contr,min15_compl_filter_desig_angle_estim,yong16_high_speed_vertic_posit_stage}. Sensor fusion was also shown to be a key technology to improve the performance of active vibration isolation systems~\cite{tjepkema12_sensor_fusion_activ_vibrat_isolat_precis_equip}. -Emblematic examples are the isolation stages of gravitational wave detectors~\cite{collette15_sensor_fusion_method_high_perfor,heijningen18_low} such as the ones used at the LIGO~\cite{hua05_low_ligo,hua04_polyp_fir_compl_filter_contr_system} and at the VIRGO~\cite{lucia18_low_frequen_optim_perfor_advan}. \par +Emblematic examples are the isolation stages of gravitational wave detectors~\cite{collette15_sensor_fusion_method_high_perfor,heijningen18_low} such as the ones used at the LIGO~\cite{hua05_low_ligo,hua04_polyp_fir_compl_filter_contr_system} and at the Virgo~\cite{lucia18_low_frequen_optim_perfor_advan}. \par There are mainly two ways to perform sensor fusion: either using a set of complementary filters~\cite{anderson53_instr_approac_system_steer_comput} or using Kalman filtering~\cite{brown72_integ_navig_system_kalman_filter,odry18_kalman_filter_mobil_robot_attit_estim}. For sensor fusion applications, both methods are sharing many relationships~\cite{brown72_integ_navig_system_kalman_filter,higgins75_compar_compl_kalman_filter,robert12_introd_random_signal_applied_kalman,becker15_compl_filter_desig_three_frequen_bands}. However, for Kalman filtering, assumptions must be made about the probabilistic character of the sensor noises~\cite{robert12_introd_random_signal_applied_kalman} whereas it is not the case with complementary filters. @@ -481,7 +481,7 @@ A more complex real life example is taken up in the next section. \section{Application: Design of Complementary Filters used in the Active Vibration Isolation System at the LIGO} \label{sec:application_ligo} -Sensor fusion using complementary filters are widely used in the active vibration isolation systems at gravitational wave detectors, such as at the LIGO~\cite{matichard15_seism_isolat_advan_ligo,hua05_low_ligo}, the VIRGO~\cite{lucia18_low_frequen_optim_perfor_advan,heijningen18_low} and the KAGRA \cite[Chap. 5]{sekiguchi16_study_low_frequen_vibrat_isolat_system}. +Sensor fusion using complementary filters are widely used in the active vibration isolation systems at gravitational wave detectors, such as at the LIGO~\cite{matichard15_seism_isolat_advan_ligo,hua05_low_ligo}, the Virgo~\cite{lucia18_low_frequen_optim_perfor_advan,heijningen18_low} and the KAGRA \cite[Chap. 5]{sekiguchi16_study_low_frequen_vibrat_isolat_system}. In the first isolation stage at the LIGO, two sets of complementary filters are used to form a super sensor that is incorporated in a feedback loop~\cite{hua04_low_ligo}. A set of complementary filters (\(L_2,H_2\)) is first used to fuse a seismometer and a geophone. diff --git a/journal/ref.bib b/journal/ref.bib index 30dc825..991e410 100644 --- a/journal/ref.bib +++ b/journal/ref.bib @@ -14,7 +14,7 @@ doi = {10.1117/12.552518}, school = {stanford university}, title = {Low frequency vibration isolation and alignment system for - advanced {LIGO}}, + {Advanced LIGO}}, year = 2005, } @@ -31,7 +31,7 @@ @article{matichard15_seism_isolat_advan_ligo, author = {Matichard, F and Lantz, B and Mittleman, R and Mason, K and Kissel, J and others}, - title = {Seismic Isolation of Advanced {LIGO}: Review of Strategy, + title = {Seismic Isolation of {Advanced LIGO}: Review of Strategy, Instrumentation and Performance}, journal = {Classical and Quantum Gravity}, volume = 32, @@ -285,7 +285,7 @@ Giaime, Joseph A and Hammond, Giles Dominic and Hardham, C and Hennessy, Mike and How, Jonathan P and Lantz, Brian T and Macinnis, M and others}, - title = {Low-frequency active vibration isolation for advanced {LIGO}}, + title = {Low-frequency active vibration isolation for {Advanced LIGO}}, booktitle = {Gravitational Wave and Particle Astrophysics Detectors}, year = 2004, volume = 5500, @@ -303,10 +303,10 @@ } @phdthesis{lucia18_low_frequen_optim_perfor_advan, - author = {Trozzo Lucia}, + author = {L. Trozzo}, school = {University of Siena}, - title = {Low Frequency Optimization and Performance of Advanced - Virgo Seismic Isolation System}, + title = {Low Frequency Optimization and Performance of {Advanced + Virgo} Seismic Isolation System}, year = 2018, }