From aba472be8929ae83de28230cd098164480ba29ec Mon Sep 17 00:00:00 2001 From: asantamaria <asantamaria@iri.upc.edu> Date: Tue, 22 Aug 2017 13:17:14 +0200 Subject: [PATCH] added first version of active search --- doc/images/tesselationExample.png | Bin 0 -> 29304 bytes doc/images/tesselationGrid.png | Bin 0 -> 18896 bytes src/CMakeLists.txt | 7 +- .../activesearch/alg_activesearch.cpp | 95 +++++ .../activesearch/alg_activesearch.h | 330 ++++++++++++++++++ .../alg_activesearch_load_yaml.cpp | 50 +++ src/algorithms/algorithm_base.h | 4 +- src/examples/CMakeLists.txt | 6 +- src/examples/test_algorithm_activesearch.cpp | 167 +++++++++ ...pp => test_algorithm_opticalflowpyrlk.cpp} | 7 +- src/examples/yaml/ACTIVESEARCH.yaml | 17 + 11 files changed, 673 insertions(+), 10 deletions(-) create mode 100644 doc/images/tesselationExample.png create mode 100644 doc/images/tesselationGrid.png create mode 100644 src/algorithms/activesearch/alg_activesearch.cpp create mode 100644 src/algorithms/activesearch/alg_activesearch.h create mode 100644 src/algorithms/activesearch/alg_activesearch_load_yaml.cpp create mode 100644 src/examples/test_algorithm_activesearch.cpp rename src/examples/{test_algorithm.cpp => test_algorithm_opticalflowpyrlk.cpp} (92%) create mode 100644 src/examples/yaml/ACTIVESEARCH.yaml diff --git a/doc/images/tesselationExample.png b/doc/images/tesselationExample.png new file mode 100644 index 0000000000000000000000000000000000000000..21801400c2257327da79b6d338e09889cc37ca20 GIT binary patch literal 29304 zcmV*6Ky$x|P)<h;3K|Lk000e1NJLTq00E5v00B7&1^@s6Z2Ww(000+BX+uL$Nkc;* zP;zf(X>4Tx0C)k_S!Y-jOSA6T%#w44CFh)T&N+j`B_nCc*#HO%iil)T6p*N*5=B%L z1qCFEVnPrR5JW^l5Rg0Q@tkv?@7{m+&z<Mlo_DLNt9z$<s=M9=fOVNjB84Dm00<3> zj<hw>B{;hrAYfktB;W?LKoKB_K2ancYio1xzg`>90JN_hm)Q69zncBe3R+*kC?5bI zYgp~=8|4!U;}!rAx;~`HXaJC#uzVygngnCiBN#J8!V`q?>wO&X8?WwT?_U^=`ekEp zs|Ns@1ppickr){O03|;xC&c;$z<wz80l?6O`38jnK<^A=O`kxbF958?u#I?VX!t&U z0ApeAfB6ac7bX(D|Mp2F2K<eG{}#9>Jg|OH6e)z5@axO}WfK|_13#}{7aiPEgoRk{ z$A$_1Ecg=j&EPu+{P{@=`4uzR24eFIbFhPD`&c~8+tTVUtmz+VXbWSwFJw)M*4@YO zGlTp|(bo2V$qyw4>RG~=0+#3dMH&3@bKxM-bU$WPu>3}3jI9HV`C;506>DbzV|YbD z?-K*<o&K`<=If{bmmGoc4>C0QZG)f-iZ<CFivf-~Z*aKT{v6=35Q@MIgaAJf10vxz z3`hWTpa=B-Y9)X_Ai|nhSQiC?!3P+J!ac&_KH)GX{AH*29}>e~{Q}^A|J@@2_`tR? ze~<ZlGy**4AD=<M7sh|&M0kvS-Tv&NJcCmI8S{_DzQ13yk#3Y;l-~I>CzK#c7NvmF zMQNZkQ7V9dVnVS236wla8Kr~LLaD(z6<`S7!=3~F%yHj;=nena;_w7j$3R%S|J%L) zXo2G|>!82C|HlkK(9m-G@}C(15bYNazXt&5g_9B@g8~Ah2|Dnb>nBby3G<N@C&)_6 zNbjrv-(34S>-S{8+WyK#5YzBKGS^`M<T(Mr*8L-kas%MnH2?%s|B(qb0>Ca0Kn=nt zCNlQ7Kh*vW4e)>(UZJdj3-AG9_&LkMnNbaB0|Q_Rtbjdm0Uq%B353^RB!~ma;3&ud z*&r7bfKpHasz4pM0&akI&;{;+UN8Vgzyx>&rokLo0;^yHY(o%)g(x9<h!x_71Rya; z7E*>ZAp^)9vV~kAFDL*CgJPg$C=JSn@}N@aB2*7uhi*dmpvTZC^a7fJ7NIrh2Lg$p zKrkRU5rPOQgfc=0VTQ0rxFZ4(5r{-YIwA*AjHpChMzkUBAqEf=h-t(UVjT&P<VYqY z4^jfDjMPI~ArBz^k&(zl$Sh<bvJ!a}*@^5!jw7d$%g8Mh21Sn|pd{dRYK(F~`Jf_D zhfpU`rKmbo8|ndS95sVlMg2rmqB+oFXjQZc+8G^yjzedl3(+;`R`f&k1bPm=fx%#y zFhUq5j4{Rqa}bk+If*%kX~J}4MlrLPbu1RkiWS3ZV6Cu3Y&13lTY|ley^9^g&S5uk z<T!4eEY1*j02hW!#TDWjaCdNHxOv<*84Z~LnJSqTnJ-xa*=e$BvQDyLvN^JCJS|=b zuYtG6AH*NU7vrzuAL3u*zmk)a^OCEO+mHv5A0aO$zefIue1?3Bf{sFzLXX0YB8DP| zqL!kY;swQ5N=iyWN^Qyml+lzql=YN7lv9+ORP<DmR3=ouREMd`sM@H;s8*=Ss0FBX zs6D6?sSBx_sfVeTX~<{<Y4m7_G>2%)Xl~Lxr&*_^r<I|#q79`zN!vjCn0A2<OD9NY zK<7)BPFF?OLpMv0r01j8qxYdtqpzlaNI%DbVGv?4VK~T;&2WWbh+&nHj!~Y`kujdJ zjPW+(G!v3ZkjaE8gy}TX4W?&IKbSe0b(sB_PcUC)9%KH_!p@@2;>VK3(!?^s@`IJY zYQP%IdWN-~^$i<}O_a@sEspIx+e5Zxb_RAeb|3aE_8aUkIS?G89JU-u990~H9P6B1 zoJO1^&JxZZ&SfqpE^V%ZT=`sgxE8tTxiz_ixbwO1a(^T+5_AY5gkr)2!YU63k1<aS z&qbago^4(sUI*S(-X`8Re0V-(K7YP^zWaQu{M`H&{K@={{4WK_1e65=1&Rdv1ilLj z2|5X83bqR_2(bv62qg+N3cVJl6xI@s5Uvz{DuNMF76}$PCo&?65LFNj6fF}S7DI?B ziUo<C6B`vri>ru-iC2n0m!Oc)k%*OOl$eoZk~Ei0m28(>k>Zy+AazEnPwJ<%f^>*< zwe%|)IvG=$RGChhHCa(vAK5b5r*f2XhH{7HI^@2{i^==RpO>FhpjWU^$W*wmu&bz~ z7^!$waZ!n1$xEqR>A5n4vbFL_<wq)L6<w9XDqSjDs*0-7sy9?Wt4XK@t2L-Cs0*q4 zsaLDdYVc|hH7;sQYZ5fQG%sjQYw>6iwJNpVY4d6OYS(Db>xk$a)VZv)qART%soSdi zT~AdnMXy^Qsc)p8r9WgqXW(E^YA|KUYZzeIXt-*mV3cIkZHzWHGd^Sd+=SD_$E4n5 z)l|te#q@z0g_)gMso6VoQS(UiP78#Esl{1~*Omg7VV3Pyz{<qxtkoN9A?pb1TQ(RQ zYnw8g54JM4$+mrV^md+h4R#y$y7s5-UpWXn#5mk{q;_<3tatqGr0<mHH0><qoZ>v- z!sZg>(s2NL!0|xMfpu2{*8<l$HwCv0w`cA`?(y!AJ=i=#JZ^hZdlEgHz0h7xUJYJ5 zL@Qz?aoyX*yWD%#N6)9oXUSLFH{W;BPs1<IZ{A<SKhJ+5KqKI6z+#|QU}4~LkX}$} z(3gY82P+PK54H-f4c-fJ3TX<(g?feF45JGR4to&J6&@cxN)jVwkftM4BMKr`Bh4af zqM#`EsE%j`_;5HFBOG%qW+qlMwk&ot&M~ezo;E%#elS5a;Y7lGqG95tBvg`5(*0zf z<fF+mDLN?^4?%~#58Xe^cR2m<+!3QAbw}|>gO5H*l}^n|{g&pG)|t+kek6VNn9;Gu z462OCjOUqZnH9%T$AgX!osc_Ga^h!}Z&rV{bar9(&Pkt>{ikG36`$J6@y{7Lt$6x; zE+#iD_t_b(Gj(~id5L+mXD!aQ<rDI=^1l^$7CbJLFRUoS7sV7!7n>J%l<=42mh6-U zmQIxEmNk`gl%FX7e$MyY*m>>qS1ULxvMYXE2)ZzN(dc4ZrC?=Y6{ae->O-|-b>AiB zOZ7EuHQ6=0wPCf>b+&a6>Xqv28#o$HHzFEi8W%6ST^_ljf2HH9#MO&UOikI>z_plb zAFq2|e|E#{#{FjH=Bq6NE#<BBtyyhQTSD7vdtm!?hfBxkP1BqAJJmW{Z%N*&?c(Vw zyUlp}^d0=2^gDZZ6Ys8fle(Ah1>Bpx?|J`Kk8{t{2R09e9-2LT)N9z=)2GvS_mRe< zuE(m6Z}uzqw+|=|v<)f@wmwmK(mJF#)HbX%+%cju(mAR=dS^^)?B2Nkc<)n_rvnp~ z6Qj=@o=rY?dp<qsJGt;8<i(elu`jn@9eIs@o%M$1P2m*BROMTdw@uTE)3;~zX9nNd zy?Zt5GrROY>iy1K`UlDnh4b9=^$W5KU5kc`BOhHq&MlFawwE(L(S9ml5ngHeto?a# z)p>R9OXQcmwUb|2zt(<}|8{@fdVOjmbYuH_)+Xy_-InrJ-w&rB^V@Md*q!2Ck=@Rp zW<Ou=h3@T<h>^r!X8|PKBK-XUSgVC=2p70^7y$s^{kx|4WdWBf5NrcuG4KE?K)4_o zk?&F0(MK@OSQQ)_8JcVhzd-($;wj}5sv+v9G%skU=|0hKFp@KIGRv}<uzIk?v7hB= z<{af(C$R9S^Lp@Q@Lv{~5Zo0O5V00b61yZmDzPgiBJC)1RQ8J8qykb=QprjAs7j;i zQ?(rpeoZs22<;-BF5S2KNCQDbW1}GBY?B7lezQdjoTZ?Zfwi|yifyS~yZxxcCnqvz zUKfo6POc=kZ1-x9F3$<C6>p3Wm#>1Kxxa5fV&Ive+JjxeV<GdQJK;1W{s^T=vnY>f zQcPNGeq42YYr@0Cr%CUUzohIPrZ~cSls{EEO)cH<n01CrruXsS6H!^o*=Z-UPvzwl zpRUNQKGTqQ4Sol16?7LqEE*^tE*USKEPGQveeV7Fg^HyMD;L))*Q>r)|G2bM^Ro`r zLk-AA)MeBa^i^~d`WpH=>ISkI(E?ibT6f#F+c!HlZmxBHzO~dffBW5?H+Lty$L<Z> zf6&wQpygp>@5R2tN7;`L^+ybN586G^9g-gA8o`fz8=V^K9l!E4ZzA@Y%X76!_Q@YF zCSSI^I{iB6jm{M7)au*b>C&07ciOY`vvcp;<}yDx%?r<OE%YyzeGFaFTBci`|J1c| z_H*#6;TPUDur~em{<qrooQ=5e-kVliT0dmA`FGfMsefYkpuN3+#h(uO2$hcx!sud! zaAde|WUuf8<hLlAC>yA1sH<pdX&dR<=zACj8Q(Ciut2P=Y~t*C9PXS+T&3LYgh`(7 ze4PAR0$zfdLYIX{MK;8^#PuYCCC^HAO3%n*<izFe6yg=jmF_9es$$iI)J-*lG>>c5 zYWM2Q!i5rtfwG~EQHb#|lk=u+W<%y5EOxEvtwn8gY@O|b?T<PXJ6?9W=ltAd*%i1k zx{G>fdfIyh5>vd-_+0Y6>G#BcI^bIn_8?oZM2L2%ZJ2L(EGZ+R5YD;n(S0$`V&BKD zC4fZAB(`Ls6uCp1hfR*yAN5ELOp8cQJeHnuDl`9h`HAYRE7>h4Z=LGN=|4T1J9*}9 z-iNa*a1PuqLKfpo=t^132<1ZOB+n~WXk0M3XkKYw<y!51>0nJ{Z9?7A`r{3`jm4KM zuGC(=+SGRK&h_3KL(R`xrdmI=t+a1+{OrWsqUvJ0&2>lcu0*%OJ@xzgJ*E$A9y;~9 z_jx}GcpTDC8i*c@dy+7eG@Lw=Jeo9?Fdp|bdV=&U__^Pt=L@HoR<8_StG$t);(yCL zO*Xwfvjpdq{`Vbomp+`CPhJREwEn2P#I=lCUidVy()c-RHQ<ZR8r#|`9QV9+&kd1{ z)$h%lQCremOFu4dJHR=kb=UhR)z6MS_h0en1ZfZfx`NO}j3Wb3I8-It0KJ4M#2UgU z*jh4IJQIG5{4|9nB^~7$RROgdd{SJa?V>wQ@5CU?u))~Rbe=ha#ekKQb&ai`y^155 z(}7EZn~ZywaGR%yH<HhSUr+!NcqP~(lp`D<q9e*Gx+>Nueoi7v(nN|&YE`;NrcgFW zPF0><{*^+bVxp3XGOO~uN~>zBnzcHY`m#o+W~P>-Hox|o&OO~+JuiJ}17I*}SY;Gz zY;3}A@)<r=l$eKG7+P{zu3J5^uCX~{>uRT9Ph-FA@W}C!Q>ru3Me6|PfgRUZZnxaa zJd!=#yflbh#69mBpB~?8zYPDt0MkH;AnKsCgA>6WA*G>*!o0%uNJ0_#h)<CtQO(gs zF$uBmaa!?w37CYBiBFQQB^RW`A96daeMA7>hd-xIq}@udJa!@@BGcu#?g{ZM=4?dv z%E{NK`g1x?*X0(U$;wMPOUm~xa4EDdGAh<8Q7M%xlPniKCv;w*Lg0eHMZrpuDv4^D zOUgA`wMKO|^==J8jd7PVt`uKwY`S}G;`&N6u7$r<x6P+LqoeWWc<0Y9@!Rfq^6n1Y z!`?UUIrU(&SGX_b(TjeKfx0JbLxm%3qxIvuPv@T%P1?QWc(wLs;%(nd@9g;8$^z3z zv*lAOFTcoq&EEi<!9U*Z`tQg8H@CFECtwr-D1x~VCzu<ts{|mk6sDDSV1B6t0BcGB z>{S55r34^6I$+wx|L^?+0oR=5FkizC1mRj!8?G@uK?Fz#C7=neAE&@3Tq#OH7Emx; z4PJ(Zpfv<DLIdH8I0;vF(?|-W3NjE`gnWej0at3?sAAMG8jIFJC!;&y+DicwgK5Y7 zz-nQSVV~fba75e{xOy@rD<fOM>%tY%D!B=HH3dZBMbS;kOPN8rOl3!Pn_8H<fCfbq zO*2pHN;^ttO80<Xi~cr)2E!dj9mYN;OQz?{-png3DXdhi7uXco2H1Vse{z&@s&KyG zO6DeT_Y%T+*m!z*WB3I5UhwA#7z*qO-V#a>RuSG2xhHyD%vPLB{FB5T$um*`(poZX zGMloKa-H%O3Ym%#N*>A<Dq5;?YNG0V8r+&(T3p(^I>NfrdTRQ{2F`{DjSd+Xm^7IV zz`GpYQqs!aI>n~W_L)7#LB%oLsn+@30fMWSTb28wr>s{BalnVem*m$MKnP3@dL678 zQWZuHXSUf$`>01TTCq3dl@i*M)RONXvN-(uXjmF5{d|Uc=G2MwY^jrPbI#^k<dL5p zE+{SvC{ZtCD*t|d>cU`UclE8BuDXW}PcJVu;jT+GJGW-F_jDq=Ozz}&&-JK3Eb3eD zcOL8-<{Ql#Uwd|7vilYPn~b*~-s!$C`|xE^Zz=E7+f|V@(zmvatu3wXwB3il^2hJ^ z)5A=j5Rd~pz#3+JqCh4nhgpPvFbj5Jo<ISzgQB3bP%AVE?I3s&CUA{kig<unMRFp| zkV(i!<Qo(%$^eywx`A3k3!uHw=g_Y(Y#3L#)}Fx%VZ*R(IDj+9mEac0RLHW)rtu2+ zQ}}suUGfU@pA?=H_u#tsB;`7lC)EJ83UxgVJI!fY6m2r?23-W*3VkU3XNE{PSEVpv zm~xppm>XHtSe~#FS--R8unV*Ia1c2GXC;>b*AjOgL5(oaQ^2dwyTRAM?<K$@FfN!Y zWF<@^JT6iqN)!_p+Y)~yQ7RcGr7z7Xy(RNpwq34N{;)!TqK%S<vbYMDDxDfm9jI?> zY-w(4{m}lYi`JvkXN6BViblr9ZYGhYS!S2a?_11RA+3dMENo-!D(nXxww)xLy<7?p z47riH8+fF9_7cgwt$d1n-}y@hqy!Ef<PDAs84eQ-&mhf38b&q5(8R{a%_dkUb|=fF z)E(wNQj*G?R(OmxqvE*miI!}=lM^|?xu`Q2&+6pQ6`m?qDOoHlKWAS-b77>ipxUEG ztai8lX=Brs(@jy=-J8u?HQN<#%HEQ@t$Nq+o>NcQ!|cAM$1eu)LwY0WV?z`C&(mM5 zzCJkhb>`H2<qsQ+{madto7Wz%f8UbbN&OXnBp?IyFkdGPbBxY_1anDspbsp<43RwK z4jqSMJ_o<?W{6ZoC(Pf-!E8+nauengGg0Gc0-A*G!LVRRm;tOX_BeJ4XM?*5b0hh9 z3_cnEjhsZjN)bh|MVUrLLv@i_k$Q+Gh!zK*L~Q7O(AO|HGLkXgWr}B(XI^7zVU1_g zV#l#R<*4F}<}%{u<X$I?@Lc9S$rs9RBcLqEBSayzCHzt3wdknWWAPq|yOMXMx}_h= z49Y%}dnf-z5vjzYEUjXy8laY`-k>q8xv9;kW2>8{*RH=}C}QMqTx0UlOxirrqSuPi zI>6?(9fN(i!xK1er(M3ey16~{knyY_QhOiq`Q{ho|2EJ$Xf)V7<Z+l$cyEMZWM8yt z%ut+N{EI}dq=l4-L)%9(Q|Z$#9FxkpbKLU8R5s}pBBvl%_)OPXoBR)jNyRiJHD&7M zqvr!Eb}!~t2~~I1Sl7O<k8Q+VF1aew)N$SL##oD6>-+ZLj?bMDx4zsCyR*>kd2hT& z|H1WM_P(Q!-t}t_ln#C#(izSfc{avA?)tQNV(2+yQs#yG%d}UuulwH2PVG$7&hWmI znN@#pFlY9`V%~DWY|-eW_LB0l#3#ZE-Dl|Ym(@349<McjJ@+krJ$S?VyTT^t7HVtZ z$KZC;PVR2_Ps=^Q{rE@u%gOHV9sr`yh3Utsy}hr3Fn^f`!0ws7z3s}qz1=Ff*ns~7 z=??k*f5pBr2Ht1NJ@!$zU_Q_O{omJr0Vo22j53V?*8l(j32;bRa{vHbMgRa@MgbI* zLqPxlAOJ~3K~#7F?R^V)RMpk?8%Q7#5)8=200II+RDysg-UTlOQ9@O`LIp1btv{Dm zMNs%EUMSkHCDu})YWOT*@zS4A{h_UZpeTGcRDwzr&`==?C`f`Nkc9c)b<XT$GMUMB zCUa(HttThvoPAyU-DmcD_Fnt!9pj=vMG;U06ahs*5vT(K&FT<FgD3)ufFhs>C;}D& zN*Z;15Wd)6$AD27OtJA@>57YQT*Ote@sEwuePTTod~7_he;?4HGY}Kw>7&<M-v<U> zsyO&UK(X<K3ZS_7#zkBe8~@le<E#9WFZuTmwNDk&0Ewdg>A>c9fLreeZn_6Jty8E~ zbVmM<O#q(Ga)vHW(?YTFT_cK%Z(PJxv8fpw|GuI3sWv<{X%=v>0Qf!6!@2;Ie+q0_ z2dKF`^p&L=C<2OrA`of>s*y&w-oVH^fmrrIPL=~3-U9Bu8o0Ix@ISu=jvNZL%9_6- zpa>`eia@9ks74y1$OFFu%1Ia^<z;p+<C{Ie{29O%9Il!^7Wnu>+h6}H0*Zhlpa?_{ z0zOIOjGn*^B!;*aRe>ly0hE;jYgPbbF9tG31a{c!U>@-H8V-+<sLBol@2&zqKUh@< zy;B4f0YyL&2ps}GNki1R|G(I9v+?#NCr<*Wv<4n}E^xD~N=O((Q`w6l5$$EZziA6R zm=3)EtuGB~cPrTqJo+SXwDzDHLJ?2|6aht`jtKZAjjm?`=|8L9+mLj{-m&-}K<{(w zn2T4ivJ#+#hcEeaF))B0AMkh}2e_wSWhjZ`hsP=dnX6;KTeHoLo)rN_KoL*`0)c>E z(vZZZw;?tSZyA0oaNb4Uu64WrzkK^YF9A58$GM|{e{&d0dKu@lmmxBa?(=Q|zU1|S zLBLB}fluCGSWX9)IOp-af^J4OEwTN?|B7)8%k+-&L%j3T><~kkhxRrMT6$Il6aht` zaR}5%8lAhdw=teQ3%6}UKKL8|F!1-+fw@!aZ_oJOXMtIK`=7-g#QJ>T`mab5R{}-M z=jCOvWKk}^4|&YKzG@*bsw*&a48t<d_YU!+7aHTAB$)`Qe7>A^l?m*TtOowM3%H1J zZ@UR7<qWWj^nMM*+sWyFML-d72m)uiUlJXnj;c`)1ZpG=N#}kx@-p>S<gq`{j=_e1 z@!x@!OMn?a^KbXCB6$Ae(9)l#16O79!@+_4h_D;*tECLUhW;xKNu{}mp$75Xi=-eO zg+4U4w|BSr-;*T{v)40%r19nseznf)weuK=_jeBml9}F=7kIxO*hyo3MS6WndURG% z1QY@1AaKtk&e26RsUHHhl12&};~Vb=e*OgC;_tT;KmT%I$-A^W(5^9sAODrG(XUn5 zknd-pv>ufq<GTZA@^<^)N?V0B8peOtkP_n{5qN^VkS|C$qdIWdDGA8hKzHW(J_!eN z=@!Ey4!DMTiSV+*_F<)Gx#~YfKoL*`6oJMfP^XvhM}BM1mh#LjmDcJ%zYpAZEimyX zz#sm`_Pvc~2;bKSX`nVzl5Y<)_3KU-Ui%W6dRxlg#|9UB8he0$u_rNY0&p377iTfT zQ4V3P%CLJp6860hfi4N|l%%_&q6jDgihv>zFa&DtX-G~oz1Ek=*W*GQ@|F1bD}e8$ zE8$PrW4JRNSUQBpZ~69u!&NWO;+YM55Bm9&GNySq!@RV`9)dcbLsud#*qyk9LsMT) z0{*vv4(w9A{tFTc^Qh?fuqU!j{8SVHMZk>!e|UM@X-Et~AH0ZfecO3d6Y7$e$2Y)0 z>Quy45m`<aH#j0VQs<zh9~`X%k!qbtqo#~Hb_QPA$VNP#M(s=Zjm1~~L*JT3gEk3Z z9EX|CWp7{s?Hh|$({5p+(9ekQR7(;=JZ%zU{Qlh#dzRE?43jp8J&hC&N6qGxT$04y z)Inxs&0ETFyMQy;e~LW?XvZ`%=*JF-iXxy02m(dhfK>j7K4YWWAYm){{4{p&ch^7m z&R6*3UEw1h(@ag})W?SE<Ri(V&wxQJ=e8jI`Yjx(bFk9)Q6#50v3YI0v}wzpz*=DH zjr>xa#PKKEI(&Y^kqrEr&+ISx_Fr5|V|KU7nN6)Z1oa91EFl3ccU9{3(`QH!d`=H> zXe%xWIGyR3aMM<@m&5c*NHED{A4avk_*MG`=260z7vtvk(vu<(JqXChO|kq1m|uG0 zAPxoR10Btes5U701G0WZ=2yfY)l2jFwl=P|h?V~Gm9Tu=bqe1WalF#q^NPQJwD5G0 ze0d*l6!-L#Yx<Yb&!2Zf8#@Bqg(_lnF^=-lnBY@iiE}6HoxW!}-q?Hkh4FQxdm5(f zr?j%|8iOta?xI1urnRZsvzKs&2Vs~rYC5+jAxNAG!te;3%=Aw8AdCtxK2Y%ucpfU^ z>xClV6a?gBrvCleb*2GN=Dy1@bGchPzM+pB20Y7;pv3*zEkMD04A1L?1ioqU{QMd_ zY}p>pD4D=1hVHletxUtD7sHR5w(t$i$SO{ijA5%NX^7uejw>%BS#a1FQ%F9~Z{-CW z-6edZGii*Y4l>Qr5{942c+$YeP`^bhxa>nqx1YZg<1UwSNWKz}WnRUyi{y!}yq`uL z#9zv_gk?2i=}R}!Wx}Jz)g6!dOFW~?syDcdH^aqA*@&sET_lpX%Erx<M=-X$yfFUn zPweo9wceO%=ezwduzOuq>=$|H=#CZEC)p$1$jtHPs^RDIa+m9gNiNr*Q(Uh1RyRBy zr^O5eHcl&3p!88}d?Oi9Tnck!cLE+mU9KJbT&|4^?fcBl%)HR$O1E||aJd$~?sDz= z)b>lC=W-o5<Z`XewC@)(4U<j|-S_baDrwHk=5@ZyHHT@-{k~F{E7zU=!i~Job-Bh$ z`g3^A{3OiOET+wTmrb<qS2N5mhE1d2Q06C+nNGhn=1;$kS$4RgOt+d8xm^7iPjt$s z{*wOd)b$X{xlq!lUd*aOr?eT&YZucWTj{rJnH?{aW!gm@hcXRGN76K98A?2MF<q-H z^vh(qbKMwZ(p@r<IAf`!;I`SFrljGG@#ie7gx$e7uTQnFv#FD-*iL_J7WHSm)y&h< zLkdtg3FvcEw=C-wilwfUKq@<brpwho&gFXh)ynk4`6`gugfkAlbXRPAq2h7{6qjN) zYSQ?uWkYIYUK$%|X`tC4Nh7&1*Y>d*N%~1!XC#viuY@(}NJF?#8g({!g(QhWHg<;? zU-DSCgAJ+>2pc>j`_{T#GqM<;;Uuq2$%|<SNz7QyYo<}kV<U-08u>$PKqal$^Gym& zmyNuoZ>5uc<v|+7EZ<r-`j$Qi7@q_I>d7M7jX%ppXG==m`bt?$`P^QvL+(bsf;1S9 z!7Rygm)nGuxalMaDZ^&Q+sDib8F&ZFWi|K`k#3$C8ytrUPAunYQ%-ldbmFjsZl%Oa z2u?7rqXhK1sbik?2^O&7Devyql&67Tuy_0DSm4zK-oART2q*##Lf}-7sfK-rR3`B~ zF9&z8yho3!<NemCkoVA#q(%Eg5{G({I84))hGUVUBm8=IJuv4Y`(FO1aCqmL@f`1^ zO~VqNi@vajs^*Sx4;L|=%Q-s*ojGd7G(`H);AuqmOO~OnnxwRl>jZJ;i=_{3AY<lw zg>F}Aq?9GME#X&PqVGLS-$(+-%c}_E8T-d^<}tKyrLVy9)pz{Iai|@fnbVnu(<JJh z%5r&1HH_)Gp81OG`PM$mhh;Lu=PTln+T9-hdiG7*{^AZ1?x&!Uh5+-kshs)rGL7F; zfQ3(mJ5;K_BH%CtzS-%e5F0&_QZY8OVPqR#A(s8w`d8tE<NJ6$Pu`R5nJJ}*fYod? zZ=~Ho_$X(0W!OC&r=PjVHclJY#7_3r!W=t|Ql@ht7s#GSdxSgC4<#lvztE)(#LY{U zj$T)9KFTF0dec}>T|}=PSq#9k)F@;FM!~u~k5GnCJe)i)g#fAunWNwzw1tqQR_`E& z%-eOq9m24U+`2I~ao6T^u0I>PNkeE4k+RbEaij}vLyTWYd^av*d8PxeZ8tH({v1@& z5Oidv-1&bcZcjICy{fQBSPP&C1QLOV*fWrCEW}_tX+3cBfa<CD5I6c6<98m7zEgNV z?b@ngFC5K#HfU=|4)IKL<#Sx8hi})?pk6(LjsFCGExU(A!Q=6(_ycC|!Avaqb}!x8 zqhK&hBAFDhu@ocwFUSh%{K2<oE1pij%b%f<nSKRZc)ytQ;N`Wx201Q%?bz!{WIQ&s zk5UH_nZe~(sMJqI-#cx6<tyJgB(8c1GcLB)q2M&kB*tVKi+*9J@WUbhDYfs&I+$sb z;jP3E3_D~Hu_wtUrtn*Ufp*-)IR9ycn~&#MRxroXKHA`NU6AB*P3N$t41I>11QVNZ z!@^r%#l{;E`#zw!m@#zw{o5W#7elx8TjDsn*cL2<ty>pk{5IjI%9uUHq+^9;WXsTO zz0=!iF3YvUi9vUsjQtxwUgx@v0?XWn-+8n>6mz`Z7~!q347ZDMMGm=*<|YlxNM4H} zTx=&MyqxEVjf7{hg;=&4rqAM9V#<;$4m@VK?F*jq#m_rlme^K2^bp%ldKS~5v1Fnt zFVE>Lr)4xxWAaSL6XW7sTQ4hbM%_VtC10@z8RNInLk#zfURMd|b5qw$>k*Entdw9Y z``K2P>#`0m*Zm{;tR1V&Hyo~ljZHXU;7>=z#y2tn#igQwWW!y^#!Y7T)F#CYlb_G9 z9xaxq@z}&a#Ps&@DE9F-Co_P2<s}Aj&*1xL!1fhR$VdDRxrxN&S<|_c`br#QJTB(F zFWr2laloG*z5@Gkv(hk_v5;bvudhfQ3HE|(Jtd&eO+E9jOK{J@x}bNGJlCs|HLDi9 zRk!u3g07fmpa?h}0W(C;4_uMRkCHfgQ>(-$1KoHSk?kz+8|%LF`02}6u_o97j^lTb zkMP`*Y&?vE`|SMzR0@Y;d<ut-eE9A9l&e{v<1+#c)+}6#M*-tKuKh%*aF&n28Y@-O zAWn%MSk*P{y4F!$({oU`pxZRO%PQ+tE!{*L5r+(GU9ZYJvfe8K(TYGqC)%ttfv(2F z+yEfttL(FG=GcD@Fwor)7ZlKd%<=aHoZ>inGL4Um8<21vqsRcxvSqY8yoapc_=vM| zB@N#+LaR6$fd)`$<IB^5F9Wh()%_#w^{W2#Pp(Ycqw(o#3W`9)A`r_T%_D5o8zZsZ zI4$)a8qqn+NJv7?(qrh`mj>JOxc9+rG(fj)Ob#{2U9p1uBkTw6xPz;e`UqNF4E*Om z_BvEPyw%AW2!w<*r0B9<)r;@cVZAD6;G|k80^y0kzI`@12tmkk$|*KU3>X01emjSy z{0u(@y0lw3dur=e;8(x0m)r1_zUXq8Ll7#CMj&LQA*h)3svhN%FtT3Ng17m?bCIL* zQsoo@XCtt8ZzV~5{yEU9l}!?ZI1ThpPDwrW)TY;8bdq7IM;@_B%bx~98ee_o&sW_Y zf<Q<~LokSI0X+9A@YK(Nd#~m<9d81c4tEGTs)i!qbOby}LV6OdTiYaY{`p*Tnnv$X zkc1~@$Byxvh39R;>e9v2r~0*!#t%PKcUKojArNxX5Hy<696$Sye7^Zj&x515TEQ#8 zbsSz*Q3MnLF9g2+x@u1%kp^wCNemuL10lQYZQ6K62=VQrMSu*2U3g*biU`5!@Ko(M zg*b#GNW)+-`6&{|Y4&<mvK!#-4;UXkD*{o7fXp`XCW)jZn<Or{fa~=9%y}f?Sxqt< z>iBWs)>}PY{a>q1!#_3-@^KVt2rjZ-l^DrqbCSoA@3}(3EC->bswe`%M8FqGNY6of z5<`Y?j{0~mq({;S2idlb>n?GM<P%TWLkV@!S@ybU$z?`<{L!|>)FqCFX>tS{M;d~F ztXGx9;mu#(!rm$upL}$_{|OgOj+kbs2m}fNAqkmcYGyWl@r7-J5R$lvGeRbC-6a`* z3OA7t1YdlS>rnlo{`#BVh8XZw9EpG<NkfntdIP^R`nzr8kfExlviPB>W3_c;-KtVf zL%^FPWOkD@l0p&}U(C(!@2-C%>$EbSI00m4auGj%MA)MNA72@}Mr}rOr<FU3^l@}E zd_Je)E_oN2_!EAtIvV)x>--j}5g)`xQN3zWoq&LlME;(tJqh`EPDtXCOCp0Lym6Ku zM(^Ih6<2us)a_npP}L2m5QcIzX$UTT&$E{ux%V0_o1^Phh44|)tZPNUB8mJ8k{Gpl zC-4%d1PV!9b{XfC-)m2atu>)A4X|bn@EymKXU=SR0nG4J!&A_hK}Nvwq#@-$qnAw_ zLK?DO)glgCoz~GthO07GK;(H{;TWh3R}4;`tjsG=$wovv(Ta^1NyO|VNt|{XTz&fR zW1~y4`d{~9@Duz3+fz=S^t_9lYuaUE{@?$D>%acXHhkAsPJhpg82ZJ|Q(OEdMlL2M zrXokr>k3B)U6hxXqj~e@LDMKKNJC2Rxn9-Si(Od!4_C!{RVVp?8Nx7UgP8Zdaj)<B z<jGi+mX_5wT2zFO8f>H?%g;wM8mwbzkBHf|3obWFxH#P862rBJ0rz-wd1gu@{~0)R z2+jZS2b|2vfJ;6utV5aoJu@*h{5NBI<#k9Sn811ygWNZh56URU&-)xqed?2TAHK9) zcinZIM=^n)$PB1Y)H;O?8`6-{BqiHqG3}-p{Or;gEL>fVGtWNBmy*i2Hz5gOT5oAA zd6b`Qh7)4sQ>lhP#A3teD#lF`LKMxu{L&(cva`;@NzQIMDc?zWkc0({$1@VKxOhGR zWo2mb``_cl)mNij4BsB*_YCHJ?GKw6;u*&wu2B5e3<(%r<z5a$7SgVGx*2x~?;Y2q z5sLG--+s&WC>P-Ote3O+tWKIVi7Uw5$rbNA1Vo*%Aq{Et<%mx>iRV{x*z{*Hm@u?C zGXHNm&OX1)^7HI*7$Py}#-eqyJdWQf9dB$h5g9*Xupvn_+p`DFIm{HxUW7#w?c1ZQ zS1**4B+BKBYgzEmQ)+l7A_f=F@_Jpj<eE~Yv_q7L&87w{yrRJ$HZchaX#V~8A^WYF z#^0QU7`%x>+@vj%w#5ts$-N2ZUDi<Cn-}LV__FipqmSaDhaSRDe)1FCci(-u{`%_? zC*Q<0P-0lvih=sa*t>o@tAQy6F`B3-ham+4L^>C9CpX7iOPe8~&o5EhVw0tTM{#`K zH&o4%{09-={e^}`kLY;eicK><QytjQb>Y#uqy(q)e3<jP_i?%Ky*vu|q1lN>j17k? zE_Hoa&nYTG|H4A7IQ?{#yN7P-8n-5Feoo^4;~8h*q=y}$uIUCMoFNR8mgjwI4!_-h z|NUG^o-=ePLLlY{6lpyBORnws-zN9^pHov)aZJ|y@)ARUt=OBp>#n=}4_}q!VOl~e z`sNlQqn~^yRD|W@lCWXK4lKHZlRDVQR)|x#86-)#%F4=I?6SIwxy>X2aTOI6xegyb z>^iT@Z-w;qPz>}L&}TrOfqFgz!-o%d9XQ~%C2e7jYbevntk`qD(6w+VxJ`SG#TAYp zo}!i(U&>i}7;+|oSchd~Jna==uM91>!r5n^?ZvHuyFK>R{WzruS5;K0hJaJC5wm4S z_8wBWxsM#3q$j~{x%4CoI1KU)&nN2nhFq=iH8?mGml{;AlUD+Nhhw-;$Ge;a@E`Bt zx;mu*)P&8AOnVd4-Y-jvI7sq1LyRx0mT~#zmtVGqqhi@Jl_4t6-iBN^Mt<`!@X&J& zPvP8i&-K~UkQB|Z)TmLTI0O7%TyxDe6+>(lIr<s7;^c8l98Vger4WU&O?VQ9TpLmm zpf+#bY?mqs|LvhWj^C+NL%=cEurb5#-4)+SfQD;pe`@(oLJk@Qt)1$sl^{3{j=?36 z#glR2dGo+GG(4D|-r$sL#O4JK1It=lvb?*4<jO0r#D)zU(4s{Pbn4W}ItyV)Lo7rg z1Y*deF>0U65ut{BOmFh^98MEA(N`uG$dFnFXSm(Q9U4!qWyyT+70x{A+wr7f2t!Ci zdKz*FSx7HKfG}iY0Mm06A;9OrA_lW|lP5_GAKv675f#Y)<u5kzHi9r}#mV6x_(_{2 z?T@ixNE;+%A@1hzZn$eLbASVw-cB;HlpclHFvQkTC*dsdLI!LZXO`b6!GZO3;Hn%V zjwcOK(o}>eX$WDI(0Hk&AqYi@2uXPMBxIMj&YkUE#E230@)aC%QV9|Q(v$!AV|(S= zAZ1}CD{sLL-i%4w9GOXJNW+{(_fV6C=pe>=+A91HKS`jPB(W=K)4+mP-@=W(e}cTl z_hCBk)7{+-N0Ww-g@~NJ2}64t7BrFlPd@PSCJ9-{PKKKJ@sc&%<gjt?h^-LW*JJr| zU^<%*+Jk~B($kQ=%^fO28q)q4dxqT0AwyN>1vixBv6EBA6C1TUa(mpddJblf7=sRT zSqW)Fk;YSEOWWAtXwr~pMk>Qpgj^f)kV+G_H^d@gRl*zrZ<5H*2MkGYZn;Gg9G_Kb zVgzLRCO;IjwvUybU9cq#hNont!x+HDreWHtU_(=v^czPQb`Cv<&q8V8gE-92y$o@A z5{2h=UG!=4FoqQ5EEOY%*fgZlggA_UlPjhfHXZ>X34TOZv7dyH1n1{lJ&7B*HomNu zswB~PoEwvZ^!ESAS%ZTH0kT_QQ;5S!!?Z#20&A{?K;$}-+N0RSWVH8D-;O2?&kB$V zkwYp<NJGwQ2=`H=PImGp34X+du3c@C;1ob52?ti5Y}qRHcFQeIr(772hP<#E(vWt_ zv{mNZbYL}oPL3uGfkP_6Ft}8bsVF&nPPuJ*=rqS#5nykyN>4&aqFXncBxLptm*>!) zM6D=>F2KM44ZQz8@T*@n?bexMr;=|S_-#$l?I#UDSK1`g{+Ra2v{B+0>BQl9(vUj! zt_aTxi}VU=z%w1XC_)mlJh+(xE5=hp5<mT^O%jb*^Col!i^81u?+3E7fN9f0vzwsl zZ3Nrs-hir=HpjfMnh_95CR*B}upkXnAyO$q80Ks$t^xI~Q38R0n<UQ3&*P$HJY*k) z?%nO4#Emy9NdyApP^7^LCYV1T7&E4EUzs&VcS9N}DUC@bXmQO@xgiQUo9m!u>X3xP zf;3EJn2Pc23OWSG=uksQf~#g$^d$D~wMgP*45z@}bQ93Cr>s*D9auziDgX2*oA_nr z+E9suztxa-$n-Ge+Puhy3aX~gH!MiQfMQ6)v(k*8o`ZmZym5Mygsg8M-$~qja|KE4 z_?&suNhJah;O=`kd>CBHHME5>J&n+6;Wu~h7guSMqD8`rG`uS%nkve>>^0tVmYFq8 za_Leo=E&Z_z=6Qyj{}1T*BFhT*jKel$Ol%2BxLrC^dz{xj!J|g@Y!d;KmK7GzBz?7 zq#~qIP!LiTYRcTSJ?30f*hur!79i4URNi=G-KF*GL2I2|I{ERv^%n3y|6_lkWy&lh z;n|araA%%r_att;)xK7eFla^gnHwKr(IVjIKL=N{3w>FIG(P<_^m^2wfj`0!cR51& zRd@|*H`>ebq+u$?RD`c{qfv;Y1YUfxDq%?3q~{?raiX<4X~i^HPm<uumsU?A)Nh!! zrux4r5s>dN<csf-Bkiy5n=)^Ic{$vsVbnJ5liZuO$@oR;+3}<islW%VTyMNl*~9S8 zgfq|Z`s={O7X!E7ZucaVB)n0IzI*9`a`!!4d#!WI?r@uibAW1kO_Vg69`oSlSkpML zd@v<b7-fkL6-7W0Sho(4k1l7;auy!a)8K|je0^`)4-1`zf+MwIE20un2wZrfe`B>w znRV}N98qP~!ZZv4SwUEq37I_Ed0*zs{@}7w>yaZ3E1Sj&8UjihL4!6#3D1~eFTvoQ zy?j4$>7_s@=iYnkr1zmiKpJD&A@9EXfcCzYr3;;E6*)_9L#F?#I0^wJ4M$Z+0F;qw zv48s;(5@Z0b~Cr%;yOMeLKd&RW?$?7=te+nTYvZiFmx!;ub*?VP`idB>qc!Fj;xkO zDkTO@S!Vph5AC_`(lE%k8Y+r_ARt?|wrT~)%EC?-Aq@xb;OAuMf~lL5MljJ0W!l{G z5!YN3N)0q`CnE5{2f+62_6K}U)X!5V?P++T;h1Y3o^ni8s3M9$kP(mtC0~8jHhhP} zUg~Po2=W8bu;i6A8it$3QUt;afggVao_`*=>n?lRbvS6R!&8o|B_$0<R*Nd72m}xT z`PM<UY?bNQ=bs<IgoB)v4o?Lc;D#lyHjO&q_Q}71Eo*_hht(l)6bAWnhhpQa+MR;S zpaXwJhvJi1eFj%1CU-bN{<8~^-_^M_ahMKIIiaem=@7&9QurET%ME7(X?z)Rj8HYR zZoS#9Pc~kX;{?~sDdop%AAYnNyN}+Fb56@b`t{?3qWb1{fc<-elC&nO2m}>@RBXld zScVty1d74Q5e|}$U4S81*R0U!4<_Po7rub<<}GRt*aW_7R<Fjl+|lsVHf?ah1s8DZ zm^kZ*<)WT)H0KoeW?Xea4r{SO`9QG&vh{t=n&*p{Q`ZWu($f&acsi>Ar3_M(lP6u= z)4q&HDOc?-!Leh<tmEj>qbMva<j0~%u<4i20!M+?Jr`2<@j=mK{FOFov~C0R=^vCV z!aC9X&$$;<pRg`MP?QffE^Ur~gB!kHg4JDqiL>LKLl~cL<4=~~7Ix()MYXm!kG93& zY5c&iW5=M>QGF8Kbmcs3-M;~CQc`jIL-Uc8l!Vrtr6@-Nx228e<boD0T3AFP1kymx zz+)4fOBjMz71F2=fYBJX?|HzIccZc7HJ4c)@>+lVH75~~z9Ch^ZQlR@AOJ~3K~%*> zmc4(1hVPNL;F?RjISOBS3H#TEI;v!bs2o5Vb<{ZGgQ+1r;;UE-7?lVtUTiOeFH6US zv#2rLakZzU;kYVMg%p8CBJlRx_CkMm+|kH{91vZn-8!HWs;Z!*Q56rpQv|{r0om>D z9nReNKmW&Gdn~+5WA7rUixxVn7L_y{RUN8~BG3o~zWdJJvh{%noUvtVBZ^WHRfnf4 zpm55Cl7>@iLiJDt8iRmr*(zUr%OZrmdp9N}XT??0a8`X#8zl|LRf#I32s8kJEY6kh z+7-CEmJc5qkaS3+$m}Uy={BTQ;GMaWhBqR5uLy)20%F71yB8QgKHM~~U0*eX*N#yW zhtTgdni%k?&5lb(;7K0b_Brf@yO3fQkD8XiXsSZ)+-YyRC3}3!?u?NpvXR&+r+{Dc zf?v}`TIs^FJW3j2S#4pc^c?EU4Fw}cWbvI<Auf$<*(y8mO`Hgvc3Oy)31arfuA!t6 zM3C#Aw30^MvC@!@L4YrcmT)N+g!57?Tml1Zt6|Dx6)FV(WXsmy|K7G?oO4bH6llUM zja{Rgv%YL^LIfO@m6C>|sw23{9N7LI_IDhNejVb2n^8@vb14?QTT=js`&~)>>be9* zFpIK$IlzVUaQ*ec%*7dLo1TU<w268x)DRxwk3d=8BE$?^1g=Dc{W*WdroCkmUg9Vy z-lk1}e8hL(eNj+Y7EwvVQT42((InMThV5@WiT+&nGA5=!o?o}e%C+p^7R(&h-&#g_ z+>@)owLGzV)sq-KZW_i}?xV3}OMxA>;62P7ZHMhYdK%s<APtmk$4rKunb9AkpWljt zElV)EzvK~v(bJY-Pe~pgPMrj-pM?9CY)4#+WF%R((ju&RHp2=R(|;UR<w?#!-tw6k zH*Kb!2Hn@~5voyf6astr?EK{~z?3QAvZ{_!EojOpX#@@0h9s<{(Gb+?8>eIkri_`5 zS^32%-ZcwTu9|@DB?Xw<VF)mF1&WLF(eI6POkSRcW8dT;d+9bz$a)|7YX@P%kOvqR z$a(2OynXQtD0P+McZ25R!KI&3i&AV}Iv=0)d=meCdkU5gnShI5*auf>J_gO7fEAy0 z#Dm#07;Xk0OHakIy*b$RV=4B$@eoE#Z-*U)#n?5kE7DUZBaa=B(xS~+I=>jR-#LJ_ zH|OH2MIY6-q#8b$2$(HfCrkiNKRuW^IV){r({PqnQLBxTMpIOMY#e$spD~X=hqv~0 zMgD<o^ef(tH|W;sL_XeqcLzF*OhfK&T3C-4B5mGpFsxf5x(%C(%(QH*-&uq~Q&(Z} z`K4I7{3Wbl#CGzMDPIAT5EeW#3O(C*!u|QVIKMO(%U@c7_4JV-oqKfR&RU(&IoX)6 zj$+M|i!o>O4D?NmN6+Cin8#9V-c^FqBT1OCaW)3FZ;y+|+yl)0(CYm#8I@2YAbX`P zSO8pf5ioFIsI_U@{I#diw0PH&o03K?00qQ5wl6aG?LhXWz|v_$(7nS0_}4dIBNsza z^8Frs@x>k-@Ad#@UD=6!2THs(hGb>klK)%MBC+zJ?ieyH3k4@m$2Iqi<<OKB5=lcl zAwfcx?Z)ixX_%I^5hqSh!|1VTRuC&xu9Z}ke{DGR-4=E+kcQR?HmM}Euw5-aBJCKe zaR=5@s1SIAZwj)j&rLUlN`uhltE3S+?0w{{q~QaUKs%KctipvOR^hUYSy;L40Q&|@ zk>B!cQpS3u-I#&NlQVGMk=4lFU5L1rNyuIK0S<~11)DG_cPIw;j>VVU9wF<mf5znT zqtN2hrTE?u%lB@(L-BW5hjgs?^D`JfsyDt`noI6slV5843pLo%8ofsXFRlF8x-Hwi z7E^ND;mnRo4+bH}opW-yN&Ob!{`+m>aGZicR){;};rMZP&pOC?Y4WuqprlbNM2!y+ zn{p8zNSltfd{gWn(*bkGF2of*&%>W+Zy3?RHe3&x-UE-_)X&b0hU*SuH19kD6SIDU zf$@pBC^H=sdf8!LK9-K$DRZ!u+uYjYR6Hkk#XafSNR@ao_u;#A+Daedw9FnCkxHgu z!lJEK+A?G{?BN^?nSD9K^+#N<F_`gsCiAeo*{ku8v?_NV9(vhi2z>Jmw~Jf?zF7dR zS~Xd=j>|~x8jh(FeZ6dws%VdlKf7?tZwpY2M5MHrbtr)Lqn<%&;lnsw7Ry<MiPlii zR@yX%X6E1@5A;S+aRQPP?JBx#@=6qrJB+e~L?p+@qGUk{Vri3DxsDsVTatpwD@!r% zU=b47oo5jIsice(Y+|vpv>3eNM)ldSN5hk`$W@G|7&bOOxgw8&kF0e2lY{zAhKnDv z9XfOfXxrAENAOPv4<1B&+Ei4CfNa^CnF-u}J8;Gs(G7_iqS9w0<j;YUMiW*+a{Ck{ zd*v6Km<-P*G7jHkJ?#~-u_OnZTPk8C$`8ja@g$7uLh@<v5k{CpB-WR0BA$)2uQco^ z=Kt^c)p3g!EpREf4IMmqFiEKPI3{jg&qt3Q#fKk$i1*)rADueUpfAHPD)!zO{n-1w z_+mxFkcJ`=as-q#LS6+`73ep48)R7{r^xKtv$1X4HmgB@_0?A~V88&2KvGjv8>QmO zlP4?E-LYebmF{h~-G;`d+bE@+6!~vlsEcHZ2Oo5jR-w~HNJE!M4;}s<IV)*+fT5Sw z5!l49kBf_o(Wg%z^y<}%i;Wgy;lhRd{r-S1UAo}ur=PAKsTP;<<Hw_{tPCkBDY)Q* z3vkw1XQ5fMW?bESFOD5MhCzb{)sjMldvD(ktXl{C>Q|AmWosa+OKcjt?Q0;AY;+n* z8jXg~*$Lcp&pnX0`i4r69z86wF!5z5>9fy1<22zfET@8k0_42UiC+BPtVYRm`|`Xh zcNJ~AtttrP#EBCXQ8#SZfEQkP!RmQT7UOTE%i6WTGtU5@d}7-zZoL%{<G1YbEo)1u z*a%2Z<EyU}4F@2gq~U-{2&M`zz4TI}lNVwTdxP{kI&|oO#~**(YVZ>h60Ee%y;WFM zQ5QBUz3G<jkdg*zP`bMj*dQP!-7T;|x<NodQaU9hg-tg|cSv_PXMx}MJ?HM6i~r)k z)IHZ+W6m+=T;qL5KtTf*`lcz*t3Bo?8r>qa__vl=Qu1m!5VOQt{^L>edYpU-xLgp3 z$iQSvYz!m$9-}r^oE$F(z_oR!qX@@2We*E+$%*|i<#2H?iV{N*O^K5_;=tG?&Dd-i zzrMcK$`NteP;lOx=Ais5tDpdVq8E3shW#E)W4m#cg??jZu)Wj3>ONR#>%fG2_cEg$ zEZB7Zvcd#`vKa%P1coO_a6?vQ%w9p0Pd#>}BP7S|aPclb(1->}c4*mIjVerkM#dN5 zy59D6f_ooMGt(!!GBaa>>U!)~`x51dGwOdRcHl<F7&2Q1;rVV_uuARZ$8F5ds)LHc zbTYfN9V+cUm{+zIXt$W+9ayf2tg9#WeeG%b-V+kr?d}|th6QSsK4u=f#cG7e(ZN<F z+V&^;v&<3>Qt)e0<5V&Gi8jVC<`OZ(;)C#xi}Iy};~z3=9W6u3H&#qbU$f*i7BomG z2&;R@kQ?M>_e#DbJDcdolDvA420^oQOeD?xPL!=S2C~%L=A@<yHsUJBhG{X-<4k(g zKJ9VM41OU|p9epwqwTuJ22T_8n@*!5yjz5ve(zh8VS5M2NKc%JM#}zM54YH2marrQ zCYU>4@zi<3`mI>V57J{&21~wdVvR4v*!2l59Uc*mnR(q6UqgBbE|B{dET{-m;YQ-6 zQ81DFbxV_aYL;;>U4AIvu_6>}&B%WXXXQMj4jTQ*OnDLYLXi*T!U31sI+TUNkdCgX z5K=%KmfmG80AJr`$#eYhU>6_=wJQ!@M~QS!>oMhs{w7q<8A6lmWO4AO<#&d8znxLZ zd0xFnxFlti$_Q94WrZ=M5Z`u(!BRp$Qd-U<?+q1ohr*aw44AP})tz;{a}K3<IWa~G z1G69oqI%GD`hs7>N=gh;z$dIvR|DEQ()f5hHNX4Wf135C5i@F8EQ&)&A_xuYN<~e9 zaf53rz;H$&wm=r`f5A<vDmgoDFTwKzhliRHVmlj2)uAy~y&N*0GAJIj)bqX6#q73` z?LH&!?&e3xIJQa}gvLrnwgY=2OYTpoKNjSSK`2=K@ntOO%GZ1t`47bVA1;4hTIS9D zo_o1UrgZo%+<z#;QDR~1D1oNvZQ>7X7$b<iINl+bxllBb5VWi^G^EZZFL3)S-(bc{ zY4PH%wAf4ggy>t1cVYf<QZg%310)gtm+a#m2MP4YLQ&lU<nLY%pRxlb`%PNN%JA#r zJ)#e>7tF*JpPFLN`@+t3$oPd&@iX*xk2KsT^_?wWyQ0}`Z@FU8G$<4(Zm9OGQIT`X zTQOqNh<i{KAGva2cw#Wa*N32L53Il8H`2iAoU;qG-U^gwbLuj<tzc5BXD$|wv@icH zGie0*_5oYR&I^(-bM>XwAe{CjW`&U>LB5<iSy}{h+r=FhawbFji6IFKT+4JWX1w#K zD2PfmFS*hJlU<fcEsqg<BtuF)yVQdR9rk3kkS$X&(51pD7T+)V)XC3CTy>Sc7{*ET z<tIZ58%1uPH?XtgBhqhyjM(&W{}M#+`6&fZUEez<Bfouvdm0rkx=2uo#_Kq0EDlf) zu|cUdp?CrUMFU->+&EQM?x?elwGxo=ID!TLAw>Le!k>;z7hUH<EqSm+BNYvXE|NY_ zk4p_6F94CM%@*QtI(Wg<X3`iVu7$i&7UzDmjYwA*OC05|=G;*^p7Jyg9&lT}i&znT zY-nbxnCAF^a-i~6zI5ylh{2MZXEikwNGt^S5P&XwQX7VHFrM$A!Bb*zz|lJ~S(>ID z`Jb9k8Ml<VzyrV)^MCTk&!I20!w2UuNdH2l&*P{+xp%rpl;K<hFOo1OOD5F`w)8VX zcu3GJqEz|pa;Vf%J?E<JFhTMLR$3HMZQNgZ=XxE_r)vS5$ZifC%!z`fG8XDD$I5)} z$t4sU4P4FGDx9K{3~YfGU;6vxi{@{?NN%d_%?nfXV^ADZ(5{QHeBHxbQMWpZ=gY{t ziZBW{-92si#`K0)K)Lp9NeHd|;U{cXzT#(rg0iFJPmb=a@3E^Hx&$~E-&<Uw$ept9 z=MjFm^Pz`mUNr~&weHpo^5KmZwTgl(@$@2u+`G_o8h&{@-|rCW!3cMFukty&Sw7~d zb}v8xLAKv7WZ*9c;r&5UnA*#IU1&`gxVY7|kov<%R*VDPsL_L}ZWiK5I-{35VF}*V zr1u2q9W|QJ)lNR>vFB**8+P@+M-$jXu8O>v()?hsU&MnoV#nN7l&=?mJu_Twt0+c1 zA?@YqTNxv+%!m8RyEHk<MSVeHj`G*Y?=fxBFPo#l{b@06r|iqZh1PW?o|7GR+#!Lv z5(&cdNfmw(WUrz&B@d5}d#VX8=RSF&2}S2DIsOVadVH8ikDa6H8c~e;YA_|6%aR%X zVi93a$ZElac|kFSaHFBK47GoJteecQJ;&Ako&x{L%shA~=)t)ra-%XfciSR+_0+iG zqOB;H-cSmsYYy*DDO-n!x%$ggvgs7wSYK27&6cuhK2CyiH;S-rl8oO6Cu-hKAA~0< z(!<r&oa~n;Ne7m@rjyIvF(k>PO-|biStc|9$izIPOcb^^TwU4yekfYxD>N?`DF=xF z1^tVSX_*aXG3l<YVxO;0Acpk;yr&t><D~ZoH#TUHpKDhtMW2eEA#_21%Jc#AKIU|s z^aFi)ckZC+(K*q~G>0Qk9kNlx7yX9T;$Sdq73KjnKIZ~FRT12?==JIoox_EQ8_++) zeRf+n#qIbfVK3gM4?T)%-)V)jHgDGjp~*h6HOxk7!F*$&z60#F-7StWe6R0{)bfak zF$QjHlb~z>Ce-x+oF=kxG1tu*HMn6v@{ad})O^<$f8-On7bqgW(=*eV|L8ARpKF+$ zq4{`$BaSKBAivW6RUI^MO)8%&IYo4sAUy8f#yM+80xW`c3ev(RN4d}SpAcVhf$?uQ zd-PvJMZOjXyieO6W@l0j6s!KhSSXI7ccRFfk9T_Mcg}5)F#odoMyYLQJRnbyH|s}3 z>3W2`)>;(T?P}F21yn`A3?nVF|JPZ^`R0)=rr<G|-}s~obbH4A(~($>s0-!B{E+~e z1<hsmhcnS;q;fCgqix3eiY2(FjPUcNqwc5$ol1ffK4aQAkEvYhTetALHgTEu(-`4q zMnMrUoK)5AIo;HEQ*UVlH_QBq9=5rL{u89JisO7~_KFTMgE^iBl-oavIU{xb_CnJZ z434NLDkkhrCX7;*cfwv4#!ltjB}>f+odp-&lHOnJlDR=X@k3qvo4GW}5>$u7l3tH* zXJtErM11y}rfaoEN4J=7!c$JB-69`bAMjsB{~fMJlTm^(FRd6i=s@&sHW;qHSGzfO ze9fN2d~_`6T<{v)Nh5}`7LhsF^?q7vXEqwzurS4naKkxdw&yTlO3Y@v(jO$9-a6wF za9eJwD+r5eKAUk^#k%MlJ9FZDVxP4c+?GZnOub{4OS62PBH^h^rrI**=(XD?Yb{_L zS0H5*v|ccD$?q&MQ6tliVM!Q6Y8T8RDZgvudnrD5@#E@SckM2mrIZsU{z}|SE%vBo zr`-ToZDVrl>%;{?dt8wU;;u*^&smt%lSdZwja9Ieg3xK%^ss)PFh$`b*U~tFX|%ko z#O(LPG(5TfFB<QRYJLjuA(meW#zfnq;i5sHr!`b!E%1!NkuBd0t#%gL3bwD7Z;7^2 zyX)?tEC%W8O6&m^Ej71zvU+{E?BL;Lx~ehZ)?=C2uRSCk4pmcfH-jG8V{B%GM`(6= zNbeMkMai2sirZM2w5on`mb4a!Swz6x;AC2vOd-S-O(KLvP2afhq>dEn*|>xn0N^L} z*QbTtSRW@2sc*#LoV=frE<Q;de95@$^|P}U7q<ff{-|!^y4&VPK5Xb=`hC9BQ8qXO zY=TuQ#v~fQ=|4uRTAFeQUz*(=V|BO27ns^J_^uJC6x2o3whZT)aY8sTkvPfUa)xM` zq~U{iqh5y!(j+Pzy*-zK$TNBC>2>n%w>Yq}%=Y#>6*c>H*EXVky+dG`zLP#%>=zlL zo6@RnT_7?Gl%<-s;2oj)+AVE#vLxnae!HJBG+aMXShLOFczwt@r@TiPu9Ht<gHpmR zR=t|(9VISJvuE$n?oCbSs+Y6hX!Y1yd0e>aJ70!cKON4va!M#)>t}%)QZn>iIGAec zZ)E}4xl{{VGj1%c=UtZW2mJ2qJXwW?m!^dqC4+XlT{Ek2hsJwz+F+>m(xm)#lFgo@ zn<6w$7)e#D3Q7LWeD4*F1B5E5%H`?`>Go%fT_&;&b6WKF(z?lzzcmv{m*-!?yEfsW z7%f|YD&Wu&X+!qsGCHypX<|-I>Ap<P3z5}@zFmwhut|v)0<x7`b%S*?%cS*^VerKr zcm3wW1GV1u67m}}O&)$5BA?cY*zH_NjR7&Ao)D4dkgoCiR>vAqzen!-xe2R1&gS)b zzx!9$nzjQp*6z-SpPSb7<8UeXO`^kMN(XjdLz`uD$q)BLmiO%z%h$JiUX3|cMxs!1 z32cW!eMv#kc)f@J{-+0{k5dSf6|Ad7pXho&!uM`@?cAu3h76|>?)1xOIMJu-ct*j( zU%~(2W`uR{Ilk*LtKZ(|*gDD<+*}WycEcy}AKTwDDCW#K9BoKx1GDnJO8BnBdjfV| z)07RzyWp>~Bz&a7=+*Zc>r5G<$`Zs>f~Fkco4(azd37%@atdQ3w5e{;Qb>?*SH+Vi zz)rBH99=N@_8^|guxKsSi$>A_Cyq%_r?Z-WfM-sL)`8shy=`lrQLI3sUElUrg#(uX zw8<iqgySsd`l_Y>8=(%zfT^_g*ki{G*6@rCNvmSljNr%OWo7Q*zbThrLyQEL2xZ`V z_cT8l_GExaoZM7nJGs@AkhJ2(q+)`4sANd23Z^-uEW%&+@<ha}M<s|v^Wv<k*DkE1 z8OZ2Qw<~{@DwnG>YoKNN{3rR>(9xLYyG~0dHDs-IL(A78S)bw?=L1|W7yUnHQ$suH zvVE96-sD5iJ2JTpE`GX*V(oHHYR<oM?e&+>mm6b0!Zvz<WW<NNB)dw(6wK1o*Ggx5 z5i+-Q3{Jb+amXokX`L}C9ZBLXbZ3{>kVGZw7`Sdz<6*j^giP1c)yhd_H&PX`G?scV zg+HK8NhTs}yKT$$fA~CBIxW2!OoO?SG2UmtEWx4{H(mcL?{InP`&Cl|>5a+StM;!S zy@$KnACh(IzTb9R7{A>0RD-dP2)3L1`XVLsm5a{@I(v$x9o@r^M!Mt%J)7UJFv!oa z1eKZuhmLef$Gk=`*w#&;(}Iqe#~#M%IFL8nofkO^y=&L~I`H!+c;n{lwLsCv>*^O7 z-^<2-Fyk<-r3sypj)GDOII~f=d)j%FCBh0gsxsulY~7p9#*gBaPTisL5^@+Nza!l8 zO+};`wT$dqZKt@3=`qQdiiC-UR5XUL@vOYp;X_u|zDA{y2%%zAyox$;xgr&M<TvWX zeYu*%Y6JR4xo^$TV6k>sXf5P&sj5Wl854{=yxE1r5ppb<GpTFbi?&8X#2lIG{L^Xl z5jzyXPdZvwu&iViSh`=nSxvQgpwiZme=!Tx2G?))Xzg=!a7h@vQaDWfl3{y|Xb=wV zq|cJRR9H7~j$U%DOrmWc77Q33Xi~_s-oFYmTG3d`7Za1sVN}6zC$A*rW6v?kP*m*g zt|Nh1NI8U=*T|gW%{xr23ncJnh=Z#KymJMDqGKp4_j-lY2)wCtoAe_5X|k#dU@;2m ze>@driLFy(t!1aIF4*(F((7RnDFA5-dYt6CJSEeWM4(PIMB^IlWSv2gy{|^C_C!2I zhxlKkZcjyHhHV^{Z}46F`bte;*E>t7RDpIS>RN`&QCZP<so^C(Iy>K#tlW6%ReSQZ zEVNP-)LI6j0tuY)I`xRe(bLm9yB`;u_)~eT&B-{-2U(Y}{%kP~4=%Zm&~J`7WHDJG zKQTR~9;ts$*JTNkf3N%!rwzp^TlE0QnO#25GX@k?REPB*Y@3_E)dr6g9~YdcTqr2; z?z@aE<E5fohT2~-1lx*bM>`Ca$7iG=;!!cE;t?!WKz0u;#}iL0Tpkf?WMw<|?}q3Y z8Ieo%>uRImSviR)D5C!pi=%ED#fshRDh#6fJGxQ``BT*Ol+r3mDTfUpld-TVljpKe z0NEvQ6Ms(qzv;&!7M<xBcy6T_Z1-E+<wL(sEobR^f-@kMp)nxobzpiS{ZsWvn7<t$ zy7wB6r#35cG>55w;yg)e|7#~7_hmLe7~fM8W3o$&gUQ!t0!4^VMSCs)n~|Vi?ET=) zTJ-oPjio*$RFZgJ6yC@WJ%VIYqZDXuCWakq56qda_nWQh#V1;|_ks>yhr-6g89Tl^ z&$WL*^jc*M#My?fWax?aAmQpZ^Rqv*{MYkNVsQhHenzyV!97@9@Pn-qKh(9jzSO38 zu1@?M1RSCme_Q9O)~nz{VTy)V9A(Wfw`PBi>VAAGDa-=>Apqd?4fUkP$}7bx^hHiQ zI&~Yfabv0o5g&g#`BXX7LFs28%mms>m}9c%+mh;8{hJ(cycXbh2gkq!i-?K%sd3N^ z6Z*EO7zo!>V>JL^hp+h)TPtad_`kwufT8SA$M^xS`NRg`1=zm@Ahm&&+_}b+U9`bG zjnBf-KzJ6o9RU!>zL0T3{b?-Bf5=T0d|;#)Hi!UNYI|P3hVt`RQh)zA2prQ}&FAWB zN&^7;f|66!gr6&M?Oy;%4FCUM`reD_NhCB<*d;$q?<1I_4$QW!w{0^PNUXMsvj6l$ z-SqjVyCJCQ1jYXG-hZzSIED6ieD=5dZL6JPyOV|i=URZ!O@OM7q$Cb$2m&HW_VYH( zHpOGBIB)MZ(!{)*n(><VvmXHg+c06R#*g441f(fpqlFQBMCcZ&o1;otj1LX%D4Ml_ zrLzHRZF(^}2a+<L<6T~tT|!9s&UCB|9Hszb;1xfx$<HGg;FSai{T%xG1ER?8Z;OGh zlOP>nJGTOX^22#YR3Pu}{G=~NQw}92Nozs*J+O&cD8RO8E{~*zN&y-3u%72da-@~% zs5k^jEc#j#!r1{M?ToPK(!f2rGHp3vFD_;-bAP+^l$HaB1S{<MV(j+nOnteJxucZ3 zsGJ%L2mUHGFB>F*4a+ISW1Zw9p<4O#N0?X^7C;qc#r$!8R$C}6u$1hqxSyfV76&AM z(?yC=6>k44YDEA9mUv5Njs2`P5FQ|3z&o|^f0xEHo$*<2lA7tsTCo^_iC}|s^K+i> zUp#>M#zsgik$g5JN)=!tAk$j$pS>~wcRdXXlWzFGyE^O1IPl_9qkt`9fK>%7M+&!w z2J6KO;^wEm0J&O@zge@EKi<pNSqN4X{CGHSr<X#G>E<MNu{+nJXr>`gyYc@H&cO#n zmnpp_@Pr$7kes^dW1d%EO&s1KVfJpPxF?h`BNDZI!6C)XOQ&;rz%A5h>vr`g&x#la z5c@N-7xI$Tee=l8>FK6nS_YjVk#G}cig10cPT~W-XUKOk`2#j=4Y0YIQl?s<5oLv? zk4S|broCa+@Vt)zFC(uA&yPk^58XuRvXZ9UzUIduxk{%?Zx%sss2|q3D*V&elAS#8 z^QAw)POwr1klY>J-KF*!T$B{=)P;g2?;6;MkyvJ}LSi5Lx%?gE-G58DZLm|C(z7d^ zn20ZB`XU^3@!cX*dKExlgq;u;h*Apv={**5VK+I!7nbHKy_|xZccuPVwc)WN_8tLn zKOYw}`bxBfcURsJPkWABa0w6#p_`-s&S#+@^d+bm7k1vHwPA`q@^G;qc8oB@3nbhn zd}eZ1<XUv29>O8u*&=`v;=pidGOXE0AMgfhaoFgOMBwyMdcE^JQZCf~L+5Nq0CW!g z-~Ger4qC+CdvC0<seuANKbFU?dzR71=ASv(r?WWRS$-SUEIWiPPd4!Dpl=%F+LBIk zHvor03Oi{It1dN-?nl)4uuu*=J8DPa1*@wbAdA+R*73Gw{_;rLD}6=C=pjr7z9IO} zFTAS=^t;uDwq|~XXrz~1hnssiWY?dsxzI4+whJS;aHczJSbTnTIN~5rZ|PhQ!MN^G zUVEdn{LDf6`*=au7RuPx9@$2<({AJ7;f6!6U%VgwW7B!g&WM$rmxQn;Mwmx5g$FMU z$E5E5efuW64M<W^?V^L&S;y&>8~K3LL+N$Tfng}}<6Z=hXu|iJCQyS;^tI@j_>i(w z%-1IJc;4vWQx+dgyL%>5dNxf~wM_JU(To!RQr|8s3GrJ|xT5Qa?lrruPP)})T{dtL zMhPxvWCVK<7=_0wI4oMs*(Cmtl=%!zP0(olAtv?)jN;65sMCSRZ=P`uiW%q!#|I1~ zwtc<~2-lzWW!{>~NM5smGpa)nc`UHI<g}dj{*o&Wz6~j~F*qdrxUgDYUad+b&YA*V z>VFiviPw+eiVSMc^d3fz=z!dIqKk<li@je;H`|?cgE>;e(lxzj9*B~NarpXSATmbF z9#jff0YRU2Fbw=)b+=K9@aqBT>Wb`_y0qN@aoFJRmALBHO6^@+>zhq$xT}!Dt0$I( z1}VL)ksBBHbiBf5lhwfSmxNqLR~q^6Hmx?s_v6=PBXG*^`B4_OPK`6a>lnm@#BKP{ zQJNc6;^inQuV^CjK3c#)_7*>{GQV2QhIeg?L;P+xHdB9}?mTeMO@g5PJ`~#gNV%U2 zwxtkkNeq1m*Lym^YUjc)ls}a6E-sAZeCs~K)RlCNjjyvYI&c{+Wjo~@D~<L!|JhCX zZrfqb7K_9}8QiG$e%cBx$6D4Uv#XJmo+`Ha;5t%TjkfZ8ldvzEj%+Kxvc^rGO00!D ziM@?)$v_9kv_{eTquk>0Ce69uhBh+iMZHJm#}KVS+ww8T7Dd`KNqt<kjLy|Vl>uCz z5yrP;Re9r<;{M(q2J&E0Mpw_mLN_Opd1vpVdo<_4HF+tMm#jdx^BL(Nh&z0v9Vu^5 zOB|YMo-VyJSsweSK%=?R4LOsY#cmNWLAr2V+W&MH?Zxf<3FC*%J&44B=C<}U%<iv6 zY(zwg(CN&Fo|uSmY1XAxi=)W76Phov%r-0Rqr8TlX}4k-Gb#>2XcRu%v0J{41POnL zjpu^X7ZQC<x>FiL)lC-=Fz|rxae5-%#_GB@e%g6vw_7{_9bv%G+=y4^g-Nv;H{nkE zdcqB_gATvtb?GHx^q9U#Fz381cu+fb+BKOf{zJ9#DaEm)$St4c6yj#(kGSH+Ah@S5 zMOhv>TnhD4YUHutA^I#gP`jx&%g<VKpP}Dj#Avj+l+v;CwL$Dq-n%h;Mf(G<Ej|Km z%zpoV&Ge8y_3Jk0D2edj^&c^p1`c!K^vyAc+GYl{KaA>Mc)SVm8v8gtef0X=VfSX9 zX{W>7l<f3-(QoWpvE<>b!JrO{q~o-|Mkf*9CjJNZwSkN%C9+L&n46p5EZG%>WXI$d z5$0LjJvTv6;{P|(VCG~L`KpN4Yn~c$+f-qOox&p>_)?^4rkrnk=DHB69QTwwI<R5r zA=8<z`xol8x(S5jcbOQl5e=!S!WUs=Gmr;}QC@S-O?yv(t|R9-U5NGx4vv=xQLnWW zMN_$0dS?d?5vSR29F@W1X%BpUe{IAB+05)OsLkySzNj+MB_o=0opC2LhlVgWuQO#_ zZ*V?x^G>e6lJ8Ork<R7F4+Y;b<y&rdoLz?`B`GgIqQBx8_<oY?R<0oEg+%#R;XCZ9 z$mvmzQp08Ji(Y@2PIEi6MWyz9HReq1UJDWb$W8CC{x|3CXNs;dg>X(B4`=Rpx{`85 zTk~B!hkJ>LORNBCAgvH?;gr&U|AKy~qpbLPK&?#XTOUD;v-5O45$i0c=S2gyxM*th zVdMOlfMY}19os+DZAQg%ZP<joUrl?~{O?Kg1@c};E+(4x)UyZrx*aW(WJkR2ZY*xg z^r{1sN|^|R5E&AWgb|w#2M+0^$a|T)c@^1)TqHh>HFAtaUnefM&1gi>c+h)o(foR; zHkECX1-C@cJoY9Q&cnQ=deWgw2f*Cr)wVPBbg5`FOv3HRaF-4}N`@J8)g@A}Vd z4GMkfXv=suOqwS?_~HGm?}Juhlvl96%iAPIK>0`-p{zvqW=9Q<^8XlsjVQ?`XSx!c z<WHR<0&Junpv!bE2;3`^*f@<ZPD{UA#g9Flu;&H_vHs?1bvuC{<(Y8Di+@PB-I?T& zhx)@|4O5ub^-+9CJX@P_zwvZB-6ZcnS`V~35q{5-hIp^;{)aXW(#rnjk9t&o+y(6` zKWgZfd^Q4!k=dnvxZMOfPu3>(ShpF!u`uVWt4~<N&fY?S#}2CuMj^w$Pcy5nQ!B9Q zoD2WXoWv#yNU{QF3#?yP<2JC$uHLN#*8jZBST=nH<t@X_zwDD3u7wu7{1UDx7Jcxw zr(r7t#Y_6N>K{S5Cnswz%zR?OfS(q-wi6V#8WyRnlUA@CRm?}KTMTt#uz6`CFI^hI zBbES!`BbA(9=%EI-zcInX<}c+)!KCCG0{5~)N6E07_7VH;<SLE^$b$?BtskPM3=e5 zYcI#R{s*z0v=|6XohV*i4Zhv}f_rt8GcoOR=9})M_j0TpfU_=dTE){Y7PDwezc9h7 zhKqHxrpvgyT;cf*?o3>{KKiE5=A*OLFT+AR9lo)VS5{U|jy8EXaa&Ig2$&2A*r78l zPF@DI9!o~_3aAzRpbI*7PI5XE@`$70Omws*!;wJN<`Z_^OK~9;_QVSg4(4-Q2l3v= zh0i<eH~c3JR@L|R@nkJO&7ai*k?y)enfuX+6+H23FC6;_mVxzwYT=fy)y^-TtC97H zqRUbwtt>sH^!81bq8!sFzmFnuad9KRT&#CxmfggUkB>i&W(P~szN2g>L67aG=SybK z71ET+p+bm&@SPa<lf=e@CS>dT#mlBO%CS511n9dYE0ABcvwQC9CQD6OJoP1d^8Lv{ z+{My5=ta*cmKy8tljnU6m(}?MHzcz}fPel$dW~NE>|7_%0M7MaO#hcg&!FuF@U>Mb zYPd`O`Oi^+x0xXNLeKjdi0uUW`xgB7S3aZwvO60+pg3<`n&SS2uSx-qH+)S2ws!;% z(i4WGh4s3G)iy*<A9RWS)QC>V?!QZs2b1ZI(__H>S`?7s3P$SN6G93A*`o^>xRF)# zHv4?v^I|*UpQ<9_jkxV6tvo#U8+o!d;+oV7Kp=S103e9b(_ejAbE#m^6SO-D{tc_F z-HZVXnDN^en(+WY2%w~jiHB=th9MUnx?g{T#AvEG39JbmGg1xAz(*0_x?gXSk!dcE zwzBWDAHH{I5*!cl9l<Q1{!#V>J$OI{Uau*IQT!Z%p0?h9gvFjbO^}SD6tF$1F9II! z1HN(<U-EwF&K}L#&q4kZy1@c~o`#E`dtqXH)JK2R+pNn2<~9Qasb<D)xwrzRj4->$ zg1!yiI2cW9eq|~FHog9fzGdI{v31Z#<(yusDYSop$O!N~+UpK9oSwHjl@vfA^dSaZ z|1gjD6ECpzOMUQp?8m1hs#H?PL6v8n9Z3QZ$a`0b=bwIifw3PQdw9g3Cu)Df{YRyO z6aV7e^iPoQ>7{T1ZKI#x_b~Gjhvssl14Dlx0h@UKLOxIcDL~n9EAr3g!#-)S1}?zp zV2zkbWannnl%{L~Smt)5R<9B$t;1ptpK8d<%|$ytPQ0vjBI=X#pSC)nPkQQ*)kXoN z1ssJ${1PM%xZ8WhF6>p1w{}5CM990ZE5t8==-}WGN_0=V=)UH07X9VzyLY6JL2tN_ z(ipSpUw5(xE4U$izq+Ak-POhQsOkEy<&|qt_C>z3T!CMYADb%E_TKT=yYNemzrQmv zkuvQ~f^}mNW`rFL2CUV2!VQ{mF&)cRvG4EYx3zX#>HFZkM^^b#(9uxZ^TH;A=oE)s zCoM32;;W9R8`?w7IeXtl2{cR4!gObspU?gy1urY<=k0AwV|i?q8Xs)W$du1D!Pb@5 zQ&P@8UaVD4)vETngk}iVBjAy5ZA~`ohyBtQ>F*Q>$B3*^`lu0=Ux!kwXu~(%?z==5 ztN_dgRa^7-an9qk2(XL<vJ@sl3vmk}Q3oO{T^`U+=~&1_&@dv>`c2N}ay7P*R*+T| zD(MAH)T3LO5e(Y<Py{XFHm<{x&Up;&QKaePy5Kue{MMF70cQ!KGB0*qT#EwmOk*Sa zT~vz};e7!}0u|8VsM>g}fzqe-B-hPgTC`ZTo|aGhtHCG|t0NsPS=v+&%_fMftyawr zEBxQ2u|@P6$i?jjoeSO;KrfP_$yC%bqW|PxUzFZHPE)%n7lairGQh@00=W{BuJ=-s z+k5$QF%sw0)f6VOi)qm&eXEPpDSn_lRsuP<52NqBaUCv3xJ&#Voi1qR)$ZV09F3C5 zm%58R_$VF1<X4!4q2Tff_w_CfT)p8iu|H3^o=3A3;q}kjq)#Ciwz#UU_JY2*Tkxm} z6nq^!5f2^yd8S8WUX7~Ge3-+F&isxGG3WI$bD<^zRg-%geA+~{<=z;9QW0W9m-u)& z%fJ^PP)sKtA{<XgZ^a~!s`tg;qP6{0?_EvF)hhyrEtXeuJLvKJ>)gp53Jx+`F_4A8 zjU!*-vR|Gsvr(Q=kmJG>XAd^opNgHSL!34U&UVTE>!=i2OA-ncqo@uX$TB2c1nS9u zO}vdN6MitjCW2(r@U~ao(;$WL;#e7*GUcSNGyULi6m^Y8ae?$kvRX}eMd+0aYs8z& z<NG5E4CWb5G-k^wKKCq8<rlJ;^q%g=mBN~jv{@{6G%W|b7@cKHx>8Q01|cbVM&C{K z(+FAnb%Iv)GupqKV3_|@l*;y917k(~s&-5J@ZxCm(;dy<J@>ItysV$37;h<gvG>(& z;v4c*y|c~kM7;Nm>Ybp@xx>BfxewjPHaE1!HZO-x&{F5g8GA7pLyEqGf<OTk#2K6l ztU5fv=kR8PNuu7e78@dV(q?2l)ivdCvG`Nx>cr}7*dS8ej~w%FD5!Yr1%7vpXojtt zf1Il~CwO8!1|_aTBPCrEFQAPQD?B54A;Q8P*qNVHl(x|pf_RGS*e`RzJ?&T!4Xx9S z!O{+mNsCM$2)Kg}j$THsehu53Z*PV&k9@J0iMM4}+2J@cydFZsJ6AiC#5IoLk)bHK z-U^Lf!_~^Asvv>9fnrnGkuZgag;n@AZ%XS2l;+~yq46&vvfZ~?&D9h#^0pjn8F;Tz zm#st$G8Uw(jbb3{x|2W1lk$sz-x6mmhym5z<1rHo<*g5%-Z$F?Ta%?}nz65xGlKXH z(s@iV3OO4-{h%fNBHhhPMnWu{y#~x!to<yN05-|Eehjaa#Lxyi4GBPpcTHECmIiYu z>eE^%ud;9Je34id1!_DCLBHRZ<*h2vp@3^Xj>7ufUi>Wj8-PUa-von4K#-H03;#+f zy&~%5#IeV5hB}fea`<5e&cev4I5Fy>ILzNf9Us%Fpq`MnDDzurWDN%dq*rTWi5*x2 zRDL6C%~9`%$4I5|+LjP#>8j3DE1uYt7bWvK8jvay!&Ykqg@T<5X$x7X(E0m;G&~Kp zLL9@pAimu8VD2E6ff<x}f4{hVA$Fz-Y!#a#Q)<?y#*pb@ksefu)iulzz98;(y;%ZF z3x!6SbL&;YcJ}?m$GDJIWt|vcRpXM}P^w(YkWec0<9I_zy*#EkpjIhUPvSU3H8r)X zDR(Tntm!L6&HwI}rA5;mAr<PVW~lfqp%Vz4(t#C{p7R!vG-K~|B4cE0dR|p7XQwMw zA_H~eY>_5&|9ej&2>UwlXLi6OG$+gHo;!t+veYJPfMpbG{cl#g+?8ZrsiAyns$X(m zrl5|qce!LIv_y=_9={WJW$*zFz$~ygXv9#xa7seK$8MKVP>`^+EYxgpe3^YzW5Hcp z5Sw0->G@UQSV*3Xfj1*`dMZrbMIj;_w~`|lJWDbjkJZl3R+L@Wl!fBZ^Rqtf$dQ-r zSAKe|I|k|)d!nU`IMV=@|F&h1E+qQH$G$U8po@tcb+&NivT)1E+Nhz9i#LAqP6{Ex zD-%;<Z*kas0<b43gs;hqc=gk2AOBeAQ-R^zA+@g%zV+|%UWH=VTpJAd4DetKDgtrR z$U=o)X<OnJYz-SA76ZC+bJ^iC^+X*Ta*qqCVa~2j&)bi!b~r_2rsK6udK?R>dr}Vf z9fZ4SC8;bYG~LWKkP6kbKc;@;C`v!m%fs0DBEb%dM*>8eU?1LpyBd%o0i;rHXV0ay z5P#|T+Khcg!2oy3ueW2V8|X@bm4?D$2ZQM9kQ|1~A;lz43}ILo!J@_tkAv9J;{SJ+ zy8-vQH)99xI3e@l4_%7w;rcZ)(=iZ@Fsy(We~l!Kg!H@<iFFnz6<8=F2Ue0iF1a#O zqQ66R$`MdytlIMdRf2ZK1oMYyj2)zSjc<e=Y|M@i8^|}wBmpI{g)a)RunQ%`N#B9g z0;tW?Re=IwG~FH_N{KPN#xHU%@OWuEDWT!3t0c+xZCzOY7LBH!xDc5E-9FRyTQ-OF zs~~jCl(5+9LZJrUTq=u7Yr?7+CrZf^?u&S_d=g|!;80tu{<~`6cp>aU$iEG|Ze%?M znhm6X1>sxDn5D?YTMuVIVjxslcI`(k$@S;>9zA~tRr3U3!xUW<pps7y$5$SIx@$g^ zbXzkYdXOjjDgs_Wk<0rU7Ekemlt<RIM_j5Kjd3DDn)~MIeeeHkz~f0*Q+ht|0TuL- zH(KBva+ElT&=!*jmWqwH|GQU!o>*6lrmwoYRz!GOQw0q?G?oA3(_21yb$U%aZ~l38 z#DG^<Q6{2k`0N4Z10G<G<?-bFGY=24esU#$-Dr6Hb0uE_WAi)pot`_XQkFcqlDj+V zfB(6XPn1{cOS?hE=ZPw$09Uf4SUK^ZEBXKNQc6hpm1YUm!*o#opb%EWy-Gqfg55<M z+#8=@b=*cIYS~Q7bsw}d;o$=9KZ1brDPaD$Kq>Exp`sblgceO1nuIWJ?!&t~@3XB@ zI3u5P{9iy;_{azDjAsZECPs!evQjAEpZoTCaVqdOfW&9A&qxRGTDbs^#|0t>_Z&C? zIDgA6t|#UnaOG7jfPp(`iHOfGG=Po9fPN4E+0CIR3`E3){v0#_!ehDLho6FEz;j~X zfPp`_G>D&lYakjWE1vnU<11c(01cM@W()?hj)q3i9Dr%MqbcgRnHU)>hfKSpB)z?x zTo&>bNt!NZt%1mGgLsV=&$Gy(Xz=Y(+&H+n4p&ED@2iy%O)KN^e8ov01qHNEE!XCB z5lmZKk<YIWcuUC=QozZ{iA6xLNf6+5*gRF9>G`}IAaj*8pu*tq1(59zwJi9C$i@9v z*V|#K)lW2)2UT9=j3vzQIZ(LF2j5<0W;;=Y-xWcmURpN_bY*yI0rG&HV}qHoM-k`{ zUaO@cF?vB4`h_U+bqI6%*?*ej0>LV<EU{XW*Fgi4nVFepn9C{@dERB}z2<z;Jk+%# z@l64EJGV}C`QmCz^FSz(7tMlB8@ce!bGU~ARs;*wzrx5EV#7E6(^`9PuC6>=2M!5C zLP<%fv$;Gxu8B@Y=9~YNTjJ_QMq><ie3*LG>2ag#f;WQpKdZ)r!Gv$GEugB-{lFjx z$(4Yzvm&0%(zYD_Oz9AlUz!WfZbH|I_)hpn138Cx{?9XFg?US9<3*g^G_#{w7R-%7 zE(E%k%Hi$R=Oiv{6(l)@G)HcCXGm;ktAZLEx_{!U-Ty;du&`i7;?Lck1V%PK)l?9N z>#TAmM|-jxw|bNYv?lffEFw7ET|7lN#^-p{J}7AaD4~##i&nGTj%yaP>#S0N$`%8{ z3?^hae}Z&jVx^!UA7$~D8NSE6(_^5@<R;TfpBat#zXrz&3KF)Ip>lF^icLz=N<Oh7 zd_J2%VJRgRcJ65CiM7Um{=<_J>w-7ZGeI-(KPJtM6Yz+MV&I-H`q^^<=z9G^R#Hi# J?49xF{|8;`kZ%A0 literal 0 HcmV?d00001 diff --git a/doc/images/tesselationGrid.png b/doc/images/tesselationGrid.png new file mode 100644 index 0000000000000000000000000000000000000000..6ce0ee485b9b308bbcc836827ea9de3f71eb05f8 GIT binary patch literal 18896 zcmeIaQ*>p)w=W#qwr!(ht7F?8chs?M+v(Ws*jC3*M;+Vd-AVt?IQRV<-+j3EjIs7Z z?NzI4)?AanIrk1zkduIi!GQq+0)m&46jcHO0*(TF6rdmg--tkx5CH*UN?C}AC`gHj z5Ggp=ezUMP1p;yfp4T+fK$EFITmBvI7)-vO={Uh1&`(ZkF|rH<1rr`>qJS1#l7}Lu zf(8yw7H%Yrgp4FADufivb8wmIm~+1Fb3d82zntw-URLf=>|WNk4({ueCP@-RWDW!c z9})eFeTJN?4Ufn`AB1cZ;+G&YQ7FG68NQX3<3*Pq?j+E!YP^)C``fPc*GU*H{Vpm( zAR;S`(^K4CRG4iHAf#6Qf?;!@*o+!aGKO?8<2h0+R?%M6z{_JVybIjwGvjl8tf(hb zW7~2_pmcB`Jn%RyY#<<V0XoY5C$MgBs$1ahsh~%QAYj3BIAhry!8}Hwih-IaS99|R zes3US!*icU=m$&+l#cg1DKdzcD4%zXW=h@~Zv`~!z|@Q5v3FBc928K2YJE)nRi(fn zjkBR(=r_6tl-u!59>|wRDVRre)+Gp|Ygosc0Mlm;-rs*_ju<Nr2`$WFjhh)wWWyKM zhAvnGUZkcM?tDfN$svw3ds9J@_n*bEUoIX05>5XQ_@+|xB5o2{2^#rmxQdT%Wd>yI zDkgM}3VNFLR`<AfT>PGb=XZg~!uDtrkP92mu}l~N9L;TO*S!3Ow^@r6-VGNBB!C1; z@|Ocn638(Z^5)%=+*g38H(f&m6WqTC5-X&aulFziUnYtd)8czVls^4ekWKIRHB#V* z?nY*xm`4vb5Z*n)9HDN^L7yE_*GR9mx05{H0Ht^A>zIMgF8%tr?}>Ve>yNCIiA2&E zg}yUl7P)MSG`_w@!nMHAiM-O-g#k0dr9EHJ4WQQhvO(W|+i&<Bz`rkTU4W?$32L;x z|LJ%)_uFPX&h6@YAA%E{nsRu)%E1OgI>kinf&@AUN<bkA00GH|{&JCVOex0%eI18J zIV?RUd0zg$q4Ugrd~edV@-`}lG<zYyvS<jj!wuZM@Jbrt0J^!t4g{w39tX`1(#?P< z2;WOd!E)CmF!u%zf8mRx_p`zm_5#@OfqSesO%dK46x0Bw!k5nunBd|GI&jf|EEIHP zlWG`fnFtfjS4jc5Oc1{r<h!r)6!@8+9utswH<&AUQ4gL5NVx&1wFpvR5EZd|DEd|~ z7!k81hEgyaG3YQ}U62kDBNY%#fKjdq75r%MRj$$ji3*}_$amrU+@UF|8<2EAjGycW zkRV~Xf=onQ$o^2&fh@w1hQT!w&7>em$f2U5Uq~PoBTGbe#D5W=5knxMhFbKa7~w<5 z?nqLQAu55CM$U<75hD`4#>EYb4AYqq*P*zfFNd9oM3WZBjSbVIL5RhW<dSVVSw=bt zKIW1rgeXud7M6Tb!94=D#A6T85_=U@&f}gkGFSZty@HMJ8!cj%3p>qzpmrm21Lls< z;n!o3$|Q_M9?3ii?ue8S$Q{NVV$u7<Se&^&)rk&cEF8Vh(x9?N)S8V0RNKEJvLgb& zZ^Dq%euwLOCfx**3lL_o%mDfpj-890akJYh2R3{fP#R4ddK#l9$g)M9Wu07|VI6m! z4KKYORJC8z5Me*z0Ooe}7U7ouc2|5t!gYLk?4MZ3I6`7ZiZ1cDP|sk`2>)o|MB$j; zD3V~3ki6u)j4kj1b!lm7A!(aY(8TFP?!?<fzC`8343!&|cohnjPL+_-(UO&t!4kES zkrMC{rxLyr)Kb{ejM9+O6pLYVFpC@WP__{BQS<LJyv3Tug~drf8}i@_`K0~`KBqeS zA;?flnhEa9q{yj@EsAi;#1!o3FXrF>l+Uj&@{tr6@-;*<Vlw0!jvcoPG7L3Lii=jv zo{nSAu8q9p*k(hGmx>(;qz*h}n}shj`znyio#vGKnkvtvZ%SbNHqkqgH@-C#KQPnJ z(sLD%)NR>~^V>bRA?p?%H(WMxJeWTEl=qKFqp355GrKe49o-%72_rWWULIZ}o)n%c zUI5;4rcS1KCfNka1i^$U0|kRkRe;W;j!QE}%V7(9OF;8=%ejLO857D6DB?h)UeI0| zM93IsdALcjY~pMiJz_n`b(M9Bbxu!Lep7y*PS4I~&+?bm7mSyYm#7y;@L{k-a1iit zuqUuiSQzLp(BEJhVHu!mq0*qpVc20hFpAI?P@<8l(7b8<krgqZF}NvNNWzg+k>Jtj zQSmWELaD>S!?I9u(Iq0vLqo#U!v#`uW2fRmlgtxS<4oe&Qobfxr|86OrF<cB5-*cK zD*omFBlJhZ4K4N;Y~?hSG`Tc(O+`(@WtQb)heU^9hmhSbJDtN3BhMo_DLE>+>NN_( zD$Uh3<s}7JrLL;>iedUAGNn3)Ei3_b4OddT5;W%Mx=DRBx-{Bl;bqfhoGJq4QYGN! z>?QOiu%*wX>gAJ`7nVru43=N5P}n-y7|rFZF3ph}Xd2cnt}UJ|7tJS+Ll5xgS}N-E zXe)6`|CBxJ%BwPpRcjNL+7{##3YZ9*Z5k+8GMNrrVHqr$H(M4MB$~|GL~?Mj@*BUb z2F~^m6|9@gtBkM>-1uP!rDbWBF7r4lUQ1rf!3GAE4roR!CW9m!CZDQxkM<u`%;Q=S zA445$b70|GWw2x*HnKF1G<rGTIA=S<JLBH9-W8wfb3^g;bFI1*IR9|0aCtc^I8WSX zJEhyr+~VEr-9y=#u?O<#z7!nyb5aLg1X<+f1aE?CGCLDKyF82GTjCYw{qCXUfx?&3 zF{LM`2d~E?N^3-JzxJoUZPjC5VzPMHY()JJ=Z4F)d>`Ya?Ye%Sz@!hcAnB%cw^6Q| zf#Ie@te%3!R7$`gBrSqBVk1ly{5q^Xj0Ma<kbBTNTs|x;f(OzBIsqZSGzzI{Y!Rsn zDGZ4k={Sj1UM|&rNwhh{gB=$ZDV`(UHwrb@MmlHK3Fa#9Yz}@pZzdtOzoNaKd4^ks zWGeCUGp{DuJNYB5)L*DE5rsn95!VEMyKuc;tRm-=I!Rng1jO)5CCQD;ye97@E~HJQ zd6RkFrB2rl!j(l1NA9OtrlAih4h2LbX2tR-I&7IG+p76A{%|(Dn~=&$cr7-OL(f?m zyPU-SQtA5S$w0wLcfX8X$?D|x%Pq{U-i&^taToQf{;CVz6@!-{k+FiR7jpnj7Bd6W zCGjF<C;221B3Ul+Fqu8cF)@eg-FxDGb;<4YTAh}QwvN_Y2db6HZTi4SyU@VY5WR?| z@kpe>*lB}&SLT-8nthqQt--Ki)~V46@hIu&7?%r2j*|q(7zdNFO_!i0&!=#^tbaa2 z_N(mjVs}HT6=Ai9u8(uy$>8g^-Mw;p5c=K9n@UX^1#RI!mL3MrKSE@N4rFID?s?|F zXxmDymz{1M(zK1cKHuQKWZm$1xS*a_xg@lvoP)1sbey}rtHwA^&O&k^1{0+7dCpla z{h2MEb-9ll%&!_Mi?O|<y>CLKLihzw1J}r}=sNSIfNPGUf^C&Kk7b6tjNjettJTSR z^PB6razZgEYBFm6Bwn^Gzl#Sy-yu)S1<B4pN=hwblRmy%ly}x`@Sf~nv>}yU?YJh{ z+3?vWvl}yrjpYsx_G8wz8~&r$37sd%QN*C@)$3c>I#}I^;0RJeU4aSuxa7tZR{m*k z8hvCvi^qyrH@kT^r!$$Il#Yf3ea)@~?}T?}pEY+cOVS_4&mIdD$KTeHuhcdWiwUp) z*z+u9ZRy6XP2_wN&_544JB<sC>(@5}CR-(Wc?H%e53+*7Y<9*12HJRSD|#>mPN7Az z1~q2B0HP6WCc=(8S6rNyV4BcC$HT==w}$N!d`sS<-bF1VJ}L&-hT=p1X<MdvR(qDQ zy3`u?nt|G!UGqJ>T`kI*i;`>i3)2gIj8vpmbPS|g3|W*mGJcv`T5^(8Do2&O@YMnk zB}JtXsmBuO{3GQCe0XnJ)=H>Tk}QU`rYA{IAra$zqHXI0h*|G?G1$z)+$6t-CfRVC zQr))Kl0HgG7BfA}C{Vd7vL@?$%NT}9=xTJ^NW_$C|D2XCTS24EwA_^d-5z7&>z_P- z*~s3Y3MK2QV$7AkYFT^hIC)z2_VpJ_oyT_-oZuI6^|~htsD6>9U%N<k3YIClzuZJs z)Ers1psbhQ-9^i(8m{D@jdg+dhWw8H1CyT)hl+p!n=y;dLz&;}?Da7yybzWabCG7T zyi@(X`Iqu_J;Iy_{IK|(=^Xlv>JQW>oTo2G5OvtoAX;}?E`J*TaPn|`uX6J4Y?x#B zd0xxnbfIh;I)h*JnVP-L-8D2mAL#6rSdKxwm4I_aT;ZF>k0;o}>51{g&B8gx@48~M zCAmu*b8L{h8zE$*&~I%2^7`c<Y69(e?�@;Y?bM%S!!Cz2o_xJr)73ao9;)`Wwb% zFF4&9ts?F0_qhg#JI9N%znJUnH5{+EZ%r$=nNiK>>ltd=S^_4Gk1k8jr9MXYxg8xo z133tfP)RWey$41?b|%fXCe;_9J55ySRJbLUCFSOp7PyBi2Od}u*c(ihHG<)>|j; zcxRbXI)$sJfl}G$o|!LNAbS33!YM%#VjTU={kLTP<df1)(v!nRJN>(-Xj|xGNoUlP zWk&OXbF%FIY%x~zmiI@$dxKCy{OW#2DTZD$CtMfVd=uFtq;sx3s1&pLVV!mAe;V@z zjYrS9b^m46bQ!&wxi!Mw+<nCZ97VZ(@<HdS<M6)U=gX!x(jAd@jcBNIiA}1lD;+gy z9@0e1>9iN~q|wGuCAqKYC|RS)PqY@1^bpzTCBj=t_^%O?^1R8}_B5>Tu8(0y6|~Fr zl(tz%BZmTf={#<GHtQkvv~~!q48j=OG6&L{V^F)o`(%@}YUVmj&E|R*`ny`EUWW#< zCN{!r(@)#*El~q0v{iF@9V>IrrxUm9xFX||T18p{zbb#({&+d*U(Gu*Y8`iyxTq3s z>gBtAe&Q=av?lm++j{RN0ON``BaahG5a|)ml*pFJ8_%20mdr$O&*tubL-)|KQLqnN zoT63@r(fO?m%ZgP#w*Lw!B-EqpZIo<@0h}$nI^9XrSa>nq`ki(QcLfFt!Czd+lqZN zVgl#S3Ju<UyVE3mT-5yec!r|(jJAvb!*`x;?tG`kXUh4>Xuf;>%Y85lB;Z#6DX~aW zPz5d(@)e+tK<a*+Udkmc7Ij`!@5mRnAaZj2T`X1L23$`io&3{09=Y>=qTe+=OOrpc zx%~%xlM<YwDswu*Nz-QiWc)KCJ(BG*ibFZ$C;jymfKXh?mq-#x@rd$BM`>DV7AU>t zNTg>(!(}X0TP2DVI3(p%vvuB*g!7H^;=@?VuiB-{Ez8myzN9?Wi5Gt>O3vr#x=Wwn zThrlV^s)doV6xOOMYFQD?)`p+&99o68lsJ*<)Bfe%~xAh%WZUQS!^J0+;ws>d>OA? z#yTx_47Bfuzo`DbX?2V@o^P$4EmWtU?;JaaV8gNt?!6mmRxpg{D7PU1wFjTL@MlPG zd+)8JF!DpFd>pCb__P=^bnThF)}<$OYG6#5e0&u-W<fup15bsQDA0O{K(SoGnsha< zGYKil>~8ZJ^d;S%wdmdG$YSx6;}^NbU5uQOQ5kNT(R+G3>R${jGz&C})Je1nbeMdU zKlE~iH~qWCdZgoM!YtEdM%lF*#GN>9zrzbf*{85pq3PRLw|dxYJxQ}x1$C)5jTiIf z=U?u*?oTpY{yB5GIzw3PTjxCEzBQ;*Z)V(`j_{81Yo)Xn2_m)5zuj(?-_JppUFq0y zZMgfkWwY`1{N!PE6l7MSepi{6v1{0Kd;=(>M_~Q*>3ry&v*}kCr2!f+m46EoWh!oM zIT%qrT2Yd-fp;kwkd-vhi7HTUAs7+w7lW07_w!d!;Kj@xCjCAHu)kW$xK+%2LST}R z<hjsQgvH$yI@Hoa7z#`}@Or_)x$I!9!k~X6xS-XYXC+|@U|5he13d#@fYmM6ON{22 za3U6jDD!zXS7d-u11yGKx-AW340AS9?L=46o84K+X-Tzvj(gE(Xj-7%*uh2ricQOo z$fOCid^MGD%Fs$}^@R{am1UJbmgbQ*lfL)Wi+JM*fat>)#~<XUkSW^#z5bh)j=Dmn z0&@v%31Z1>*~cN#zRa$3Q+kiW?jEJrZOPT$&G-TAQSV|O3Jc~0dL2a?Ney`&xrNeQ z=eJs>%u7N){s^lQ@QfCVHA_Zaf<x97I!szPvNXXc28~vkgJs^4C(dfRj+V|!`<LF< z8F$O4L5PW9$S4GC%v82Cgzh%Hb)}in^)CZO6L)ERgF*w$tECH%Yry{HDcVWTdBqj# zF!tRceJMMOn1WDW5!~@Bu$9oo5zHEwC7(+qn6+v-x06s-66WI13G29Cp3F6HXkI;o z-S_6gZmdfaY#k@|w07ch@ENEWV5k%sjcMM!y^negOkHj()2muT;N=+^+5c9_u<<l~ zYq<YG()RUi$7S-iM1bx{``m@#1ZNL-97o^7*`V$j{XHYHiE#a}lkxe&RMD3DeD;rO z9Q#aDV^k%5<>{}rkgltt^M~2Kuh~ekpkKOizIo;eQc#6x26{8WIvZvS0@nrP<<=|E zvO%8`xgt~8Md^pT`f7+dkc`Jl4U6yKqS226B?Y$nX$iyTc1-P0RWn&*sW&fV_rMtT z)!1<PvJNF}VhpPc*Nqfi!aOKHxH;%v!Cd8_Wux71w%(KNeo3O72Wo(E;)CZIx6nef zVL64+1zUshqIxw8mgG_|WHL;3NaIe$)tuAFu0hf4s^6`Ft-rTJ;s!P&PXS&;aUl?j zVwzO$E}v}cS?oS1?<QelNToF)|4PcsphO~#hKl(I#fMCh2AiUgah$q?{930`w@T%$ zM%XtmDn7Wb8774$3r7yW#Ma8%;=ASR)LUu9LjH_dkUUqo_~77oint=pZz;dewSQ|6 zR|Z+rSsS0woIv6?W|3uf=)fStwZygHx)ylcp^UUwdWL&$Jjy^6fFa=Tz(>PkBcQ@h z5{--Dh_m-k^$!k>Q)o~+Q;rqR;qUP)#gAsHxkj?*Z#7?pC0lnSpmO)Ll(1q>9_>t` zXVa)rYtqpQnCS~S8Wz;$uH=hXI5)jn7|$ln$|^3dWsGaF|3Uvgx*p$F>ArZ8tt`LP zchXrj#0YUI;!b*KEQ0znIjdjM$<KA`v^VFmi-;TXIyRn@u}Dzuvk|{0^@{F8Mg)4p zq)!R&7&+k!<|yTqi%8cC7stDUfFEMa57~A3EmTF6^7qDYkF_kU1L__stKzmfNP6H| zPmO_ODsTnU9hWn-SwI`KBk7|n96$DZRHkUBG_qK&Sem$)h*O+aoO>umC<bJXsP%#O z)FKlXOFZ+sX}DRknVzwQVU*dd(NL|xZpl{U9?RC=jvi&?h0Y$|2<E6QX)fI`sVMO; z@)yOOhJ%tN<4c0Dws@Dq)<VKjxKaKD56X8j<kJ2;x=N!moZRHn^<3tPmeTQJ<MNE6 z#i}kNDx)1+9pkI%uPeOMd5huGU%K}i@dLb?`%t;Q7wdMmw6^fw#9meqDR~fHKJH=e zymL%(o3kzP?%1a)y4nUBn&|(s_^pts`w}5qT@aKQD)L8yMfq@+0WV>+*>Ux%Uxnik zjii`J{o(7PX2uI*kkkMhIlJl6ZvH+fjWgq(DzvK0yvAZ@b&O5FE<+2=%$)t&#h^xn z_VM%GX7pqha#QCjLR<aRshe%P#Cy-J-d$&Y%=)brVj_Y!!3_2Uwj0MF|Jl<oSEB=N zBc;8d=doYv4b{)<zvWK9&Q4r(CZwtQ*fN>j^%$)?^jJ?H<ln_#)%8>M!mODcT&K)) z>JwD6`|#b<d-1vUJ+<{x>|<D3{T<qLo|QZM{vbcb*Yo-y@i;bec1{44Hx+a^RNB!| z?a$<s2HAa<*3oe~+R?E_)4b&S4tj|2`06vz$qcb>Hp+LW*kr&z`S!Z8=mSgz;upyh z1+)mj)}ierH64L~bi_YCz-WmgLO?)7KvJSYDz3n%S<sni^UG&B*@pz*eqces)8q%! zBMWJwf#N6vKk8JWV4|!h34xklfxk8%r6AR(j>%3n_q)RHDGxR86c6^7Qv_$^`flM6 z)rkhnLnzoU1!>}aJKW7~3TB;SThmxmUGr*v;=V2~EmJmgn|$Z<%C1tDfi$A=^#hai z{kZUus2pO0rTz^hyGaD64J3$wQ2>|SSE2*1!Vu_ceNVo;+L93Rr~DG1(>F}H122C8 zoNasBVG3@AI^v8Suy{gh7_`6B_AwOl08mm+bUVK&rp9ijmy81a`;Xf^%elyEY8@a* z3+`pI5D^ZOgiOIJ(iz(sXzsLZt?OkRUbSb9FLx}6(6K<!v4RNDJECkWc{Zt+PjKz^ zr&)L$S*i(u8$i$n(nOSb$k9Q>V?jP0NVhqQ#g?qI|8DJs3d)I+qC{zo+nkvH3e2aU zNzUJ^*ZMdHk)I0GjyRi9R^Bx7X?OaX71QhOVd>KEUu1$HfQy?XR(U#>Bf^gN_OG(5 z!7V7aaJQw@Cu=fMhQrMJ1@eiXo7$aVF}zxbOV_?id^|)J3ksr%3Ez^9tn70sw8ik| zgK4ZhTl65ZO8Ln;ra-|3v3&B?_?kl!9EXz^^$7_Xi7hs0R+jirwi8)_)iL0Bkb0nU zR7ofz9&@aQwo1M&G&Z3W8It{D!*^(tS9X#+R(NC>x2fuiX1J!9c6?<B*`C;S5eJ#T zhJWWKBad~OUcx=&oyAgNhFbKz#eqBE2AWO&+9yl78pv<z6^}Ct`tP7El(B$w=Ob{# zVBWW39an3_5awJN>v*lsA1JRt`m8(@ug>Oz|C<Z>vFyQGJ|~TbiZeEO^?nww9ZPn( zM-^=}?!Ot|P=;I0$$FVgH__;p>%e&K&Ln@N+W)&3jKn0&Hh71hv4^ILs<z0`jjJtc zeex%Gsu25Q!X_(-cxZ5T-cw(&KtZ6(?Ti(Z=g8y3WvjfwhFmF9R&zdr%j+Qq-O_xS zaWlmyKSN0k0?k8ZH*g;4>8kWIe}tDc<fB(+opY~l1lo}X`<Qa;;kvLGY%EZd(8x<@ zoO|`O0>e75#;A_t>&Hfie7w$FJ7_FY9J3of(!Xe`u9LXUyVF3U1hgW{q`<#5+684r zfy+14%&?C!4CkK|Q*T0GTN|r2TR}qSwl7l{#5a+X_G^(qCi#4|yg;V4A|SJI{kr>r zt<q{A*x8&a1|x^Wh1NZtp8a8X+chl>{wYu>r9q@QS!Q0iXP|ZtAL@O|cwZSUf8SHR zchE)^4<QNMt9)ue@iPMC)7&Oy@y2okpV@<Y?RFoeVn-ozwO$?mA4WI)S(nM)Sn~6- z#&b1)X0ZBnhVJ7O&+rB{H?!@c`-qDG9}k=|wyy$zRdw7=CSB$JBY+BEk;w3u;xA~b zv&rDi`tZBe`af`kRc{7!x4S}i1BSK-oMe!zxbz0_qyKgU8>`HhRGn_qy64rJrqI}6 z@~^%@E7;?ByA4oSf3cy`YHAPFJX-F3M78a)=s@gKkFFlne+H4S)D7=@g6{)LIu3l( z7s-j$==;<tG9>wKs%c4HF?OI3jB`#N6TZH6(7HN(uNqrK9B16UdwkzC+vr~8F4&)2 zhywi?JQ-89=a?X73|t_@>}%-m%^>@>ir}Tay+<y5IECq<vsOg$sh{QwCPHvQpVWb= z{f@c}Tg0pcSjQVe4Sttn`k=jPEjR%Fw9Tbyvw}58sG<&((B=_{vD2QWO6!K*lef11 z4zsHTbE5=~q+)cf^U<lok}z=n=yyNNfH?8tudw$I5$Q$(5zn%KLVZS4on*`kyga(P z$*N{lP$!*352K&V*V29IbU)y>jFl(gs3M3ap7bUx7hn|W>Lj(>`s>h;2T4;t(KBy( z9PRQ3=T>KYeizjI#Avqt5~ZCa=#<>#v{(#@UNE8ZBL87Qf~2~etLN{iX}JHCC}NgB zmzyyoXm1?m{`JZk>?p>!JzZktTWSjNk=1?C^Pvb9aPmn}-<MAf%>GHkYu7WI`jTS% zH&o&d|A>z;C!0_FbZ)zy#p!;IM?_)ynst@q!k}(Wm5%ens^{iHVNv2-DyzKC4V*9K zS}AOl6qO2x{FK%`STI~_UaWE0i1y;qu!kVXuEi8sZ;;i!S6F1DDWCG#KvO?AuIjhq z{3`6Fj4J#y`@pe;p&F)1_m3a0b{`cj{k*E3?58<2AT`J|54B+tbGpl!Y$Uj;*GnNd ziYn_<e%;eAySk}opKsJ$MO97GqO-m{3U{4UB8^*$+Tr-HbUTs0_KR$RfnMP>W+ocU zpT@5n7%UPE{uyl^7=zXSpTRWy_~xf!R|JjOkN?aP<OCTEhbS&-eGKFIuqy^&z{X+$ z$|~YFh&VzYWUWDA(c&19f6NH~oDu5?A|4^mFRS=^f0X$=qX9ZhQEN3!&WssGp0Q)? z;6+v#g5u2uzw`^?U`?Wfmt%g2!kS)}A1U=Zke=l58*=oS0Cv~X6QrkSlRIWokKusQ z=1%Py+>eO?UB2ZBS&pzUj~+4UJL4|mjD9b1qlh{6*Brj|7uTYg0y3OhEW5+o?wC0x zH!i!qprf?#Yrk0v=QAGmM!naS=~bT0LjK~;YoV?T!WH0SApx`Hx@=3i2216rHcREd zTUA9TdlVV{?H!vW<j=Jmv9eqIq_w>iDMA)%mVaj~x(_b932v2^1O&ieirqmsHJtYB z!`vx1tCNt1IRovu8dwmMzvL8^AS@sdkbam{0VGLHwt{I#LiO`V)oQcwGR{6N0fC<r zg8V9uPwXzBi`l@LTEL!WC>JN52FKq<;=DI?RSxTo<vv{MMY^C?-t_%gx6OTNPRr6g zilB{3+sjK3Oy$wyW%7$=*qBSAx#tY_{+_^jo`Xn?+oyK%fLMeDs|bkcrVx0SJx%jG z^NiRykJw^+?aLYLSTqi^adL|#QE+|xq%Q`#fF8#W=;E4^r+R1T51Bu+&a3D^Qx~BR zS%;ULV&PC!%ZG4xTY*?a1Y<e;1~O=DjYfWlEp^DsFK#7oT7H_gvHiM-mc)st3G*RT z&BUA1Yz14iU9fQY4f<R=a61A%;MP^36ft1yK7XRzp06FkEd^k|$9=>gRZz_s@W^5M zB^ETAkSZ1_YhXIBd`|mwbrbuXRXnJ!u2rqGG)GlonRyOFF)g2`g9WJ)6Nx;ySKUn3 z2MP!-X)8NBWGg>G+MJN^ZSzLteUsUIs`s<WX=P1QRyliM=KLgzd^5~E6nSwK?qv@> zw?M)rVd53zV(8v4+WD645z1hCY%Qj^Ziu=}#hjV`;{jLsZr178`=7(txM9Mo0#1#! zqO!W<VBd5>Ja3N6b8`+Q6Ztjqic8|-40%81n}qn9<z0&HLCo7UZ#?UFlfwF~*+%y6 z=KHla=qDa`Z)fIxiG?=)U?q_@SRDmUt5E(BFy-4Ke#uew1g{)$w)!HU6102xQBV|P zL41_O$60;Bc!>b<0a8Xl2meq2qq81wPWkuc*BENSZ345(IO;MPNS*L~4&Hq_P~{JH z0>~Ecjx!H*0gcdCnw73Rf^U3oTs*b2?Nfb6;JqfO@s{0Z#eer4(@gsaaJHTXmTDJ4 zNi$VPR6g`t5G34}boyb}u-k<Zq_G4O`RKGzdv<%VjNnQ8e&^PFJ;kfWy|9&c_VU&T zv?zP`Bn{z1K(hY{2)!SuXuj$M#YthA7dLRTMT`5DcE^>Dt-gS@PhO4m><aztn=zzy zLrG~#I%%ue7%RTOE93s;54-R(T>div?XiLjP(Q*3gTaqL8r;jC6IpJ*X*GCXMaki2 zx8p74au6z2!??YkKQ|HT{tOh(V*icvbEcfPIB@P*gkYVKi)ZeVR+U?+h)%w<d|j;2 z+i@Gs=^HR!h%%$vbC&9lSAoqm1a|UdI=gThZ0y00s7`PXN+XXr;E{6zaw5B$pz|3n z)umW#?pWaM3Yp1%Bg1lc30j(nZ1|3IPI;5@4a`+Fteci*lGR7GYo|E7@8J}N9PhT< zm7FEDIJ&gByYfm7e4QW7^~8eOkcU7`J~n0ao6f=A*P+5APVPCU7}1Zj(jefue7ssB zSTMTj2fMg3Dfl<Xd0UekMcd+zO@!MwD)ie6klXHToR@iGxVJBM4{^^;&o?S&Fw#g< zV4ijG52K=+=g0Gtoixm)UYY~t7u&dgq$;hB^Y}_f^dzol_(>b-3o8MNL*;EH6)KQ# z^#Ze(N3fC)Um=T25OxN!KFrE68?X$U8ZQi7I{{VGyz`A~?oU~<5-aotHl>giHG$)m zhDv8f56Q3|af89wiNjKC7vMJQMGgw}64!=E70em^q`IW;gFbiq>5<+8`|cbyl0#Pf z%sgrF=*tf|ip2(5NO@({&#h8@$#)c7+;}^lTOOKTBxs+J$DMPEbn){@%$+N{mwB~i zFT4#%m>07_#zn`Vjyqf{o1jwp`K|*MIHw8QSZ11w=Y&Cn7;^h{=$uf+O!z2rx3UXZ z;E`TI5qhb(qzEx9T7&KzeXGgxq*|yLXsGDwJZ6-bDt4T8T4sa$Fa@(EW3%Ll_X#IO zE%Y#LyK1Pub{xc8M45Vc*VQO3*LnN(+)}y3k8`*NX<kM6J1)Y9wI8<}E3Cfv1nTmZ z<t1+Wc<p&OEpp8Ay3ftBx+etGDNTOwYO{*3^NoSUb2kfJIgk~KLAv>RI<9T*Hj+Nv zHGj=3R2KQ<s^js^guI)K6YH)Q#&~iLn~x8%#J%}^BcAmm&5pANer*?Nf?NEYvN9`d zlShw0ZzH%E_@*BYu}pSipm&2!y%h60$i_i}b(!~b`-i76<vx*;sPRmJ)WKRCXIXi9 z<!(Q!{LOC^8Mdb~KV1g3jVg%U=R0(-26v-p5bU}6G~pJlTPP=j8e934hz)6LxehOV z3a88MG^wK#Vd$a0>-k?vUy8nUR@nV)1@ew+jV5-{LEl#R&4uGCG9h}0?@||U_9gpK z0<J4{{_*QjBu`$d&Z}%@+a{j%t0!&e+_sFa8P1i^Ohh96pES5QbLy3qyDXzZO+_{o zCtdM$RWml+WwOCza7^mum{*EL>t9G{rJS`#Hinh=r~ZZ;$!tY|Yi%yuJ&=U`9rx?7 z$Whp=-!?~StNRT@N{_)Agcn;PbYA$aC@#|6qw~HdQDe4ND!{yWZ@BUlyS_K#;g<A& zlO?^lG_#ZllTJdx*4;cj%6qD#KBB(rI4agxb=xh&gDwp4O?ZDJd>CKyKwamXPouf? z!G4e#GkTbQp_0qO>ds<~uv;xR+qyX7k9AR$#b=vGfmGCj)?#C!$5K*fLlO4yXtOTD zCz%dCW|^tt4Zt3uVAlCBd8(<6@UxCEkfdT@I72~f3q7!*jU#@)N&WS09Orj7A1>#H zCo4nyh489fAAD%Cb96#!I*tNp*>?4&W7)h@P`XMv7n9d#!Y<!IgbYd##S*35MIB72 zgV|d?v=+FR<<mcP`r6iszLa6Q!9Rs#Yox{Lg<1GIAB2f%I7^)|<)i6?NuHXPWR3DX z0L!3MNBAsTwHzvWV6c(ge0Alp>bf0;#cUMFz8<P!-ZN#D#ff^{_U-(=woP<5q|Xoz zr=T;75}L%Fg#rQ*daG@=E5q)otl`c+5#9~A{Ek1;@G7XwiM@)t4YrAmf7!_>w!qaf zuQ1O20+Xt4ZT|qu`*mN<(?$o<`Naw|*;pL`UTsljCvr?ygIKJYzCQ{$ZB&LMY{b1= z4Cc<Xi#cO_Yajx(jZLQZegE->HlA_S!LPle1%FcDotK7|_Hvie2OAw7T{p=2$FRzC zh<RMww@S@C1kbWJ6ZJ_7CvpxakqIVIs)8Q(HEII<61dk{mggxE#DS@wB91eG$x)nJ zvA|vSs`rBl97f#&kj*rG{l85TmT2|t;^JSIY`>Z+Yimao6;WjVn2~ZlsU50dg@4>* z2&my7n=&H^>p0`TOKAdmzdqa99d34AwRg}~x+gOGy0S3F&<6KMekHog^PEw)o^?&x z<?OUAaqUz=FEp6HZ)#SV1P`PrKfh1wN5qcNB+q$!zDs@US!S#9wF`qkyUo4b(X2B= z*>7pJI;qIze{U(O)Os)qH-ua!Ik%tyOak8t&Hd)UFf1rEm>{;nwO`(f{0zKwNUiL6 z92DVdO=(cv0||GRB5PJ&0DrbverRVT!G4%HMB{a^Hp;SM{+i2;6o&*)i|xdp7}Nd6 z34w>|`=@6sf*N)u)7LVcLI<AKx4z{f2Z=AwwTXL*;oifNp{nY*)4DlcWE4xFL0WMs zn{GYStx#$)<+K$Qk(VFQ5*5KfhEr4s7Z3OH8v!^g9pqSMiH$U}USc!MFT<K4aI4W> z1ILjBl>^BZvqgqSbEQ~3&K0mcdj0yCi+{_DT{h3k=N%E6hwGRaMMdh-BtkU~*u_K= z{-y)`UIv&b+YW4brs8mFLA6CF{GIL7@-L}Hat^+y5^rI-p`yk!ZfltA&PQS$r#E$? zn7>h7TF}9nmL;*6dBk!^lp_srn#WeyTXgn1RatXkGi&dre^UO=%j*vWrcgH14tKBO z2A6O`wx3MYG_)-~MMq)k1kx;%oogTCE{O1J7vtGWUx`}a#W737h+K{94N_oNn0k)x z&Z*$YO2rp)R<+xbPrUaaINhtr-*z{)1|5w|6lh{b$;xKb{5V!RCi#kD)h1SsYl9p? zGp%rL8sES*)QUs^t!-t$x~I%pE{~wVSvO(dvPN@hVD>1O>#g_%PXF?Qn8T+NJIC`) z>M1BSy&m%-&fO;RPmR;LT9a66(grER3{{=f*^11YGc*MOAG4c3TJ&6WDfnEWpKORt zFQUK;vHyS&m`Ik3T$o0v?}1o{D_Kx{7Nrx<!a+Kq*G9{hLihe`cHh?+K<OLbv+;?d zF}}cR@g&p9B9qqqB&_$_<G#6tg^;5o>(P3JGUJp@Z|*(h&u`w}&-l}?1Klqa`BZbh zW8Kei$*{R6-Smd`Nfc;D#yxUWDZ6;UIh(Dt@*OiZ&zKUQeZAmq!qn*pkb9Dz)TWmA zvjaGQX0ADQu&+N%8vlOQeyPy}HyYoYZcm3}s1;V~`z<3qy_lwEnv)`LK^4JqTPaG? zUyg}Wjr`-3wzC;%0Kq>Z=iF}Cl3rVbts`n4$Rl|@$uNA3YxUEO8~R$kl-9J*C<>l2 zfH|GMeW?a)Ee_b4`jMp3WGucPQg`sk!hw%s<cqhxW>D?&z5W&FLjMG;%obTkXsH(6 zFJCNcepx`CbAb&^kKvA22OB(@?k?zfCCt`uX&;X8<lpRMW4y7Z@`Sbc8M#=1un;T{ z5kCrB5$3b#`Y?=5;!aP+_H(233kB^^t$9}g)ZsT#t{MM{aXfai=|Xrr3Vh{`bI%mQ z7#&5j+-EWjeQ+shbla%sB7a!j(T)-&dH~xa6}B?Y1wcqwkz?&87+dN|IN5O-aW3-= zToS!_I3F}<0nzABl$scUI(&g+Vta46!HD4C>|x-rtsm|mdrA`Ff^7{EFUQ1;1>-N5 zk=lxqY71)o2gr2*L5te4b`AWOz()g#ks?5S2ARHK@u<1xWM+M&|H6BFFe)5-!Dg9H zC?PfzFp(RVfcQVapWt*TXs_-+Aa@5aab)=d>LaCy6~qG$7so?7i~r0p$S4J;B`1O< zJ~A0!uq`MBwk=<~lz&@u1136f@Ctvfk0jxnhYT6m_HXMEXa%<T`oGBk*5?9H%tf?H z{<*cJoNr!SCVb4tUO>?O$hqeIz<(M3ZT<iI(fkZn`CVp-YQAA#oRxR~g0?VOPbqY6 zU*Q5N(JdWq5qNd@`n$!jl80las4!m@TY$CJ>G01EIA&Mfqbna#PC2T)R1Tx#w_ri( zyxJpO1Gxqdkz_=bUEaSHnHp@q8WmW8;3zK+sW(nZv={y0;Gib??#>ykaia{T4SC#& zm9KU(&GcEhqaJCg5qs%r<gHnyDP2y9JJFo&6nI@b(enHC?D^`^#9-D?ybO!Kw<+r} z%H~|yh}P@CC2mL=8+!Mk(|V+~lCkpXX-TgN*pbKIG!HpZDl}mQ`kw|!0|Nt4WloBU z$RL)~FePwxu6Sy6;Jt<=t-K6NhW6B|Jf3h&O@oFbkWy;Z5CbhPcw~O33S=lEWGF(~ z+PI=bT3T8ag)J@EQBhIm=H}hn_2&@jepf=*NnXk}E%o)>PCCWriT5B7F_P62Sv0m0 zqbtgl5?N__*hPCQ1^mBWA=~P!ZmArk#HVbn7fS0-_aHWxgw6&tIQ2pJ*FoSEylksH zt8a+fxxW*rS!?_2sfyS=o;)ikD*6Fp7%OU8+5kFoa`I%?P+u?>BAlxMK6QU$C}k9Q z#MI~3&}|8j7)H73ux%B{jiN8jiI;7Jh3dhjEd>c}uXT!*o#u1DSsQ{5mv_N2OP6h} zcVw38$PbxKtA3+Wje4S2{Hx1!GXg~2`;zR>Cy=E{Bse_`Vxsmr@Fo_VzKEb)M)*vb zKmf9hK~xX(F}Q#Q4iDDn9dZ6t4-9|~)B&Zbehhr|0fUz<sm6d0ji?g@NH+Lz#V54( z0;uo)#gwnZ|1>n<01coM{D%H9*u(@ugYR;UKLspJ<VS~;Yi<I}QvNXj?g9)R1`@0x zKh(j*51_HIv6|o?gNGm>u|CFZ6P^!61Z0y`IGVn{X+8%3cOJ=q`TKBDk*r3m6LlG; zDo@!LoCp$(2x6Qt5|k(+lxQM6==_gJY~|0N+Ska_Ziw9PCRy7K_^M}2(=*K7Ro;Tf ziVC+pztK`EoaHOgNv(>lx5Am6T>vhcUzVX&vGfKzYv`r{Q>P%x+jx6mU#i(;i&d3J z+)HlWtfRAz;Qpi?@@FuTa1WCIPonUM2xK@2MHQ96JV9g-G2r2x#h!BS`ltH&!?Yp@ zlcB$m%J`Ig=DnrW+>BIs2tz9?D@Fhns?FkY8RUO^*w(PIv9VsRGp$|J^Ps=$g=OrA zK)}ONu8aZXL3_Ek-;Oq&%m$1osVowlTK`sHPKF!u8`a%uU*v5R%$QvT<p8Lx{BJ6q z7Br6>%LdblMS}Y==roW#jw_pt%GD~vZI)_?8%@WOm$LahjOgj<nQu)wFrvf<;3h?? zp5cp&XJ)nbNIpc)pG<I?O0ubb`cpIs0Fv&$ni}yT!7%YaV);aC6XKsyr8_{%4(*K@ zK6HRe5+K3;KX`CIf$2kl)h_rCRG3Qyu=OxAr`BvAa#{dz191eL|FhBmXVxKoEu1f3 z2j5dEjwiKW*#%YQmM7!KAjFM;K#D*oMWGi1<g8pu;gZcD`1tq>GwbWo-FMomjkwb~ zoU^FUTQ(1QuTc$I)FxRBIZ-02*dMOrFy&IrdGF{UyzuN#qIAxPDvxsIA|O#I-dJ+I zbVFfCd-%#ezaCBVRV=p?GrMJ9G*FlrY#s#$LQzL2GFJ-|9i2EgInS3x7$YX&Cc_S^ zj3EZrSwly99)S?_ccs0?l~DXVqtfoGeJog`NSq{GoFo`P^pX4kv1o1b-rhGh4vxH^ zKY?w>J8{h+Mj^-5?_?WxCnomh?7ORMx_!52_pkS4x`(t?ume=F^ZfW>ND96`%hGsJ z=uf%zh^-uPa_D(x2g*?JE89kAgrNUP&Q+gB(&2yRS5Z;PZEs)I8XX-al97@MlcO5~ z^ar32fAj*76Oa)3!G#9jWavN#+tJ7S!2JtR$TkB3CKZ~z=wC!s1F+;v$R7q3P`d|& zQC#8!f&V6=!QXubE=0gX0ELKFWG4O$R)E9;ygp}3b~6L~gR#n2$<IItQ0OQ8H0FRp zKfncGp`-qbqyWo9{?GFN8*>0%5ag)eg1Sd_X_^1-LTi@NYA0D^fwaHYwXH0i>B%ri zAw5NnYt6sn%0oD-3P18n<>QNvnT7Ru0D1)eD??LZEOVYqH>V5NpbIu93zuZ`m*n!d zBa4iVMzFH5wH>-+Z3M;jZH;q^9&=J@=@XS1DnVFBH0OgU^A1Vdil1j+&}8;lxi?Vs zN)PH!viSX-aWmije!+UJZQhHRB54yD(;~il39l5VAQ{6V1@cv#Ut4Yn*U+s-fu{&V zA?NmK)P1%7njjEI7z$hjD%c-@5Q#Hki2`D<a1b#lG;ksxy%2!zz7oJ8mkVfTC0qP= z%BINLIT+Uy7m5PY@T;Lg2rRF#a6l_PBV&4LDKt13;<vXkNBaE6XX4gb|If9}<$#PO z0~Ke2&;52xrDj<6__NWs(<M9D4Kq*biNR3?3IXT5SmNZdU}#r3CcK)q>GqwmIaJZx zxLz)ch>abHFR289;4i{+@F4udHn&gz0jYQ!t-1c0y#RiQjZEWDOG)%^?=j0IB_)Dx zZY>RoNl8M<@$o1!yrP02z~^Gd6R6#*Jgc0Py2RtH=FPTc0y433e3u7d(F1HbRKH8r zvR>3&Fk~QIp2YfeWtYcJ!7*zWY^{sd+umiCwizHL47ZaEOsdUNb(Q0m_olymI2Iz{ z1c2s!y@Ig$)Wbf2Gl6qK8viU-{m-5MliUH5Wue%x(QHCV<sg1lqr<PfdZ&ARDcW1D zS-{7%NFiKy?HP!%-f`?VVDL*6HJ)$l&a}Lb(&Sf2x(=*&Fq#H3sx=bLD-0DEhLVs* zhE~u0Pk@HJ<Xuj4BF?rre$&rymV7FIE2X>dJWH_IA@;njytd+2(e_fhkTQP5)6`l% zUTwYVfL!KnN}5C4v@J_x9F;(@qT<0*#)R(or02SqjfxBtgrp)vg(Yfl&+Koh3y8{m zwZcKlL^XL4Xp=$yGO6eiV4l=8Czj#!L!@eN*X#t4m8*juuks+ocRN#|gb|~J<yKb0 zlOO>?q=lK))yTm?aW*!#&B|pL&~d=;I3GtlOB+|0mvvqV?5<5VvB-nJ3m+?Xm$0rS z8hl%qJ@$PISX<tc;Lx)JkZDry25R?C=x2Cz@<&$Ew8CC#GGlwsaG;U=K$RUVWwH|y z0CZ9@GBzdG)v*jKDJdC+g+NIXf#wklr_=TziBpFKqm~4vihsE)MO}lnY(SY!cY7nH z%33{5ffuN|^pCi~=&G})#d>=tOPkIkVf~%^nL-GHJbvUFh<bTZA3n+e3gCc<qjNvo zfMO8=Ug~j^eC5*}1p_h)tSzYY&lKnXnunP5*g*VzgGM_S{)mZ*P5lYJDEK%#xUk?L z!86Hr9GFKYp&}L|j(6Y$$7kWW6J=O?ddjkDo=+cNunmB)T*dej_?XZORbN+=muhFK z=NY;rh6d>=*>X63CU!g^EN{c9ChuEXP+Q19Lg+nO8UlVL+;r+ikP~`d08@NTjpfU? zW@oSJfVY)z9>Yj@+FO~nV*YE^QY+wOm9<Zl(O=R2ucsh&^xN)eQZnl;mx_+gDmhLN z5!TE@I}-dU0NE7lC&7XvL4)~OQSsZsL|GXXgI>$T#_U-iy=vg9lHr&iW2x2Ls!?{* z)ypIStCmuyz?|?7cXk=i{orS^hWT{trDU=&lp8V578E_+7S<^hGW>GaMK5eEpfCY| zd?#08#Pg?Q=r97jMiQ0BaHQ=`s!4UxHxw3p^(5{bRD;9PMks34|3m0bcocz|3B7TI z-m$&JAOUh@c~$G^Q?I1}Luag4$8aa1$a-jS_aZ`k15|Lu?pNcB&7Nf2@2!Uvcn-tF z<~BA1>m44`)6;&9jU0h}hN2{K441=!qyidTKq8cK_Nq6%WkIQ2rg`NoUw?n)2<zDT zPrc}mwFzAuk(H(>9^i+@Q~45zR`m7QERDmhd3n8BCUEgT-04sZr-hBt_hx+nIGJ-^ zMT5w%{QgoBC=xuB{+=z3*Gj>Dk?1f3$uC}manSm~Dv&gcj4rA_n07kh)FDi$SeMBl z{;N2pZhf;vY=zBZ5*gPJ8bqK^CKo77I7>q6Ms~nFwsQHms~r^{8LlIA``2|p=hGkG zwZ9IH#Cbo9%UO0s4h1{hTQan0pj#M8l}{a&8J5Y4_WyH;6ziK5CF~b5Dw8;S@y^)M z#z@w}!omeEGJhla>LAkEp;u(9c60OJPhUkc^RmeBWx#~+sDv`9ucJ_GbCPGwu;V$D zg}%gSPeTbMDQ?NRvdr^y?#9a+x&NoA;V>W`!%Bj}vM{V=ja(aRFP39J+}M9uE4dKA zer|?6vvDQP)=%uvn43kWEL+EDWuFa(43|cN5{I<k*wmzLAxd{Mmsr7y+#d{hNo_Q( z5)m}pNG_($EJPvz5n$H%Po=BHC|^v(hkJfIYKattL|;LS$%NjXgsM)`y+ab)8{tFH z;IJPWl>6gQMA^tfugOA5J?>5{bb(TCNC+<d<EW=&qg6CDEtccXy5eh|S!ODE5*`7S zS>69$e-h240>Gm?F*H%8eu}W!UgArX;z?<;n2B?<zr$(ffQU}+MM@f}{Cj%C1lk<T zw8s8&5WgeRU-Ok2lU;X%YfvgH+WMKJ9}a@WvI@V<6py<WQ2*e5IhAah+7o0YFn1+d zI~m*~&FHTJ*w4h1A@2kk4)^XRSl@L91|?yog72fyCp`63YQn@mG;+nY3p5PMh5%e( zeV5emyY*O@IiRanYWW*7<a6F63EB6)V0c%t*Vlu7!k=YjWyifunWmDlzrL6peck5l z5`f^;>(o^5Yls3`UNTtX7#V;d{m!axZPJ}WWW`-RRgdey@pJpc4F4%ow0Pc(vNf~$ zn_C^vKD5CZel}W?8ZkF=)DJpz%3%_@_}nZi5~AWT6q9smIs8*)B@!CvpEb7HX}nMS z^Ys1%X~rtdfFc+EPJC)a{~kb7^?#YxRS)4KgJHY;5;X6LZP5F{FHn6g&DPIM0#esV z+rp*j9~u%G3>CmX7HdO~7SE|HWE~rE4hM-slmVJ13D-il7cax{G*-HlJ0RTHKy|Ra z#p97lkZ@x^vb`9~nMaH;C5MKL1dKGFVxuWDg;`;uO_*Rqk5iXbN^}b-{mzo(I!Jvb za;t#$5+B5c)`Gfu$iLqX5+0(oO%jOSNisKw9M!_35b76E(f9ws){;T|7D<AeM|>x_ z_P^6$1cHkU<T2%+yKP!IL?l%l^XK}(iRWCj?+j2W^-8C@tgB2B6^^5tpdErpG!Hdr z#*_55Mlc~^25>^AYT4XyYT?kU<1NqCYqk2gonZGDYQxr;mMCe$R%``n3C|f7j?`wX ztP#^E<I)GkW7tRN88S)T=y66?6Hy^)r$c%2WxwhD9Y~Nh%1MwhVnC7M7+KzVklYXQ z!^7=4`hepB{vrSX#tFQL(tZrPp#TUT>8c><M?Oa+0El0&`s9Cowjf6WAVsELc(zXr z@c+4oz^?m#pQ!;1j`983nr<AA^GQvksjiyZ0=mymoP@KD<0nok+0EzIjur*Taqt{} zGed>+twBW+@U>d`%#|o=G+Jve@_#m1`hu;Z4v!s{@_RjP0{$FFNKw(~IwB%MQCIgg z8$sq@ihUqXgrxja`?rwu9u^xy_4C2|mV884m<FAL?|`2m^Y>0qo>6x^=kD6CVso8{ zcmgcSh5Ifo06p^g6o+_O2M!05@8dl&W>@CYd2v!AO7v~{`JfJi#LB97L}u6t6tV!k zuH7JJ@3W~g))x{CW%Arr!LBF*DoD8bto2wyNol5%hTFdNZ!m7_z_N>)fHyn{DuLjB zp`#O-P6%Nn)MoPj@yST@zmCpWU%~b8;9p`@+QV5YcT?r+jF}#1Eq5i~|Gb98fMbmJ zZX9b>J`}0KUS6wqx9FcQ(&M}QLk$qfVi)t69d;z)2zZ;zRUH<XZC2R80!D0Mj@_=@ z9y&)ahkFq#Jl|ctDeM&C2UW=@XH|j!b*gql?=+Z)o`0LjTxJ&|H^bNWTvQS8bbNEA zOK(AHo1dGj7l7=1?{HMm+Qj!8|7^9V1iQKuZvLLhJGIzkt1n7|42R1BSFTptQp5=L z_bCPG%M;1(=lOg4=L-26%ku;*@%aZ3OWv;^2#DZ%j>Ry($`j6DBsMGFBb$a3#GXDF zIhMfC1y4COR&)t7@+YD#$krdW@pjil?4qn|KY)(lwSv1Zn?Q>z&DDA96ylZ|`Js1` zMCZlD{Xtj12P-t|F92Nt)A=6}uOJ9E*(U5X62Hg&+38(fFC2^UCz`2H1LjGJ1AZ*T zk@@j-9pG`sK(jdX@_)c5B>?z*U=RfQSYk7?`yG;~6GEYO4=fFU76D0#$%$488~Fb( Dk{Iiq literal 0 HcmV?d00001 diff --git a/src/CMakeLists.txt b/src/CMakeLists.txt index 175def8..9bc4cc8 100644 --- a/src/CMakeLists.txt +++ b/src/CMakeLists.txt @@ -65,7 +65,9 @@ SET(sources matchers/flannbased/matcher_flannbased_load_yaml.cpp algorithms/algorithm_base.cpp algorithms/opticalflowpyrlk/alg_opticalflowpyrlk.cpp - algorithms/opticalflowpyrlk/alg_opticalflowpyrlk_load_yaml.cpp) + algorithms/opticalflowpyrlk/alg_opticalflowpyrlk_load_yaml.cpp + algorithms/activesearch/alg_activesearch.cpp + algorithms/activesearch/alg_activesearch_load_yaml.cpp ) # application header files SET(headers @@ -110,7 +112,8 @@ SET(headers matchers/flannbased/matcher_flannbased.h algorithms/algorithm_factory.h algorithms/algorithm_base.h - algorithms/opticalflowpyrlk/alg_opticalflowpyrlk.h) + algorithms/opticalflowpyrlk/alg_opticalflowpyrlk.h + algorithms/activesearch/alg_activesearch.h) # locate the necessary dependencies FIND_PACKAGE(Eigen3 REQUIRED) diff --git a/src/algorithms/activesearch/alg_activesearch.cpp b/src/algorithms/activesearch/alg_activesearch.cpp new file mode 100644 index 0000000..cdf5da2 --- /dev/null +++ b/src/algorithms/activesearch/alg_activesearch.cpp @@ -0,0 +1,95 @@ +#include "alg_activesearch.h" + +namespace vision_utils { + +AlgorithmACTIVESEARCH::AlgorithmACTIVESEARCH(void): initialized_(false) +{} + +AlgorithmACTIVESEARCH::~AlgorithmACTIVESEARCH(void) +{} + +void AlgorithmACTIVESEARCH::resizeImage(unsigned int _img_size_h, unsigned int _img_size_v) +{ + img_size_(0) = _img_size_h; + img_size_(1) = _img_size_v; + cell_size_(0) = _img_size_h / (grid_size_(0) - 1); + cell_size_(1) = _img_size_v / (grid_size_(1) - 1); + offset_ = -cell_size_ / 2; + renew(); +} + +// Functions to fill in cells +bool AlgorithmACTIVESEARCH::pickEmptyCell(Eigen::Vector2i & _cell) { + int k = 0; + Eigen::Vector2i cell0; + for (int i = 1; i < grid_size_(0) - 1; i++) { + for (int j = 1; j < grid_size_(1) - 1; j++) { + cell0(0) = i; + cell0(1) = j; + if (projections_count_(i, j) == 0) { + empty_cells_tile_tmp_(0,k) = i; //may be done in a better way + empty_cells_tile_tmp_(1,k) = j; + k++; + } + } + } + if (k > 0) { // number of empty inner cells + int idx = rand() % k; // between 0 and k-1 + _cell(0) = empty_cells_tile_tmp_(0, idx); + _cell(1) = empty_cells_tile_tmp_(1, idx); + return true; + } + else + return false; +} + +/* + * Get cell origin (exact pixel) + */ +Eigen::Vector2i AlgorithmACTIVESEARCH::cellOrigin(const Eigen::Vector2i & _cell) { + Eigen::Vector2i cell0; + cell0(0) = offset_(0) + cell_size_(0) * _cell(0); + cell0(1) = offset_(1) + cell_size_(1) * _cell(1); + return cell0; +} + +void AlgorithmACTIVESEARCH::cell2roi(const Eigen::Vector2i & _cell, cv::Rect & _roi) { + roi_coordinates_ = cellOrigin(_cell); + roi_coordinates_(0) += separation_; + roi_coordinates_(1) += separation_; + Eigen::Vector2i roi_size = cell_size_; + roi_size(0) -= 2 * separation_; + roi_size(1) -= 2 * separation_; + _roi.x = roi_coordinates_(0); + _roi.y = roi_coordinates_(1); + _roi.width = roi_size(0); + _roi.height = roi_size(1); +} + +bool AlgorithmACTIVESEARCH::pickRoi(cv::Rect & _roi) { + + Eigen::Vector2i cell; + if (pickEmptyCell(cell)) { + cell2roi(cell, _roi); + return true; + } + else + return false; +} + +void AlgorithmACTIVESEARCH::blockCell(const cv::Rect & _roi) +{ + Eigen::Vector2i pix; + pix(0) = _roi.x+_roi.width/2; + pix(1) = _roi.y+_roi.height/2; + Eigen::Vector2i cell = coords2cell(pix(0), pix(1)); + projections_count_(cell(0), cell(1)) = -1; +} + +} /* namespace vision_utils */ + +// Register in the AlgorithmsFactory +namespace vision_utils +{ +VU_REGISTER_ALGORITHM("ACTIVESEARCH", AlgorithmACTIVESEARCH); +} /* namespace vision_utils */ diff --git a/src/algorithms/activesearch/alg_activesearch.h b/src/algorithms/activesearch/alg_activesearch.h new file mode 100644 index 0000000..c10b3f8 --- /dev/null +++ b/src/algorithms/activesearch/alg_activesearch.h @@ -0,0 +1,330 @@ +#ifndef _ALGORITHM_ACTIVESEARCH_H_ +#define _ALGORITHM_ACTIVESEARCH_H_ + +#include "../algorithm_base.h" +#include "../algorithm_factory.h" + +// yaml-cpp library +#ifdef USING_YAML + #include <yaml-cpp/yaml.h> +#endif + +//OpenCV +#include <opencv2/core.hpp> +#include "opencv2/features2d.hpp" + +// Eigen +#include <eigen3/Eigen/Dense> +#include <eigen3/Eigen/Geometry> +#include <eigen3/Eigen/Sparse> + +namespace vision_utils { + +// Create all pointers +VU_PTR_TYPEDEFS(AlgorithmACTIVESEARCH); +VU_PTR_TYPEDEFS(AlgorithmParamsACTIVESEARCH); + +/** \brief Class parameters + * + */ +struct AlgorithmParamsACTIVESEARCH: public AlgorithmParamsBase +{ + int img_size_h = 0; // Horizontal image size, in pixels. + int img_size_v = 0; // Vertical image size. + int n_cells_h = 0; // Horizontal number of cells per image width. + int n_cells_v = 0; // Vertical number of cells per image height. + int margin = 0; // Minimum separation between existing and new points. + int separation = 0; // Minimum separation to the edge of the image. +}; + +/** + * \brief Active search tesselation grid. + * + * \author jsola, dinesh + * + * This class implements a tesselation grid for achieving active search + * behavior in landmark initialization. + * + * The grid defines a set of cells in the image. + * The idea is to count the number of projected landmarks per grid cell, + * and use one randomly chosen cell that is empty + * for feature detection and landmark initialization. + * This guarantees that the system will automatically populate all the + * regions of the image. + * + * The feature density can be controlled by + * adjusting the grid's number of cells. + * Typically, use grids of 5x5 to 18x12 cells. Try to make reasonably square cells. + * + * This class implements a few interesting features: + * - The grid can be randomly re-positioned at each frame to avoid dead zones at the cell edges. + * - Only the inner cells are activated for feature detection to avoid reaching the image edges. + * - The region of interest (ROI) associated with a particular cell is shrinked with a parametrizable amount + * to guarantee a minimum 'separation' between existing and new features. + * - The region of interest is ensured to lie at a distance from the image boundaries, defined by the parameter 'margin'. + * + * The blue and green grids in the figure below represent the grid + * at two different offsets, corresponding to two different frames. + * + * \image html tesselationGrid.png "The tesselation grid used for active feature detection and initialization" + * + * This second figure shows a typical situation that we use to explain the basic mechanism. + * + * \image html tesselationExample.png "A typical configuration of the tesselation grid" + * + * Observe the figure and use the following facts as an operation guide: + * - The grid is offset by a fraction of a cell size. + * - use renew() at each frame to clear the grid and set a new offset. + * - Projected landmarks are represented by red dots. + * - After projection, use hitCell() to add a new dot to the grid. + * - Cells with projected landmarks inside are 'occupied'. + * - Only the inner cells (thick blue rectangle) are considered for Region of Interest (ROI) extraction. + * - One cell is chosen randomly among those that are empty. + * - Use pickRoi() to obtain an empty ROI for initialization. + * - The ROI is smaller than the cell to guarantee a minimum feature separation. + * - Use the optional 'separation' parameter at construction time to control this separation. + * - A new feature is to be be searched inside this ROI. + * - If there is no feature found in this ROI, call blockCell() function to avoid searching in this area again. + * - If you need to search more than one feature per frame, proceed like this: + * - Call pickRoi(). + * - Try to detect a Feature in ROI. + * - If successful detection + * - add the detected pixel with hitCell(). + * - Else + * - block the cell with blockCell(). + * - Repeat these steps for each feature to be searched. + * + * An example can be found in src/examples/test_algorithm_activesearch.cpp. To run it execute ./test_algorithm_activesearch from the bin directory. + */ +class AlgorithmACTIVESEARCH : public AlgorithmBase { + + public: + + AlgorithmACTIVESEARCH(); + virtual ~AlgorithmACTIVESEARCH(void); + + // Factory method + static AlgorithmBasePtr create(const std::string& _unique_name, const ParamsBasePtr _params); + + AlgorithmParamsBasePtr getParams(void); + + void setParams(AlgorithmParamsACTIVESEARCHPtr _params); + + /** + * \brief Re-set the image size + * \param _img_size_h horizontal image size, in pixels. + * \param _img_size_v vertical image size. + */ + void resizeImage(unsigned int _img_size_h, unsigned int _img_size_v); + + /** \brief Clear grid. + * + * Sets all cell counters to zero. + */ + void clear(); + + /** + * \brief Clear grid and position it at a new random location. + * + * Sets all cell counters to zero and sets a new random grid position. + */ + void renew(); + + /** + * \brief Add a projected pixel to the grid. + * \param _x the x-coordinate of the pixel to add. + * \param _y the y-coordinate of the pixel to add. + */ + template<typename Scalar> + void hitCell(const Scalar _x, const Scalar _y); + + /** + * \brief Add a projected pixel to the grid. + * \param _pix the pixel to add as an Eigen 2-vector with any Scalar type. + */ + template<typename Scalar> + void hitCell(const Eigen::Matrix<Scalar, 2, 1>& _pix); + + /** + * \brief Add a projected pixel to the grid. + * \param _pix the pixel to add as a cv::KeyPoint. + */ + void hitCell(const cv::KeyPoint& _pix); + + /** + * \brief Get ROI of a random empty cell. + * \param _roi the resulting ROI + * \return true if ROI exists. + */ + bool pickRoi(cv::Rect & _roi); + + /** + * \brief Call this after pickRoi if no point was found in the roi + * in order to avoid searching again in it. + * \param _roi the ROI where nothing was found + */ + void blockCell(const cv::Rect & _roi); + + private: + + bool initialized_; + + AlgorithmParamsACTIVESEARCHPtr params_ptr_; + + cv::Mat frame_gray_prev_; + + void defineAlgorithm(const ParamsBasePtr _params); + + // TODO: parameters? + Eigen::Vector2i img_size_; + Eigen::Vector2i grid_size_; + Eigen::Vector2i cell_size_; + Eigen::Vector2i offset_; + Eigen::Vector2i roi_coordinates_; + Eigen::MatrixXi projections_count_; + Eigen::MatrixXi empty_cells_tile_tmp_; + int separation_; + int margin_; + + /** + * \brief Get cell corresponding to pixel + */ + template<typename Scalar> + Eigen::Vector2i coords2cell(const Scalar _x, const Scalar _y); + + /** + * \brief Get cell origin (exact pixel) + */ + Eigen::Vector2i cellOrigin(const Eigen::Vector2i & _cell); + + /** + * \brief Get cell center (can be decimal if size of cell is an odd number of pixels) + */ + Eigen::Vector2i cellCenter(const Eigen::Vector2i& _cell); + + /** + * \brief Get one random empty cell + */ + bool pickEmptyCell(Eigen::Vector2i & _cell); + + /** + * \brief Get the region of interest, reduced by a margin. + */ + void cell2roi(const Eigen::Vector2i & _cell, cv::Rect& _roi); +}; + +/* + * brief Retrieve object parameters + */ +inline AlgorithmParamsBasePtr AlgorithmACTIVESEARCH::getParams(void) { return params_ptr_; } + +/* + * brief Set object parameters + */ +inline void AlgorithmACTIVESEARCH::setParams(AlgorithmParamsACTIVESEARCHPtr _params) +{ params_ptr_ = _params; } + +/* + * brief Define detector + */ +inline void AlgorithmACTIVESEARCH::defineAlgorithm(const ParamsBasePtr _params) +{ + params_base_ptr_ = std::static_pointer_cast<AlgorithmParamsBase>(_params); + params_ptr_ = std::static_pointer_cast<AlgorithmParamsACTIVESEARCH>(_params); + + if (params_ptr_->img_size_h == 0 || params_ptr_->img_size_v == 0) + std::cerr << "[Algorithm ACTIVESEARCH]: Image size not set. Current values [" << params_ptr_->img_size_h << "x" << params_ptr_->img_size_v << "]" << std::endl; + if (params_ptr_->n_cells_h == 0 || params_ptr_->n_cells_v == 0) + std::cerr << "[Algorithm ACTIVESEARCH]: number of cells not set. Current values [" << params_ptr_->n_cells_h << "," << params_ptr_->n_cells_v << "]" << std::endl; + + // Set internal variables + separation_ = params_ptr_->separation; + margin_ = params_ptr_->margin; + projections_count_.resize(params_ptr_->n_cells_h + 1, params_ptr_->n_cells_v + 1); + empty_cells_tile_tmp_.resize(2, (params_ptr_->n_cells_h + 1) * (params_ptr_->n_cells_v + 1)); + img_size_(0) = params_ptr_->img_size_h; + img_size_(1) = params_ptr_->img_size_v; + grid_size_(0) = params_ptr_->n_cells_h + 1; + grid_size_(1) = params_ptr_->n_cells_v + 1; + cell_size_(0) = params_ptr_->img_size_h / params_ptr_->n_cells_h; + cell_size_(1) = params_ptr_->img_size_v / params_ptr_->n_cells_v; + offset_ = -cell_size_ / 2; + renew(); +} + +/* + * brief Create object in factory + */ +inline AlgorithmBasePtr AlgorithmACTIVESEARCH::create(const std::string& _unique_name, const ParamsBasePtr _params) +{ + AlgorithmACTIVESEARCHPtr mat_ptr = std::make_shared<AlgorithmACTIVESEARCH>(); + mat_ptr->setName(_unique_name); + mat_ptr->defineAlgorithm(_params); + return mat_ptr; +} + +inline void AlgorithmACTIVESEARCH::clear() +{ + projections_count_.setZero(); +} + +inline void AlgorithmACTIVESEARCH::renew() +{ + offset_(0) = -(margin_ + rand() % (cell_size_(0) - 2 * margin_)); // from -margin to -(cellSize(0)-margin) + offset_(1) = -(margin_ + rand() % (cell_size_(1) - 2 * margin_)); // from -margin to -(cellSize(0)-margin) + clear(); +} + +inline void AlgorithmACTIVESEARCH::hitCell(const cv::KeyPoint& _pix) +{ + hitCell(_pix.pt.x, _pix.pt.y); +} + +/** + * \brief Add a projected pixel to the grid. + * \param _pix the pixel to add as an Eigen 2-vector. + */ +template<typename Scalar> +inline void AlgorithmACTIVESEARCH::hitCell(const Eigen::Matrix<Scalar, 2, 1>& _pix) +{ + hitCell(_pix(0), _pix(1)); +} + +/** + * \brief Add a projected pixel to the grid. + * \param _x the x-coordinate of the pixel to add. + * \param _y the y-coordinate of the pixel to add. + */ +template<typename Scalar> +inline void AlgorithmACTIVESEARCH::hitCell(const Scalar _x, const Scalar _y) +{ + Eigen::Vector2i cell = coords2cell(_x, _y); + if (cell(0) < 0 || cell(1) < 0 || cell(0) >= grid_size_(0) || cell(1) >= grid_size_(1)) + return; + + if (projections_count_(cell(0), cell(1)) == -1) + projections_count_(cell(0), cell(1)) = 0; + + projections_count_(cell(0), cell(1))++; +} + +/** + * Get cell corresponding to pixel + */ +template<typename Scalar> +inline Eigen::Vector2i AlgorithmACTIVESEARCH::coords2cell(const Scalar _x, const Scalar _y) +{ + Eigen::Vector2i cell; + cell(0) = (_x - offset_(0)) / cell_size_(0); + cell(1) = (_y - offset_(1)) / cell_size_(1); + return cell; +} + +inline Eigen::Vector2i AlgorithmACTIVESEARCH::cellCenter(const Eigen::Vector2i& _cell) +{ + return cellOrigin(_cell) + cell_size_ / 2; +} + +} /* namespace vision_utils */ + +#endif /* _ALGORITHM_ACTIVESEARCH_H_ */ diff --git a/src/algorithms/activesearch/alg_activesearch_load_yaml.cpp b/src/algorithms/activesearch/alg_activesearch_load_yaml.cpp new file mode 100644 index 0000000..48a24dd --- /dev/null +++ b/src/algorithms/activesearch/alg_activesearch_load_yaml.cpp @@ -0,0 +1,50 @@ +#include "alg_activesearch.h" + +#ifdef USING_YAML + +// yaml-cpp library +#include <yaml-cpp/yaml.h> + +namespace vision_utils +{ + +namespace +{ + +static ParamsBasePtr createParamsACTIVESEARCHAlgorithm(const std::string & _filename_dot_yaml) +{ + AlgorithmParamsACTIVESEARCHPtr params_ptr = std::make_shared<AlgorithmParamsACTIVESEARCH>(); + + using std::string; + using YAML::Node; + Node yaml_params = YAML::LoadFile(_filename_dot_yaml); + if (!yaml_params.IsNull()) + { + Node d_yaml = yaml_params["algorithm"]; + if(d_yaml["type"].as<string>() == "ACTIVESEARCH") + { + params_ptr->img_size_h = d_yaml["img size h"].as<int>(); + params_ptr->img_size_v = d_yaml["img size v"].as<int>(); + params_ptr->n_cells_h = d_yaml["num cells h"].as<int>(); + params_ptr->n_cells_v = d_yaml["num cells v"].as<int>(); + params_ptr->margin = d_yaml["margin"].as<int>(); + params_ptr->separation = d_yaml["separation"].as<int>(); + }else + { + std::cerr << "Bad configuration file. Wrong type " << d_yaml["type"].as<string>() << std::endl; + return nullptr; + } + } + + return params_ptr; +} + +// Register in the SensorFactory +const bool registered_matACTIVESEARCH_params = ParamsFactory::get().registerCreator("ACTIVESEARCH ALG", createParamsACTIVESEARCHAlgorithm); + +} /* namespace [unnamed] */ + +} /* namespace vision_utils */ + +#endif /* IF USING_YAML */ + diff --git a/src/algorithms/algorithm_base.h b/src/algorithms/algorithm_base.h index 13a4af2..28f3c3a 100644 --- a/src/algorithms/algorithm_base.h +++ b/src/algorithms/algorithm_base.h @@ -43,7 +43,9 @@ class AlgorithmBase : public VUBase, public std::enable_shared_from_this<Algorit virtual AlgorithmParamsBasePtr getParams(void) = 0; - virtual void compute(cv::Mat& frame, PointVector& pts_prev, PointVector& pts_matched_in_prev, PointVector& pts_matched_in_frame) = 0; + virtual void setParams(AlgorithmParamsBasePtr _params) {}; + + virtual void compute(cv::Mat& frame, PointVector& pts_prev, PointVector& pts_matched_in_prev, PointVector& pts_matched_in_frame) {}; // Factory method static AlgorithmBasePtr create(const std::string& _unique_name, const ParamsBasePtr _params); diff --git a/src/examples/CMakeLists.txt b/src/examples/CMakeLists.txt index 96db02d..27850e0 100644 --- a/src/examples/CMakeLists.txt +++ b/src/examples/CMakeLists.txt @@ -4,7 +4,8 @@ ADD_EXECUTABLE(test_sensor test_sensor.cpp) ADD_EXECUTABLE(test_detector test_detector.cpp) ADD_EXECUTABLE(test_descriptor test_descriptor.cpp) ADD_EXECUTABLE(test_matcher test_matcher.cpp) -ADD_EXECUTABLE(test_algorithm test_algorithm.cpp) +ADD_EXECUTABLE(test_algorithm_opticalflowpyrlk test_algorithm_opticalflowpyrlk.cpp) +ADD_EXECUTABLE(test_algorithm_activesearch test_algorithm_activesearch.cpp) # link necessary libraries TARGET_LINK_LIBRARIES(test_factories ${PROJECT_NAME}) @@ -12,4 +13,5 @@ TARGET_LINK_LIBRARIES(test_sensor vision_utils) TARGET_LINK_LIBRARIES(test_detector vision_utils) TARGET_LINK_LIBRARIES(test_descriptor vision_utils) TARGET_LINK_LIBRARIES(test_matcher vision_utils) -TARGET_LINK_LIBRARIES(test_algorithm vision_utils) \ No newline at end of file +TARGET_LINK_LIBRARIES(test_algorithm_opticalflowpyrlk vision_utils) +TARGET_LINK_LIBRARIES(test_algorithm_activesearch vision_utils) \ No newline at end of file diff --git a/src/examples/test_algorithm_activesearch.cpp b/src/examples/test_algorithm_activesearch.cpp new file mode 100644 index 0000000..495282f --- /dev/null +++ b/src/examples/test_algorithm_activesearch.cpp @@ -0,0 +1,167 @@ +#include <iostream> +#include <iomanip> +#include <cstdlib> + +#include "../vision_utils.h" + +// Sensors +#include "../sensors/usb_cam/usb_cam.h" + +// Detectors +#include "../detectors/orb/detector_orb.h" +#include "../detectors/fast/detector_fast.h" +#include "../detectors/sift/detector_sift.h" +#include "../detectors/surf/detector_surf.h" +#include "../detectors/brisk/detector_brisk.h" +#include "../detectors/mser/detector_mser.h" +#include "../detectors/gftt/detector_gftt.h" +#include "../detectors/harris/detector_harris.h" +#include "../detectors/sbd/detector_sbd.h" +#include "../detectors/kaze/detector_kaze.h" +#include "../detectors/akaze/detector_akaze.h" +#include "../detectors/agast/detector_agast.h" + +// Algorithms +#include "../algorithms/activesearch/alg_activesearch.h" + +#define MIN_GOOD_POINTS 10 + +int main(void) +{ + using namespace vision_utils; + using std::shared_ptr; + using std::make_shared; + using std::static_pointer_cast; + + // YAML file with parameters + std::string path_yaml_file = "/src/examples/yaml"; + + // Root dir path + std::string vu_root = _VU_ROOT_DIR; + + // Setup camera sensor by default + SensorBasePtr sen_b_ptr = setupSensor("USB_CAM", "CAMERA_test", vu_root + path_yaml_file + "/FAST.yaml"); // Any YAML with sensor type setup correctly + SensorCameraPtr sen_ptr = std::static_pointer_cast<SensorCamera>(sen_b_ptr); + + std::string def_detector = "HARRIS"; + std::cout << std::endl << "Which DETECTOR do you want to test? Type one of the registered names [default: " << def_detector << "]: "; + std::string det_name = readFromUser(def_detector); + + DetectorBasePtr det_ptr = setupDetector(det_name, det_name + " detector", vu_root + path_yaml_file + "/" + det_name + ".yaml"); + + if (det_name.compare("ORB") == 0) + det_ptr = std::static_pointer_cast<DetectorORB>(det_ptr); + else if (det_name.compare("FAST") == 0) + det_ptr = std::static_pointer_cast<DetectorFAST>(det_ptr); + else if (det_name.compare("SIFT") == 0) + det_ptr = std::static_pointer_cast<DetectorSIFT>(det_ptr); + else if (det_name.compare("SURF") == 0) + det_ptr = std::static_pointer_cast<DetectorSURF>(det_ptr); + else if (det_name.compare("BRISK") == 0) + det_ptr = std::static_pointer_cast<DetectorBRISK>(det_ptr); + else if (det_name.compare("MSER") == 0) + det_ptr = std::static_pointer_cast<DetectorMSER>(det_ptr); + else if (det_name.compare("GFTT") == 0) + det_ptr = std::static_pointer_cast<DetectorGFTT>(det_ptr); + else if (det_name.compare("HARRIS") == 0) + det_ptr = std::static_pointer_cast<DetectorHARRIS>(det_ptr); + else if (det_name.compare("SBD") == 0) + det_ptr = std::static_pointer_cast<DetectorSBD>(det_ptr); + else if (det_name.compare("KAZE") == 0) + det_ptr = std::static_pointer_cast<DetectorKAZE>(det_ptr); + else if (det_name.compare("AKAZE") == 0) + det_ptr = std::static_pointer_cast<DetectorAKAZE>(det_ptr); + else if (det_name.compare("AGAST") == 0) + det_ptr = std::static_pointer_cast<DetectorAGAST>(det_ptr); + + std::cout << "\n================ ACTIVE SEARCH TEST =================" << std::endl; + + std::string def_alg = "ACTIVESEARCH"; + + AlgorithmBasePtr alg_ptr = setupAlgorithm(def_alg, def_alg + " matcher", vu_root + path_yaml_file + "/" + def_alg + ".yaml"); + + alg_ptr = std::static_pointer_cast<AlgorithmACTIVESEARCH>(alg_ptr); + + std::cout << std::endl << "... Testing " << det_ptr->getName() << " with " << alg_ptr->getName() << " ..." << std::endl; + + // Open camera + sen_ptr->open(0); + + cv::startWindowThread(); + cv::namedWindow(alg_ptr->getName(), cv::WINDOW_NORMAL); + + // The following line is used to remove the OpenCV "init done" from the terminal + std::cout << "\e[A" << " " << std::endl; + + // get first frame to set parameters correctly + cv::Mat frame; + sen_ptr->getFrame(frame); + AlgorithmParamsACTIVESEARCHPtr params_ptr; + params_ptr = std::static_pointer_cast<AlgorithmParamsACTIVESEARCH>(alg_ptr->getParams()); + params_ptr->img_size_h = frame.cols; + params_ptr->img_size_v = frame.rows; + alg_ptr->setParams(params_ptr); + std::cout << "[" << params_ptr->img_size_h << "x" << params_ptr->img_size_v << "]" << std::endl; + + +// PointVector pts_prev; +// +// bool full_frame_detection = true; + + for (int nframe = 0; nframe < 1000; ++nframe) + { + cv::Mat frame_cur; + + // Get frame + sen_ptr->getFrame(frame_cur); + +// if (full_frame_detection) +// { +// full_frame_detection = false; +// +// // Detect features in current frame (and draw them in current frame). +// KeyPointVector kpts = det_ptr->detect(frame_cur); +// +// // Check new features +// PointVector pts = KPToP(kpts); +// PointVector pts_new; +// PointVector pts_existing; +// whoHasMoved(pts_prev, pts, pts_existing, pts_new); +// +// // Add new features to the buffer +// pts_prev.insert(pts_prev.end(),pts_new.begin(),pts_new.end()); +// std::sort(pts_prev.begin(), pts_prev.end(), LessPoints); +// +// // Colour new feature detections +// for (int ii = 0; ii < pts_new.size(); ++ii) +// cv::circle(frame_cur, pts_new[ii], 5, cv::Scalar(0,255,0), -1); +// // Colour matched features +// for (int ii = 0; ii < pts_existing.size(); ++ii) +// cv::circle(frame_cur, pts_existing[ii], 5, cv::Scalar(128,128,0), -1); +// } +// else +// { +// // Algorithm call +// PointVector pts_matched_in_frame; +// PointVector pts_matched_in_prev; +// alg_ptr->compute(frame_cur, pts_prev, pts_matched_in_prev, pts_matched_in_frame); +// +// // If number of tracked features decreases, rescan and get new features +// if (pts_matched_in_frame.size() < MIN_GOOD_POINTS) +// full_frame_detection = true; +// +// // Draw optical flow +// for(size_t ii=0; ii<pts_matched_in_frame.size(); ii++) +// cv::line(frame_cur, pts_matched_in_prev[ii], pts_matched_in_frame[ii], cv::Scalar(255,0,0), 10); +// +// // update features found in current frame +// pts_prev = pts_matched_in_frame; +// std::sort(pts_prev.begin(), pts_prev.end(), LessPoints); +// } + + cv::imshow(alg_ptr->getName(), frame_cur); + cv::waitKey(1); + } + + cv::destroyAllWindows(); +} diff --git a/src/examples/test_algorithm.cpp b/src/examples/test_algorithm_opticalflowpyrlk.cpp similarity index 92% rename from src/examples/test_algorithm.cpp rename to src/examples/test_algorithm_opticalflowpyrlk.cpp index 37457d2..585066f 100644 --- a/src/examples/test_algorithm.cpp +++ b/src/examples/test_algorithm_opticalflowpyrlk.cpp @@ -77,13 +77,10 @@ int main(void) std::cout << "\n================ OPTICAL FLOW TEST =================" << std::endl; std::string def_alg = "OPTFLOWPYRLK"; - std::cout << std::endl << "Which MATCHER do you want to test? Type one of the registered names [default: " << def_alg << "]: "; - std::string alg_name = readFromUser(def_alg); - AlgorithmBasePtr alg_ptr = setupAlgorithm(alg_name, alg_name + " matcher", vu_root + path_yaml_file + "/" + alg_name + ".yaml"); + AlgorithmBasePtr alg_ptr = setupAlgorithm(def_alg, def_alg + " matcher", vu_root + path_yaml_file + "/" + def_alg + ".yaml"); - if (alg_name.compare("OPTFLOWPYRLK") == 0) - alg_ptr = std::static_pointer_cast<AlgorithmOPTFLOWPYRLK>(alg_ptr); + alg_ptr = std::static_pointer_cast<AlgorithmOPTFLOWPYRLK>(alg_ptr); std::cout << std::endl << "... Testing " << det_ptr->getName() << " with " << alg_ptr->getName() << " ..." << std::endl; diff --git a/src/examples/yaml/ACTIVESEARCH.yaml b/src/examples/yaml/ACTIVESEARCH.yaml new file mode 100644 index 0000000..9b7954f --- /dev/null +++ b/src/examples/yaml/ACTIVESEARCH.yaml @@ -0,0 +1,17 @@ +sensor: + type: "USB_CAM" +detector: + type: "HARRIS" + maxCorners: 1000 + qualityLevel: 0.01 + minDistance: 10 + blockSize: 3 + k: 0.04 +algorithm: + type: "ACTIVESEARCH" + img size h: 640 + img size v: 480 + num cells h: 10 + num cells v: 10 + margin: 0 + separation: 0 \ No newline at end of file -- GitLab