From 90961926d7dfe7e0dfd0f2337d5e7a632d773a66 Mon Sep 17 00:00:00 2001
From: swittl <simon.wittl@th-deg.de>
Date: Mon, 28 Oct 2024 11:31:25 +0100
Subject: [PATCH] first push
---
.gitignore | 4 ++
CyXTraX/__init__.py | 15 +++++
CyXTraX/artist_bridge.py | 92 ++++++++++++++++++++++++++++
CyXTraX/data/cone.aRTist | Bin 0 -> 1918 bytes
CyXTraX/data/cylinder.aRTist | Bin 0 -> 1743 bytes
CyXTraX/data/olaf_6.stl | Bin 0 -> 182484 bytes
CyXTraX/generate_atlas.py | 58 ++++++++++++++++++
CyXTraX/load_maps.py | 23 +++++++
CyXTraX/mapping.py | 115 +++++++++++++++++++++++++++++++++++
CyXTraX/mesh_object.py | 23 +++++++
README.md | 94 ++--------------------------
scripts/01_aRTist_bridge.py | 50 +++++++++++++++
scripts/02_generate_atlas.py | 22 +++++++
scripts/03_load_atlas.py | 23 +++++++
setup.cfg | 18 ++++++
setup.py | 3 +
16 files changed, 452 insertions(+), 88 deletions(-)
create mode 100644 .gitignore
create mode 100644 CyXTraX/__init__.py
create mode 100644 CyXTraX/artist_bridge.py
create mode 100644 CyXTraX/data/cone.aRTist
create mode 100644 CyXTraX/data/cylinder.aRTist
create mode 100644 CyXTraX/data/olaf_6.stl
create mode 100644 CyXTraX/generate_atlas.py
create mode 100644 CyXTraX/load_maps.py
create mode 100644 CyXTraX/mapping.py
create mode 100644 CyXTraX/mesh_object.py
create mode 100644 scripts/01_aRTist_bridge.py
create mode 100644 scripts/02_generate_atlas.py
create mode 100644 scripts/03_load_atlas.py
create mode 100644 setup.cfg
create mode 100644 setup.py
diff --git a/.gitignore b/.gitignore
new file mode 100644
index 0000000..a6b232a
--- /dev/null
+++ b/.gitignore
@@ -0,0 +1,4 @@
+/build/
+/temp/
+__pycache__
+*.egg-info
\ No newline at end of file
diff --git a/CyXTraX/__init__.py b/CyXTraX/__init__.py
new file mode 100644
index 0000000..d45a132
--- /dev/null
+++ b/CyXTraX/__init__.py
@@ -0,0 +1,15 @@
+from .artist_bridge import CylindricalProjection, load_mesh
+from .mesh_object import MeshObject
+from .generate_atlas import record_atlas
+from .load_maps import load_atlas
+
+import importlib.resources
+
+
+with importlib.resources.path('CyXTraX.data', 'olaf_6.stl') as file_path:
+ olaf_stl = file_path
+ import numpy as np
+
+ olaf_mesh = MeshObject(olaf_stl,
+ position_mm=np.array([1, 2, 3]),
+ orientation_quat=np.array([0., 0., 0., 1.]))
\ No newline at end of file
diff --git a/CyXTraX/artist_bridge.py b/CyXTraX/artist_bridge.py
new file mode 100644
index 0000000..a8a9123
--- /dev/null
+++ b/CyXTraX/artist_bridge.py
@@ -0,0 +1,92 @@
+from artistlib import API, utility, SAVEMODES
+from pathlib import Path
+import numpy as np
+import os
+from scipy.spatial.transform import Rotation
+import importlib.resources
+
+from .mesh_object import MeshObject
+
+# !!!!!!!!!!!!!!!!!!!
+# Read carefully
+#!!!!!!!!!!!!!!!!?!!!
+
+
+# This Script assumes that the detector has a size of 6283.19 mm / 2 pi
+# and the curvature is set to the y axis with center on the source.
+
+# The source Position has to be on [0, 0, 0]
+# The detector has to be on Position [1000,0, 0] and orientation [0, 90, 0] / curvature y
+
+# The Script assumes this geometry and only moves all the objects!
+# Preview must be reseted!!!
+
+class CylindricalProjection:
+ def __init__(self, api: API = API()) -> None:
+ self.api = api
+ self.x_resolution_mm = np.pi
+ self.y_resolution_mm = np.pi
+ self.x_px = 2000
+ self.y_px = 2000
+ self.detector_radius_mm = 1000.
+ self.objects = list()
+
+ def translate(self, position: np.ndarray):
+ # for obj in self.objects:
+ x_d, y_d, z_d = 1000, 0, 0
+ x_s, y_s, z_s = 0, 0, 0
+ x_p, y_p, z_p = position[0], position[1], position[2]
+
+ self.api.rc.send(f'::PartList::Invoke {str("S")} SetPosition {str(x_s + x_p)} {str(y_s + y_p)} {str(z_s + z_p)};')
+ self.api.rc.send(f'::PartList::Invoke {str("S")} SetRefPos {str(x_s + x_p)} {str(y_s + y_p)} {str(z_s + z_p)};')
+
+ self.api.rc.send(f'::PartList::Invoke {str("D")} SetRefPos {str(x_s + x_p)} {str(y_s + y_p)} {str(z_s + z_p)};')
+ self.api.rc.send(f'::PartList::Invoke {str("D")} SetPosition {str(x_d + x_p)} {str(y_d + y_p)} {str(z_d + z_p)};')
+ self.api.rc.send(f'::XDetector::SetDownCurvedView;')
+
+
+ def compute_projection(self, position: np.ndarray, temp_file_path: Path = Path('temp.tiff'), output_full_ray_projection: bool = True) -> np.ndarray:
+ self.translate(position)
+ self.api.save_image(temp_file_path, save_projection_geometry=False, save_mode=SAVEMODES.FLOAT_TIFF)
+ image = utility.load_projection(temp_file_path, load_projection_geometry=False)[0]
+ os.remove(temp_file_path)
+
+ if output_full_ray_projection:
+ return self.convert_rays(image)
+ else:
+ return image
+
+ def convert_rays(self, image):
+ inter_image = np.roll(image, 1000, 0)
+ image += np.flip(inter_image, 1)
+ return image
+
+ def set_cylindrical_mode(self):
+ with importlib.resources.path('CyXTraX.data', 'cylinder.aRTist') as file_path:
+ print(f'Load: {file_path}')
+ self.api.load_project(file_path)
+ self.x_px = 2000.
+ self.y_px = 2000.
+ self.x_resolution_mm = np.pi
+ self.y_resolution_mm = np.pi
+ self.detector_radius_mm = 1000.
+
+ def set_cone_mode(self):
+ with importlib.resources.path('CyXTraX.data', 'cone.aRTist') as file_path:
+ print(f'Load: {file_path}')
+ self.api.load_project(file_path)
+
+
+
+
+def load_mesh(cylindrical_projection: CylindricalProjection, object_list: list[MeshObject]) -> bool:
+ cylindrical_projection.set_cylindrical_mode()
+
+ for mesh_object in object_list:
+
+ object_id = cylindrical_projection.api.load_part(mesh_object.object_path)
+ cylindrical_projection.api.translate(object_id, mesh_object.position_mm[0], mesh_object.position_mm[1], mesh_object.position_mm[2])
+ cylindrical_projection.api.rotate_from_quat(object_id, mesh_object.orientation_quat)
+
+ return True
+
diff --git a/CyXTraX/data/cone.aRTist b/CyXTraX/data/cone.aRTist
new file mode 100644
index 0000000000000000000000000000000000000000..7c92c2eaf561c53c85465f8c73db5686732c2ff8
GIT binary patch
literal 1918
zcma)-dpy)xAIE<#qgg}p?6^K^N240!Hg;y*!?>m~qH&An8KY<tja(O_n3O^kiXLfQ
zCn0MXMlO#n(TrgXQMrWKlGu1GE0+m7w%z*U*?qm9?|Gfqd7X1U@7Fo!^T+47<7H%_
z000I9jSjxY8fFxI&wu~`3jqLm006=$r_UNjo(MV{)alKOcx|%lAd#GWA_8?S@t12Q
zI;XN(4;C$Q@bYlnw`gdwD3Dh^(0&gW!m`Tg@Z5R7?zpV*-T4jSMiSnqY2YpOS^7Mk
z=r&eJ(qF31cgBCO5^bc}G(rd=XVy3HZ5*K|#?l9vbXGL6ncElaT+lI27+SX@a-up1
zRhB;#miAdX$Kox|`837^1ut=fgN!uuHfg9H`nj#<sQ$N`<^^<v<6gyjEwVPT6Po_b
z&+!iNqvHrZDf)>LE3dXKRF(RgA9`2CggmmY5!OicI0Z=sztDof@7tawwUjX25&z0l
zX{N8DZ2cxFjAINH;#j`2TJ{Z3mczY#UVi{Yy{oahv=`RTSc=k#uFMTcQ-osu&hS*=
zPF1nMK>AW@_syZWhVEGln{QUVZX{Sf5NLM>-c|niwBUz<noYl^*<r6~xkz02@u4+s
znZ4y*FKam(t%Kwig`RL*R+0@tm(qL$%gp6_mR-#K)uqmn&P_SAk7A`Uo1!drRX#Ku
z_ib6e@X-Z3>t6_7;?9SCbMFnPFz2s&5Nq1h3?4N%BxhCxL<TxJt%wh2I>>9GGm<}q
z(a5@UaR=hqx8-b>^-5TFU$^_ltApBI_FRTkSg~CegKwy~E4Az6eI9DAY1DC6j(Kfq
z3KnPHUCPfdY>Kmpe@IJ<1J!PILUQN(l`^&e8RC{0R->zr(C#O!7G)HbDKJyQvM9_c
zJ!!H{X5HMS$6|14<2Y%2*@Z0ym`)8!rC4IXF(jig`d5X}P}-}|_f^ymP&#W;TJPFE
zv|kCjLca!lr(0v64l^5Ra0zx8dDy;%5W159GL!Lk2pp5+0}m7*eQ%w(i5@B1Snm&R
z#cQcvunw7xR2sf6@E1eI!bV+3*yJQC%C;6<&l~Y7R9#Lax)wvl=U)-5g31eS#;~4M
zilXNC)E64lPm*D6Meu~g^}B5J<CLpc%b>>ph^)Se%VznWF3(cH$t5lJZhf;TYoG=>
z-GuS=aMTIDo-Z4$;PG!eXL%Qg{GXulcOnBEH4HGW$y7=&EFN`=kw(nULMX~5=lEXv
z?<BIfL8L|BY}eWmV>50I9jZ6T!;QjQr<K`?Au?Q{bho{_HY|cDmMU8}w}7`*w5f7I
zbdtC_0zP`zj}_^0d~IfDm;Y8VJV6<eVZ@)Tfm~i8atJM3y_o_Er`b}Jd7J<A(By~3
zW<y$V>Gfz34`<C06?N38(%7{>e|<$6Mems+azUEWfoFtlTkg=C=l?m6sIwTOCYUDr
zKLo#l93s7#OL{V&v0M`P@%<UHr*brMvhwk(pKAty@`;mc>|;*QJ9}QLcVK6kVlE}M
zYV<ngDTP(KN+@c(#Rku<_fIzz#-0S+k_0<ntfcob!>xMNc2FP_>5PxqhJwrBm2!hh
zM{m%u5^rv|ved-{e*r16Y3AYv(c5)n)~wTiku;Nb75|}OcuKUyoJ+nnvAg1Rm+c$D
zrWE&@p?D&InLj;mfA0>d%wE-He82n=tt(y})fdHeMlE?q&A`H_@5BlLJFJdqrk>}3
z_6l`=V(=g#?iQ`i;aQb3&IhPJuD$$$;Y-UB0Cc1E-naEx!_`&%82G{AIN0O?yLOwx
zTQepbSJF3D1AdJH)0=6d2|Z)T!5_Ky=!}lh1<&OPU27%jzMGVy_pT^83#E9qXA8R?
z9d@x-OpVR(FcxK;L+xW|)GNk|yp{OR&EM_5=ZTcl>#2n?V3+LFk!B{ci@DnD+$ZWh
zW3)Qt{#wimQX$-*Rk1%j2Dgybd~jgi$#D+T>QUu#+dyUkmR4S0?6YEm*^jv;rP_L7
zWUtKh7U2FzBm=)kvLuB+mn4}1cS=N*(cepUNF+pBbi|3E&!SJ!9@B+j;U$}HlCJkf
zWQIjK;!tQD$`oyiF*kR@U`^1*_NHi@nLP?^;^biXT|^|&d!e_Fe?)-j?jz{)T)j~k
ztUK3D@tFt8Q|KcMk41V4k0J~Q-b@|!UcmPd8l?56r=)e!91Paf)Yx<ojWuq-7@MLc
zi;}lVNvp~H{f8sDn51ApWAgU+<@@=<{k1FGY5>>;Y~gYwCjH^jpGW;KvhC8R+C~~B
k$Y;;~g#OyKPtf-LM!!J+kB9DfklZ#bD`{(zk@oxSH$@Q~#{d8T
literal 0
HcmV?d00001
diff --git a/CyXTraX/data/cylinder.aRTist b/CyXTraX/data/cylinder.aRTist
new file mode 100644
index 0000000000000000000000000000000000000000..d7ccf000d51c5376729ecdacfa9ee3d1a420b7e7
GIT binary patch
literal 1743
zcmWIWW@Zs#VBp|j(5VcMOv(TBd=Co)11}c?0}wIfCT8a86(^_WrJfD*E>>0-uruD@
zlw7od=ge($Tb=aDTi0$%t)A=MbyZ>J8;-QbjZcGK?q8nxr6hONrtK#tWwY(t|L@o0
z%pF2r5^;}9Ti%ssr#+gdKRap_Yu1(!JtyUdN6us#Pt*>qckB@mT9}v}w==uf;XHfy
z-siX2RAM^aR(NS_|H<%qPwUS7A3ia~>OD_PA6;_ENcwU=L1XKzedbJSm+KdO^X1I`
zXqwOP{9pCyUJi{8x7isw*@^sb6w8?3EsK4$)vQQSbiMSuxDQp$)BW^Vmw8T{u;old
z*U?j)sfr@+*7h6?=b!Jk!{cP?{3ErKj#mavd#rLc$?DAGIzyF0Gr8P`!!6S7*A~}J
zR6AKPNA$3b|6?1ur)zY&l)OL9nCAMg@xX^ie-3gbpMN{MWSyP-OY_j|xSPcufx0TH
z^)j}S|Gr%hx*?u;-h0xLDPPvC{w{ezN1H?Q@nfF2yW91eKNTJE;JV*rC?(DuCU9(0
z&I~?xXR8IZClZ9;y6|86n!RP6O<(o#HPau@j5CwEXOp-8+j{Q^mp$*UzYWrw@NN1r
z4sDrh#b@3xyR<8QuZ};vNYuAS!ZVgLpV|`2&AsJ6E9<HU4bD|90b%o!*5*t2Jjra!
z|I6I1`A$0D>+xAdpVWJX34bk4seQ`+A$iZ~WydaU%Pnk7+m{JBJe2F-6B{13J-#-u
z>bAAv{EY5v>6$D@-2d0+o~z&$ebN4>zwKbOawGru5B~A<{<~!FnQec=XWPnyPx=em
zcU?QsQTFZO^UwEL%o=}bcd)Was3mUMFC_Nzr;QWy-Bny!SqsAN9++&Bl=ma3Wc~E3
zPBy338rqy#KcV4%*zUK#f;-dKbX1$@f0zFw?^dS~_)`4e3+M8r#7}`6KK)mnG5yu&
zyUm;Q@0je%SKM@&{hRKa<XPgbUXs42O<GT^y{vrVndIDy9Ieyb^&h-+^VU9QqtkfH
zX+?HyZl=uR-?~BzKUS>h{a6rdut|DXcJu#<?(8?t=7p~de)dXmsfx)n>0)d7-)nSL
z*ni&ndGupkqkhz<tDLSZT=GRf9~C#wHGa%Dt4I7vRoU_T_L?tiw%<F*vqO~Q7iY(=
zqKB>Djzr2RnpR%))9#*e!dCzFdC@}KFClDe*R0;=$m8%m%jZFdz*W6X-v#tqzR3MP
z>Gx={!|t=P5}&RvUMCx0lEmP@K<<72vT$APh<V>uH*D~7Kk<F>wa!!SGo;qN35vSk
zURz|XFze9iDYIJ(eg~hqZSNWSa{J)`p(EmMD`viUpSAh4^YvVf&kNp9yK?X%_s*^L
z422i3KCtdt;i_v_x^_+Xmh+3^+g%=hXp7K`xbx%GzE|%U>J%58Jhak|iQl3sy0QK_
zSMUv$m#%YSHwf|TFP|Ye>%FN}BJ%^`zt46*f5yP}MEQ&K`seeX?pb~M!3jT`?@WJO
zUi7aPp2GKXvC7@O?|bKcYxMtXc5$&t{yw&2Qf7yb+^_PR{Uvk1Cs)*$u9jt)Z1o=w
zruRIp<V&{n5&U6vrf3J-7msP_s-4_kwUa+4%FSZ%<?K9Ytj4)m^`!Uxd!P242<GPs
zUpysElP7gcil<<@XV^5)X_wA^+xPfnTUcWHhN69UHO0lYSU)I`XuCaquH?;%Tl-X{
z4*vDKdvk+a>fSS#o2)%G3>p5j2SCd>25=b%EY{$F703!mEGp5*t%Mz@BsjmcC^;3b
z2vSy-9n;w84V3f)Vr_&X3sVDEM+0L=19M|@6H7}s6H_x|LuYeiM+;{IV>34wy^5mZ
z;II#8&!6@`r*Z1Mug*!IFas0QfD?Yamjew#E`(nwDAx+P5Tv1Z>B-Zeun*p6e77=b
zKYz-kX}sOU)ZE<A{ED%u;Z_qvb7P<&Lx49UlQ=W(5*+9-AW&e~)(D~z<vBth_D}(-
zW?&Ft`0w}xSRX*c2xb+g7EmZb07wfvFkInb2GfK+)F7Hb?$beNA|?z2yjj^mnpuGG
LCy+h>ic1Cnk&Dp5
literal 0
HcmV?d00001
diff --git a/CyXTraX/data/olaf_6.stl b/CyXTraX/data/olaf_6.stl
new file mode 100644
index 0000000000000000000000000000000000000000..908979f538d74e1d0809e0c0c5bdbc77ba3246a1
GIT binary patch
literal 182484
zcmb@vceq_u)&0LHB1qBbYotn15Gm3ULgb#k@3}~jsDK~|y@-eq0=ywA0s%A>F@g|5
zBs7(xAY7y)$vNkq3xp1dbP|yw0i_cH>MQX3jJ3-hdz1Y0eSSXAeLUpMamHMG&aviP
zyR03;|G)q1T_uE7{&HdK8*41s_T9#bouOYJF}Tt5^9J<a7d*ek;31=)7`S}J-=r7v
zdB~>vX~Yt&Mu7M@9MhIsF!SFvj_+l^+qsC}b-wB^5B*owEC2Tpel_8|*8bb{&$^?=
z@${to2i|^W-GLc!INuPu`f0=xtTc`vygjY8;&Q`U7Z}F}A8a*v^oq9*eC~>)2F{1W
z`9gU9%~{EpUz=!0HUuH6R)q)R*uOow_54F)TBp=F<{x*#z&!rM;d}+Q`u}yWif{Gr
z#pmk(>zfT2!z9L{m>1t>hy@RiY3+5}<mi=KxmQ($My=5o|C;^(h{1P=g+~8jjy0AU
z)|$x?zKD^#mGc$9mU%V&3;na+JN&#>J?2-e_?ynb<d!)6mkM<K{}=C72#UuiMlP5s
zjnws_-}oXL5r^~Df4lha5lgTV$Jduyu<bB7>iYOSnp4h7AI?|&X8r&5TSHzAeXD=}
z`8NGM=b(}Lu|gZY^N5G^QRf@_g`rO)mS820UoA7NvHhZH6@BO^qK#hVl{lPl=vS+s
zMl8Wf9H%@mrtwNA%funnh(WZ`D}5*)&NpOJ{Z#d8=sR(|+?w1t9Sxpl8q_!#L+|I&
zh&Y^Y2wnYD^=jxlas2n~X^roZ@yd4p=5Lek?^Vu<!}*46s-LP}4SgpLw9)t|8mY$|
zZS*pqcf{d*L+I+K5lgTV$DV(=u<@TY7R;>c18wx`R}^#StN!xPe~nmzl{nBw<E|H9
z*w*!-831kcDrcn+=PQ2oq?Au575^o~;;b^AFK(vRH-F9I(W_(SzZ5G@<4ari_sk}O
zZslH85#J}~t1LIH;q{})E9I<?mH$$KuK)kyy~@UMKznlI?+=Y>c>T~fdA>#t<yFB~
z|Lx+xv%LDSmE!T@+%b*)W=(F?^?^2|^LVVBulUXS|5v>l`c@ohqalsd^?^2em9vV6
z^9`Y^pGGXfN*rjTM;fW?L%oi2RvgY({pF$miuI%P%KtqCw6U!;QrCy9l#PnH^9`Y^
zpGGXfO5;EqGo_KbKF~%-*=Y46^68gP(yVimJ{AAVYU+?z#lL5HMLkcKRL_;OAJzxj
zm_=9CkUlh{6;^ex#FKH~xBRfyD+{N6NFQY3z)|FtdLRF#ArGscvU*<d4SlC$_4cwO
zQ~eG7YY0oBkGzwK97Ep;F>213)?ZPT6y;X#RTc3Nozhnx`LG7nw<%|Jto)avR~@{v
zN}yj{yjSuzgnOy8uZ0iKC3$BV<LrDx9#%h%h}~r2)zEk1SfAeNR%W=mK72jkd_(B!
zrx8WJTV`MOy*NI)-^BK&J^ekmy?ac@*8p=K->)mHZ+~&)!3k^c-Q91?^?DEFa|lO)
zSa+-bo*Qcjw{j2ZuRwfbbpNbdK&YQm?UN>Q#DF-QFN769oVH2-teG{0Te-I(EVJdo
z?aB8|ZN0hZn2vY%2W`1?_x+VG8#ria9FZ@CcegmWJ$KI3)`A+sy>f3uxb2bM+B;8~
z*E;j<V>-SDxc%fWcej0R;ed{h^M!CLh(kb}QA4<udmF-oTaRu}I(m5f{)NXB*&f#%
z3#PBt{X2gehx3In6~yNz4{zU7L%5ZD8^Rh`pAJe0@9rCnVUV$i!}&sZ6~xXXN4FoU
zA>69E7CHa8d9D2(-K{;1k!#M-D5bmih{d;wiDK@2A$<Axd96=_(6MraTe-I(r1wv4
ztxMcrBZ9sLIAF{51`Z`>#o>G*d;!Gzi$S=RbBC}ddAkzXD1Ar|)QmU|eSQ#)h{O3p
z7!G225Pz;A+{(QTVW-jkJ-dLAK70)TqPG}_^M$bUCjC7tf{;EOQMGb!d+(px_<)Rm
zgJb1u0JPD2Kpu}O->Y+`Hlz{h!x3&(<=*?ad5sZ7Q2J2kLN79)`bq?Qi@b8a5Y{_>
zUSk~)Kd&L&%DF@M-^s(11IXRG$wXfRppB;Vk>^#(7sB_C9-iC_LW**Pd*$ASu*cS;
zlS6Vv`5J&RG^KOt!}&tE7{o{rcfsKZw{mXIGd#LmatRudK70*;HkvUT;cz~>a1cA9
z5$VGbZspt|pp8Ttkv@D4Ku$H4jXGA&7sAFMzPK2KTRC?K|J-+CvJ4qNi%j(H9&I$0
zjpA^=5H{IwV)D7I`e#WWj&LjYHUzYhNF#auDAyc0qE{Tw#|jd!URl4trwxa%2i(fJ
zskMy!#QUZ;qz_*MXynR9#oYNqSPR6y4^3^nR71FxbJJx#vRkqg883Y(Qp!5%LlIOq
zio^LrcpSv3Q|2|4C5~___cny>$;5roh~{=*1LT?5)Uk5D5Do{i5(v}BhrDv`5YR><
zjflh70BA$HYnnbHANl~X42XMc2=~gp4FPR5q>;RSl-)hr=+)>o4(AI2Z8W5jynd7j
zw{mY|{b)!d(uc1B@)<62I3M#Vv#2yuN4ORBqvWP9A#bG-=|g%@)+y%ttZdXcoG%2l
z(Ibsm{fG#+a&Ke(m?e$m^`l&KNaM;zaX4QHXk(T%B7M}?1BF*Y40$_+-s*FVT>4Oj
zEV}y-vU(n4Ugbj@EonshaD-b`xu0b|e430ui){DqKJNe`hx3In?}T}+vq30J9N|{Z
z9l}?cr8gxLt-p$E4l+^wRXkSC$6P`;wkwDS$X#DGxRrCqYS5+{B(p*I8X)h;;(X<N
zAxxsTTJHP9BZnj0%DI_+sG6#t%64A^sI!dONX(rNuRsh(RX?jC+{(G>^y!ptW=?#W
z2zp2Foxbyha2)-@%MVU%yj4TEm3vDa{Ppbi^ub|`HSV6!@$UYS6@MSsA5TwOd+?>^
zW!>X9J{nc6O&)IzVyiocH9l8ExRrYt!fqgD4i0Or1;V@g7qj0T(C<252>XF}?Dk<T
zdF2SVa&JSprqgLZcG~>b?zd0qcz3_kH!tW;dwUW6@7jaT7s6#Ark*;#wNDM<R?Z#5
zmquOJo_*1%_O+b}9q;bfn>($0_S~(z8i(_RumOmBFCNvttcGwa_cny@yg8$N)@{eN
zFKbQccz3^f<M8fg{E5T)LYM;LsXrgr{$mZ{R_-n5!ukD!L5U+*QCxFOXAJ8y7I8RV
zhBySor8R_Gxwj##%&0eF^w%(Q^iZxj&L9>Vy*QjNggG49oj_b!L%5ZD8^WK7-i<_Y
zPa>GlYF*8>Cv1H5zzgJ-V(xq)Y)*z;0%GqP!mZp}t`FWhto3!W@v~%NK8HsReO5MV
z9L~qTlDAUTYBhvgxwj$o-#)A%j<rGXO}XX((W`6}hx0{_#y4`W98t9zS|PF8K5c&E
zB1XOkBlqqeZS;Pfc<5MF`4SMn1|fYo!mX;@&s;Pr`8pB&0}&)r<r?4@bGL35dF6Z|
z%mVQZ5LeX@Zspt|9DUny$=+n*6=b67g=T2Q9BniUeK=nT^Fizc;&(NKTe-I(9Qx*r
z<fmvv+3wvvV`zek@o+xY&LD0B@%tLWt=!uXHXL<b@+mZOH6vF)6|;_UppE9w(THO1
ze01R;9z!EKx{h!w_cjEykw_!bN8T03H3!;gGM~pe(fLBS5X7cvL@{@STRC@}sgpOz
z_|?f>@9rO7@%M54p}cawSaFg^^SRv-Zspt|eDl|{lT)Zs(qR5gtl?;*sch7-a=y$f
zbQr6eBizcl(I6vV{-^UB-(ut%VX@{w8%<?nUO!5{5Pk>Zuv6zZbgUfVUO6}C;fSC#
zq6qpLAdg^E<8VH@B@pXhJStKBaD-bqcL?7i6R#oTHMi#(7kjEa6Pp@`^U=M6kVZ@&
z5#d(O9l|YG?*&Tw@a{fWR2+x%g>Vmu9YI`PL%5ZD8v@!$q!H;O_EY(q18wvw8+ELl
zF9fuaNFzGBj&LjYHr9`ZG?Ld3&AV~r>UETjI#$jX0@`RuBZ_%G_eF$TId=$Xqd_;5
zNFSP=y}MU7io^LrKpPF}N21=x5pLz&Tob0Zl18KtzA4unXk$RxsAHvHD)_P*C9lNc
zh^m!)`_Z4yZ_#5V(ua5VnyZwJ8b_7y5M~}}M90bzZdK);$LvNmh!M;yeOz;>w^BB0
z9L~r3h#7Q85SJ5kN4S-9hwu|-={?9q&D5$L%DTLMbT_4jh{O3pcoM|d^K5s7TRAu9
zuc$WgL=Nxn)sZP1#X;tmd?Ad>5r4odN4S-9vr46Eu8Ky~UwL=G3$<L?C=Tc2Sb_N1
zVi0cS+#yV(`mavL??(i^qxVkV`9k;uh%=ZI)n7Tnt=!uXUb=RL&X?91+xXJa<2&Bn
zFW)mSu0I~RbCtnWe))L!zh__4`xid5qW{whouk$o+gQ7Xa4YvPgxBGi@Y%7gjX<cM
zQti|6`2@fF!P%D#h{O3fs|Mm%YmIGfUPHK*dmF+IuPxj8+;gk8A3AJ&hxw^obNuk}
z`?~))c-`*9a5x`xJBY17JWxZpm3te)AFp53?oL0Uz1w%jcf7kl^Tb=a+rEBeSL1NL
z3~?TaaW#Znxwp*i4?Wi2^Thu4?FWqScz1v8M=t1I!=J|Cd?CyP@ih>4)evsw-coU}
zz7y0S2=DGS2BE~^d>P`7#UR|uy~Qg=eJ983YmD5x`+pLP$;3pt<a{AK2jUtK->f0r
z%DqK{MDJ@v@qQwxS&f>})jX)2TAAF^IGhi!o?W$l5QxDV!mZrf5Y}FMZ0l{Z@e5?4
zX6Itfq0b*dBZ|56h42iBS3xM-9pP5)E#2#9$2Jy%_#z0tDc2kzdZj9HI3MR3KrG6=
zazxe2y(QBVW&DGT+`Idi4qmsZc<5MF`No473*z1y!mX;@7fnAQxt0h@A9>a9O6Nqd
zsl3uSoG*l5g6IeFts26u+}jWiI<Y^wDbK{bdx;!qqgm*~`9kOgaT^Hf!x3)f-qI~S
z^jPu@G|~l;cVUs^j1zBZZi7a2telUl;UK;aLW**PTRC?KS6sg+SxCn3!pOb5|H0$;
zHH-Pm`QluXyb5CH8p5rd8`ZwHZ2EUHUNf9`_h_T}J{r+DoG*mELEHr5t{TFv+*{71
zUb{m2X)=Cea@V_iv=P@Iin;UAk%71^uLh2AEB7`8w2?|9(ua5VXrrlY6o>P1y%xlF
z)Tm8s2)A<X5Xkg&e==U#?%h4wXet{uUpZffkVfhVw{mZ}eooAfC*!4$JW{b2K^skF
zqsHNUbjTns2BA685pLz)avqLM+?$M_J(TUSH_S7!shB%ohS(j1=_4ZC%DF>88;LX`
zeRy||Hk!&t(?|Rq!gvsW0wE4ZxL3J8A~!K-)Y6Fbp?*qPr?~`eG?k6wa6V!Vf(o5T
zAC7P<=MDjFB+`iV;oUvj=v6j~!}&r$8;LX`eK^9coI3=x(U3-@56#X)Yk*ikV$7Y7
zEFsgSk)nPS2)A<XSU+0QNM1k6H3!-lP&O*B^2(5ZMjNV8sq~>*q&(Hpty;Oa`(?E@
zRkkz2a?PQ+O4%rVRQYxXAyw58ZdK(zk=acdzbg^+?jCKZFNwZ8UkK;t2<gKSZsp#F
zun)8J^<?55nLf%jK!4U>#h5!EXAMBml@<CZ5pLz&v3|6r5$VIbd&baJHY%^2FN7BL
zX%`Uc`W)d_=%eHg;VP=;yJWnw-Mjm|vy2?h7sBsB{2PQ6<p{TO?hp>3`aeL%-$Bg1
zqbH}z^?>uSjso#75O>!QZspvpoiE#@v;2`$8p%FA9q;Z>zj4_t8%Hnf9x!F;!9|;P
zdza;N2tNWb;;1Q&9cu`;at~Q^f>`c|DXr~5cy}+vzc%d-EYI&cUkImyI2S}(L%5ZD
z8^ZEeZO~b!KW*Q+OHaqU`_pE>(>>}t+jT#ZkE7t@9xo8fEe7FUxwqW=vCQh7Ykzr3
z`}rMuI*OVyP};tA{inNI^_<=PC>+km*;){{f_S!ua4YvVgxQxb*SVnkRQtWYo{o3-
zhh06hdnkVzhx2jwCx~evmfEh0a4YvVgr8#l0H_b($W_#pCNzed8H+fak9dHX1LD6m
zgj>0{A>7BP_ha<0GIH<kClL#cUdPJ$SUZFG5s0U22)A-?>0Yl%+n*$gcMx;UYSfH4
z4#jwXa!dK@d?CC8;@?+p(0QbWa4YwgCwPvW(mIxG+<{EYXXhAmeO5MV9L^WQMj)nw
zP_{e5t=!uXjz3~bV+x3`g5aBS%>klU*(iNDUkKZQI39#F=!mM-&<ZJpkM*a?G90UW
z7`b=%Xrp;5@ledGe8ZQ3aH}f!-M_pfnM?%d6G0+Xt~rQcv&bvw3*jOVlR&&sL%5ZD
z%YEV9r;-QB#E+8g-rXO5^~|QSG0&@#k2^U)+yml2<gO#!EB7{pxtA}OP9Wo#hQqsi
z#?VwY>R35n2zP_{K8WRO2)A<X5Wc(2>ggOZegPx*?jCJ4HD75Q&KJT25FFt+`#8d_
z+*|J0x@v>;Tr&Q4V(#5N+UR=+jYuEP$KALfIw0<=A>7Kj$=%B~Nna%6mAl^Eqm4dg
zqsHNU+%F138fn%LZsp!`g%^Z0A`b8F(MF%LF|QvbU*?rGB7Hc*y>f3uc>l@`(q+l`
z`x&`+_h>`$Xe#E;$6ezbt5G2Ct|8pYxkEr3ieMsrFe8_1jy!^q!}+)`1H{ik{G*0&
zEB7{p`^m%~lJUx2@9xn?Q`x9-I3L{-2x-Lh5xW<+a&JS>@1wy~9NyjMii#Z07s9b1
z<k>Pcgj>0{Tn%T`(unlo-MwO=Y&4xmJ~RkI8c}t2gnQ-QhJZE_X+-+)?jCLQDjTH_
z=L-RCB+`iV;Rv^KZl0B)w~|JbyV8TQu2>H!8#NB+i}j-+jYuDka4YAI^`j+?<n^Om
zbD)g@Wuw)P$QSEJOJ3#mqeO_W@XERIn@%Y(@&_2XclV6Et8C2chw+UBA&p2Mj&QG3
z2c?nhh}^;KM#q%q^`on-(>Eisz9e!uUkJa*5z<HA8$^U#Id|5{C=&yj=-s_)hq6)Q
za6X=3V_v@x#Q#x09N|{Zo%L|1<Qi1xRjfJkjx2IGAK6ZQIthe2GDo<TbH}-)a~m1|
zH%9K=J=*9ViAI#W&c{`K>g;R~f3G3j%DF>Wj(*``GG09cQ7YFQ<W#vHaJ~@Kja>?&
zTSK^&dmF;%j~m<h!$s#b7M{LFhs-P29ACTnGx76Dv%1}vHvRPAuqWsD>T?LIgBW)4
zIgS6;5N_oja-JH*?=L#1B@TU~43s8>_|%j02lTtn$9>@-u38Mjt=wC>*Q2-VoORCu
z?N^UqqodI(1C>h$e|+V^s}I<}Ta2UN<9s5BnIK-MA>1qXHiTnljOg6E&dm0mlh){X
zcmMV~?{&Yu!WG>U7>Dy^h)X~`QbV|vd&`{z&#%>)-3*<+iEDI}eZ`vN(lZ|G>Q5Zb
z7s9?E=7Q*{A>7KnWo>}<&Y*Szp?-==6OYw*7{emQ@?$uhFN8Zm8~|dw8p5sITh`}{
z`dE(F1B~3e`)i5CcZi8NoR7I3#9<(wsv+FUy+u((?;4``Pa>$CqGrT3hhn^d+)~V)
zFNEDe+yFu`cZ6HHw;?Pb3pXMg|3U8NvvXW?=<`@KA`a(c{=WE})><Imts&gXxp`s)
zHEjxF0SLY+*Bl^vH4eqx`M4()#mcLH)(};zp%oJ60PZ;;`BOe}v{9}(&_?qh;*pQT
z_?`f9BZwzBR*rD5+*_{Qtur%`M(!tqdDV}3g*KWC$t#V+`9ipP2?)1x?hpo=Aw7gl
zOvyy=?#aZavQc^EeB9#&;(K}SI>N2oTka5geyy}i#_s?}-i5_+Fovddu4CnV+(Qon
zo=#FjxRrC$Gt3x~-blvZ&9PzzDAydsqEGXc#^HRN{QxltVr~uLR_<+_;nLH{c<IBt
zdvdCe`8>|M&PP0s-Y(q}#4|O7Te-J%OUI2(UqU0wM0GC8K=o8;BYxJgaz37@0`W66
zB7Hc*t=wDg902hI87~fI<Z{jNBI8gtDzBU`^GX`2BizclxeCC@FCgPJGkABOkG!dI
zI9~{;DV+=AuQh~Qxwo+zq|%7=;oW^6!KTLHd>KL-sUzIVz2*L0GO<m@Yo=yKF4r7+
zCN?dvB3}pxgOIB7ic=!oEB7`8w2?|9(nqeGrfO}jsCcZLkGUO$c=Gz8vjY*~R_-nL
z4l-(KMEcMSpsZ6~<ugN{={)j<fHo3oMEY=qd*$ASfHo4kfK>YM?q0o)vQaU2K6YCG
zA&p2Mj&LjI=KhjwrcRZK(uT55<Irbiqd1%|1hmnRMy!5Bgj+dx2xy}vjfg{)MOml*
zN)<=hX!RrV#rmPD_929OrO!j&p3AJz&&Xe7<czRfbEtAE8#Rt9AKGY3BhrT>+{(QT
z;aX-l>PPwz5%lg}bFQ*c9L^WQbs&BS;=vlit(-gNZb!L0Ceuf`=C~~DuOf%@g>XCb
zx)8adN`zavw;>!zwE-0s<=wqHvID7w8i(_7-5CTuR{B-Ea)eu{zmi6zv=AmzXJ?Y}
zPcU+Ikj0uK?<}LvoiA3L&L2SBS3|g!b92Rn>aV`xbt0%vxVQ%4oxbz&Oaq8tfOxZp
za4YwgTJg?-o$a5xyz%~uueQCrU-<g^as9FSiPO5>jXpQH?U`@(>N8gqL7e{N<&8&c
z2)A+%IS&Wo8&6%{62iOt95JBZb-oZ*1aab%m$#N$sfutb_m-zjSKF&|%%jJ)S6}!_
z+q?U{cU*n&tX(H}R~s5f<YPB05NCn-WDVh7xwj#lb@Gm#ZSJ4bUgnus+TPut^8CjJ
zpPg}2SL1NL3^5kOTYsw}+{(S>delance<;8yz}!%UuiQtmurq&{`k-CP5g<&`7*>^
zAbwm!xRra$6%(w_gy#?t-rZ{qLg`pJANPTRI2MGwa)evCw;^oFsAq7zUS#Cn-A^YL
zdlD0k!}-|lAxHe9hHxwA))Q67wvQl+8xldyYSfHDX}fXen*)lXIGm3a3y5PutXxC5
zm3zyX39@h{vhgi)H=ms&hd#fFM#SNK%y1x90r5r+;a2W#2oIqqW&9!#d{eFgdeB@k
z4(DT41Y%_n3u}m~)zAtlgi9YiHu+yhz78YLt7RO=*}G0|7V)U^9S7nl5L?y|Zsp$c
zwA1}_l3x<@<<Uo8_2Zg@2sV{hI#$jX!ss0FA7btZw{mYeTf6$l)6>YrA7$BIt~t;~
zv(Sh0(ffcn9mHuhgj+dx2-7#(Je3Be5AW_ze*R;9Pooj(!};hLK->Z1hc$#-<-Jy1
zS-^4fj%iB9FJ$C-HyOVn7JZtpbgY~&g!4dfyxP)-Bizcp4FPSWza!%}CW7AGlT&>}
zG0ua|$5}NHKLN2?4dGVK&2wJw9GL!vj8`UlcaJvW`a|Pzz6>!4Li%uoTe-J9g9XR!
zWc;fjyt@|yRrRg`hx26|mx7R2j&LjI=2=Td{zo!i`p8EZ`vi?#*;tID;L8x_gV5~b
z2=~glxn4^Ir4i}FyL+_Jr)(65^JNHWMEY=qTRAs9HJK=lm_A}Jl4oM%a6WoW5YmY0
zBO=_&xtZaxmPVwHTsiSrp^c`p(ex4dLO>g-G$MUC!oBjfh@K6flSm`dhj;gAL-SRi
zj+OIebBQz}eK^9coSP?jh@LbeeRy}D&v22$`8bOMLK=}i9N|{Z9Rk|WOr0tdH9Hq;
zj(nz$9L~o&Do0rTh~2$gId|k}eH@N|!lBAitU2<E6FHoZYJmO<uN+aea&IqZR+z-d
z*J9+}-J^||vx-NRZ#sx0L42-;aH}f!hRkl~62Yai`cbYq(1!YwI1c9vLEV@Ti-<X~
zEfH?z-g0!ArGH8$o=hfscb|7MP0c0F$ND2joK!=&m2=beQEjBC$soMDXAFbNMjb2X
z<Bk*%c$QABA>7Kn4dHC+><(z;ZAR|heco9XIxqRSRs>>45bx9w?v;DXbxx}P(PaEE
zBIq4GGhDeIaK2bSI?{;r;Rv^KZ&{!J@1)L`7hcy`(06g$yZa5-+bpg>zB_h<?o}f<
z8QlKO_j~o3J5oTL58{Iw!mZpx?z053<-+S)Lge4X+Lt2+4&is5kE@m-E(Rg59N|{(
zE!Vu>KeY4InWwg&+3ezWKHAtPZ2#DLgA-pnseACyI3gc&B8WRcyjw%KSI!;6vs>)Z
zK^*PP*S)x{s3`-LOaDFhlY{&IvD4K!oR4P#K<ofw!y3Y^+*>-!_j)?7oU&SH?bR-B
zdv`zVns>YU6NmF<h}S@@RztXzd&|=`Sc^x9Tt#usu|H$@KgJ>s=i}a55HEm`SB`Ki
z_m<~z81++(etSml-TjfoLZeqMIbR4{gZMj$jcW+Ea&Df~BYJ~GaS;*BXSG4aV~022
zA5auE4(DUP3=jh#KBytw%DrWa!ndw#y-GGdK_=$2bL7zH_0fpN;d~*iO5VNzLfP&J
zw{mZJD&?)~8q(qGAo3~`Ii$HO^Km#I=aA9gD<B@PA*xnGE2I$qaptLsGJYN-S3gA!
zj~s`;c2aXY;vs!h`Tmw8UZ^45%Dv^8>ZP7ZUL=B>5<wzWt~t;~Q`x9v<$NKG%n|F<
z5N_q(a^~)o)ly~R#$=**_he$T(1-Jd@D_;YK&)CrxRra$eDz*W`VtvG5`=g6XroWr
znCqkD<LH8T0mMc%gnQ-Ql8IaFk#0}Mk74BMr;0TP+UPq1jp$f8UkHB%@f8qKlq1~A
zy^V8XdJh@@A`xUJE7u(4R3A}{^Pux_b?N;>)4zgvzlLxt=Z^IweGiSiNG5uBzrlK&
z#r21dmGi~9JzWEhyir5Am3x~Z{)<Ll2BB|IPh4~4h+c6xU*^>?H1bjn;a2W#tOltx
zB7JCt$~w&@XroWrD1A6zhIkglb2Wrpxwnji2udT;2QzZH=0F>L%0|W9`7(qwB7Hc*
zt=!uX&_*ham_Fh-@=R>%SUDfh-GPusOdk>9R?f}Q#abGXKD@g}8%<@S=_B&t6$oiW
z`f!AM<!cf4@nY1{i1gvzJ=$n08+ELlF9fuaNF&mRBizclIhR57=o?b$!@K)@hKu8H
zJ~RkI8j(I6;a1Mg9@yF3o=P8iHRw}b>9ev?$IAKGosztz|4F3}N4S-9vzNE@2jVpl
z-rY;%%0_WGU*u>}aZ>5S5ml>LKXPvuFe^|$()o<syL<Ik%Er8Y7$5bceLIM}ew3>e
z_o~XhKC|0PL~tY#%qx9l*PN?tRLq?(){piA5bM_vZsp#x(=D^KGJY7D=-s_)$H%E5
zd0v%#JmHr0KIwWjgnQ-QhVV4i<_$7_D-i0uiZusg7*saq`Y8F}03lU<yoPYE+*@iC
zb#@yvemzF+-96eEEaofc3qjr3&LH|~2)A-?S*y`6JW0l@zhWjU*BsvIJ0H(Mf}p2P
zo~j|-%DqK{uN>Fe|D{_Sn|0P~dv|~G@!6XqK6dDy1BXx9eDKegF5O&#&zwC6@naC%
z)evsw9&%@9#?k5n;oW_X=<0W!FGHO8(ygt{YY4Y;Z|ULIJGyh|oN4X1u3EFL(H3it
z?>x7`;6v-3+_iB;z6^0Bh!<)I_sYFx_lw^h&^fd7QhSFB)@&<k)Qs4R-1z0;gKr&i
zNBmvqW3NdN=YrU}hHxwQHiR?w{Yqz_scUy8O<%LE?9-p}b(u%s>n_Eg#^HRNNd>VR
zh~L)`Zsp$cln~ZugE|p}R7A}f6vslwa6MyDJe-eRRY05r;)oi;t=wDgnq<@({mP8o
zyZbrB;vizO5**H_{bXKhZv<l98p5sITb?B&dR^lFJ0j@a{cV>n-93z))i|7wd+0&T
z0&#B*;a2V~XPwEyJILSflZl$0i#3NnPemi*a6ZnLym)Kteh~j|R}pUI-g3_!YLZ5z
zkG#snct~>-@^LsHJ9)l%Yh!lql_RQF?(O&HOiONK<g5L0&9-;<XrrlkD6gt~C*}z0
z!x3&(<vywNQnEP_+?5D=cTWVH%0_WGANPfW*b>D4HH2F^H_v}eT{}IEOiZ(EA6f%E
z`d%}j55>dz*nJPgFG1W`L%5ZD8#S2jOU6$I;oUvj=qq#%SINh_Aol%ADpfUV2)9BX
zL*8mGNiQVhr4MF+a?L?3`il9g!WYm+dI^ZtY6!P-Z@Keqy`xjjiHcy}6~{FPIn_rL
z<FRr+t`&jU55zq+gj>0{v3{iMp^-0<yWZXZ<oK=P`a>~yzF0rfJ`lY%gj>0{B}bzP
z!n^w%(NxTxFV+unG~|^d+{(S>sb5AejYuEf-RC23+BhO#hLA>dbRFSdIXAoS5kYB0
z`d~&b*Bof0PuVCA=i@3e2x&z6aD-bqH`hPNL}|qIk*xu;e2u<4AJ<qwNF%0?h;S?C
z=8j*isZ*)+p?<1ZbD)i;vQhEKGe198Xd{(Iqz^~9Rj!ZBTXtf5DUn8`5AW{tnW0a|
zs=^o0Mk0+!AC7RV3U@#oi8K;Bfb3jP?A7xbE{?<bn5jWXBhp7c_eF$TIX8EblZEtH
zsWMTsbFt>oXJw<t;d~(|Z|Sj8>BA9j<=l~@B@XF>Z^|`*G_Gun@i4x~(W2r+uN>iC
zIXBO%F)PqxrP7CY_j%<k#$kL@K>RsBx{h$Ks@$3n74we}K~;KXoz??rLw!l)a6ax&
z1R;$`AC7P<=VoszX6av$iN}$N>R5_3hdLQ$qvGLwG21(5g7|I?;a1MgZlP40{mA&U
zKrn`K&7m5oY*aj)F9d$Svo{EJWR7qv_cnyhsG2_~<JDQJgDln@>MS3mri;V**u53R
z86Za05N_q%S%*x=l)Ol`dq?k`zVnqmb@FBn;a2V~XINf6p>x@qGaC;LA3oB%d-Ybz
zz>N>xv$x>WXU?~Scpiklb%a|LpR>#t$E9n|Y~8=>KG)~pboGrst1Hv*I$y@|EQkm3
z$`Nkm++43+{s*1EAA3f-b>f#tdUwxEJ@^1U?-%oN6ntn9#KRzda^40vIKsVhZmRQ6
z2X$tj{IX_+lk%v=aTK4$;e6C65RZe{dBd-bEUYwgw{mXofxKZ%r+MhQo%27QjPUM0
z&%`*6;+Mq7J&+($5Q>K*+^YC2Zmv{geKV-(aO5hAeKgvL>ks{|^JR#eK>XxihgJ}7
z<=m{#8TBoU9;Z(7?!Ne}m^&X=mvY2TSAJ|nVWm}yTe-K~txojD5XE^<eqy9|_xgM&
zxh0O`m$KZw1H@NA=vX<zt%}d$X8)_TX10{^%0!`xHHSV^KO%?oEg9if&dqxo)|%PK
zIDR!UuQD+n`do~|`8ZbSukgwdRV(-Q&STF=H1el7R(Z7?)Lg62-$o-kK2^TIfRIL>
zeDbUc!mX;@cc1)n(kA9tKm3uAM5<hK=<^NCy_>+{d~}u|q^fIb2)A-?Im>zIy6M$q
zqDJc7eeqej>wLTy0K}CbCjILlBMK|&z^&Zd5cvIcdo=R7lRkfwclV4zePW!iieD06
z)L>d@&=D1@%8F&?P6wrP$@s63?RhsD&n@UP^&@gr_yTH5rIG(m+3+SuxK)Ka5byME
zGX9{WmK*8aJ!|tmWuRj2d~r@pvpMmKtsLQ2&K>JV`W!R0!r|S0@mU<smmw|(G4=9K
zjx4N{iEicGu^Ob($b-XwRg5FnSG1vZR8t(qFX?l}p{j;cC*{8t5pGp{7B_jtv64ol
zk9>r&7eO0Z4@M5>%Z`;aB7Hc*t(-dqw2?|9(g!ngx#q|t*pxn;FIJpX8j(I6;a1L_
z<*qcs$P0bMk>{D%RK7Z22xudfMob?O;a1MgwTVoFsq~?KN?B)$>WeYR*E^X$&_+s~
zl2<1=!o4a!>zCLAf>BE&(nmf6#Qsygl}0R%T)`Qi>PI4tNFR=Hs|t5O8=4bS=_8-X
z;yBPoQ`smE=L-RCXiiLja^4kNIl`@+o2yG?A^lGl4~DO-%Q7zeEDq-j0c}Kt^Z}qm
zxRrCq`k@+?N*{bvt^xE}*%;%I|5ElDZL|un98vH&BJ&n)sB)%vU%mdwyh6rtFxsxN
zQHUzv91z=qkUkvYR#k4a(UwNg$4#1Nl~-cdXJunvKT5u`XHd*laY}@H<=i}Z&%AyO
znRwyH+p)h=)`>&scJ*2L>U<$=3*rV4C;sanH|G6Lxk@R#a_+K+OOK=C<egWcsJQ;n
zv2wl)(ThfouOZyZxjDB))tp1dzt70MyGI-9cMF}De5?&X`~}1xnM)kO7)mSmw(O}B
zWkucz$8i*&H4f*?dg@2e$bDo<VI@7gm2+p;KTc-FJ*=~4;@y35bz#ET4Vqf}e1-D^
z`kY;xSQo@sYY4Y0xU)Rg6_Zmq<MnkA`lhITLTtx+{V;yl`M8b)0*;ova)evCxA9zG
zM^`iMziQ1SA8qW_wc64O_|SYD1t0GQ<ZQ6k=v*<WAlxhG=BY+jkh+@j%?s8{@<>Gv
ztr~SLL>$h?GYud(<JDfNhHxwQmU{<T32Tjh%Jel8@9y)}Z{%=3c4-Eov$g#-gj>0{
zJfn#9&p?d=;obcKjNv`TqBuKWcEv>3S*F$yZsp!`9xgjG(OH9$7ahH_fwM#L{wr}f
zUkJYip;i5pGpY!;a&I{^LG%WQ;x$B2ImKK%D31KB%AoSf`Ph*Zgszx8QbV|vdrPe)
z3w3Q`B$=4c&XGf(IhPfWmGkiq1rW;LWornxa&Ng;4K)=AX5@0sA<c0HHIBpiB1gPp
zvK1VTs9MF^m3xae5?v#l&&X9FsTuJboj3Uk8d2_6`8c|Xj_x082)C+oPv`8)#zb&0
zV$O_Qt~qpdLD$&C;d~+d7=*5WtY1U8m3te)UYsT3{8i`2WMba2L=CRO)dkLFbro~x
zV;4aXr-HbxhHxwQHYzG!|M(q<ybFuHdU16@x$At{S?2>$mGt3=idDY;$T~94Mx|$w
z@go>HGeEiKAQpYaeC2%cj92<o5L|t&Al%Biv+Ez46CWUg-reVC<>S2TeCz|s`O5u4
zEUY2i%DpYGG_<xR6N`>s8Z53ZC?3w2UCr1RjeMboa4YvVR)Z8BHl!%;?uAe`io^M`
zYb@xnA=Ns<t=!vq<}Q^+qz~`zHF9O6a@YB|SB+yOjYuDka4YBLTplr(Mx>8CQgO|p
z2r3(+28}Pf{vnM>AC7RZoSVD6Ia4c*m_A}JqVv4UMjb2XV@Dtm(unCJBHYTkv+EyJ
zomBeBmD8t~7grag59i~FTz+3yW~2{CRIDm17POH_BhrU=_h>`&RiBQP^KlP72<k^F
zeK^9coI3=?Th~9N5AW`osWtER=~y`*J5zv=Mx+l%xRrC$k&%VE(ja|kc2?HuSmiTy
z<ZwQYE(l%CkUkvYR?eN(DC$QleW;^X)`<gc#A_iMhx27MO4nGV4@Xq3+}nxF3cCLB
zC?ogoUUQY=A&x3vUY*mtew3>Qx2ke$KGgM(u|$wam1_=NT~IbE=FW$IAawm>%^Je3
z+*@`KV3y|Ik8}c==-qvBbwT5BzF0pxr-Qh>hHxwQHtWcg@%w{7Ic0aR-b&eM{Z-Us
zcGt`PAih^axL59Nb~WQSWc<B!OL+$wRl7sh??w*i3t<cD?CBsrs3F|Sy^UAHJKrSZ
zUnAz;(R-)we6fD0r%q0-A>7Kn&7RHJn>!Bw+jnu|-Thlb&!0DS@7L(;`E&AE&wq>r
z@saJS2)A+%dENlTT<%eP1`hA;Z|8S)2O!V<G@UQ=D&-zwdF2SVa&OrwojcFMFHUVQ
zxB111clX>KJ9uXH{CU&)sGK04;2y=tH>)Dt%DJ<9?l$8d#nJ0toOpMyyXO9uJqr@=
zt1tMt6OeoE8X&$=L%3J&EgaklsHZKs+c5F&p1Wg<=g;GP^#xyscpbzBHH3TR-iDy3
zMRbqiXb|4r=XWvcxq$dAuQ>HNg#Un;2SQ#s!o6~D*)^L{|Bcaa#K^t7*V7fc15hz{
zJ`@E)PuT2KL%5Z5qYt8|D88%+W^--q)$_Yb<2am;nCA%Y*{&em%Dv?tFS1bi`&TkC
zpPl2HL!Xt6IzG-9!m=P1fLK^VxRra$oQRqVg!(CJUmS-$T@mX?!56~PAmr84a5%!f
zs(Sl0cV+5+&es^ZclY*mMYGDMJ3qO9xATu0!mX;@Be=ViyKy_8CxYJH6T#wHkf!s6
zkbqFkx2hrB%Dv5==8#4{K_+^4pWi*+&GV|{TQb7Ea&NOc;@>3Wc|Iyp-BYDe%|KDQ
zXFqZ{ALrq?qxc;V^2!lz<=(=<UCUo1<5y$ks^6-Zs;_9HPuZxvaz5r=5MKtt^HCLq
zTRAsfICo!n$@qT~LG@G0I@Q<wsRTWPSK*6uVydU;=GGAIm3tfK_EZ{~PbPYIPp-!M
zuf*Ye*|Qn%Qw`p(A>7Kn<@yJAs_#t3zX?L$DC=w-dKN_WqvYc$ArK#M=e<<x2)A-?
z<Gs%5Gi1E<p%E5q4tu(y7)RvGp8ueJBnxVYd>qlrxfur$lt!cvMXFeH<Pp?U7pfn{
zvnOaCZKTqO^x+7%a&Kcbh);7U6TQ37Gf~gi7I_u<mW*((oSUZtxKmyC&`Tfcr}VA#
zkt?cC9L~qS=OCmJdF2SVa&JRG8;LX`eV~nU&4D&z&gxh>A8)|`A&p2Mj&LjI<~e$z
zCyhuS-rX}(Yu^t&d(d<~?x6=k#Yv?PN4S-9$NHg}I+Z^189Daq`Ai)-oR6KL^Cw26
z4@bC_bI1CjijxX~Hp(?eUU4Fa^Kt(!`jbYa4@Xq3vb)8kwD{ypyUWPmWaOHwly&;X
zo~|hBV9xg#^T?ebR6iWyURAlbVs?9h2yQ|IHRmd?G$(4#RW>Rf&KG0e<{8^`n;OEc
zoIC4%c=k0NP9}PHZ%<cL_;_Ct^SYj@TD^vFuiV?LBhwR4DG2ZG8H4)7s5|G&o)~!-
z#HVWrw{maU^Ngyw7a6}bBlqrpA9^cgqmGsH#U~CseIRzNA>7Kn&7P^gmW=-|5%i88
zZItT)=gSZ?K`2Wc;a2W#2<)5CnVsz|lBmB?2FmxZ=lk!ncU!Y~l3ky(-K;*vj#+<Q
z@zn~#t=vQ2JpjjCc3fbOt;D-~A+!U+6n@wFGOtbr!Hy>tgj>0{?3B*4`s`oQ{>;Kx
z67TM{*DE_xbj9I(+|dYP0f?<@2)A-?LpX;W1U6^Kg@6C;l_ZZ^SI0_w)@la?aX26E
z#{sbsh}UZfw{mZ}qm7*qX0zi0yDcQl&gGgT-~AwRIA7eqq;qGs_f7@jR_-nPEnz(s
z9`<KQaup3~9L4TZinH_4$$-$#TtBHH+{(RW=3>;^qv7Mvypo`Ya?MffE)|az$4p*k
zd9@t~?J)K($I21zm3te)Vc8BFoh^u9KC8uXD5ta+mVVdyxPLcCDCUlEE9VYj4EwgN
zO#Z%1Cg!tq<j`kjqd1(8z1u)2<KL|z+{(FGv7jdHnDq}3-ra-fWiOM+;e6}}lq2Mo
zBdS(ID<ranT`9Cj+;B$j-900ZIh&8e__R9=`<HY+UqiT8Rql=1T|*hKY$sB2{>$?!
zYJ|PX6pxaR=c97O8#RPmxwq^cz52&fW#acUeZ+OL`W-Sc=Bzl#n39j48U%aqbdIPY
z+{(R;yOG4*td0kPa>_Lax!OF4oYk>%zH-;0&JSt`x5|61xO#~@<E5LC@$7|^cz4gt
z(5Gw^hx5f9PtqO`?{lmi;a2W#wnNiFWIQ{VC3#mI*BshiO4%q5=gSZiKuA%Ja4YAI
z`zNMtX6lcw6stkG1}Juy5{L7}YLH4J(uc1H+{(GL{jYvPjbeYF!mHTb+wM|L=gYi$
z7sPUSmH$>mxRrCq{S#AZMEdaVJ|B5g<8VICp0l^lA3(7GRR!Tz&dvKK*r|v5p&e+7
zyo!B-BB*S%c*M`Kex%Zf^pVFlBHSzYHp^XR>(ulS*Bp5!YM-i@R|Ve?;^-lLM1*_g
z+^oF#eeEG6edNlCYYwy#e_P{lKCXFz(9T@aha=p|*CJf)W7JfGRQgaqrL5CjlFtl%
z;&8qY&_+VVNu>`*xRrDBE>ohXeb1y1X5@0sp&3ruC=Tc2`N15KSDX^zR?f{mcVr=Z
zjHc3uW@lxc#-Y#3MsYY_tRD^PM{4yWBHYTorDs5Yp8)X^`p7FwfhhLEs_^mdGxR5~
zqz^~9SNc5UE!t>PKhnH@6i<^E^|cs>@likG{#SYZC=r@@q!Eq0%Dok{8$DLaZf`|j
z630QVc9o6l88kO4`-P76qfNz0^ZHS)QrybdBC&pGC(rcPWTJQX+P(94)DY`S;^%Dt
zt9<Xa3c|f2uS#yN%1~|QknyTGD5qR=Fov$OQNBB0tRJ13AbwRtxRra$6%*?0C^CKp
z^x@sTc9-g>#)`xFGQ_SRH2XNht=!vePuXHlWF{-u9Ny_WAMd-PU-%jt**vQsCBm(o
zJA^Or7J)6;ar))c*GRp)pO(K1g!h;<zda?sv!i!uK4)+9_y@aMzga`Lm3x>Wc3{W%
zw?KGzFT}suuX-td*ZJa(@16J9)mmOT!mZrf5UywE;fvXE`ftatk&>U~n&W5ejjT61
z$am-CPFWDwg7|w4;a1MgTc_D`S$p>1Gii-fQ7hIQxAFdrk-SGk9L^WQ2_Upb^<6cD
zTe-LF5S#4;-5EV`jg;9r?!PNd9Les_yg?vxIA6Bg^xe&%cMnt$Zsp!`KLgf!n+NYu
zNOKi+m9JYf2B8#l=VKRI5bS}icZ^gJZsp$cJ{CrOBBN*j_B5Xj1~rbOiG@b5W959=
z8$%X{a4YvVgnpvOZtR_>h@fUQYQ~_(@gH(ZQ51*sajyXgz5j##+bamSa&LL-0$I2=
z8u>pmQ8RL}=Fn$$M~@uN#~O!B*Sk4hsv+FUz2%LbsEPgCJFkFHm0@;{<3Jm|#W<WV
zgwLSZ&w}_*4N<jnZ*OE*Ze=|0H&8svH2~UZo}R~}%EylH@%<l9*AQ;y-trCwb{C&a
z%z4|5vLw!b>UYpa^Uvg!;^BPUa|eQ5{yXg7UO~8(drR-r4C$U^;@B+P%Qc7It|5JB
zE^)pg#EvzDTe-K~0l{wXcaia50fBPLH3wtRZ^v;sUxwfv2OV)Z!maXNtE(fDz1cxC
z1MkPtd{wSFh=q<(<ZwQ+gnjO>UJSymoI8Y@*nPiE#y?2}y}Ku;`ZVt<9?qA&)8beV
z$`VJom3zw$2fRgKA2h<dol@`ai??ei9?nP40OAD@(uX75%DLk^s?tmHisRk=f!R1>
zU2?t*aU}><97nj7b2FE4tbRqtYvlO|V-3$o-qf*jzU)}>#;GK)ABtc^xRrCW<2w=L
z9S6xAVy;Lj>oi~G5sVzp$Nmf;c#~41EOCTexwo7PCwHY0(??u$<eAu1Je)6k+sz;d
z-iV<35mye%QMYn#&PHJ^jYuEf-RFvm<8VH%ROblk!x3)f-0{s0dY4eDk>@i&>^bw9
zA#yk$^_=4*jnolt<=m_dh@LbeeP||A)+rwO4A*D#RpcWcAoNZPWr-u)EBBV`A7r64
zB7JCfF4i3ROx>q=I3Ih3fsjUu`cWX<%DH*NDr(|AZmBvlzA3wVX<XT;aX25({NxDf
z!x2?0_ZDr$_kZN|qg(@Mu2MFNqsphZvFTVz9~_es;Z{{{z5hcQALm3{bEpGQHs<jt
z`LZ{L6fw_xgScvNubeyUeRw}ds!qnc`-xe974_kKynBFoU8>^!I2D9jIX4{u^@)xw
zB2Z4b=Eys;I1c9{OF+<LMXwy;R_-l(<5D%1@tRA#yH{_eH_+%<IbVj*+lP4nM+M<l
z&K>8(&Q<vg=N&yWwcg<u=UwN^dg@<*c!`)h!mZrf_@>y-#S`P3VpH$#SL&J9RW{DO
zbCteN^X|i?r(e?h7d~^og!e8tvp1$!5N_ojvho6P^Ta8wZTIO(^^G!6x)WlV>6Z-X
zcb$(;AH<Q_TXic4w{mZJD+zC2e1bRkwsz@BHCkn$Vld_L`?{|mylz(<&KJV}fcV>D
z5N_q((lhV|$$NO?<`X;gq<PdL$GN=YvdO#9#o>HIh{tOPw{mZ@H({R18#n*m*OPj8
zf5S&E=w8pCIGm4%40(6v&p^CeL%5ZDO9ZjjJEVD=dFtJLA7eO&v53R@LijZZy)j*0
zIl`^nTUHYsryp><c=LPe-MwOQ95GQ|IUn6pj^O?06@*(kH)qv|-p7dIUBukGdvbN)
zQ{<L7oR520K&%9UH>Ot*Zsp#x#yDb1>r}F_m)y-~=eXw3=M~Y2IGhi~km)}Gv3(8U
zR_-lNLGY&LALj_PQLZ^a^cLfAKCbYh*wb>a98t9zS|Ra1FWyQzoRPN~d0s7J?!I>L
zx=qDHF|YEylD+G<GrNXxt19;(Z&y8zm`fjd)sGx#qgmvY^U=wG_#ucVYY4Y;Z#jp|
zTVfv}6Bm(*-rb{(W}y$~i@BTV{j$=BBizcp4dF4~cKcm2e(CLc(!2|c<M<`-<CV_m
zaIBn<9r!_<4?<o!!mXT}v**0$_aQP~*`9Zkk%L$?UqvIzE9c{F9w6Evq$o$Ym2-2J
zleZ7wPsS^P%w*-7gPdwA8^z&#WHg9{AO>p)w{q@S4bn5vNP}$m?jCK#o<|(c$Js&N
zyL=WJkv<&ZR_-n5cR`$rMx;US?sG)1IGm5X2eB%Mq=s-S=VnK1Mt(jS|7S*?kFc57
zk0!H2Q~Gc|uF8OT3IuP@t{~jXy=DDD1od8T=_8L+<UkwBtH?oRDEU}(g18^VV>N_Z
zxwj#pjZ_*jeZ-zB&&0^#e6ixB%<qZmBO=_&z2%G-){ny@eRy}DG6rR%=_B%S9S4Lo
zBCj0bUin%CMKNk+y!4UJ0P$F%4b4|g#oYN=P2>pPe_lbjm2>m1c%mnbNFVu37JK!4
zhKn4|7XsQyq!H=E5pLz)()*BwbWExAq1jnkr`{)@sUwH;ab=XbT^h0a5fN_X+?-WI
zf6|Eb!8hfaLmF2$>R34+d&Gl~Mx+l%RIS|G8<-WSA8B4c$~8b<Ib)owd}yOBjYuCH
zlM>-pRqnaWZWHsEdv~unS9zszI3MShKu|wY)elFwm2<N<6|?k1WFl`fPrbV*6T8<@
zL&V{H!~=xhfBs$#;a1Kaf;y$s$@up{cz3TFsB9F6^JR$XAjIJaw{mZ}--oJs9~sZP
z;$y#CtU1(KDjUV&d|aCVF#v)$rdJSd<=otbM8EJaGF~xPCtR#Kywi8SI45>q12IrT
zxRtL(*q{2^6*?z=c5Gw4qsOP--Opd~_i_F4*GX&lz23a6d+NqV_v$nEa)LN>t+9;_
zY6!P-4@3AHh%46_+u90*clWRIyQgh@^uRoR*ZJ7d8pMpX#<o6JL%5ZD%NrzLTedU#
z*;U*395z1HXqADg8RNcrLHE437Y)dF=VK=X5c5GiT0^*%dmF-Uu3yxCdin|N-M%wE
z_3nPdxzoDm&)vE!4(DS}SrAVy2H{rjEoX`zdaS*}iT&+64j7+$cfW1p@a~uS6NmHh
z<O_(egP2`IxRra$^UheeLEQ<$yL*j6C~-I+cjJPXwHSn3xwj#_%c!qs^k3&#<+DLo
znwUW>G<tD3A8Th2*MazE4dGVqE%zZ4y}gLyY$BM?YLVk~a_Tj5OB~J@!oPFGoi&78
zxwj!q`Rv%1vQhfTXXnVF&-2iTIGk??q1<(ZTe-JX=e5T+-hyL85PVavInYM0vT;5f
z&KEfv@=6?zs9FuJka!pNv#Tc3$lZ+GyL+_JTg0QvcQ1%<fw-@RaH}f!GQT__Ih_cO
z|IYYS{gh@lY4F!`w{9x0G!EyZ%b0#bay5uuY6!P-Z@J6j#Qx-dGEv!{cQ26xZ8W#Z
z_2GOxoeAPkAnvIl+{(SBFL~&(<S;V+HW1$3e{t@#<^{RVoiBtZL2L#>s&#~0xwokH
z`bCN6l5a3_@9v33vzV`(kGIQ%_{=X)Xzx%%xRrDBtp01urprCIYWof%n0Li-&B09F
zR5mJioi9VYheqlMw{maM`L!#gJ2F$N8hCe)He!#Zm^&Z$w}YVKwAQI1+{(RW&wdbp
zrbc}Mgm?F7qp55Zhx5f6l}-Y&ehuMP?k&4CbF6M8<CX2q$mJSf_beW<_Bme&yMy@C
zVi0cS-qMQ_LFUC|91%nt<(dO+G?k4Shx2h}8^m)URB;^PR_<+<yV8j1Bld=9qp56E
z%$<+xIOMK0V)}>(w{mXYy^QrkprjA)?$Jh5*(eU@W9JqSPl1p=9N|{EJ|Z`F#4~DX
zMEdaVKA#yPhx3K-S!SPGL437_a4Y8y0c|AGh_Zd@Y~GE18`|hqHfkKsN4G@07l&{w
z_m<uF=#-=p=_9WOkt3g}W6Yf|1a)K5i1gtIw{mZp8BminlGl%N%^{5|8<khiNB0VX
zJ}#9$_^lFAwTks4_x2WMh4CCK^;h2ApHDqkHtJYa`R)Qi{YdYxA>69U{ZVGO^NFB3
z0PpV6hWe5?4(H=J7G}_k7lUvs_m-}YSz4LcB@?~7S0|%v)Hs}vJ!`VwC+4mr+{(S>
ztQys3A2NOb4%NV7&4D)5Cq@qE;~rrU;}(N(E9d6Q4^{Jb`OKgWvRHGdmMa?-bLZoV
zB8Y_`q$o$Ym3te)&Q$*|p^<xupm+3Wqg)R-UkKj-@fi@BMIGT*&dtt7zn<N`=k{TZ
zweFsfdUyZ)<NL++#{pm5xbK9u_wF9P<$AsP%uEg9hQ%P<%01+*h#;Q1eOT+$aG-~B
z&5<Jp^t;Z-`zAnqa`v#+nl*%5xwo9z?sVF(oie}m&D$rW`DkN({o$57cUNBJvVm_8
zjU)2$6a<KGojSkuts26;a&MVSMqSq)dGV<BWt|BrdMMW%kDUDF?rzU5j2zBK-wtA{
z#UR|uy$xZlH)pgzGV{3h?^_d6@9y84zE)R%I#$ldSv3$>+;&|1k{ZIT+}jXF@%vu^
zbrA^f?ssJj4>J}WE9Ya=M1Fr?@zfA*<=%3g^P*90js6cDEAQ?Vi^;^~KscO__YZQM
zz6#><8p5sITi!ZN^sXX`dl5m+YQ>sk;+E?TOd)4=teg*y)8@CX1u?#ca4YwgD_G>M
zGXAq<Vm>=}HAm_5GiXE{&KCmOXzdO{igJWoxwqW?iJA@pu__4f?m_ex<8VICS>7?M
zA+N;Yh^p1l3Q5o1&2K!%$aiDpd9{q=VC20&ARdZ&mGAmf=QnQ4kCh|bs>*%*#iNof
ziJ<hMeo9%Vc@S;%Dz6lC=i`2A5XoW?ZspwEv31*V$@j>_-;s&l-IIw;WurKpkE=@{
zrh<?@9N|{(E%)5LIV0JYjQ;}&^;5-~gE2IBK_iN}^M&wj5LbbaSB`Ki=jJY?QP(Ah
zlJQq?te64HH3!;gDjPKp=i`oe5X*zOx`uEo=jQr(r<05$<M$xu-rY09HIG0e;&8qU
z@vX%m+{(RW<iDPsyx1MqT7yjV?jCK#^@om?^KpkMh+hs4Ypq^GxRra$U92GfuRE;q
zKOnrj&k?;Ehx2iU9>fAvB@Rcpm3zz1svN67kny{6th~EN8;VD7UO#l)A|G7_h-W}Z
zQI2r0oSXBBMDQjuUUQ=QDP^7Xp$IA)b*!9^-IG8Z1R}3ECBm(on`=d6;&;jT-;#;m
z-RGIuEb^-0;~EQyBSBnRL%3J&EjyB8EsaPY-rb`O>AY$Bh<scP2O*6}AC7RZ+*`7Q
zQA;DzM?M3@zD==EHfkKs7XsQyq!H=E5pLz&^ggG}Z%8B32QzZH=1`0kMR7PEeaR9K
zZspwU1Wn#bBhp7+4a)U^vQfv%`FJ18;IM`?qTF?aTe-Jf;T;&(Iv2z+5UMQ7I^`9J
z0gXesn^%U&$30%uC}~6-j;LD2`jLD4HD-m|Iabn#clV5ZK-nmcD&K2V=f^<oQbV{^
zmHQTEx0Q*Y^x@ro-T_44oiBugLG*&iD^9seaVz(h-K?0UPa_kh5AW_5X8l!+x${w@
zKpX=?`f!9>xwrJYRGX1xyt*at?$Jhf=UnH`$Fs;FE(Rg59N|{ZjRvWjbWHJB<sD>P
zbEvai35_VPoG*k1h`(MmD)}Qv*AZ^z+<GHhr~OGXUj3DK^yE~z9&o-8eh=af5bCcS
z;a1Mgu0#7xY)?x2dv1C6nAE%banH{iP&UqZev7{8qn;SJW5wU3=kqy)l|lT^=>DER
z*AQ;y9#U&Tj2_)TOC0$(v4-b}jyRkz^J=rjAl%Bm<%-D`2e*6gpW1q3(J`ra_a8a>
z{((<DyY7IEBl7Vkvn>y9KXvca)~hvyd*$Bpo}Wi|Yp;I%yw=a(J|^|<{?{wsI`FM4
zj*1-4$CLaZHdqY8t=wDoVcB|gd&kMc+jlNJCS`Um*Bt*n?t%gRX&lbS(={M&0daQ?
z;a2V~SE{kz0hBnryVn?g#aJ$Y!}++=3B-1bLAaHB%az+7&1+r3=+ER>X*Qr{bTyy<
z53yL5m?)Q=FNDjEo7XxV#LsF7w{maUjpTu;t#ydv%S6z-`#V?sP3LTKR>#Wu==ATK
z+WG<rWr-u)%DoL?{x<!ywjh6%yZP+gmA3VHH#8y+=i>?%h|L#+a4YBLT;DeRJzIg0
zqP)A8=8ADRU-YV{@X8TYtDzMV*PS1j+Bgr5EM(-2uv~MXjou<2RldjWo!U5zV<mk!
z!mX;@um5OXV}r+bYoAI4y}KWG#ZkS=M&*_B(eEBVuQBhDEK3~WR?f{nsgs8%pCuFT
z%Cfy&bNt_N7xexo*N5|QP6k9P&)vM=iJ9nD&P^vXYIJhp(V3#uPZev9b632zcQ-Vm
z+;u*7v)X!e@(u`j<p{TO?hqC|vRiW1lzFYwIac1?Gc)vlI-jqck89hH?v^~9&n1p<
zE9d6fzAX<<W;4URNd%e6$~6Z$)vIjOIGm5?CP3`UEc$v4;a1Mg8;JLtnEW8k=63Jy
zcYl6fud-44a6Zmi?mIC#aAg0inPj^o+{(S>EG!&r=M^XaCawX5P&SIg`7*C&jLd43
zBizcl>EX6GIGKh<qz~`z(MB(`L$8jN^Rc%&$LiGkrZ#k}9N|{Z%{xXO+b#LTvGW=y
z6LUqXSOcJqUS*^7;e5QKh%EWmlzELaY6!P-Zmt55iPDJaBlc8iqgUCeW958Y1t8m3
z198WYq9Ver+*@{UJ$iT|ji};ycaJuDm5n-9&d2pA5YmXea)evu`iR`DSQxc5B7JBE
zP}Zq_pp9N-qsHNUtRR_7q!CqTN4S-9vs(`{oHQbRcz2IBH1Ec7I9~{8qalqb=8kYH
z=jQAny_Ga#^&?*ca0LyGXdKSR+L^qSMx+lZs92@Am2)#+(J8$QVipMR?$Jg^;}C~_
zso>*%cU$()l2@uYj;LA<akHl^wN@I*>qoif$SY?&R#m=@nMao4SgF=>OiF}XRk>F=
zc3x|HG;#_N^zME<bFQ*c`fxthgUq1Pi1gtIw{mXw#rgj5_9-JrC-;$w-rcKqC>xc#
z&KJUPdY=`@_Q4v$t(==3ny5BseSdh;1%YzPH3wrDP&R5D&d1qW5Pw4>Qj{ay%Dtse
zq*FSddaG>r?jCKZ-;LvNKJLE)aUL2ut%h(b=ceDK`i~>yUn7FdWaXN}JALPi^`pH7
zh`fH32)A-?x#Iqp3tMZBT(IpsjT2Mv?nm?bIAtUAY~S`XXWsJ0qMvQ8&)n$*;+dTm
zY`dt2a4Ywar)xkwF=D|?ad>w>hTnaC(a&zy?>b+`vB$^-GvBQt+{(RWZl8Ew>(^iC
zpLIKke6$17-9=NcyJf-|D>0T4ea@F5&INHt4dGVqEr><aS{E)etaTv>MXgwK>~a4k
zx4i%6E|J6eLij3(UxB!whHxwQHiVU4dLTLZWB+d7GJQ<y-95cCzkk{-AAGP?pY!oH
zkM`u&QxA=4odSm=+{(S>8vT?{bS}K}*UA1XY@d2}ztmg1+%lfAWN{wi<0+xJV_N^7
zIk|Nj9FB0WoSXALjQTVkFA(``Fd%IoPAoKfaa8$Ee*Y7#^WT{kkCh|bs>*!~(HleD
z6?64d%(XG*<m$~|A-BZgeB1%CY5%PE4nMCoSVOp#bBAyUS*ZL~?n({CnnRznY#h<&
ze4JY%)As|R+;xOoxwrHT=&#^VKV{?CTB^$8F~k=+w%s@9aD;p1-0au)h5nv1@?+)Q
z{lqg?8dbzY8Xw~0SoIXg$`Ni=<vwtkVT}*o&SLJ}J=z#m<dySrC&&ArXk3`b+!1c&
z+`7|gOk)xmuWa}3eyMQUDE{KHa=s8=pF5_p_pA>i+{(Fmx~4t3vGtrxQQqA%hEavi
zosYLO4NPu4`@n}0Zsp!`AHbq%jX&gb2{S;s<{%cM7=0Xv^YJ#1x2HAsUM`zU9N|{(
zEj`1;^BRw#5oNpj5M`Y-NQ_4*8<khiM-K<${dGT#a4YAIGj(GvG$MU?caJtkDFemf
zd>P{Dofgbg{cwa^xwo-?G#(w1)hK<Vtkc}iOg-w&ye>IkhS&>LNrR4XE9Yi?&dC3X
zMl_drcb|`(vBdgO@^QxN@bel!LnAu6j&QG>n?pg&_d_Fj{V3NQiePwam&oCKAzbjz
zw1za2*N+n6R?f{+LgemzG$MVdpHkN8Sml`*Ih>Dp%$nSgM(PN+qJEUzTt%KUrXh_8
z;oW_%sL0`b>?8cZn1(b`N4S-H%l;)CCuu~-%Dek~W{4cl$2gW7){sVYtQ_H1?k&%G
z5pQWEuLdzwq>p@tiyY3!y@Tua_edjkgj+c`_cM@%bWDw+ehjSv;(Q)CoR70n^j1Ph
zAJR{Oa4TPn@B}Cohx!pY&_+jDr<kLSo0W}LKl0CLBO;^^M~E-eh&sI?ZuZ$`9+5_*
z5AW_ZS1B9CQRPD$v!oFnD@V9hl^boeq!Gp3yL+^uz9e!uU#uUki_nN-?g+PXZg%Wt
zUVkpnUGMJ6#9Ne&@mLvO2)jHorghJ4S?)T*y>f2uNo`GToqz__EqQmJcVv;H!WSO9
zZE|a;M=}jM!mTRYLA}+pXk-Q)c?TKS9K=HXZsc&jIG40e%;yqExRrCq`q7$~&xzjA
zt4B~a>R35nhWHl<#oQ5Y<=i1`_|S!|Pk&~?wjVT3?C^H_!R&9J!F?~xKEJ<k@4@c*
zdkp^k`SFALtgQ6^4I!#lL;ucsI5<A_nFTZJ970T9bC3AD&gZ|c&qKdewHo?&-tBbQ
zd95|q>YsIAZ5#)ky~p4$MrRz(w`7j0m2<53#<bS{%MWW^TH|=;CwmN@wR^_ld`sr2
zS~<r|I;9zRjcJ`x<2d_>JqGos7&zaOIjUC9QN;{>qcI4T#d+x8Ig5!$@=RX2m3!rJ
zsxtUDiiJl1VGbQx9bp|^w{kxJC7u~2-imuY=E^BWQ5^nD|1VyljhWI&J$L1avQZq)
zw`7j0l^<QSv8^;x=aA-#aX8<SIjUC9!LjOD+_C!K+2hN2RQakbL%&tE8v1v3(pi33
zqllnK>Zq&6mwDxURR%{?t(=2QY!sR3h$GJ)U+Tm8mdsJLat?mKu{aH`Gi7|KbLU$!
zN7c$XI8Kel%_SEdF}|FyoNvh-RV(Kp-i^i0iRpL8m-DXkEt#Wg<s4|EA&t~?_tJyL
z$66wNINy>vs#eZ{HX71MokNJCE;-+lIjUC9QN;{>qme5cKRjRMN4H0_k6XD{o}rvi
z%wpcC*AGQd*(eVGrK(T<TUDzeZ@E60<!+;1Ka`336NmFHnWJju9JvNR#Gx?=B@X9X
zGDp?QIr6#WL&r+7(CEeCd`sr2S~*8PCpPN!Lph}=io^Mq%u%&+4)s>jNIiEIO=Y7v
zoNvh-RV(LEZzYY?IaG0qaX8<SIjUC9QN;{>qxxFJ!!xw3&hhBFRo$x$wGXpJ-T{0l
z=IQ{7yz*bF`sD0Yt08Zfj8G?2=)-@h?0sUkyOn!YMW{YWBXtd`w<>h*zx4m&mFlcC
zQlCrITNU$_^DUX9YURgD{em=7*N1wNV%~MWC394*oa3eaC$?YyYJbm+caQ0KcmI0N
zp0R%1`oP|UmyOzU@LT9!pE=KP;KcSp8}|2{UPHK*dl<q+2TW`qxnuvVUx3gz%0TH(
zh^xt8{jT%z&KJBI&x-rEHH2Haw_KYTd2oBVho-h>E;^>8(JBMQ(Z9i-gI7+@IGm4l
zEr`F~J+*aL4dGVqEmt!h->tp>3G-TqynReZQBwwrWBV2M96aKhjKleOYY*4FM*VzV
zYiteSR_-l#{Ei;o-jHjAmn=M{<K6ul&+IX%KOHOQ<0*(Qk8ZC&X?XjZ8p5sITkc%u
zisHT7jBft{gm?EEgHYmdKHeF|?;pGE==O{n!mZp}?!P*5Uh8{L?AHDkBUdamN*yc3
zLZcUl^KoYT^m(mwAKtBfR1M))?k%<U!Ktm=M;zR~od|k&ubfg8#o>HB3vkcW*0`+?
zZr@ZxxRra$o6Nr2KT8^sJ~Yb|YYu%@Hj2af*j0YR{#nw9^x+7%a&OrSX~+H^X+-*v
zu9bB<R?=KC4(AI2ZS+VZ(uX6eR_-m@Xh<W{hvugu9^<EtFXK_=LmLfgMEY=qTUEKy
zMnf8rJ~YoN>okt(yN)mO%K1V-8x3hh`f!9>IXCA^xUMaYNFS;f$~tkpu=e;;AI`@P
zv6F@;(unlo2)A->u7+>Jog!#N`cT~~bUyQ)aiz|kk9|3}9i2!c(uX75%DK7w^9k-P
zK_k+K>US|;?fT@na=vmtcK&*Jw?rC|J{;jz&dtvCBMwfYMu?#LkYe85=k9Ukyz6}2
z5x@1pi8La8IKr)*o0(z%iHS5KeRy|2?DlbG{cyfmKN4v~`f!9>Id`ldi8La8=$oSU
z2~pG~=Zp0tkw&BsN4S-9$NG^-BhrUPSk!QhT-hj{J726Hi8La8IKr)*JJydx8j(H}
zsiGHA1eJ~AaK2bS5@|&GaD-bqcdQ?YG*alJ=&6*6`V)uq#rly*BZWQ+gj+dxtRIOq
zB7JBK`nJ%yP~vdDSU(bJMEY=qTZy;ej`gD<jYuDgg+?h3#X_SOhmKjn7wbnu8j(I6
z;a1Ka>qkQxQT<R(DO%!CPAQ7waK2bS8q$dL;Rv^K?pQy1q>-Y2DC@+bnOfN>4(E&Y
zqlfyD6!oJ(xRrCq`Z0_8kw_oXcrgxDoMIf#7wg9?>PI4dIHGDbv_j%8J$ftZM<RV_
zt}5c8%2~vt%BS9n`jJQ<j&Q3gw|Xn;M<RV_&Mopv9Y7Ih=L@PIE$T-ieK^9coIBQ!
zHuWRw!`A@nWD0#aU#uT(>POUvBizclWBq7TKcYUoyH{^j=-l~Y{b*A^qCOnqR?Z#k
zN1OW5kUmt)i}^~uRWV;VU#uT(>PJKRaD-bqcdQ?6>PJKR@Qz+RNipv_U#uT(>PJKR
zaD-bqH?t2<vF>wod{VFD-TfElW$TaO)Ak-5{qUZHKf7}Lpgwc|?yqLI-+XXbW8b?c
zR1j|E9&%OYyxHyF-ZZSW3kZFq43zGK_`$e6;_o^i&pYGQigSjwR;nS~%DtsVwL9&5
zc}D7t+b481T4kU(?)&qegD>`H9L|>^eh1>r8p5sITkgHt`nvWoo~wSQGohoXDFelE
z;~9GnF1_rYaU9N<AvReI!mZrf5FUSnCuMo!>9*E{j(7L}-FeSJ{V4{{$KE15yR^-%
z$F(1=A>7Kn<$fQ2|EKWG0iiMITa80w5K0`*$Gryp{%_#9v4(If_m=&T8TELMmyVTU
zp;3xMvC!zn;e1>r1Mv+I&(;ub<=%4tF43Do+!b@>l%gdL<&>f*4(H>}We}Gx2H{rj
zEq8L<Jgg;+@T6W`O%!Vm&D6?9aX25lay>Y#C5=cQvMLa6<=*m4^-aSX(unjST^Hk!
z=8ADRA9uvh8P<?Sqz^|_t%kUn6L~638qu-R{8YqawQr6u<5A^98xbLWIKr)}+-M__
zMig_+vqfIL|9|7kymCIC3|j)it(==3fp|(!8j(I!FO+qP`Lcf+SL(z0LO>gdG$MUC
z!mZp}_LO~NMk0+!AF6wW&c8lsT&Z*CV`n^`Pn1Ta4@bC_bJGEAeO)4r=vb+K7xUG~
zG2_bl%K3PX1%xyreK^9coSQRu?M@<%DCX)zih1|0O~#e;uJeU}HWF#1j&LjIj`bsv
zMx+n#?w|a8w*HVloG;doL>iGk9N|{(ZLA-OG$MWIo1*p!QPd^pi}fRsMx+l%xRrCq
zYLG}H8o5ST)NqYl*;tID;L8xwi1gtI_sY3rHAtip#axjpdJ#oX*(eU@%Mj8?9pP5a
z9qUITjTHJQdMah2{>0&Yv3?}dNTH7c;a2W#tRIOqB7JBK`c@j$7=#js^Tqm+NF&mR
zBiu^71$V3;i8P{-D;651I1~$wUL4LB>qjDuNFR=HE9Z{&qalqb=E^BWOB{-^q9_jM
zi}j-+jnolt<=nA;G^ihX(mPgzV$GqMTG=QL=Zp2DLH$UI`cWX<%Ds*CqecBlqz`Gl
z7>6oOF%IX8^`k}oNTd%(RIP?qNU?sjs2_3UnyZR<sB#wZsPd_|5<>cLgj-d))mu?N
zV$3z?7I~!(pvWudi}ho%2)A<XSU=j-kEjpt?$yZ@`f$EjKibrfs1HZDm3tfON1OT)
z^`RPA=v=*3p>yYp^`lMwi287ZTRC^EA8qPKJXWgZ#eAjSs+g~wFV>GX^&{%T5pLz&
zv3|6vA8}4pPg2ag>Pd=u(D`EhXj4Bvgm5e8=Kgke@jPtxv5ixY9^diqe#Vh|#Wngi
zo9sQf#ld?Go_6N=L4A(9A$D$FacpCB4dGVqA#ZSjW2rU9wzdPIZ<K-3oe)3z-#y~*
zIv@A-fq3&{V_Rp{5N_q(@?6nt%XW7B`>JjB3hQXJ%0O|<-EObJX)k0P&c|~_Af9J;
z$`@(~w{mYe!@>?``<#11`%~W;-%-?*f#O)W!d`=o-S>>+a6YO5i09cOuCIo0EBBVW
z^yWO)ZcOQK|LcJ99q;bnzID$*{V4{{#}gwUwmhl7J++2#EB7{pkDSooo<8@n_Afzb
z4Ek2%&=`aghx2iF8+*{*{orHm<qxbP+{(QTVZk{kv=3qQJ&atj&?v>BSZMU(a6aA(
z!0vk!ZdlY_sfKVX_cnx6o>{d$pD3O}1eH^YmN=ACo0D7Oa6aDK_x!4D?Yqp5y%mI8
zxwq_$x9Ql{{$%5i$i#YguWS^D^YL{5N@H6W{%M8I{gbK)w{mYoKpPEdMEa1fi*ZPE
z%0_XByWpcneQaz)8j(I6QMDRcA+dMc-&akf5$QwoQxT8vuRgwvN0kq4B+`iV;Rv^?
za-)qz8j(IU&lY)g`&Hx0ymCHvtz$nxX+-*Pgj+c`yLe3LPoxp)L-j&gr<kv_!?;o(
z&KCmONTd<z!x3)f+##TiL>iGkRQC#<ul&yLrOus?EP3#;L>iGk9N|{Z&AiL5k<y6t
zq556SSDVe+y_~O{k9}BfSd>U3(uX75%DLI4;I(B_X+-)^A5zS_b7t&b&b!VR_aslH
z5$VGbZspwEGr=yN(unlo-TimZ+P$nF&PNaTrxj9ZMEY=qTRC^EAE`7VedwE__C3Tn
zin`={v3{h|i1gtIw{q@SKT>H#`p^iA8m^Hm8>K<#i}fRwMx+l%xRrCq`jJQ@(uX2d
z^dgF&vQZq)7wbnNjYuDka4YAI^&^o+3Vjqkl`>I(;&8rLKN4xA&_{uAE9Z{&Baud=
z4~;?JN`o4MP~vdDSU(bJMEY=qTZy;ej`bsvMx+nLLZcLiVxiHC!}((UNTd<z!x3)f
z+_8QnQ6of9Ii+ZcLph~v6o>Q0`jJpSQt87HZspvuel(<!qJ9)>4$aicMsYY_tRD?&
zq^KVS!mXS;){hqTBb7d+@nRgRILbzGIA5$EE$T-qeK?|OHMByC^`lMwNTm<WRYg2h
zIg5By`P5rcKT_$#5pGrGR&Pc9i2Bf+TjZ5GfFiG)FV>GX^&{%T5pLz&v3|6vA5kCP
z-K&!+^x=H5ezd6{Q6G+QE9Z{&qfPyY`cMrlbgtg2(7E%)`q8F-M145It(-g7k2duq
z>O-}>n6K1Z74wzz#rn~venfpZ!mXS;gk`Dz>oO-kO9a)E6!Wfnl49O<zU-YA)Q?2^
zaD-bqcL<-ne3Q;gKbq3$-@B*d-Thl{W$TYeZ`f;a{DQp(fBS>+gZj)`4aCNW^ERj&
z!mZpxcC-fZ(_^Q!#)9zfUWolq-y{C6^YPRch<zqZY28*sxRrYwLW_5$?AM>R*V(0~
zqtPk@rMov;dky~m(|bk^=L_L25FfohZQr<a72#IyEjuXl2AA!oU($YWhn|k2rVJFv
z^FQ2caIK3n4(H<v7Kqc%y`<e!L%5ZD%ey(QSgtdbH{M*?*VFOt{*h1bHK;$u!1>ti
z4aDa<Pqnw%u8MFg_m(&04LsF;l6T~63_@ejw;G4WAe1<qFNA;b`~QJw23|SBt=wB?
zIF8dE%dXzpiIFQ78l^ZC3yoeJ&c_?_&bg%huccS-tn$?=!mXT}9f64VF+}m(#N4}k
z<&>f*4(AJ@_v*C$3W(Ke2)A-?xpUx0Q(Dr9^il8bm5t(Xz7W=#I;C|3?*{x&4dGVq
zE$2&)ozjp-qz~y@S*K$q%@yMicfrS5wFy%i(unloh^m#ZSkOiyjYuDwpNe>_`Ok4>
zJgR(XBaud=4@bCFl^bm&(unk-dA7)_=f6F!%q!<(_kwdTNu&|!!x3)f+##TiL>iGk
zR4<fuI#!$fbN3;A7$5t?b)HJ35$VGb?v-=XC-Sx_X+-)^-79py)=zgYb?$s2pp8@-
zkv<&ZR?bbFmR&uSMx+nb?_$3CMBnb^eC2#-aOu@kX+-*Pgj+c`dq?y3FKI;j@b3P~
z)pjrEUFQn{ZKTqO^x+7%a_(3^QfWl`@a}$<m9zDS^x=H$PtAMJq!H=E5pLz&v3{h|
zi1gvzy%0rRa=utUQfWl`aD-bqcdQ?&G$MUyghdV4$d!%7I10X4KT>H#`f!AM<=nA;
zq|%7=p-2_Ih$5(L6o>Q0`jJW_(uX75%DH3xNTrcNA4N~4Ow^w^oG;doR2nJtQ6SvP
zxnuoEq!H;uW6-zKpvEASIGiumk3<@gJ{;jz;w`vi{Yaz{=|i#5D8->zX!PQ6zF0pJ
zX+-*Pgj+dxtRIOqB7JBkE7lx}v7#sr=Zp0tkw&BsN4S-9$NJHrexybHDApXBsg;f5
zaK2bS8q!EnKMI6fId`ldE$T-qeRy}Tic^fk`C|QOQ9n}Y!x2@hp%oJE^kyERex%Ze
z=Bgqds+>hUs(k9Ls2{2H;Rv^?a;vwZex%Ze=G-E$)BzNE<$SSzw5cDd^x+7%a_(3^
z+SHG<=&y>5QzujE!}wzTXj4DZqQ5E-?v-=L`q80&B+`d!V4-vMR)x-;FV>F^^&^o!
z9N|{J7K!zvL;Xmk57qKwzEW>h%va7A>qm$Bkw_noa4YBLd9^=o(0Q0Su|ovCqgPK-
z%!AIyJJ#tJj=e5TzDSlh!mXT}9UYDx+j;P}=QRH3XKQr4yZ_Yw*_!#2U)XE#s$DoA
z{_XLD`pj+%ARfB(oJO;Ta4Ywa9UVa2{M&O{2?%|o43zGKSm~lY;_o^i@c{AnOV4S2
zqK0rQ_m=zNj^3^_`K|-nM@?CyqtPk@#c}L<dk<dM*fVlCAG^wfSmmAr+MCr7Zsp#F
za31fx|L%G-+Z#+;qob%P1I2OBU-laO>Ju4<^ZkFEy$Q4}Rdw~hHJut0Lt;=G8^k`w
zE+C@s?2-#2JR&OHdq7kO5@Ul6sPO}V7g3{PhoVsgG)C-wid~3ss^IczY!RbigBt9j
zD0Uz*|6%u^zg1_|o%_IrF}^p(9izspHP76;)|z|n>UQZVRv|ug$MgGxj%yJf<*hB{
z559AU;rKnu@TI5jF!=0#-@iR<l|OReeEO!|x9>1~N{F{CAw0@kTgqvQzN>iNCWO!K
z=>sS@oKIi75#r0@`RWqFqrA1c&$-k2{g>ac%kVJuoLtaLjup9}_i#9$j_#e#?~i)p
zF2kSd=sLoqytP_wTKPU*S^O(ykeVV}a8OfZ5f10m_wGJ%M8BU9$1Nc|%DHvzRdtI-
z&<CUQ+Bv}HRjLsV=hJs(RJUjZeK^9SytSoB8#6S5KG1dc1I=YW_@(623Qf`<8bKe9
zXrsKf(ndui=!5Yo^YN(FL#BMRe9}fmBk02s9@TP78x@U&KI|Oeyn7rn)s^#U&Tz-`
zD;f!X7{a5RTTjaFQ93k&KJW{wj-3D7SqD%0a6UabAOspgACB-STmg;nOIk(nJ9p^N
zNa(}PInMgagD0IkpH>(9_8mGj68bQNM>)3;Z`h?nBj^MF&hhF+Uq5IXubfZMg}-r^
z4vnA>M|c!2oug7I(ng0y&<FF79Cr_S-$B#3>wNlNkq~GEeK^9SoICuZLnG+JXZJsP
z^FeXV6@55g_(z9E(1#;D%DKZoIy8bl*puD|B7MpE!aq7Rf<7GKQO+Iy(V-FaK@Zcz
z={ePy{Ybv>j}DEX4@Y>boICuZLnG+JXZK`~YJ|i2!aq7Rf<7GKQO+Iy(V>ylN6x9J
zME>A#zVMF@jif#j;Ze>V{?VZk^x?C6`Tz<J=L`Sn&<Of)gh#=f+~FS;ji3*5K`&7!
zxuExOIA8cjMI-3L5gz5-;U5)^pbu(_Y{5Z}$s!!i7yePv2>Nh@M>%)+$4vgwrGMl(
z2O~Ar2#529f6U|`UHV5NJj%JlKj!j}4t=2U><5mM{cyhUkGcG#Lm!T4qt?!lbe5`d
zME=pC55}s@2hN%KX!)31$v-;u;RugvxtUwZKf2I|=2lc4`CtZ+b>)2FAHDpe3w;>E
zqntbZV~~G@K74l1OeXc=eBmF1{3G<?2#;dS&DMs04DydrA9l{c+$weMeBmF1{A1LI
zAw0^t!#@W3M@1iad5%}it#Z6_zVMGh{!!6~BRtB_MZ!M@`A0<`%t><GWloagp!0=)
z4DyePJ{;ju&aIU^uRC(M_8aHS_Ph0udY|2w_mAfvpL*G$tFt#BzWR`ZHm$N<%C$n=
zPl&rLAw0@ksJ8)y_~kdwneQwFd#FHk2jb)J-aP)U@#(685O){ikxK|~Rnyw?s_hON
zj=1djez@l9p59V{F-9G8`05FV9~vADUpYyL)4y<h|GaBkghzR6^#t7u_8lJanM?bx
zTzPd*)~G-@?zHXUtEc_$VZm{p@#&32A>RAhOZ)fzSBvl{Z*3`O-~0Z<d$zyLaKdF*
z_ddIS{u>Wl<&PXVUnv&}aVsHSeR+%UC~vKvj#l&&g?f<?^ntzfgFb+Q!};{J8zFup
z#M9I&M|hOCR?lqfIK4u>|EPLSF6bp3<bvMA;e1;6T*s*h@m(EVM|hNT>*-JB`zU4c
zIAxHUB3p1!Q)Cei=PTvY%J)_w{!LZl2#<1ZRlBMXji8U^**(<=hx6$>aYCRG^x+7P
za_&;3jTst2ALu&!f#$Lw{8IAiX%QjN2>NhD8|7y#(ndui=!5Yo^YOyB>Y3+7K3YC$
zqoNV?;RugvxuuPYM$iZ2Y}VD&e|hj!SI$?8v{BIr`f!9tId>`2Mu$ew2Yx}-k@J(D
zaPXuL=PO0p=+FrIaD+!Wx9)A<`~DpoK_B>D>inBm9W?3O`7{d`0*#;#M|hNTYX$xn
z?AxIc^nripc(u=~51Pg+=hNCELZA`!;Rug%ZaL0&hjnNKeJ~HnarffKA2f}-&ZqCO
z2!Te>ha)`7xx+s?G=e^Sc3&Pno`0YZ=L`Sn&<Of)ghx4d_(z9E&<A_c`#?;-WPITt
z9U4I&j__7Fcfmh8G=e_pVR|?{ry8>#6CeK3p%L`q2ya!xjem4#1bvXHoQse_su2$7
z3;*cQ2>Nh@M>%)+M~6mIA33L@68VF}`NBUsG?MyAghx4d_(z9E&<A~BFB+r|px|)6
z@Q)6SpbtlQ6uikD{!!5g`XCqd5)N`f@8NL1@Q;c{(1#;D%DKZoDjGo_)D+o*gB+7Z
zIGiv1qoNV?;Rug%?(mP9{G&_%$a9Vusa3OZIA8dOzCv?km;R9mk8<wt4}JIS$PRs=
z@$3hVll^eM@Q=CtqeCB#XrtE7kV;{0g+|Z^V^!t@=gfSxe9W!T2>Nh@N44C{t<VVi
zV9d?BVg`_P<$U2Ez5Js?ACB-S=MMiE<R76ApWQQ)Nqsn9_{Sjs2z@xhqntbZV~~G@
zKJdWQIdiMjx$}j84DyfAha)`7xx+sO`A0<`czKRjczKRj&KLeM$UiFjaD+!WclgI3
z|ETDLIZ2MY%t>+_biVM9LH<$Eha)`7x%G|1XC6De<yYs<ezfP=z0d9+^PTbh;~%fq
zyTe~NeD%W{H?6X*GcO_D_rK@PzPHjMJjz?BGe{vW{MEVh9}B@ADiGa)*zd!e<9D4;
zbtS}U|9kHII<?9XO;mHnGW_OA!#OW|QU9F>oZZt~DiDsxo_+Z0{Aq^<hx3*4Ss~v0
ziWl`?yI+g&DCgESvIlG){^RCX_OIXR?4GPqfpA>4hu&q`>9F8%KHcpTVs^7D`@h<`
zMR=6AR#)zRy<zy`3%48oWV^F_pWSb{@UT_>$bs`|^bz7RA^vO$;Zfe&Quuv%{!|G1
zz+U=6A3(w3d?Ny$WrRn0YxQ2Lj?;bB`+KS9<bqzpK`!V$9L}e2jOaLBD8$QkbRFSQ
z-dcUFK>5DSZ=W<=r3_M2WD5>ziY&t6eBw~Pk5}HWRh2lxqnumsd0lhv{Qk!tJM6iK
zs=YY}*u45^X#@`E)3rX;EgA`Z7{a5RTW)ZT^e4pEqz`nR{Xlcs54e+0*RQ0%JE~O=
zR;wJ*M)?_wv{BIr`e1y@eC&DJAyYnDK53((5%l2*k7~K4jfzIl2jgtk)#tu(@KjgM
zr>Ajmc4b8)=)(~n<=jz;9kq)-@C&MroNwt4p7h~-y59c6?K)}~eK^9SoI5nwp%JPb
z-%Fj}u0Lqfx$}hvZxG>8&aE%3K45c)M$iZTo#WNN9(m9-UO8VW(ng0y(1#;D%3G`N
zZU6R39U4I&%tLbAJ>WhEP2;Zf>FVHbp46cc^x+7Pa_&;3jSh{3KI|Oe{X2~3ALzsR
z!aq7R68bQNM>%&X(ng0y&<A_c`#_{GIiKE;R;$nm`fx-OW!(B|sCte@&<8zC52xo;
zWA-EY!aq7Rf<7GKt#a=0j}DEX4>FZ=5i&?M!r^@39~~M&ACB-S=bmb}OMT>=Dk@PL
zfy4R6+C?L&k3@Ksb59!V&<A~BFB+r|px|)65g`p;58+YpCU^KpMI-2gT+mB6$OXNJ
z!}-ELDjGo_j_@dNZTLq;Bj|&gB3p2fW3mW`^M!v@G=e@H;Ze>V{xL%%lYcDEIT)#_
zMmU@={9`8n=qCSIAUw*s!$0O|1bv|K><5mM{cyhUk2xAaAC73F<R4ROnOmU|^ubt_
z`M^0dA1xnqEBQx<J{;juEjM#3`A3I77<03(m;q#6IbZlkFaPM!ha)`7xyPA|{3G<?
zv-_CIEcD@g;U5FFoAaxzQIB%&X$}|Fjt8dBnOmjKoo}4O$#JfS@F;I>_{SjssOSSP
z&+&@6RgPEA7ydEGKPvigghzR6!#@W3M@1jZNpjp}PLkuG^M!v5@{fu>9N|&Ut@pgH
zK4JLJtInH!<oq3bpWWZ~*W>xeBOkeW^%hqjzPiian^xJ@x*tM3_L}o%Z@sWZc$BwL
z_XmY|-&N<$PZxqcR3N$o@tJ?$9KY*)8ry|9<(l*6Z&#}v;ZfdNt#fym<A(cOdTM{a
zx9`}~TPhHaqdt82>R12a(BN>sQf|KUal-+ZoZ7$o-7UhSoLl*L%Mrtd3%}og?+rWl
zWQ_`h<6jRse6_#tVZq^i()qs}F+Aa-@AtoWU5oH2Z>`Q!AGg=A_4#)i-g5emz0dA{
zN9WP}kpt%|<^7M{YuNLYI}MLMtwng0w^rZI)bDREp2LNp5A3BM^Z^tc&ZlbE@4rqw
z?^dfE;ZfdNjk`Kd_tNny|GZ;QF6bp3<bvMA;e2{>S;y&<I$l51(RGAJIk(nGQNCA|
z_ZKLG)D+o*gPJ0Xa5!HnPgcI~qrC5=DshBIId>`NsBYh>YW(g+JI49P+Bv}HRjLsV
z=TpzWciwz1#H&>8em>w)&Mn75e?p)Sbe;V`bJ-7mDfzVa=l9N=oh!suYLz3}C_iJ7
zHYyrHAB<0#kGsEE_bC?nX!)d#ibl|fBRs0*mNqIHK_85>SywBcKX|Gu=hL?{FZzB(
zBk02s9_8G+!t(q(b!a5?Vdosb)-~AFhx6&3xKr-bp^?ysAw0^tOOZA@G=e_xz0~<*
z-+9oabLT5X+UU>-`f!9t`MF3b(ng0y&<FmV<JDOk51Pg+=PO0p=+FrIaD+!Ww{(7&
z<2p2gKA4B(xck3<cF;8LI-jlx?|fW`M$m^NJj%ICkv2Lsf<An9fA!7B^AGgle0tMD
z2sDB|9N|&U9sbdw5%j^H^ga;jOU|b^EyRIF(1#;D%DKZoIy8bl=wW&|J*OJ8AITT~
z(V-Fa;RtV)bBBL)Xas$bsho?DL8=iB=L`Sn&<Of)ghx4d_(z9ECVedCR8%5=a5!K1
zM~6lxeJl_j<=o*P9U4I&^ntyMyYvAR9L^X1(V-Fa;RuiNbCK|mibl`}xuBPDkPCVb
zhx3JhR5XG<9N|&U9sW_#2>PI=$QB&rm@LBKeBmDzji3)lc$9O8f6UNG`bVB~Fj7;E
za5!K1#|(|6e<Z@AoICtuF8}Dz2O7_Q;5gY2=L`RsqY?Dsh&F2N3`rw{#u52Phdvmq
zG9Nf+=A-3fZYBTd(1#;Ds^w;GCI9Hq2V-v56*GXWE9VRU=;a?B`f!9tId}NSApaQk
z!8P4Qwc{PB59bU27~~(LJ`CYe&K>?S$Uj0Kcwp+BxmD`i`NBU2`A6u(5gz5|BH<r{
z{G*}|ygbJ%ygbJ%=L`QB<R2A%IKrcxJN#pie^m6roFvCx<|H`|I$!w5ApfZ7!x0|k
z+@*Z$M<)#{cRYWlk*oOZKDR&kfpHE5AKU9AvQ>^iT@zw1iY-J@f$Yg`5YE>iN=AY3
zsN5bo9(jxB3@0D|a?SY47J5qsvM0B*RdA2}Fht3y<yJYj?$Pgm)NuR@e;6EOjS2*j
z+x%|+(#T<ml2Oa8N(BF}Z_9pY@8OMqzmwh^C|jsL{@9b-{BHiz$YF?*QOm7L#K^(#
zgMx!Tu$Mi#&F|(fjU0w38MWN1M2s99CwdPDxuBQq$!&f&e`(|}M9HY-RwZKOAm3yW
z4)RL2*pu7*ZvN89VTh7Z%dJYp$RTaa(Fh#NvwJF#-_2hdISe7t+F64~<@UreTW}O~
zmi@@>@ps2xnmBF{k$jFAIn?vII97%6Df5xrAX>il5FVA=BZo3r7dbDCvsqWU4TtmT
zddWHnkIL<lLzUPqYPaAQR2_SA8xH5wia_fiJSw+G4*h;d#VdUU-=pf-liP4OU#Rv5
z5gwJ>BZrPtN5#W|e^Yhr$!$2CFGimmM0iwgj~vQ(N5#XzJcO!aPj18Ee0suk9fU{a
z_Q;{S?WlM-e0HDPa5!I#?Kg<<sN8Ni!lT$j)nz}}hQs+L4iH-$;ZeCgax8jY=?6We
z>e!Rp*^lI#junV48MWLh=bkbc9At{B1CiUA54hQG5L+^8xmC_RRbod!s5t&W<Tk&X
zzcg|zsvX3Zj9PA0BKZGFQNckU*vp>W=6Ca#MvjFBgCnDsTa}2BV=*!W2gie6vM0Cs
z-TbAI!w_3CYPnU37&#UrTyT(Avc;a<=6Ca#Mh-)4$*AR4C1T`YZiPn3$MWo+3gmb5
zmqrdlY;lA~<@UsZM+FBO&wk|g_`Boz&ctzph~#s`$iX;*M(71&RpukNL9~49Av`L#
zM-IkA#zi<7bF;2;8xH3S|F}VfN9Fd&!MIMv!-02Db?nJ)IGiu4{RR;pmD?i+{zS#Y
zfd^7`?8$97oNt=-T@T?=xjk~=vs64BcsW(ap4^7R`KH<H^$;GF+arha&73+A%t@#^
z_T)Ak&KG0i4I(@$w?~e}oran2G}N5Y^SgO>;p9IXug7xk!}fSL<6K>hdhLZRBF}e-
z@)nMFOdhYh8RrPW9;yxw5S#&naK7oT%v^V68iYrAYo{}>{?2dTvC>1T4i3(?I0Mdp
zB;Ry)$JN0G;jQx4PG^w)Pv5X(B~w%#9Go?B1`LPu>1#K-8q8JB-@dLzc$BwxJQKd3
zuB5)_^c^b|#~&PdcL5IP8_%{c8n4zi2#@mC>RVg-{lyjGN*~w@2Ymnqhx1KW4Y=yu
zAUw)jD^2K1?c&OIB^UG(4st>7;c!0POVx4u<MHZwgYYPCZ7Ez~x!HJyzEV?U3l4Hj
z7U6Kd=_(F)0UCryd235~l&;h+?lee;ehx6+U4X;+bnR7FgSnf5J{;ju-de@zN-cLx
zP7?x+XFt$f_Jd#I-V1!wRW<I)G>A6JTgz1m?);z+#;VK**G)1XEgx4t7I&8#ghzR6
zrz?;H`e4k>y5jmz)|K;(S2@uL`f!9td297_w5}3yCk}n!9g{v5cQ!aa&Zp<wbtP(X
zx358Xlyi?)gwcqu2v<BXb<Xv;)Op^`81umuwZ&b_2H{cO+VRRZ8qt-zikIhj#dXCT
zubfY-_2|kV8qw9k2H{cOTCK~etB+_zS0F2Mk{oxrKAPj8^G$aFxQpB%Jjz=;-f2K1
zx{_LRMvn&b?gGch`Nq2$Xas#Y!lS&k;~f(;f<D+o)xiNGeaZQzyD~HM;Rug%?(h%p
z{8aQo52-ph=sDGx{Ybv)?h;oA8-%yYxu-i=6@8E?styh^NHxOYeBmG5C9CMe5gz5O
z9q+`Uk<<rO2S?srfW!I5yM1UxcQZmCiSQ`r4*%eeB3Fb%wd{q1K7fM5`KG&;Ty<^`
z9tCf5hktNqlq=gY`p`=_$OXNJ!}-ELDjGo_WG@jO#SzcehJREvf<CAzvIPe@CW~-5
zU-$=ik+}=dAUw*s!#}uF&7FpDgFNS8q^26-aK7*l?uK{iABpfNZ*BMocicPlfyT2R
zXgvGjeA8Wh?#eWXHp<UfxGTz?AN0Xkl}DGcD)Z6uEuK!`?oxyBC~xib1Vo2E7;`gI
zjJa7?&KLgC%Rf5w;RuiN)`oxZBuD52@5qYFy9?yQ`NBVVS|s%02#<2^@re}qN9Y3&
zOr7I_sdMMkDa82HOX$N99_8HQlQr@W?u^FlHODKwJjW~N3;!78AKZOy5FX{+;U7F9
z#2wz4h3B}-oFvCV=bN5F;x2N7@F->;+1l|<i&yGNS-ojdeRiMkf;{NBL*jXKp8aBb
zdYfbRla&^c&wCN&Eu7vMnSWmh_N4cL;0eI^-NoA+BL{D#yjgFI=&h6n;jL;~%ahAI
zdB&5=mENX@^RyXH0P?#HA5TZ~bX)&}`?Uy<^45+|(49CwfnCX3&P8}?jwb-&a6VPL
zo^5~G_>^~p@F;KX+S?p@wzm50KHsi@!}+GC4ToifM|o?fCmMMozS0NwqC5Hk3J&L+
zo|@#T`v&1r-r7>CjuTJz@3QmRm0ZwEILHOPhr{{wRo(IF&p}UrHVBXM)~a^LCro)m
zp;A+13l3_EEW+V@rMzu?ij}t@8iYrAYjq`6PpU88v<UxLI|tB}K-DZ9&ZqBS>1lZ0
z=J>j*o$)*o9_6i_p7<Tm2fEIFpt<Y^zm$CX&ZwS-U%ZvlAlj%k?&-;8-h8QyPni#%
z8_j&Qd_4JyM$m^NJjz>JinPHK*u2408E3Puc>Xo(%K1u>HYys?Q>+cbqnvx~ZH}mR
zJLfoQ@h%AYaK7nj^-&*&@F;KX^hA8<1K&%X^E`6u-1(-b?n57r@F+hQ8K3M&Bj^MF
z&hd)ps&l+@zUl1&-acy(9_8HA8w$MPR+)$7xXbh5Iqo{&_!a~jK_8CrC~xifrUe?&
zljPNB_xYwb`Eb7RZ4NYoJ{;ju&ON;`GD9EiN$&$O`O@NT4*15mQqYK=x^EEPD(9Zw
zeBsIEIJ)WK^qgwUeoTD%c1=Ye9Fs(NtGu=28$4(PeUPb~i;zL85f0}I|LD*N`f!9t
zIrmh%JV{=Cc26bp2Z!@bwR=5;M>+STL7s>Y)v}jymp*`k!}%r+^3;8U@F;I>_y=!>
zb?Ad5OE2Ld7xW$u=L`Sf?K9pUXb>L75zp3!fAEGIZz#klO19u2$7B%>=L`SfEjr$U
zXb>LdtquR+O+DVU2sg-cju@#`vv4?H_y=zTcIh99@F*0JtquR+jl&Lopz+*;#<L&J
z7yiLpjl7l8Alj(4Go<OwN#1-3cg}n;R%Jd~KIT@uU4uRt-4fwZ-rDdF-T>tdo|t81
zrWkXxuADFYqnCg1mQaK6DCeH$K4Y~no~9z_`E~{QaK35ob3KGd;Skx{X$}|qzynj~
zcwp+>`KCEs=)(~n<>w;fo5k{viazl29DVTe9Iu>jdV86-&l-eBIrsF2G;g@YEIh|u
z<|H`|I$!w5AphVkx(4A<-r8xUDOL)pKD*D=rGEV7gJV4&-hF3#T+QlNT3bZxGc<^N
z-<>FL;k07bd?z8;limk{6%ar;-?%E*b>pfk4Wfx^&RE{Bm3kH{w^VwY9?skLtbmaH
zNIqR*5rWlWo^|cj4Z>UHt<_noR_a--KvT(D&P9&Wc{D2^z~Ow;Dg~@^(;z&`TdOa{
zYNZ}l(s}A-S683iv-$yla5&$z+5)TfGzgFK*6Qo0`h8XeI!y@rz+QAmA3(w3eAB8C
ztUA;nJjz=ugjVWdWuy0~=j4K3!a*+RJsi%bs~_X)7sHpYY!M#it<`sa#uYY(6O=(}
zifq9_j>#e%&Nr>%!75M<!lS&k(@G(%G!_1_b`CJEE(M44DTAt8R<jCy*ja-|d22Px
zkeXOA3w@yL><605e&7sLKYXP~8?36rs$31Cjq)=VR;55ATA89UK4m^wHzV`W^0D#-
z8bKe9@F;KXv;qw)m{rEvtSi?4$hvaAX_XsRNox=u<=o>+I%ouaFngiu$T=(Lq&}Qa
z50z+Lk;Q6n4Z@?mwd0CFXas%Wd#UrYbcHr`?tIg#L!l2xc$BwxT-gYXpbz{z$1B#A
z$??khV)Z>%PwLT!BRtBvwL+{`pFtz&gLz1fyR46s<Dm2D9=#BYRRkM^M>+Sj(iAHd
zR-fJH>QdCJ^Np)np^?ysoi%urw{}`FYlc49limj+eaZR8Rk_dz`fx-OW!$=Esh%%Z
zZi%Ct9!}4x#_UJ(g@3TRSVbQklSFu{oICu370g(HCh98ZB4m(igv0r!Rnl1Hra^d=
zbBBL)tfZ6rpz2Z|R3jYDH?8)@YCWNkM0gZ+m8}i`U`4qOeel2Rg@ZnTf&;GPn^vu3
z)u9IAQSc^r_y;TFb?Ad5OE2Ld7xW$u=L`Q}^}P;#IKrd6wc#JE@YkUaYKm;ZL5|5H
z9L^X1!775R0@WZq%DKZoSgDYerbho*oO3WzQ;l#qU-$>BA+r8egYYPCZTJT(E_Ubx
zjb}g5c=p5j!arEmkyW`GL>ndlm|Dx+3XR|&j8&Nr#;VLm%g5Y`)y2?<BRtAmJFS4p
z3T83O$hu<8&AM{FX_ZY@Nox=u<=o>+p7M{-htKZuj;vSb3;$rX(9nk?JPL=%)`ou!
z@{iC59+=9(15@YDH?4Xa`f!9tIrq4-s{Es(54=1_AG|!rE9VRUVD;9DJ{;ju&K>?S
z$UiFjU`~?bE_0F`2c0kcgH>o*MX*75lymEyxDz%G_j~>6vz=bLL+7*m{JO?p{Pn@B
zn-^>UvR%q6h1mTqr_T;~NsI6(Z=vqE3-PYkpFV$t5bR0s1HlT={I2s&UnOBx`3B)p
z-r7>GJa*6Fs?Q$LpK{_39lcEtXSHEgfX;p-UnviM)}F&rTI>FR=d}oLmA6)oqpwrE
zPAl+lbL<WsS<AUdeq93&=hKxteF0;!%726KC~vJ0`ci>bg62yC9o5GlV-%}D^9P6X
zO{-1wwSorWQQq2VMe9K;zIXJ2y>R5$HQ;c*Qch5d{y<+{xO@rWQQlg8WlUc^dGYvV
zhmKs(OE|~{y@$j3N?EM_J={%4*AX7&tyPs=azuZrzECk!2B|5s1qV4Mi*Pt!$yLnx
z3PywQDCgGLuDV4d=wo?yPc_2fe5FVmb2NfJ9N|&k+ES#A85%(!=sNqsdfV9#eku8;
zRpt3ANrPyk{ES7~sAvRzFg|5IST{WL(eg<f6^)<|M|hOCR?mKDb#XL;J{V`SuJY>|
z94qIWR{7^EF%80_oLkrW_9%QwAoPJ>q(1m^z@(4A9lz!QU*T&ve61k#;RtV)w|4p>
zLg)kEOP#aAJJrbVI^X!!1vG*_9N|&Ut#5{D)p|6dmFqkFJI5>5)z9(D`KGT&^ytG8
z9_8GnNE;m*K_ARRa@>8#;#(l}!})ZLOo+u-lp2IbId}L+heptc&+hZ<8gMvY_(z9E
z(1#;D%3E8?<@!>=rQ;W;qE+dA`E?CAoG<)?uafXps|Mjw&OLtl<cs6V^>K95!|6HI
z2#52P@>(6Mw`q0!2TD;KlSFtF{m9mae{^UBeUPb~i;zL8G4qjp;U9b@hE@I>gtyAM
z!$0^E4PO%Is5t%@6RAZ0;Bda_Yd3tYph0+)w>JEPFY0vo2Yp~K9Qkz(IGiv1gRk~<
z=)(~n<*g0>sAvRzaAfHv92^^Z4~O%GfAIC84t+SnqntbZqoNV?L6&l^&Im^}!r^@3
zAACiLuV6F?k8<wtkD2_VOaI7o4n}IK5f0}I|Cq@?y7Z4kc$BwRI@i~9(8wODb~K*-
zz;UQXIGk_#suf=)X%KDH+8L7W0*qhI;>#!D&Y2I!s?0~rCvEg-1br~NCBmb;wc#Ip
z0jxtGjJa7?jJZ@J{cyhND`R{mra^d=bBBNMC9}{6-jVviJ0^V?U-$=KOACEC!dvC7
zoxaEx`oIHI=lOLF`r&-zSKs6xp$|uRlyirF4Dt`Y3>5R-9DVTe96y{d{DZI8RrG=T
zC&HtgJN#pie^m6roFqpw<|I@j{cyhU55B_3SCkrrM>+TSrST)jFU58~yU(wSf8q;z
z{$TOFb+$|Sg}&|lmZMLcefRJU4Z@?mg?cJXh}-FF=HC~BJyalV2f-IG`CaGJH$e4m
z=WiT+;(We@@F;I>eCwAlFZTC-^oEY!rib&jMZSQU{Ybu2cGb5Ux9A(42QDGJRo+@<
zP+w~P;cMR1|NA30bUwT1tC98KLxaQlN;y<t0AGCNxIuW7w^pO5zSR8Ci@w^wZDm78
z_3_8C%CC#l59cf8uY@>Chzpkx9_6jox*v*ulu)Ze`0SoOfP%yMbmK>ey@k+MLmPxg
zd297RqrTMq`|AC}b*#t*y`&%Hg5Ja7e7YK?ub=)(NBDh92#@mCYGn%LdtYU7lQKw6
zku5mLF<FGe`E<u!`F@-b_gzAGl($ybWK_3{FOA39#M(K)_;qnOoKGuSt8UQ<`tb7s
zkMh=*B5lmj2>L+R*^m6XI2`;E)jw`a8#6S5J{-|T`56meCFNW3=!5Yo^TBshGaoIV
zv{BIr`f!9twcOH1MI)gPI|o>O*Z3`T+U|U%@C8>iqOZ7i%oC_nGU!pxt?O48eYK(y
z^nqWbKKLpv)d+|4Y1StM8bKe9@F?fjGu!&wEgA`Z*f|H^<4v7AUn$Z?hekplhVUro
zE=Aht&<Og#zjM6eyTUnMIiIeJ>N~$^1bsNdqi|_7!Z9pG+UU>-`taF3U**nm(D|fV
zA<zi=aD+!W_xPo8G@>slcRst%uZxop=L`Sn&<Of)ghx4dDgUA`HQ(i_C(a%wgwO8t
z>*8=Y-}qJbAN-#aXZq@NgYYQl9={cjM$iX6Ob@5$RAcso>SsH?4%{ui4%{HTRo>e0
zj}DEX51-wWL8_78b-wVAZt<1l2H{cO+VGDKjif$uPDLg12Z!^8e{^Uh^^pjV^45lb
zR5XG<e0EPCK*8aB;U5)^pbtlQ6uikD{!!5g`XCqdGDik_4~O%Ge^fMrJ{;ju&K>?y
z(FporB+GLSa!eNCAos}^{!!5g`f!9tId}NSZ1JV>aDzPOV5FuR;c&k2kD2_VOaDlO
zM>)4D@u?@y(FpqR**zN1emGw#(#Cx8RrUtaMy;J8>G?#BBYZ0!eK1yKJ{YU=__Tb|
zMz1euclvsEgYYPCt>%fETggAV(1)D^Fy>}mIbZlkFaPN94?iF9C~s}}M=$^Aa(<N+
zhj*ktoG<*Nmw$9Qze<EhId}NSApfYLkHyynsa@t)sdMKG{}|*SHS}Q!kHYDaJN#pi
ze^m5=m*;qem*;rpeBmF1{G*}|M|hN<i-dm+@{fu>n3GU-^uuTR&KLeM$UiFjaD+!W
zcPYob@0P<aZ?$pus^d0wKD)oi6?|)O^Xd<tb<pbWUpRF2rSIGn+j=TYh+l8Nakl3Y
z!lS%}rTkoo3;%HA{2@a4>>k8(4jDO|Pxpj{c;faO=Pz19c$BwRW66(yclh&f-@gBk
zCvWQLZF=~p|NYR_?`>V{NAhX)10nAGo!j@{TtaxOytSo#^8?rQ*Sz|q{$Wqp)RDEE
zi+t$ihpyiIcMpx@<9zy-sSvk%%}M>9O9+qh)|T>#OD^y4bMn^ymPc;ts6PG}qki?M
zL*tKpIG?UxU3z(c)=69YOO9v}9_6joD5~H8yU$(T|A!FtfxU3h2T*W0pRUR1_unPN
z&D1JKc$Bxcl<(;{U9R5Kb8<m1;UE|E9uDU#<xg~+z9a-k*AX7&+{(wdZr^XKEPhuR
zWK<(taFAoN2#51YQOft-h4|DG!lRs9&vM>s;~b5k5A?8h4zPKZYJ|i2wCarN7LA|}
zM|hOCR?mh1;l>#nK_9qG_5;mjKlr8OlRl(BG=e@H(MEY|rHzV4&<EpF=3~3ZY?|`X
z@<|&Nji3)lcvQ<RZB#Tu1{r6wu72~&L#DcNK7HHiH78XxLX|keqnulBfS$ayq7n3g
zU!*>M_x*=V`fxt2esR*)ibl|fBRtBvHHZ7$<rR&f4}33m{@2GIGU?p;bO%caG=e@H
z;Ze@5Yx*C!uA&ioj(_KPwd>x8OyiaFl_G6aG{Vt!ghx4dDbhxVM$iZI5ULJ+{KtI`
znZ{k`)6*hCpb_-p2#<2^@Q)6Spbwwjf90;chcS-3&KLgCp%L`q2#<2^@Q)6Spbz#?
zb=i;fCFcwO=+FrIaD+!Wclbw#M$iX6Ob@5$R3kd4g~=EG(V-Fa;Rug%?(mO_M#vzU
z%DD&`q#EIHzVMHVMyL`;c$9O8e^fM*`p7vImB=3)&KLes(Make5gz5-;U5)^pbz@M
zUh+X7K*8aB;U5)^pbtlQ6uikD{!!5gJtr6R5)N`f@8NL1@Q;c{IJ%DTDCZ9UsAvRz
zP*Y?J4suKu;c&k2kBUanha)`7xx+tZ@{cb4BhNX|2Gt0M^M!xR<R4x7M<P7Rxx+u^
z@{bOE`0O6X$$mIr_{UuS(V-7Vv{7qkNIK`#+)Dn@p%2EY%m>by`Dpo=Tgg8<^x+7P
zYPp$P$v+~4jJa7?%mA{koG<*Nmw!Z+IKrcxJN%=Ue}q2pj;uIlGN}*e3;*clAE6IN
zc$9O8fAsQ?&<7rvI%jT`I(NSCk6!)}`f!9tId}L+FaL<1<K;PCF}KR`%K5@Sdih5j
zT}OD7bBBKn@{fu>n3LqV%bX;~LFWtq7~~%neK^9SoLf(~z5e6<A6&Z4?AD)nO6Rlt
zpKKV<KlXXxL91`w?a<YOKE7#{ZN2mJ&mZsK{o!q9k6uD}l($eLwGcnKWSjZpgkTR9
zi0(kVXZOwVyUth2OVp~PKeEkyyCsB2d24mmpdb3T{?pgyFTeOH9lfOj;n;fjLs#GV
ztdYa{N_m427u^8DqrA1H?Ej#1`@5aBTmQ}BDIHm(0^xZ0Ee~CN+IvS1=hJ&$Lflq}
zE0+)+<*ijees)Ix%kxg?KQw<z=d=4AuH3xJA31QoQvT_uXY_ac=n4JjmJlB0t<^OE
z{r-K#bD0qIfxYyDK7fM5`J_Sp{+Gn_(ItdOd236#w|ae`dQZ>E1-*oWT+n+soUfE^
z7Km>xAw0@ktGlnt-dW1xE0l9;ifq9_O_4=7oKN?i7YM4v5gz5O)t%8xx0$06^ug%7
zb`G$4m1=~;`AU&C=4b?cIKrcxyA)|-hDOi_y3T%}x$FnOlzgQ~8#6S5J{-|Td26MO
z85*JIj8B=5O{eO8_eDNhK51iyMwSsC)pAQ46^)<|#@Vc^S0Auxsw?Nyir+$@5vs%y
z9_6i-pPzR^MI-0~zo6>K`FHQMY0`)D=?&?Ro>0*U`f!9tIk!COXJ=G2f<ExQ)cN`U
zbI7D~=hN8!(=#d>K_8CrDCaIk+NfxRp5xy+ULEk2L#FY{`9=g9K_8CrDCd@H`=O!{
z^uatN$KAjA#39qT>wF^ujgWIkc$9O8e^fMrK74k+&Bf#S2l{Zn@Q;c{(1#;D%DKZo
zDjGo_>`Ctfk-p@7;U5)^pbtlQlyirFR5U`*>0x>}J*OJ8AITT~QPIdU!dvCs;U5)^
zpbs*Ya}hE~HNxS1;U5)^P$iD=C~s}}M@1v4kDOCciTuIgeBmDzjif#j;Ze>V{!!5g
z`k)W&MT7JK6dcYM{!!5g`f!9t!JFLS9~F(zb8<m1;UE|E9uDUV{|E&7aD+!WclgH)
zji3)|ifq9_j>#e%&KLd>2y*TSk8<wtkD2_Vrhnu)2O~Ar2#529f6U|`HT@$I9_8HO
zA9MLfMIUH9`+?(RKb$Z8V=n)w=))0h)Y=(RDa@_pAC;anR%JeL&df*4$J|Q(aXo}b
zwcO0D<R76A#@wtcW&l}N&KLgC%Riz@9N|&k+VGEF{t^1{**!Cv)Q9thfAsQ?(1#;D
z%DKZodih7_0}o7{Gq*~eJ74%mFaHRAIKrcxJKUg`e?-sm@*J=5@*J<6Z$!vHD*AAQ
zM>%)6K`;N9p%3OHIqoti$#Ky6Muhw$#zaSWlymECz{elc|K&ZlZupm9Jh$`N{bSG6
z+egyK&2N3s>P2tZy!x_lZ;EYk9C1v4%)VPUymATQQQktWp?~Bt{qFbNdcjMCU=J0D
z?m#?i$C1PN^j4!<b=5vwFF0cf;ZfdNt;F`Aqx%y-dCL65*FCqRw^Se;N1e8L^_Akw
zek5Nh7w>X(f6*nU%r9O-c&og%rJVkS1N$$$<cj$-e*WB!tWkk*T>1RXtFJnE<ZwPc
z5BKE*`%_Q8Vt&vP!lS&kS_9&tyZ7HacH91pYoFWs?EbW;ZC>S%emI|2b$syd{l}lR
zZU2rXghzR6OW9e!f9gYa@82r~ePA#Bpbwzna6YXA`Rr}`mECvm-?W7AC~s{k4}94b
z^KV>sV82m4Cl~Y*4st>7;c&iE9)Hpm^GAH?!2XF#2#@mCY7M$ioHE~k*Q5K7D}&S&
z*@A-{lSMe3uaw{W+$r<DggAc*;Zfe&QlyOw&<OfK4{PTDTDM;sfy4Q9MR?z>7oZXJ
z;RuiN*6QluJ-2Q^Bj^KNXFt$f_JdzazEY%(4QK>?IHHa6)=C>QG=e@DpE4iM6kq0}
z<&!pMXas#Y!lPPlX=8>)&<Epe*3~oKwP~s==hORxCtoo`Bk02s9_8G+E_&>?6^)<|
z{DP_@=TARP>$ffR;e4e?8x@V94@Y>EbE{Pk*}b9>^nvfC&hLEUrb*|{SBkVz(Fpo*
zghx4dDbhwoBj^MF&hhH3r)`?XE9cYuwO=~0q7n4r2#<1Z&3AV_x}p*E!8|0#-R|!;
zP2;ZfX=PO*&<Of)ghx4d_(w$}=)-6C|9TYbM34G#zVMHVM$m^NJj%JlKPnnQAM8o*
z1ChSueBmDzji3)lc$9O8e^fMrKImb3I6bEt(Yf=5f7H;(62haLJN%=f5%fW(axOv!
zsYW=QFZ`pT5%l2*kMh=re^fM*`p7vImB=3)&KLes(Make5gz5-;U5)^pbz@MUNlG_
zK*8aB;U5)^pbtlQ6uikD{xL%%=!0C)OE|~{y@$j3!artc1bsNdqntbZV}?f12Q@{u
z;Gm|+A{@>a{xL%%=)(~n<=o*P8{{7~{Ugsg&<52Ahx3JhY><D{^p8Y%lyirFTp<6b
z=mU*sKX9Dvhx3JhTp<6b=))0h)Y=)6>PmAf`A0<`j8&NroHO&$^2OY0m!m8G;Rugv
zxtUwZKPviQ%+0!D29R~-eBmE+`A0<`j_@ew4*%%oAE6JQ-Qyjp59bU2=;a@w4@Y>E
zbBBNQ@{iC59+*03Zk0NBzVMG;{t^0cghx4d_(w1Qn4u56JjW}%JjW~N3;*clA2amf
z2#<2^@Q+^pF+(5BNpjp}PLkuU^M!x(@{bw%aD+!Wx7Osae?GgX634&s+4HsSdw%ZB
z`JdnJzaRSB*~#5{2ydYgXp=ouAiBFrj`QJI57BZv3%#WRZ=B=4PyXlm?eDPlg3pR$
zJw(gx9Au3Oym5{@|Mb-P5x3Z8{;%R#57BZv2i3>ljdOhDf@jU&@E?26UoMXI5G}WJ
z&<FP3I0wHE4;<?uT5jhc7xePRIXF%nFOJoEh?d(qs422_;~eChyp!|w5G}WJP_Jv}
z0PE>P+BhGL{8k+6AzE(dK-V|j5A>HDH-=bit#dFw-E=;5tTyDaS`X23I|t+JP1luj
zKFgf1hiJK-1HZUweW-S4S-a~YT5jjS_ikF}`u$mIa6LrJ?Hu^`O^;VPPO}_K)<d-1
z&cQt7rpI07dzNG3dWe?WIed11qhtFFji8V95G}WJuqVCmCON{R)<d-1&Or~;!`J6n
z9IFkCCF>ztZs#CVITu-<W0CV2W8!*<mfJa~IR0*YOkC9NjIn(^M9b|Q^ntxM&au#7
zaIA-Dxt)Vt(90X=Sd1lctky%c+|EHwk*yo&Sd57?a=sp-<#rD8oaY?t=X^1?&weY8
z^$;z$bD;5??#IHTg5$;zYprz-#;Tjn$HHslSgnU>xt)VC_onM=F=vRJuZL*4odfT<
zX?-l_K2f{tAzE(dzyohu=lByEL4)fdT5jjS%WrzT!e`M4W664mmfJa)lic*U%X|Ti
zppW$sE%(bmf9CwFfBD}J{omQigI4ug8UJ=!x$iZ5#r7X=*k$FRzujwP&sXfT!nU`G
z?KMOjwf6s|EI7_z<^VB&<|E>Foo{`PHp)4+JoTUF_qf~E3qHQok3YTF!&i2@a_``9
zzV$iUDCc<Q&rh9y;O5)R&syR*`<V}4x$i+Ehx4t^(MCDP*Hwv^T(bB4#Y-H2@ui2a
z@J9}uZ+(t7$~juh+8+7<YRvQ6|LeZ9T13lel}CB2ycSvp-$O3w{q;FGvK(QKu17ha
z|B{|IRK8DE-j{PuO_4=7{FnY;S|x3qPyH^}E-j%N;c&k7Ioc>cy6EqQIM7`7!}-?d
zXrr7%+Sq_bmXFmZ@3ilfkCv}xS^KRvYVH3^Q3hw3LC@5GzI)%PuAHxBa6}vB9IC`w
zR-z*=yXVG9AI`TvM;qlFU(?ydXKz4*fBTA!lg^!QeU3KDIdq(6H#C;~`0|a@c;$TS
zbF@*;A#Kdi$a2m<R`#35UFTb$qm6P7)$I(8EbHSx|7gF-Kb&uUjyCEBIIibWAkvqd
zZ+(t7$~juh+8%mNHC}(bS{&UC99@s{R(XY%vdH;t*+0l2)d+|GQrjl~tu|_{wWUZK
zGc>a7A5<cLa5&%k9Bq_yNE<UWvdlprK*8aB>vObG&aoIvu6L}+1-*yE`PS!Xqntz9
z*gzJS^+8RMML3*qeU3KDITmC4Y*`;<lWK&+`PS!Xqn0@EsOxdyIN1;9Tc4wia*h_W
zwg-RBe0YUU?i@$gqn2AWqL^8Jd;4wX=ltZ<*%?c7A7%j5)va}`{FmA`IeQzm*4p(E
z%w$p@{!68N?1H`LS6=e0*$dSwkMdTv2x()EMwa!#+$weMzx4moDts1=ERQA3t#Z6_
zzV$iUC_h%rt<cD_KA4l_xa)lDbF@*;@ubHc(?5LIts6e{i{}nLyWjM=J;Fbp_{&{Z
zK7G|5D~H{A-xan?d4>@0-R1fSkMb6l@-!j#-F52)9~K9Ds6ccF;%)!6cl@sNmGTrJ
zj@;$?2#@mCmhyidcy#|SpFL&%-`72Nptn>Y9B0&(l@FZv@ZfMh-BA?c@K2vI|G^T%
zqrA0Rk4(=K-S>G{%uo6Ga|g0U1;TO8=T}y~cIIBe;e5K=Bg6-V_{SxLM|o>Yx#B^)
z_n*_#RKvB;jcTV(;W++Hd&VC*a6Y}ADa0FZ0O3*IT79WW(VrKp7l+U8=>sS@oKN3N
z72<bxyFS9BXlb^#l!vI-UsLa&ubz_&dPzUX1-*yE`LxQ75N{WPqw5Hd^46B}BISEu
zWswZ}?4Ft;i*PufR`^rCzbV9bmk=K1+@(kx7oZXJ!6=jG9BfmKa5$gd&Q#r^k!6HO
zd297nqx6SH&<DDv>Nr+tF8kqpT2)2*LnF(GHp*KoZOqUJ`e1y@eC+>-jZ;2aK51iy
zMmV~T@F;I>DbmIajgUdcS*nhH-1Vsarn+)It+FWu8bKe9@F?fj8@YN?1dX5%{30vv
zMMwPgqz~sSMcSxnWEtU6&RvSMQPBwc!1q$;hu`;+lg^z_Zz$}xdqpG52#@mCmLhFb
zG=e_x?;NlG?xIId<CXK3B5hPO!qIhvM|o>iR}Vb8q7gF4JS4~6UANtT8h4#fZ{!Mr
zM$m^NJj%ICkv1wCq1t_R|IYjGKlz9Ag@05uvW)O3=MMj<XoPBKPkJAS^d;vD|EOqW
z8R1dh+VGExM$iX6Ob@5$R3kcfzVMHVMmV~T@F;I>_(w$}WROhdT!ai#jc_<$_(w$}
z=)(~n<=o*P6^&%I=bVa4<PQ$#3;(ESWEtU6&K>?y(FoN}AK05Z2L*@og@05uvW)O3
zTAHm5|Cpf>^g%A@B^=~}-owE$OTO@r85-f}I>Mv8wc#H#G(rZcDYBL0E?I=b`NBVD
zXas#Y!lRr!{9^+eNjJ!I4n}IK5e{lN`NBUo$Um-!@F?dF|F}T@QK@z`p8dda)_!-5
zFZ|;I`N#DT-YT}&TFcx@{!!5fV^!t@=gfSxe9W!nA8~XY;Zfe&Qf`0374wU~cwl|4
zGRT-qU2)FA3?TFDeELGgi?5g;r6;zYyM*v4=MMkq<sVV)ct=(oGnv$f^Of?>r)}Fm
zN6!v^Rn_hYkMh=rfAsQ?sCJ*-Gq*~eJ74%mFaNk6!lRr!{G*qD%+Lp3p5qmBs~oSK
zFZ`pIf5hnH2#@mChJW<(j~L;6M$eoi$6e<O|LEl(`Zn1()6bFGqnumsdFd_IJAG`M
z*<C;Jl)-2B4?VAjf4q912d#Yfzv{{(-@oq)+oe2Ph-Y88&FuC|2#@j>Y862t_W0N~
z^X-LT4;6^+Kz#5~d&loOpPrf&;=UJdGrz|Y!lS&kTBUgy`n&0kr^j9Vl!4w-fpFaF
z@!iTvuiYy+oKJ6`3311ld~JT<62harwOXlgw{!bz^p4cKho=l=jS7V0q=!Ci<+>lP
z1c&qK8zVxTD#W=<2#@mCY8@86TYHG!WPQ*4DTB}MZ~KfrR`??a&Zk+Q5I+{;ZA%D`
z^46B}97P`@)Z2xi5A3BM^Z^tc&R0q=#GeYWbqV27-dcV4RlVL>y??)YPA=$W94qC5
z-oxR1@(&>%DFjE?5gz5OE#=Y5_rEENhbV*86xo7<nj(vFtc*`jO}_YR^S=|~$Rz|7
zlC5&?QlyPJ8bKe+vwNx$4(F3TRJUjZeK^9SytSoB8#6S5KG1dc1I=YW_@(4i&ZR#z
zf<7G4My;J8l_G7-&<Ofqe9C-0|H=DJJ#YD>jTst2ACB;-mRs7WXoL(h&SqWx;cXu|
z)s^$<Nm(J#2>Nh@M>)6N$kkhVXas%W7gQZNzxtK?Px^4aQlyQFM$m^NJj%Jnq4yKf
z2>QVHQs>YAzmJ-9?tEJBRR}bKJ{;ju&aG$YcRRPD5%ht7=Xmu$Z+`SNUOAt>A0-4D
z;pjTTqnx`GX``YM^uatN$K4(7{FrInbv~`!A_N*iACB-S=MMj<Xas%u?EYhqc+BJ<
z&KLes(Fpo*ghx4d_(w$}=z~4!eIU}8oG<*Nq7n4r2#<2^@Q;c{&<8zC52xo;BRY4!
z@Q;c{(1#;D%DKZoDjFe!WGd$(WRPlv!}-ELDjGo_j_@ew4*#fVB=wPVDk_mbIGiv1
zqoR@2M<P7Rxx+sy8bKfQfxT#uK7fM5`NBUc8bKe9@F;kbJN%=f5%fVW=p`KFg5Ja7
zeBmDzjc{}w;Ze>V{xL%%=!2RfTX0ZQWHFAF@r8fP&<Of)gtyAM!#`&7kDC6G=NydG
zR3jYD7ydDmf7JAkM0k{QhkwlF9~FI|@$3hVll^eM@Q=CtqoNN-v{8P>!rV&!QPBru
zRptZd%zU(b%&p`f6@56uqgrm}R`QR?AY*RU6*GXWE9VRU=;a@w4@Y>EbBBNQ@{iDm
z&+eJYq&}Q4{G*qDggzYMQO+Iy(aS$VA9!HuoViu%-1)*kdih7_!x0|k+~FU+{3G;%
zm*;rJEGEY*=L`Sn<sWf$9pO>V9sbeFKW6BIIZ2MY%t><Gb-wVAUj8vdACB-S=hjNl
zTE+8>oi@(yecYzOXZPhl9~%Dg{kuGP<&S=U!^+?6wsGYuZI7$zo^|Jqv%4)JJjz>G
z%4fuJ*bW=#dkMiFDiGa)c<Uz~z5<8yX(ebOUUlb<^ZP6zJgRAJd5BgD-7>DfJJ4Gy
z5RM<dXT!=qFW58M-td*X2{C;0_Wk#tys1HWl($xIB59@2mz;J|f8-N34P=c9grmD;
zw({A7_lSNtpWf{g;?r7HZr>$@M|o>YIaaIOUZJ%=&p%?*;IsSLQQb<_A31P7z11kh
zJ%squ62harwWS=R=x+=4AtC4kd+7�oEhmGWdEuDSt)M|o>YxrJ5=ovZiz&~uyJ
zlM8wehx6$Rk~&UTjBBwr2#@mCs#VI~7G?2+${;mGw&({n^~b6$IGnGPn`t%D#|d%q
z62harwYn~Pw~h0wRE_skB`(kIsYYto`E-Rvb-V5RZ#n$Q62harwWUZKGc<xe&~^3$
z%~6e^5BXK{={uv+9~wa)j-Z8WQGUiEZB#U(^%DoiC#sHPb>0boJ>{e2lQt?EK_8Cr
zsFqvWsAvRzFwSOO{nfknkGi5C&R2@GQPBwcaD+!Wx8Akfvb7$eN_@9!7r&tD=*Q}r
zkDm15e5FVm6^)<|M|hNTm-11q6#6z*{0D`=_fqE<J>Y<-S?bF9^hIc`&-g+i-nNAB
zDCaIk+Nfv*ec<0YUj6xj2TtRa^XbaA5NJdzk~RpB^44l?Dy@ErM$iZIkQ{ff`Tq`#
znx!Akr}f)~KqKhG5gz5O)tjtZ#S@L751-wa`yMp;hx2Jw$M@f|LnG+J5gz5-;U66u
zK_Bc%@4Ht0pc+FT#uxt4p%L`q2yc~h>w3F-jz-W2JxmX$=Tsxds^P;wIy8bl9N|$7
zH~vx42>KvXITs;=R3rUxzVMHVM$m^NJj%KCg>F`YmPS$^Ij5o$`J*4s7yePvNa`aI
z9_8GnNE;Q6pbz@MUNlG_K&3vCuM}ycq7n4r2yYd<$sPVt(Fpn=7xXg6D|!!y^M!v@
zG=e@H;Ze>V{!!5g`k<!B797+R)d&Z<PrmSvibl|fBRtBv!#`$dB>f}LIbx(%&60EH
z3;&p*k@Sy5c$9O8f6V0{9r{4y*$*6tY7Bkg3=6*SkGcG#Lm!UtR<WJj8cQ^{LL=yd
zu`2U{bLR1B`IuXw5%l2*k7~J@TcHv3!I+zM#SDPDq94u|{?W@nLLZLsDCgG5pi%m1
zs@;E6CHm~1nM~HJ^M!x(@{iDmBRtBvOSxKZbEGQ%BjUgVQ|HXBs7C6_`E=hvh$q|t
z!lRr!{G*qDgg$(B&)h1<E9VRU=;a@w4@Y>Ew^rv-a{s$({atd-oFvCx<|I@j{cyhU
zk3s%X(T5{E%DGGVM}7C~9nUy%w(>U{2A|#Uy~pn1AD2DvAuIcB*?r|}TXlC-+Zs!R
z*jtEh3E@%R!cyKN#P^<Y;(P}oe0C4wu8+|@F#WFc>3c;&T=nD==X)<9nkYL%(#jP2
zQrb)3+x3?`dc#0(sX&g;{u>{%a_XyBg2VYr*-40F-q-c7UqX14x3-j1^li1R`bN&>
zkJvDfH7XEIy!ssvU3vPCyWns>saA*w3USR6!lS&krJQ@t?++K~OJ={b=Z3*&_uqK+
zY=u9LmGkN73h^r;u3c#n9_6jo3N(s7Td3{DK_A%5v7!&4;BY?O=M>`OLj2wm!lS&k
zrMyPHev*3s6ZM>2&`UVT1-*yE`AT`e5C;nJt0ja-d235KPWk@XJMTVxOBwXpJvBuZ
z;cz~kZNI1M|4e!R%o4(*ytSoB8*?-g`dFIXQ;l#qpLG7L6XypD@dWf?2#@mCYE2HQ
z35}o+bWPP!S7<K#!7nABW_?dSafU|Fha=i3KNpcUDjGo_j8B=5XCI?$y^DObe9}fm
zBk02s9@TP78x@V955`&Qihi87jph!Ex^h0Paw7y9K_8CrDCaKa_4-QF8&!#SP$l9Q
zc}{+pV-A?~;e4eWwEEThG9i8{eK^9SoVyfhqeCO;1K*?SsNJV+I%v|l^XbcSLcB!?
zTIC3ja&A4hr0;B@5%ht7=Xmvj4<9^@SI(#P6NNw{=)(~n<=px@*E{dtp%L`qv-=O9
ze8@EJI$tT$Mu$ewha)`7xx+s?G!pu-bAXTbhfMzAe5FVm9U2LJ7{a5RJN%<VBk03t
z_aM@joG<*NLnG+J5lxhFhktZv1bxuM^l*AkHKKFp3;*cQ2>Nh@M>%)+M~6nx2brSk
z&<7c$8sTuh@Q)6SpbtlQlymE8LscRgNqyv;ib~`U4(HQ)$U>lz)JGyb%DKZoDjGo_
z^ntzTgFb+Q!}+u-uv&yhXq6*83f|-n|EOpLeUJ-!2?x2L_i#90_(w$}=)(~n<=o*P
z6^)<|MzTETpr*(o9L^X1QPBwcaD+!WclgIl{?SeTF`WbKvqCk(;e6pAGc+>!#{%I|
z&K>?Smw$BV!)Nz6PWA(5NWSoox%{I;AC73F{EUUUmHeYaAB<I*51cde(eg33l7DpQ
z!x0|Uax=G*e{|@BF_*eRAIt!<uADFYqnCel=)(~n<=mxwNTc+7Rf)GyCHm~1nM~@#
z`AWIX>R0=>3z74ytWl40?(mO6{!!5f9!S+uyLe#g-1+pxmi*}-geXf0kMh=re+=@E
ziazl29Ix>59Iu=&{9}-RRP^BpkMh=r8w_V?OuR@L^cg*Kk{oxPPv1DyZ0vu9_}mi0
zqnum!&rjGmoORaev#VdS!{D>~x7>BN@Q=4{n613`sk^Prw%Kom?Na_oh+Dns^x3Un
z+9Ev4Tc~~r@$$1ypWjjl_E3T74#dSj-Yb6B`Sk9R5Ier<^!Zh4l_NaLTU*Nej@xs1
z@~4jIKli*H26{^c!tv$<AH4FCPwx>N&ZjC7V%JX|(ckRE7U5COUCJxZ*k!oj&gb`A
zkKJJ)Yg8Z{xB1TQD_VX#IGk@poGHXdmk=K1tu5sV-@e1}#y!e#<S{!8KD&R-Z8xm&
zM-H4%elEoS2ywS1ghzR6OZkYRUntbQgrE=Xr62SG6dcZ{uiXgo03n{RgzzYDZ7CmC
zuV13$b)I@oF6bp3<bvMA;e7HaAzmj0N7oS^<*n7+N7<Vzi(gdEsVTAr2Q{@-wFQUs
z>FS^m9~0s;O9+p0Zmq_yDx?~z#O2vN)d+|4>3d#6++T>VFCje2TdQ+Usp-B#d_xFy
zo&7*_*$=pruM}xxmaTF`8|7y#|MbZt>Z*EviF(fXl==9FUG|^y(enL`5X0i=I>Mt`
zZfT>U5i-a)n|1Z1-5)d6mGhM%Z3F^+IKrcxTl1?uN{2?!2Yx}-k@I7Ic;KWD=PO0p
z=+FrIaD+!Wx2~{!`wkr%K_B>D>ilOr95U(L`AU&CIy8bl9N|&UtuIuZu}g<W&<Fla
z)zOb5-m_^MubfZMo(q9SIJ%DTDCgESvg7vb&<Ofq9+Knk0k7FSjl0gLcP)iLBk02s
z9_8FIQg>(sefaGDaaU}f{KNT5kv2Lsf<7GKQO>QWMT9^j=z~4!eIU}8oNsIu8lhE=
z@F?dlMcU}l2>PIh>EZO8YDDMGSBkXJp%IR*BRtBv!wot#LI%lH&PB)|)d+|4jR-V?
zJ{;ju&K>^Ip^?-_&Z($G{@`%F@Q)6Sq&^biQO+Iy(V-FaK_A#l?a~KOa5!K1M~6nx
zha)@+-sBGdsAvRzkPCVV2f3j4a5!K1M@1tXT}OD7bBBLaG(rZcDY69zHANQTaK7-5
zibl|fBRtBv!#^0QyY!E>bAWwTs75%PFZ^SMM$$hL;Ze>V{xL@*=mU*sKX9Dv2hNau
z;UBnKw8{}}l%KIk8$B99AB<I*51cde(eg<fJsRQYI>Mt`ZfT>Je?$fubF;3P0c2e{
zU-(BL(1#;D%DKZo2Kh(m!)N!*WKti_7ydEGKSCdl@F?dF{}|*Sp$|MTb<W%>b?$uO
zAA|fO^x+7Pa_;btLH<$E2VS1z6?3Z`ubeOZV~~Hu=;H{Fa_&-IA@|?wQ%BU_D1*#N
za@_TqzVn5D4DyePJ{;ju&aLOj{^!V{eC3?k+yC%Kz0dA%{;FNWKmO=p4_mq98@sN2
z{jmL3*e+#zA)Y71MN0^e@)oMvh4|%{&zWB+gwO6l?0n8%@w?7f%3Vf|`8i7nkMh>)
z$+PVa8{T-?@%_WEy}GBjR3O^^`zs!_^7;qu5gg7}%F~57?(*aNpIp--Jjz?Ee7s=a
z;m<#FX@Bo4ukOhj6$r;sr#*P(K3C0x!}&(U&O+?CgzzYDt=@^d_x*<}x4+G>$7NUd
zKD&R;|2$-cKXTxFrMy;%vxWH362harwOY+e(a#WSs}S^oz4U`VfP%yMw6c*9?-k;q
zO9+qh)|PTR_4-vhUUyK>$pyWHgIv&iIGj&U$_nvPA?~$=@F;Jso*YoVk5CqmQqFyL
zPfd|UIGj&Uk_+)vA@*8Ac$BwRD>{7ToH-gnAIr0Qsu2$7ljB@*&iv=s9XZTy-6A~7
zTU*L^zkJT@H$pr}2wW!5InZ481McM0a~W5hGecExQmY)%M)?_wv{BIr`e1y@eEj_<
z_n-37@;zFJql9?E62hZeZfT>U5%j@0n|1ZYZ4a30%K0?<ED-3!5gz5O)l(MR-=;$&
z=mWo?>Zq%i-tyo{AI_)i3__q0^x+7Pa&C>(_r8BOSH(Y32z)Pft~DoD{&Jyn=hHi*
zLZA`!;Rug%Ze5{&!M@!_Rs0U>IsTpF)$8wi=rmqApYHq!@s~ne`(G`>qnul_rR@&u
z&<Ofq9zxY|tRA@2Vbi$le5FVm9U4I&j_@dNZTLrrM$m`P?tk<2!zTZ5K0UvC-H{y{
zK_8CrDCaIk+UU>-`d|-Lm;Fdza=vh!4vnA>M|hNThktZv1bxuM^l*AkHKKD`n0%#3
z8yy-!ACB-S=MFdM&<Oe<Q#ltQgH&UTSH?FY&<Of)gtyAM!#_GSlKRLw6_v;z9L^X1
z(V>ylM<P7RxpiKvXf%R8=mUGn#~0KGP;fY3_(z9E(1#;D3f|-{McSxn1bvVTdI<-)
zp!aY%pWgWq0*#;#M|hNThksNwf<CAzvIPe@CW~-5U-(BwBk02s9_8HOA2az!m;SMK
z4zSM()d+|4g@4S@Ncu-2Jj%JlKjvrzefaDi$H{)+49OS%F_(XI=))0hl%KIMw~~K!
z=!3B;^MP|_K3YC$qemm?!x0|Uax=G*e{|@BF*oaq89>&R^M!u|0)05bqntbZV~~G@
zK74l1OeXc=eBmF1{3G<?2#<2^@Q(qFpbtDSb<W(1Y9!~*7ydEGKSCdl@F?dF{}|8+
z`oPO`ykc&Z<CXJ;e+*~@eK^9SoV%3g$^Ad95$<Q2A^VJ;IZ2L#&ZpI)g*f!`<Lkc4
zxg$KvxiuGk=CQ-kKR<VN?4D=$KD$5WrL1+e&&r8M?z!^%ukW(*oR{tw+j_c9h`+ez
z+}VGvv<Q##7D^vNJm=@<&c7uDd#FHi0OF<3-Yb6B`AYes5ZhjJ?)*7wl_NaLTU*Ni
z{mqkx7r){~{a@Yx?4I6Ifp9$F`@60@>hgz8{YXBYs|oQ<A-;6K7U8XOZrvGuz~*7&
z&93YpvD4W-S)&5sIQApEt=#gkS#UU?-d7XijzT<V3E@%R+ETvy>kY%zFWhc;%ywt@
zKD$48_K+3+$cOXkSP4;HwB4}t62harwR$c?zyD34?kWU*U@!fk51`<1KHXOn;u0Ya
zSVDM|w^r{#s@Gps?~hZ@$pyWHgIv&iIGj&!NDJ{DAs)Vj@F;Jso?=z@KCdi(RT-qF
z$QB&rm@LBKd^);9d_;&{RV9w_C~vKvfm9X#K-KsI)$a1_{%O((9L`sYv@u5`=)(~n
z<=j%O)bt}Et`q`YXFt$f_JdzazEbWb1g*Mi3DHLR8H==0(MagS&N<G#*Q2I<w0zP=
zMI)gPLwJ<8R&$1%U0Klx`taHPFU~n&sw?Nyw`GMuBk02s9>r*y)uYkph1+##1byHa
zR2@10$JK);eK?<H0596ELnG+J5gz5-rCj#w4c+%u@h1v_@8vniEB4(y>D>ABy*nYk
zCd8wb5FX{+QtbmacW5N^Vdos5yWL^ac;$S03qpv;2=QX{VF-`%*6ObQZ=Tel5%j@4
zB*)#qzsup%xa)knvMmG}K_8CrDCgD_DbGB%3yr7}eRluqXCFTKhx3JhbZ7*9IKrcx
zyA)}oLnG*eJ(KsTAL&cZH?|6mpbtlQlyirFbZ7*9(8Kg_dQLT>b6S{u;U66uK_8Cr
zDCZ9U=+FrIAX9nHK?bQtIGiv1qeCO;!x0|k+~FS`8cBWRoQg{14-V%G|LD+2>LU>z
z<=m=v{XQB&A3nRM51`<1KCRXx1R6mfj_@dWlRNyQq7n4rvwLzu@8NL1QlyQFM$m^N
zJj%JlKPnnQAJi1tVk{xYWDySM3;(ES1bsNdqntbZgOR#R|5!T**k^@mgv0s5KW6fe
zF8w1B9_8HOA9FN<KG1mf1INjJ;0(zZ{(-A?=))0hl%KIMw~~K!p$|Lfz&SG?Egy3$
z`9~M}FoZ|7+{~@yA07Je**!CWtSjdW|LEl(9r|#DNBOx(_{Sjs2z~hMo|#PQ!}-EL
z2Kh(m!x0|k+~FSs8bKdEyJv2dI(NSCj{%LK4@Y>EbBBKn@{bz&uyYRPRykfdU--v>
zMnWHk@F?dF{}|*S6@4%#$#IuCNsfcg7ydEGKPviggh%1!IlAgh`05jeO;?>aySU%6
z_u2iKe>|RloOI&eD_=fpmzCeVdA}95HSP-WX(4`ZXb~ReE!2D4LL7Y6d2<l#p#sqz
zh?_rtulQZ(E9DbHTq4A(TIC3j^41D*m*a-pUwUf)#dq!4(_1PKj%}X3>&gZ1cvx^a
zpH>kRVrL;PSVDM|bL-CNTaFkm8@}H^_lzBTvPK2Maq0fMt?aYyEI6D`V~G%N6XKdB
zghzR6wGPW;_ZpsX%AJObPusEg+5I)&-F<~Wa^QS=QdWq^3-Qe*ghzR6b>&XcuM!Fb
zePA#Bpbwzna6WxoR*2UML8~0$QQq28{z1LoQ@#Jn8+Po;1-*oWT>PnW0*CWyojW1!
zC&Y7ftQ_G{&aGpm?4749o~H~_Q)CMcYKknv;d~l*h4_LHZ+&Nr@F?fjH$YW|7pNM4
ztJ+<j-BXQlIG^4H_}+Q*&j@kzC4@(LYxM+#)PzRR2fEIFpt<Y^+{ve&fA756st~_(
zVT)*^{EX$FFFm!cs^@Q4&l#UGAA4<j)Rd2wPui$x1bsNdqrA2HzQOSQT9v_PDT9o&
zSy$&>b-+|t&Znm)h1gYyFR4l#;Ze@5=SxnxQ-?;VMEruPBj<m-{UMV+oG)s(LnBnX
zBRtBvH7|MWUL6_%f$ycx?|<&*N$1X|Ysf;N5nAO4k8*C!uikP*hepr`{+;91k6(J&
zG+sGhDbhxVM$m^NJj%ICkv6)}h%(4LB*)!*o_Y8*?mAy7(ng0y(1#;D%DGFCHaawd
zK74lngrAM)ALzsR<by(>5%l2*k8<wtj}DEX5B8+@fk<C+K0Qq(4m5&39N|&UErfcG
zM$iX6Ob@5$RAcrd`NBUsG=e@H;jMD+QlyOzji3)Qm2(j?NHxOYe5FVm9U4I&j_@ew
z4*%%TNa`c!R8%5=a5!K12esR!J`&+k&K>^Ip%D=DfxT#uK7fM5`9=+*5nAO4kAgS3
z!#^q-K_BFTUcx~x=sg_H7yePv2>Nh@M>%)+M@1v}2Q@{u;2_6j5f0}I|EOpLeK^9S
zoICtuhDOpq*3JR;S)m%?aK7-585&9dNQ6f@clgI#{?VZiG@kvyak3vcL-K`x%+Uz?
za6}vBXDrf2k4Df3V^!t@=gfSxe9W!nA07H|gh#d9(ngO)&<A5~))g~=tSjdW|LD;O
z`f!9tId}NSApeL;^w~W#nbe2#jdLIQM^w8bJj%JlKL+_nAn?G{IdiMjx$}j8Fo%m)
zIl`lyyOcZ2H9xPB;io!QczKRj%q(-fa=!47LH<$Eha)`7xu-dGJzp7QPLks;bCMhf
zoo|{`*UOg>9_8GneCJ0e4PU*(`7@p86`$Sb_Vb^=?{xmb_WFoyl_OBsgqVwB3sF=c
zdvY6u^EHT)Q6M}jw?~c*w|LI*%i~|(-}<_@JjiGF%yp^8Y!%#NKOQN>$wJU`M<kyk
zMvja3KWaGl1wRC_g{)D5?8$8q&R5F2gxL0lKkT=8#p@oFQ6N0Zxl7sl!M%r-C+(zn
z4aydu-RCwO`AcIz3{f&_xmAf6`@!#nf`dM=mp!@7@8&O!9EK<v1;V3pd*t9aalGK*
z_~$uCZu7hOOCyINN=7ZWDiI?G`6lmhkXJItp4{el^Or^rLzIjH;ZeCga!4C<Gy=!+
z?4Anbck`D<4nqjEcGlogxjk{r790g#XFqa#{N3@FCXO3KB%dQj4)we)j#XiN%6#ND
zh?Z|Xgh%D}$e|3@Ma~Q3Y}Qq7!{K~-a(NwuN9Fd&p-Su)wOjBDs*XLm4TtlUvL3>t
za(m>^?{`$Z(pT_3s*XLm4TtlEYHtwXQMo;G=s0y$JRJBpRmYy(hQs+vSr6e+xjk|y
z-yIbX2lEiBjy<^zhx3)P9>Sw?d*o2vc2qnZKD*CtIGiuW_8UZaRBksM;Zf|N>arhf
z!{K}r2Z$|>@TlA#ITk&y^n)H!b?nLQ>__rV#|p%jj9PA$b59uz4l+g6fynL52i$Bo
zh%Fhl+$!guDzT#<R2+XGa+}}HUm7_U)ed4yMlH805&Zw8sNkRv>}5}G^Sk*=BgaC6
z!I4qRtxCkmu^1VGgX2Lj*^}G+ZvN89VTdgmwcM&ij2w#*E;z_5*<w#_^Sk*=BZncj
zWYlu25;1Zxw?ZT2V|jK@1@gQ3OCyINwm8D0a(m*yqk;pCXFqa#{N3^VW8%0$MDjUe
z<X{{@BlLo?D)W)sAX>il5FVA=BM0LlqZ}NJxmj1a4TtlEf7~F#qjG!XU|gr-;lMkn
zI`-r?9L^WjeuD^)%I%Q@f1={yzyql|_T)Ak&Nt2au7~ib+#WgbSt|Z6npyhn9$%*d
z;cz}Z*Eh~ZWAt%^N9Fd&p?ou^4g_-&s*XLm4Ttk-jb9<Q2yr{j!ZQklM>+R&r(vc$
z4K-)<{BGV|IOXiWj_a|U`>?HZ99?(ij>#VkEh5j5h)QmLX}n|d8C@g$eIeLG)v*Tz
zXTTtwZ@Mcp*Ik(g;Zfe&xN<ivuH036OVz=_*%oKO*^lIltAoSh>R^NLR(WfuGf1vL
zRx(A^!NFN0XTWecUtHxJ7FRhNghzR6*Y0NMN>Qca_=AJ9U;f~5zUgdx(AB60;Zfe&
z>5P8(mJsxTy>QS6P;fZkc-7#Q<5lMd;Zfe&>B<CGwkx@ymvE2^dJl*5#np4JesJ}?
zL3ot6R!@J9S6H}0U#Tgw1qV4Mi*Puft}TsMakvZ6AUw)jJ6$Q_PJ?vl=K$l~1vs3q
z6lr6QM$m^NJjz?ECo^@WmOCcDBLo`HexSMR2fxJgSMW_&)wnCuAlfKzt+c_FJMR2c
z#;42&*G)1XEuXZ()j{qqH3*OL)=pO-xr0?1bF;3v{*!g(eA87<?vgbKkMh>8-OWHB
zct`4kE2pUs=bNt9hCUqOQQq3=ig4%y4@{kNJuY?be0on<*T)ujEvq9u%DJa2+dcZg
z%X7Tqx?+x3&Np2>@6m@NJj%JJEA-sqt;|Vs+~xXcj)TrO-392;ha)`7xu-h~bM%oj
zdNi1K7pPa~3;*cQ2>Nh@M>+R+#{`X_5B5-XaDYf(a=z)V%nW@v!lRsfy7N=f2R%#=
zr{`2-_9OYGyGs>)IKo@y+|wN_u0TeWP<8Z!3{s77INx-atfCKcn+T6`?x}XUk{bG;
z>foRf`Gdpxrn`Myt!)q<<=m47JN$zdu@?^d016K08}C}85nAO4kD{g7+Ud?HSGHsH
zp_g!w3wjTS^G$bOxq99pJj%JlKe)ro75W%O$rc>sm@LBKeBmG5MdmI*gYYOSCtDl-
z!JTUEG=v-EIR_&()d+|4g@4TCA6@!KB0LHOWNX7exZ~cT4>X>8(0KO4`KG)2+?8n%
zZIqv}ES?<T&QG{=_K>kE^U?Axo=)KIQiJd)Z|(F11b47vmXUSEn45LweA81F+$C!e
z9>vI-t)1pRQHgj*Rvg}u`f$GKX_3%}BRmR+$ktACxX=e4n5xACQ|HbX{=ri((JDuH
zlyirF42wIXF?-GN3NO#`%K5@ShQ;022H{c8y>>SPeK04<ahEwsj)Tq@{xK}>A~y(+
zV&;*p)pKNeihO%LDXTXvs?Y9kzIgxq+Tl_0JUY*Qu|2-c@oc>@qPIC3L_Y6Dl(%qv
zW8^;L6Tg)`>3tx00uY4rjc=tqQE!ast&|4gQQq3|$+NGFPcBz_n;xETS7bkuZ+g0o
zr=uH$x5`_qt7>{u{jPdOQEw<zvX*lZo|@wcKscOFt|r8a5PC|lL3ot6Rz25~>bz~C
zH!UjF$6t(5i}%mra6YXaEyPcSxZ4uKqrA0R#awTTES`w3^ntx_&<9X(@U$U+Y)?;3
z^3;8U@F;I>DR)q>zpmatOFbtS^b(GIyMo_!zEU0`#NzFN2H{cO+ETu)e1BYTD12EN
zq^8Ig9ORfR!r^?S{L}apD{nzG2#@mCYJIrH`yYBzHvD7l9AKZQSvZ_;dK#DaKN^Hb
zd236NHhAJU5a>Glf#$LwI0Mx`=0i{1$NL|8s<=9$jq)=VX@e)1dGn<*K4m_55|e7=
zcUwN5{9L?U(;z&m<(4+;;tih4IGc5qZ&$$Ke7d$I1R6mfj_@dNZM=!pE#5?`_(kf2
zC&#m1ov##WqeCO;!x0|k+~X7RXhcuMS9~vZ&QtKIbLX3$x(|Ig!lRtK6ltSFBj^MF
z&haYWuHaZX-}v?b8hOML!lRs9&n@ZcYczsBn1|%J%k$wm4m#iX76ckWACB-SZ|(S|
z1sXvgKD*DiE9i&wO>cAX{zrrGDCZvE7(pW-*puD|B7MpErngdN=)(~n<*gmxd_g1V
zgC3@b=i3$8kK~))uHotE2H~x8?(mNejo=?-D(52kb_E>H7yi+q5%l2*kMh=re{}MX
z4*#Iy_+w0@68X!zqWakm|LD+2>LU@}D(4RW;EgJth!54W7Y=?C6u;|y;UB!U#Z&hU
z!lS&k;U5)^pbv6EFX3Qhp!aY%U-$=apYirUgYYPhc(ykDqoNV?K~0e@ILI+sgv0s5
zKX{9dw;&pXM|o?*KjQt5^p8B}V5FuR;c&j`Z9v}t2>(cgM|o?*KX~IX5NJI6fyT2R
z&Nsc)$XhAlgNbOP*3OW$42wpE#hWkT&Y6#VyMlhSe2cd;dAp`Tc$BwR&$nxIlYeyR
zgE2QV#h9CQ<$U2Ez5Ii>gc^iLF|uZBr#D$cA9zPrJKmA{aK7*l-j)r0IKrd6wc#Iw
z{3G;%2d2)MTcz%tFZ_eIenTIQ@F?dF{}|*S6@B34IbPxAR3rItzVYp4`3G;GH3*M#
z?(mO6{!!5fbCMi)nUmx==zQZ_*zynFqH7Qy<*gl8n)*krW~G&;s?Y9ob*aC(_<(8s
zUAD*7tbU=@NVNV{gUI*YiSibXD`uUbHBvyZC%rFMmjdB@v8u{&^Kn&`2H{cOTB&w?
z^S<AB?bVguQi0@vx9wR0A^X9K3u8ZcyZ<9vqv4usT7<XCTU!dNGqFmOR%ohZE$1S)
zS*(S{?>e8ZEeXL&P+A45L3ot6c3KHxm|b>trTX}zAFO`B9~{nC3M*N$qLo&&Y7idf
zt({hMU`3!xAJ_{At2lsy!}+FFBUp8)L3ot6c3PQYxTAVbF6bp3<bvMA;e4gsc3k~}
z)sq^8M|o@Y{n~Mb4OVEX)D+o*gB-Jl6dcZ{Ck%ynqVld)JQ{>Yd237glve6ltTYvF
zuyzhGt}X?K^Of@7TIucA<7!q7!lS&k`o4iy>RGIqRnc|!gH=wlA2<Wm4__(L2CJ&D
zDp!MOqx_6Ts~WL#3oCzB#;42&>t<v=TE5o^@pvIPx{mNDZ*8n#Ru?OnRmRz@E7t$W
zx^li)C9N)2Nox=u<=oRsI-w8zBK5(_IjIll(~L~(iY!)ptB&v}Z|$@qQ0N2SOP#YG
z3)M&toNronDD>e7k8<v5Wurdyv3Q%D+GSmt9Iu>jM4*w-hao)5xpn1EtIwbjGRQn6
z$6eO>$#Ky6N|81?G=e@H;Ze-~a&#S6nnEL5NvHbkK3A9GSUF!Q(ng0yv|3Mt@F?dV
zSIj~qS`jE(mEH$}YJ|i2!aq7Rf<E}IM0k|9c3SzXq7R?l({rjZ`;mOQlB#2cMmV~T
z@K!l@_(v!IsQ3q&%DD&`q#EIHzVMGu{!#G{M|hNTPb;x?_y-loA9Y0~@&||Wg@1Ht
z1pjb^M>+Sjq8uv%g=*Og2Ymnqhx1LV*0JhPgYYQl4*y^+y$*da0?<o1$OXNJ!}-EL
z0)ajp;ZYp%Y;E`lEBvuSQ;ecy3l4Hj7U6Kd@Q;c{&<8b@2#<2^@DElhWTmNagFNS8
zq^26-aK7-585%(!em>w)&K>^2ii@n66|R=8LgU#F=L`Q}RYz9kY7lMI+8NTcawIE%
zg*#_H7^^ZLEuXZ}qmj_Z;{9`uu17H<WoxGuFgx_Yn46hm%+0!TzG;=s4t+SnqZq5Q
zwaSNPO2=z8t20%Jct`33?|AA$AI=y4!D^wQ4@Y>Eb8Fow&Bz8R>RfT)fvIymFm>*H
zx{o8oWkQ^~gzzZm4mTKB887r<=Nx!>j#th%BIF;ez85oO>XdwVl(%+VA@*5X&FV**
zAu}h*QPgMp&KLeM$Uj&`ut9hfGmqpRzyERZnWxWw^pYJqpWWxzHU8q1gI6{#*8XL?
zl-=~Df_LbPAy+IRJjz?BZx#!2$(g6mzaoUs?m@5uG{5V7W2@ezFNWx=Bn`r&ytVrJ
zkk%?bTHpV;^m#jU^fo=5)rMICI{T4)y2C33Ukv&D62e>Mt<~LEecRw)#})WHvX*lZ
zR!wFFXgHj&lot!}5h4El*c}>#M|o>Y`Fnk-;G#XsP><Q6qx$$`jQaZG`-5;e-?-ZJ
zhlS8;><z-BytTTjrs$`O=N>}P2lm22A3(w3e0q;Xh-VA2=MuuBytP_qdtA$Vpy%X*
zUcx~x=sg_Hr&S#n2(91VAUw)jD?caS%Hms<K}NNltCM502#53O`qd|o=r0oDvr7n%
z^497JY}GB_$Uz^B&THoY<JUFdaK2JLt8W_|AjJEY5FX{N)#?}0-(!V%mk_v2_Jj4d
zvmbC%{qV_g*1rGIAlfKzt+Y|m2>M`r%6zbHc;=(!d%M>0{)P~5UqX14x3&~%qoNTq
z$T*vI#ro%2SI(!a&O)FOs>Bf<<*gmRByfLK;?b%^W-qA^z8t^`=^P*D)0YJHD15CT
z^x+7Pa_;eq2-m6N^+kjZ-%Fjd9)0TE`Sg5=)`ve+2z_;-L3osNk6(5`BlH~q&hd(M
z^{Gbs;d~<kjc{}w;ZfdNJ$tOx_t6OYU>-u%p}~(Vz6An@^OYiPbZ7*9IKrcxJ4Wgb
zji3*o-RIXe;BdZDq>T=ZpbtlQlyi??4B1%~|1}}lL)Fob{JI7l&NqJl1C8jb3k|}f
zoLkQbtLKZA>*MIAhtqSa5f0~@zOKUeKN^Hbd27QzIy6EC$yClo$RO32`AEL-j}DDc
zC64e`k&kSxj+H8rFA1bRs5*`nmB=3)&NqJT295mxg`E$~on?8(9~ZWT8RuBEt78d_
z$TG7@rZT#&-}}6Bq%AE2O6sf%GEj*bMTFVKGHfyYv1HBI;#`2jU~`)xGx-GBCexnx
z8!bhRhCe#kW=llGIfr(0w7Se}-SxZAcc1TiUwTTDFFnb*uX8>3d*9b{-S2nK`5yd(
z=tOvwbC)7*aH0+;B7|yr77lU%1&8x#ZWm(l)P)|xqntbZgOh<;GABE-goA5C_HZ~~
z_(vePx{mNDu6V9C{G*}~^x@MzI~?5zhx3JhR5XG<9N|&U9sbeCKic$<yysx2rW@gK
zzVMF*jii4h!lRs9_i9oTCx)O8pYCxSx)BcN3;)RTKYECM)Y2W&^yDl~J_&cudazez
zJ^Fl$r*m=sM-SmqzS{7QPX599A2F9?rPy<`ubeOZqmzF`m#|AG!lT$(bG7nX?b55&
zyKhwQ;vIR<k*8~Lt(;HqO$c$h5KK#s@F?dF|5!jH=mQT-o#TP2bLT7N9rCA#g;+g6
zc$9O88!X5_VrRh1a|gi7bN_I@2}1r6yN@G0%2zv?kl&@bftkT4dS(*3k!$6A6U1gA
zbP8V&;Ze@5ldW}%^5@^Ttyz81T<g>Qk;hj>&;HC)`g`7s=dJTt6S5FDU9qjXc!2OI
zUtuYi3USx_wv8_lf@kPJ)((ObFnL|)(-+}{xZ{d#<8=dsNBL?eCodk>8J&-9m}|+F
z4upf#7C8Yk^GH6u8?Lh&H_f-*(*uOB%2%uNQgu@E4xQ2Y+}gR8s^u)gsgYmOd7W@L
z-{h3aj|%ah1B6HUYV}!XorHPSgJ11*(s=9BJ*Qvt1BdhJj$Vl05#q-O2#@mB>b_Rd
zUldAbfVSkovs^2300oEhY2OtB9-YtGLwJ<0R%bgY>nnA=$edb`B^=a(?BQ^}Qf^;F
zaCIHwQNCJzazOR|cU6%J(o<9m4r)vl;c!0nuIl}~5MLf3Jj%KC{V4S<8bKe!>7H(c
z!};Xr>RU8|J{;juzFPf0QTjt8=mTA69(lSr9K2HU$(^M?G=e@H(U0;w7HOk$Iw|^K
zf697rZfe$}&nInEG{V(&gh%zcrHzV4s37}n_SNZ&XQ7ja^XV_}3!zhvqe~p&QO+%d
zPKkX+o%o&6T+8&5_Z)k5dUEQ+`AYeFA?_98IQ6b0Jj%ICkv3X1f<ExQ)H&zzrp}#D
zQ-%;<72=>)<p_^*?oy<U7LA|}{F|<$cR5!$_bca{AkYX`*AX7&-1-!f&UQv4RM4k;
zPUX&h(D_P{Hd-`-J{;ju&K>^Iq7n4r)BQ^pPZy`=&KLgCq7n4r2#<1Zy}P87n!h5A
zR3Uu2&(p=>aK2Ku32~cFOlMV&@F?dF|7g(&`XIyfa5ASGGY`6-$CI<-@0gqp+(Y=P
zoICuZMI%&@O64p<1?fg!*ZIOfS~NnJIKrcRwWUZKEgDIE&~@|`oyZRy&Zpnp34umZ
zABpfN=MMj<Xas$b1JA-i4xr$0KFx4?eKdkT9N|&$CU^KpMI&TREyxlMYC-mJIA8cj
zAh^1Y@F?dF|EOpLeXx_|JqI<Wif}kz_(w$}=)(~n<=o*Pjr^lc|Hyj|c51p24(AL1
zXyhMl`bQ!>%DMIXMClKWpbww!(Rk+JeEP=7RoljBgjG4BAGLIcq<c>7Bl3?HeXv($
zJ=m+V9(_J(qvNc2ou1u8cvPR8*-HM=QbG3I>?`)%>?`LB|LEi&(It-XC|_;(M~6nx
zhfnu-N9x1*!aq7Rf<7GKQO+Iyu^|7b=mQT-oikgd&Yds(V*!nz4@Y>EbB7x&$UkCd
zz{_*L!prGK`pWqx2>D0sK92Ax=MMi^kbhM4!AwHep+TSYoiF@jLH<$Eha)`7x%CD8
zdyilE=!xr_Q{KO+_33`q12txzD?k0V(X(HxqjT?FH{x+A%Z2#ZN$Z=J4-g*ZE7bev
zpFDnH*U9V0rwZZI{WtZx_n*HeUf20bIU>X#oV0#?@&MsczFM96a`3o?mwkPC_le6l
zwPc$fe*AB>qwRlkR&Y3<?ox#~<(cK(!vlmz`D#nK=2M5ekNxrXZo|7bwNx!<k?pJA
zI6C~?C^($2lpnqSaQF3(Z12t(AUw)fTgpE?`ee8My7}&*OE$H1A3yZcE3P^-eyD--
zsoRBEEyUgd!lQh(`u3Tk_X!1p9C#KEasUN~^XYSzLVQk$eFKC?`D*p9rLz8luGczc
zPA$k14r=iestFv<SIXUsh;s)BkMh-)@;|EgTUEu+se<fkR0|GziYmsnGCutkAl3Wr
zLhKzNd{w?$z4N2K{Xcc%N$SMmbWb<J;e4fBufF|(5U(B}Jjz#FinP(75%ht}WFBZP
z^ME_~N|81iG=e@H(U0;wmh#N<`VM8jN13xfWj(IH{G4e$`h3zxMI-3L5gyg&-t&>|
zwNuThAp30g)k|M-?zFF*uatX)cx<st9N|&Utuscho3CgDec%_VkMI29c~gBjUn$Z?
zMI&^(BRtBv^}Fy#pR8yE1iqI#zy6`~r#g2&eN#^eG=e@H;Ze@5zp(MC!xfF75Bxj#
ztMC8WhUtFgd=ms3;p#fVqnumQ(!t}}_8ZH)c~y`(gswv$KX>v4(|y<ZN|82NG=e@H
z;Ze>V{?W=m#`uR%_xsjfF!c}T3;$^4A7lK(5gz5-;U6s;K_5Ir*OAAG$|HTr`NBV1
zG=e@H;Ze@5-#{vJG=e_JFg={i=|*(!eBmE08bKe9@F?dlMcSxngbGrroJFW0-3W*C
z=^ZQ~&<I`P2#<2^@Q;c{QXe@}(TV)P;e6pA6^*1m65&zK9sW_#2ncfES?WO!px|)6
z@Q;c{(1#;D3f|-n|EOq$%&7%g!a*&_9uDUV{|E$E*AX7&+@(kx6^)<|dWvenK~K?*
zajlFm{G*}~^x+6!m2+#})oi61t3@Bf>7H(c!})ZMyZTl$R*OCy;Ze>V{xOz+wCKa9
zdmM*ugo9T~zVMH+{G&x5j_5}%-62Wm+7-|U`e3ihdf=Q{k3Ju>mHeYcACB;-KDV^d
zp%E&`o|}Ee1dx5@eBmD*8lg)Z;Ze>V{?W-lLLYcXb{vyT>cjcMKRWqGbh{%w%DKZo
zI{8N+@W9kLvsLQc`NBUs`A6u(5gz5-;Rc=jBQnR!bH8G?%KgguCJ6aQTwO<alyirF
zET9qe!Az3-u21^T7yhvz|ETE05gz5-rM&aDhr9pScS7^3&u?yhx_`@~XGYK7e%ytl
zFMj*X(c4d6H{!8AO(n$dK6pa&>H)%|e1-aMh7kX@?}YJ*LiltK;^SXi6R+!hddEbF
zvxQ()j_@d7t#+S<h3@x1d|-U>-p#E~_qSKwt3R+hIGj)KGzc+&(}D3@2MCXHZtYhq
z?(2T}t`*%q3!7W2MhBwpzj}N$+JEO+!Qp)J!87jb{z!=V0m7qvwc35Yx2wDIGuL*z
z#+zH8?ytG~O(TA&f%EC_ehBgBLfkz-c$BZUl=BsR3n(GTfoI8s96-V0e0rZ#i0g#7
zbAa$DUu`L;E9<qoUVlmEb`KD>P)*=)zEX}CV!aR_A0Rx+SF3lCRqwH?_?xPrPxtf`
zRfNO&O1W&uf$>8^Y#bmw%2!*8v@u2_bRxU+(mlZ0Bf1d|=PTvcgt$<MpB*4P%2%uR
zgrz?;f<ADWyyrl3nFrj-r^iB|5mx1hL6og_&xa2*k1O*FlsWrT-UFVucKwKMB#%Dd
zh8+i*PYdyu0m7sD+|ouxBj|&DHv8(-cbzlsE9Wca<w9&&thpmR%2%u3;(X@Xibl`}
zenHpKyE`uVjj2AIuM}ycq7n4r2#<1ZeNy&&yDAz%ANXGC{KC`Do$B2AN|81y8bKe9
z@F?fjw`Et{SJ4Riz`yA_@_5gK`d#nhe&u{}P9e|;`f!9tIk!GvvanFEmPRgB1%0~z
z)n(^R_g&}HJ8?oBy!pWR%>#r-Id}L+MI&^gPxpWQE9Xu9!}-ELs{CV&e>lRUe6`^p
z6^)<|pYB1VFF9ZMM@1v7$`Kyr-1;`TGDjomgACKd$vnDI`f$Edq>YM3(1#;D%DKZo
zDjGo_R4Qi?Do8h?LFdyq6jXCGLd_lFQO+IyQPD{1BWEf)ksmmmFZ`pTk<>>bJj%Jl
zKPnnQALPKZ^e#Dog2Va3KPnnQACB-Sc#}K)qoNV?K`qD<4r)R6a5!K1M@1uXtqkE&
z&aKQ<Z#05Fe7dKns3IKBr@zd6^MMA9pbtlQl&?1Yqmh5q^pB-`fU`$*BOK0G3bU1F
ztcpI^&lBNM&K>?SmVZ?A;nO{ilX>t;$rt`HmVd;m95INRa_cV}X&*r&=)<Rb_NuH$
zpHJEtqY?Ds2#@M>Gh4|&LLcn8*;h;e*;mdN{?VZkYVHV+a_;btPW}=4@adjOCiUTb
z;UAs+BlO`2k8<wtk52v(`oIHI=gd~AbLR{H=;R-v4@Y>EbBBL)@{iDmPxp9v?pMwi
z{?W-l;#xVvqntbZqeCO;!zX%XlH7NlFZ`oJBk02s9_8GnY`*xt-CvwOKlisk*xLGZ
z|JwgLEBxc|-U~;+e&nptn{GLK#AAJqOo;EUnxDIGfbb|^q29p~;=t+iyYCf(XXrq5
z2V#8n+IU^()18+PhgQw+zJGx5C||8k@Ob@Y-P|KLjGsHawI$p1@SlFo>d}g)*93?2
z>2EFz@rEzmF#gs6;ZeTYQg%MEvAgku`^PsO+S*bzI*>f>x_b5Kv^P&UoUfD{g}CD=
z2#@mBmhv}etnA*sWm&iP`K_%__rLT{`plYssDbn8uaT@+*}d*42#@mBmhyXweoQ=j
z#6b=`OCIC^3J&Mf?|p=LUWkVV2#@mB>iu(Ny+he^t)do-GcL6tdpMj=?+J6Aj)L$g
zU#-5irFwrrRs5zZNKa8M@}Q@vA{@?F$}3gxZ9;r|fbb|^tv-{gzC|PGgWcKg0Y-Eq
z9L}e|@1wp&Bk02s9_6dm>A=z-8bKfEI`cqtnFp_weEN;B^oK^!ha>t?zS?j6%MH!V
zr(M?l*!Q=#u|F-=;|nM0Z+NR7hm=R3PughE2>Nh@NA<a%x_W=}zyG|k`=DyhKAU~@
z<xi}e_LcMLy`pRPHy`-NjoqySghx5IzG<;#Sw$o01HYi_sQDAWsP~2!_2GP_NE;Q6
z3=tmXtJQCv&sbT}2>QVHQs?iyW&Kp=&ZqaCg+L?d!x0|ktJNnOpV(N@2>QUkbH6(H
z<ofA;<$U^Fh7h5V0m7r4Tl#qYWfhH}59W~EcONe2O!r;q)903iKqKhG5gz5M4gaWU
z1bz5)|CN`WGxZPW3;(ES1bsNdqntbZqoNV?!87T7AkvqdFZ`pT5%l2*k8<wtk1GFY
z@DDOf4<~cF5uH0<_(w$}=)(~n<=o*PRsPXXb1Idy2o<Co;c&k2k1GFo0m7r4JN%=f
zk?i)Ispv$0;Bdb1kBUZy2#@mBhJREvf<DNBXVD-zfP%yM!aph+K_8CrC{~)Q4gYA+
z2-k{QkR=?{g6!dNzVMF*jSLYU<=o*P4H`io^c2;CgPx*_a5!K1M}tPtha)`7R~!B@
zC;zDFA9fEgq8s6GzVMGZ`A1FvNQ6f@clgI{`A0<`Xgu@4aWW5_A^E~TcFR91`fx-)
z%I{d1tu$j*^ub=0^}so^9(_J$EBQx7ACB;-J~y+KW~`_ndv5j>6F~Nr^XYTQLTJW%
z0m7r4TRPvmtb6b4R@SH0i9X%q9dsi;=zQTHo%|#8;RuiN)#|q&Z&=ygbIr2)fH?5L
z)H$<N>fHHs?<_}q&-<6vCk6<Qa_;btPW}-)gHQK(dG1%v7yi-7KVE?FC|_;(M<@Si
z&<8U~?z_w+x$in(_(v!IXwZivJjz$A_tos5PxnWTtPTGt<)V>4)_Zr4-#Px_YyWlC
ze>c~)vk<;QBhW*hp##yKzpg)?a(ot!SqSIubC4|^cx;XhSKK;&<!k15?-R!?gmd>f
zs2Ux3Y>wLw?Hup_krT#$DvntQ=bmtU>*t>u@Bh&DmhR&R4t9pg&#)e^8gChwN7syZ
zh{F+nMV(uJiS2@)>|lEShHWi5@N9BSP#5+1^qsgz){OV>-ZH*f979B|%DMa3idvB6
z(Kwd)xK3OzuGK6=pW9jJDXMjB4(d(asrf8KpW8Xebm<;IpKc0`Os>_^V`<~FXoTLK
zh3IoT2fBXoJkVcq91F3uTIXPYq8{wEtZ-@V{`k_1=jWdM?OVrtl*cSYpW8XuXJ_@*
z65mU|`%}%|{cz`ampEo2`rOh+oGOPt@Qb7Aqo0$s(V!9ZF+}7_`rPD!M$iYompXrm
zX%LRd5A~Ne8Z?4F9N|%YZogLeckWlt;g2<69opICUNQ^O=XMU}klc6ub^Wn^-+sld
zP40=a5YFuk`W(95bNK7}<Jj99G=e^6A)LF<!87T7{<{8n$`Kwl3*p>-4l+y+_t*8u
z(`z-yUNQ^e+<guzm9vPyu0NjEykSqAg>dc(2O3Fz&~>>dX2141q>ToRq&^biQO+IO
z&|hXMA3`7Gz_ZCQK^?_E8Z?4FhKO7hyvg0aR@8zlkH)dYCv7xngljbm(dTv+dWvct
zn?u@|LnG9D7NXDX9Avt551>ys#U3)bR!ff;_x2|J<7kLJw{xKJ7tdqSqk`jDh$U|4
zV6S@ddMtWvT&r1#KDTqQ=g#V@rE9e~Gepg2A^O}qPkk!RN_}K;zM|vO2eT9Te~II5
z&#oDt^2sgDjjB0r<A{vvbCZYsqoNNyFm*n08+>D^k7wt%jDO+THO)=paKvOq6GWfe
tuN7XN`;~L}<Je1>MdMn{LiD+vgPA1vU4LDF9D5=&b?k|=5YAo7{{sm9G~fUL
literal 0
HcmV?d00001
diff --git a/CyXTraX/generate_atlas.py b/CyXTraX/generate_atlas.py
new file mode 100644
index 0000000..7506562
--- /dev/null
+++ b/CyXTraX/generate_atlas.py
@@ -0,0 +1,58 @@
+import numpy as np
+import h5py
+from .artist_bridge import CylindricalProjection
+from .mesh_object import MeshObject
+from pathlib import Path
+import json
+
+
+
+
+def to_dict(mesh_object = MeshObject) -> dict:
+ return mesh_object.as_dict()
+
+def record_atlas(cylindrical_projection: CylindricalProjection, mesh_list: list[MeshObject],
+ atlas_name: str, save_folder: Path,
+ map_bounding_box: tuple[float]=(-50., 50.), number_of_maps: int = 4) -> Path:
+
+ grid = np.linspace(map_bounding_box[0], map_bounding_box[1], number_of_maps)
+ file_index = len(list(save_folder.glob('*.h5')))
+ save_path = save_folder / f'{file_index:05}_{atlas_name}.h5'
+ file = h5py.File(save_path, 'a')
+ projection = file.require_dataset(
+ 'maps',
+ shape=(cylindrical_projection.x_px, cylindrical_projection.y_px, 0), # Initial shape with third dimension as 0
+ maxshape=(cylindrical_projection.x_px, cylindrical_projection.y_px, None),
+ dtype=np.float32)
+
+ projection_points = file.require_dataset(
+ 'positions',
+ shape=(3, 0), # Initial shape with third dimension as 0
+ maxshape=(3, None),
+ dtype=np.float32)
+
+
+ mesh_list_dict = list(map(to_dict, mesh_list))
+ print(mesh_list_dict)
+
+ file.attrs['mesh_list'] = json.dumps(mesh_list_dict)
+
+
+
+ for x in grid:
+ for y in grid:
+ for z in grid:
+ # x, y, z = 0, y, 0
+ print(x, y, z)
+ position = np.array([x, y, z])
+
+ image = cylindrical_projection.compute_projection(position, output_full_ray_projection=True)
+ current_size = projection.shape[2]
+ new_size = current_size + 1
+ projection.resize((cylindrical_projection.x_px, cylindrical_projection.y_px, new_size))
+ projection[:, :, current_size:new_size] = image[:, :, np.newaxis].astype(np.float32)
+
+ projection_points.resize((3, new_size))
+ projection_points[:, current_size:new_size] = position.reshape((3, 1))
+
+ return save_path
diff --git a/CyXTraX/load_maps.py b/CyXTraX/load_maps.py
new file mode 100644
index 0000000..8d5105c
--- /dev/null
+++ b/CyXTraX/load_maps.py
@@ -0,0 +1,23 @@
+import h5py
+from jax import numpy as jnp
+from .mesh_object import MeshObject
+import json
+from pathlib import Path
+
+
+def from_dict(mesh_dict) -> MeshObject:
+ return MeshObject.from_dict(mesh_dict)
+
+
+def load_atlas(load_path: Path
+ ) -> tuple[jnp.ndarray, jnp.ndarray, list[MeshObject]]:
+
+ with h5py.File(load_path, 'r') as f:
+ maps = f['/maps'][:]
+ points = f['/positions'][:]
+
+ mesh_object_str = f.attrs['mesh_list']
+ mesh_object_dict = json.loads(mesh_object_str)
+ mesh_object_list = list(map(from_dict, mesh_object_dict))
+ return maps, points, mesh_object_list
+
diff --git a/CyXTraX/mapping.py b/CyXTraX/mapping.py
new file mode 100644
index 0000000..cacc2f1
--- /dev/null
+++ b/CyXTraX/mapping.py
@@ -0,0 +1,115 @@
+from jax import numpy as jnp
+import dm_pix as pix
+from jax import jit, vmap
+from jax.scipy.spatial.transform import Rotation
+
+
+@jit
+def map_source_2_cylinder(source: jnp.ndarray, maps: jnp.ndarray, map_positions: jnp.ndarray, radius: float = 1000.,
+ angle_discretisation: int = 2000, pitch: float = jnp.pi) -> tuple[jnp.ndarray, jnp.ndarray]:
+ """
+ source: 1 x 3
+ maps: u x v x m
+ map_positions: 3 x m
+
+ ---
+
+ cylinder_value: m x 1
+ """
+
+ #map_positions = map_positions.T # m x 3
+
+ angles = cylindrical_angles(source, map_positions, radius, angle_discretisation, pitch)
+
+ u = jnp.expand_dims(angles[0], 0) # 1 x m
+ v = jnp.expand_dims(angles[1], 0) # 1 x m
+ w = jnp.arange(0, maps.shape[2])[jnp.newaxis, ...] #+0.0001 # 1 x m
+
+ uvw = jnp.concatenate([v, u, w], 0) # 3 x m
+ values = pix.flat_nd_linear_interpolate(maps, uvw)
+
+
+ return values, uvw # m,
+
+@jit
+def cylindrical_angles(source: jnp.ndarray, map_positions: jnp.ndarray, radius: float = 1000.,
+ angle_discretisation: int = 2000, pitch: float = jnp.pi) -> tuple[jnp.ndarray, jnp.ndarray]:
+ source = source.reshape((-1, 1)) # 1 x 3
+ directions = map_positions - source # m x 3
+ direction_norm_factor = jnp.linalg.norm(directions[:2], axis=0, keepdims=True)
+ directions = directions / direction_norm_factor # m x 3
+ intersection = directions * radius
+ sign = jnp.where(intersection[2]>=0, 1, -1)
+ intersection = intersection * sign
+
+
+ z_value = intersection[2] / pitch
+ z = z_value + (angle_discretisation) / 2.
+ beta_value = (jnp.arctan2(intersection[0], intersection[1]) / jnp.pi) * angle_discretisation / 2
+ beta = beta_value - ((angle_discretisation - 1.) / 4.)
+ angles = jnp.concatenate(
+ (jnp.expand_dims(z, 0) % 2000,
+ jnp.expand_dims(beta, 0) % 2000),
+ axis=0)
+ return angles
+
+
+@jit
+def map_geometry_2_projection(source: jnp.ndarray, detector: jnp.ndarray, detector_orientation: jnp.ndarray,
+ map_positions: jnp.ndarray,
+ projections: jnp.ndarray, pixel_pitch_mm: float = 0.139):
+ """
+ source: 1 x 3
+ detector: 1 x 3
+ detector_orientation: 1 x 4
+ map_positions: 3 x m
+ projections: u x v
+ pixel_pitch_mm: 1,
+
+ ---
+
+ projection_value: m
+ """
+
+ projection_geometry = projection_matrix(source, detector, detector_orientation) # 3 x 4
+ map_positions_homogen = jnp.ones((1, map_positions.shape[1])) # 1 x m
+ map_positions_homogen = jnp.concatenate((map_positions, map_positions_homogen), 0) # 4 x m
+ uv_px = projection_geometry @ map_positions_homogen # 3 x m
+ uv_px = uv_px / uv_px[2]
+ uv_px = uv_px / pixel_pitch_mm
+ uv_px = uv_px.at[0,:].multiply(-1)
+ uv_px = uv_px.at[0,:].add((projections.shape[-2])/ 2.)
+ # uv_px = uv_px.at[0,:].multiply(-1)
+ uv_px = uv_px.at[1,:].add((projections.shape[-2])/ 2.)
+ uv_px = uv_px[:2]
+
+ uv_px_flip = jnp.flip(uv_px, 0)
+ # uv_px_flip = uv_px
+
+ return pix.flat_nd_linear_interpolate(projections, uv_px_flip), uv_px
+
+
+@jit
+def projection_matrix(source: jnp.ndarray, detector: jnp.ndarray, detector_orientation: jnp.ndarray
+ ) -> jnp.ndarray:
+ """
+ source: 1 x 3
+ detector: 1 x 3
+ detector_orientation: 1 x 4
+
+ ---
+
+ projection_matrix = 3 x 4
+ """
+ detector_orientation = detector_orientation / detector_orientation[1, 3]
+ rotation_matrix = Rotation.from_quat(detector_orientation).as_matrix()
+ detector_horizontal_vector = -rotation_matrix[0, :, 0]
+ detector_vertical_vector = -rotation_matrix[0, :, 1]
+ print(detector_horizontal_vector)
+ p3x3 = jnp.vstack([detector_horizontal_vector,
+ detector_vertical_vector,
+ (detector - source).reshape((-1, 3))]).T
+ p3x3_inv = jnp.linalg.inv(p3x3)
+ p4 = p3x3_inv @ (-source).reshape((3, 1))
+ matrix = jnp.concatenate([p3x3_inv, p4], 1) # 3 x 4
+ return matrix
\ No newline at end of file
diff --git a/CyXTraX/mesh_object.py b/CyXTraX/mesh_object.py
new file mode 100644
index 0000000..df02eeb
--- /dev/null
+++ b/CyXTraX/mesh_object.py
@@ -0,0 +1,23 @@
+from dataclasses import dataclass
+from pathlib import Path
+import numpy as np
+
+
+@dataclass
+class MeshObject:
+ object_path: Path
+ position_mm: np.ndarray
+ orientation_quat: np.ndarray
+
+ def as_dict(self):
+ return dict(object_path=str(self.object_path),
+ object_position=self.position_mm.astype(np.float64).tolist(),
+ object_orientation=self.orientation_quat.astype(np.float64).tolist())
+
+ @classmethod
+ def from_dict(cls: 'MeshObject', data_dict):
+ item = cls(
+ Path(data_dict['object_path']),
+ np.array(data_dict['object_position']),
+ np.array(data_dict['object_orientation']))
+ return item
\ No newline at end of file
diff --git a/README.md b/README.md
index 89d45e0..f8710b3 100644
--- a/README.md
+++ b/README.md
@@ -1,93 +1,11 @@
-# CylindricalXrayTransform
+# Cy(lindrical)X(ray)Tra(nsform with Ja)X
+This repo uses aRTist to generate local cylindrical xray transform maps / atlases of meshes.
+It also defines the used `.h5` dataformata for the other projects, and I/O util function for the use with the autograd [JAX](https://github.com/jax-ml/jax).
+As simulation backend [aRTist](https://artist.bam.de/) is used with the [THD API](https://github.com/wittlsn/aRTist-PythonLib).
-## Getting started
+## Dataformat
-To make it easy for you to get started with GitLab, here's a list of recommended next steps.
+The local cylindrical xray transform maps are stored as `.h5` file.
-Already a pro? Just edit this README.md and make it your own. Want to make it easy? [Use the template at the bottom](#editing-this-readme)!
-
-## Add your files
-
-- [ ] [Create](https://docs.gitlab.com/ee/user/project/repository/web_editor.html#create-a-file) or [upload](https://docs.gitlab.com/ee/user/project/repository/web_editor.html#upload-a-file) files
-- [ ] [Add files using the command line](https://docs.gitlab.com/ee/gitlab-basics/add-file.html#add-a-file-using-the-command-line) or push an existing Git repository with the following command:
-
-```
-cd existing_repo
-git remote add origin https://mygit.th-deg.de/roboct/xraytrafo/cylindricalxraytransform.git
-git branch -M main
-git push -uf origin main
-```
-
-## Integrate with your tools
-
-- [ ] [Set up project integrations](https://mygit.th-deg.de/roboct/xraytrafo/cylindricalxraytransform/-/settings/integrations)
-
-## Collaborate with your team
-
-- [ ] [Invite team members and collaborators](https://docs.gitlab.com/ee/user/project/members/)
-- [ ] [Create a new merge request](https://docs.gitlab.com/ee/user/project/merge_requests/creating_merge_requests.html)
-- [ ] [Automatically close issues from merge requests](https://docs.gitlab.com/ee/user/project/issues/managing_issues.html#closing-issues-automatically)
-- [ ] [Enable merge request approvals](https://docs.gitlab.com/ee/user/project/merge_requests/approvals/)
-- [ ] [Set auto-merge](https://docs.gitlab.com/ee/user/project/merge_requests/merge_when_pipeline_succeeds.html)
-
-## Test and Deploy
-
-Use the built-in continuous integration in GitLab.
-
-- [ ] [Get started with GitLab CI/CD](https://docs.gitlab.com/ee/ci/quick_start/index.html)
-- [ ] [Analyze your code for known vulnerabilities with Static Application Security Testing (SAST)](https://docs.gitlab.com/ee/user/application_security/sast/)
-- [ ] [Deploy to Kubernetes, Amazon EC2, or Amazon ECS using Auto Deploy](https://docs.gitlab.com/ee/topics/autodevops/requirements.html)
-- [ ] [Use pull-based deployments for improved Kubernetes management](https://docs.gitlab.com/ee/user/clusters/agent/)
-- [ ] [Set up protected environments](https://docs.gitlab.com/ee/ci/environments/protected_environments.html)
-
-***
-
-# Editing this README
-
-When you're ready to make this README your own, just edit this file and use the handy template below (or feel free to structure it however you want - this is just a starting point!). Thanks to [makeareadme.com](https://www.makeareadme.com/) for this template.
-
-## Suggestions for a good README
-
-Every project is different, so consider which of these sections apply to yours. The sections used in the template are suggestions for most open source projects. Also keep in mind that while a README can be too long and detailed, too long is better than too short. If you think your README is too long, consider utilizing another form of documentation rather than cutting out information.
-
-## Name
-Choose a self-explaining name for your project.
-
-## Description
-Let people know what your project can do specifically. Provide context and add a link to any reference visitors might be unfamiliar with. A list of Features or a Background subsection can also be added here. If there are alternatives to your project, this is a good place to list differentiating factors.
-
-## Badges
-On some READMEs, you may see small images that convey metadata, such as whether or not all the tests are passing for the project. You can use Shields to add some to your README. Many services also have instructions for adding a badge.
-
-## Visuals
-Depending on what you are making, it can be a good idea to include screenshots or even a video (you'll frequently see GIFs rather than actual videos). Tools like ttygif can help, but check out Asciinema for a more sophisticated method.
-
-## Installation
-Within a particular ecosystem, there may be a common way of installing things, such as using Yarn, NuGet, or Homebrew. However, consider the possibility that whoever is reading your README is a novice and would like more guidance. Listing specific steps helps remove ambiguity and gets people to using your project as quickly as possible. If it only runs in a specific context like a particular programming language version or operating system or has dependencies that have to be installed manually, also add a Requirements subsection.
-
-## Usage
-Use examples liberally, and show the expected output if you can. It's helpful to have inline the smallest example of usage that you can demonstrate, while providing links to more sophisticated examples if they are too long to reasonably include in the README.
-
-## Support
-Tell people where they can go to for help. It can be any combination of an issue tracker, a chat room, an email address, etc.
-
-## Roadmap
-If you have ideas for releases in the future, it is a good idea to list them in the README.
-
-## Contributing
-State if you are open to contributions and what your requirements are for accepting them.
-
-For people who want to make changes to your project, it's helpful to have some documentation on how to get started. Perhaps there is a script that they should run or some environment variables that they need to set. Make these steps explicit. These instructions could also be useful to your future self.
-
-You can also document commands to lint the code or run tests. These steps help to ensure high code quality and reduce the likelihood that the changes inadvertently break something. Having instructions for running tests is especially helpful if it requires external setup, such as starting a Selenium server for testing in a browser.
-
-## Authors and acknowledgment
-Show your appreciation to those who have contributed to the project.
-
-## License
-For open source projects, say how it is licensed.
-
-## Project status
-If you have run out of energy or time for your project, put a note at the top of the README saying that development has slowed down or stopped completely. Someone may choose to fork your project or volunteer to step in as a maintainer or owner, allowing your project to keep going. You can also make an explicit request for maintainers.
diff --git a/scripts/01_aRTist_bridge.py b/scripts/01_aRTist_bridge.py
new file mode 100644
index 0000000..918f662
--- /dev/null
+++ b/scripts/01_aRTist_bridge.py
@@ -0,0 +1,50 @@
+from CyXTraX import CylindricalProjection, load_mesh, olaf_mesh, MeshObject
+import numpy as np
+from matplotlib import pyplot as plt
+
+
+def main():
+ # Create the cylindrical projection interface
+ cyl_proj = CylindricalProjection()
+
+ # Load the test object 'olaf'. It is automatically installed to your site-packages
+ print(f'STL Path: {olaf_mesh.object_path}')
+
+ # Obect are loaded in a dataclass: MeshObject.
+ print(f'Olaf is a MeshObject: {isinstance(olaf_mesh, MeshObject)}')
+ olaf_mesh.position_mm = np.array([42, 69, 0])
+ olaf_mesh.orientation_quat = np.array([0.1, -0.1, 0.11, 0.99])
+
+
+ # The `load_mesh` function load a list of MeshObject
+ load_mesh(cyl_proj, [olaf_mesh])
+
+ # Set local position of cylindrical projection
+ map_position = np.array([1, 2, 4])
+
+ # Generate local cylindrical projection
+ half_rays = cyl_proj.compute_projection(map_position, output_full_ray_projection=False)
+ # output_full_ray_projection=True is the same As
+ # half_rays = cyl_proj.compute_projection(position, output_full_ray_projection=False)
+ # full_rays = cyl_proj.convert_rays(half_rays)
+ full_rays = cyl_proj.compute_projection(map_position, output_full_ray_projection=True)
+
+
+ # Make some nice plots!
+ fig = plt.figure(figsize=(12, 6))
+
+ ax1 = fig.add_subplot(121)
+ ax1.imshow(half_rays)
+ ax1.set_xlabel(r'$Z$ / mm')
+ ax1.set_ylabel(r'$\alpha$ / rad')
+
+ ax2 = fig.add_subplot(122)
+ ax2.imshow(full_rays)
+ ax2.set_xlabel(r'$Z$ / mm')
+ ax2.set_ylabel(r'$\alpha$ / rad')
+
+ plt.show()
+
+
+if __name__ == '__main__':
+ main()
\ No newline at end of file
diff --git a/scripts/02_generate_atlas.py b/scripts/02_generate_atlas.py
new file mode 100644
index 0000000..9c1cfd3
--- /dev/null
+++ b/scripts/02_generate_atlas.py
@@ -0,0 +1,22 @@
+from CyXTraX import CylindricalProjection, olaf_mesh, record_atlas
+from pathlib import Path
+
+
+# Just get the temp folder to store the atlas
+FOLDER = Path(__file__)
+TEMP_FOLDER = FOLDER.parent.parent / 'temp'
+
+
+def main():
+ # Create the cylindrical projection interface
+ cyl_proj = CylindricalProjection()
+
+ # make an atlas
+ save_path = record_atlas(cyl_proj, [olaf_mesh], 'test', TEMP_FOLDER,
+ map_bounding_box=[-20, 20], number_of_maps=2)
+ print(f'Atlas saved to the path: {save_path}')
+
+
+
+if __name__ == '__main__':
+ main()
\ No newline at end of file
diff --git a/scripts/03_load_atlas.py b/scripts/03_load_atlas.py
new file mode 100644
index 0000000..2a0b349
--- /dev/null
+++ b/scripts/03_load_atlas.py
@@ -0,0 +1,23 @@
+from CyXTraX import CylindricalProjection, load_atlas
+from pathlib import Path
+
+
+# Just get the temp folder to store the atlas
+FOLDER = Path(__file__)
+TEMP_FOLDER = FOLDER.parent.parent / 'temp'
+
+
+def main():
+ files = list(TEMP_FOLDER.glob('*.h5'))
+ if not files:
+ raise FileNotFoundError('No .h5 files in temp folder. Run Example 02!')
+
+ maps, points, mesh_object_list = load_atlas(files[-1])
+
+ print(f'Loaded Map Shape: {maps.shape}')
+ print(f'Loaded Map Centre Shape: {points.shape}')
+ print(f'First Object: {mesh_object_list[0]}')
+
+
+if __name__ == '__main__':
+ main()
\ No newline at end of file
diff --git a/setup.cfg b/setup.cfg
new file mode 100644
index 0000000..d616638
--- /dev/null
+++ b/setup.cfg
@@ -0,0 +1,18 @@
+[metadata]
+name = CyXTraX
+version = 0.1.0
+description = A Python module to play with local cylindrical xray transform maps or atlases.
+author = Simon Wittl
+author_email = simon.wittl@th-deg.de / simonwittl@gmail.com
+url = https://mygit.th-deg.de/roboct/xraytrafo
+
+[options]
+packages = find:
+include_package_data = True # Ensures data files are included
+zip_safe = False
+
+[options.package_data]
+CyXTraX =
+ data/cone.aRTist
+ data/cylinder.aRTist
+ data/olaf_6.stl
\ No newline at end of file
diff --git a/setup.py b/setup.py
new file mode 100644
index 0000000..fc1f76c
--- /dev/null
+++ b/setup.py
@@ -0,0 +1,3 @@
+from setuptools import setup
+
+setup()
\ No newline at end of file
--
GitLab