From ac7c114649a0383c560d2c298e1a79cfca829eb7 Mon Sep 17 00:00:00 2001 From: Sander Roosendaal Date: Fri, 19 Apr 2024 18:17:54 +0200 Subject: [PATCH] fix --- rowers/dataroutines.py | 146 ++++++++++++++++++++++++-- rowers/interactiveplots.py | 60 +++++------ rowers/tests/testdata/testdata.tcx.gz | Bin 4001 -> 4000 bytes 3 files changed, 166 insertions(+), 40 deletions(-) diff --git a/rowers/dataroutines.py b/rowers/dataroutines.py index aed97603..8671ee65 100644 --- a/rowers/dataroutines.py +++ b/rowers/dataroutines.py @@ -1,4 +1,4 @@ -from rowers.metrics import axes, calc_trimp, rowingmetrics, dtypes, metricsgroups +from rowers.metrics import axes, calc_trimp, rowingmetrics, dtypes, metricsgroups, metricsdicts from rowers.utils import lbstoN, wavg, dologging from rowers.mytypes import otwtypes, otetypes, rowtypes import glob @@ -35,7 +35,7 @@ import polars as pl import polars.selectors as cs from polars.exceptions import ( ColumnNotFoundError, SchemaError, ComputeError, - InvalidOperationError + InvalidOperationError, ShapeError ) from rowingdata import ( @@ -178,10 +178,18 @@ columndict = { } def remove_nulls_pl(data): - data = data.lazy().fill_nan(None) + data = data.lazy().with_columns( + pl.when( + pl.all().is_infinite() + ).then(None).otherwise(pl.all()).keep_name() + ) + data = data.select(pl.all().forward_fill()) + data = data.select(pl.all().backward_fill()) + data = data.fill_nan(None) data = data.select(cs.by_dtype(pl.NUMERIC_DTYPES)).collect() data = data[[s.name for s in data if not s.is_infinite().sum()]] data = data[[s.name for s in data if not (s.null_count() == data.height)]] + if not data.is_empty(): try: data = data.drop_nulls() @@ -1566,7 +1574,134 @@ def getsmallrowdata_pl(columns, ids=[], doclean=True, workstrokesonly=True, comp return df - +def read_data(columns, ids=[], doclean=True, workstrokesonly=True, debug=False, for_chart=False, compute=True): + if ids: + csvfilenames = [ + 'media/strokedata_{id}.parquet.gz'.format(id=id) for id in ids] + else: + return pl.DataFrame() + + data = [] + columns = [c for c in columns if c != 'None'] + ['distance', 'spm', 'workoutid','workoutstate', 'driveenergy'] + columns = list(set(columns)) + + for id, f in zip(ids, csvfilenames): + if os.path.isfile(f): + df = pl.scan_parquet(f) + else: + rowdata, row = getrowdata(id=id) + try: + shutil.rmtree(f) + except: + pass + if rowdata and len(rowdata.df): + _ = dataplep(rowdata.df, id=id, + bands=True, otwpower=True, barchart=True, + polars=True) + df = pl.scan_parquet(f) + data.append(df) + + data = pl.collect_all(data) + + try: + datadf = pl.concat(data).select(columns) + except (SchemaError, ShapeError): + data = [ + df.select(columns) + for df in data] + + # float columns + floatcolumns = [] + intcolumns = [] + for c in columns: + try: + if metricsdicts[c]['numtype'] == 'float': + floatcolumns.append(c) + if metricsdicts[c]['numtype'] == 'integer': + intcolumns.append(c) + except KeyError: + pass + data = [ + df.with_columns( + cs.float().cast(pl.Float64) + ).with_columns( + cs.integer().cast(pl.Int64) + ).with_columns( + cs.by_name(intcolumns).cast(pl.Int64) + ).with_columns( + cs.by_name(floatcolumns).cast(pl.Float64) + ) + for df in data + ] + + datadf = pl.concat(data) + + + exprs = [] + + if workstrokesonly: + workoutstatesrest = [3] + exprs.append(~pl.col("workoutstate").is_in(workoutstatesrest)) + + # got data + if not doclean: + if exprs: + datadf2 = datadf.filter(exprs) + if not datadf2.is_empty(): + return datadf2 + + return datadf + + # do clean + if "spm" in datadf.columns: + exprs.append(pl.col("spm") >= 10 ) + exprs.append(pl.col("spm") <= 120) + if "pace" in datadf.columns: + exprs.append(pl.col("pace") <= 300*1000.) + exprs.append(pl.col("pace") >= 60*1000.) + if "power" in datadf.columns: + exprs.append(pl.col("power") <= 5000) + exprs.append(pl.col("power")>=20) + + if "rhythm" in datadf.columns: + exprs.append(pl.col("rhythm")>=0) + exprs.append(pl.col("rhythm")<=70) + if "efficiency" in datadf.columns: + exprs.append(pl.col("efficiency")<=200) + exprs.append(pl.col("efficiency")>=0) + if "wash" in datadf.columns: + exprs.append(pl.col("wash")>=1) + if "drivelength" in datadf.columns: + exprs.append(pl.col("drivelength")>=0.5) + if "forceratio" in datadf.columns: + exprs.append(pl.col("forceratio")>=0.2) + exprs.append(pl.col("forceratio")<=1.0) + if "drivespeed" in datadf.columns: + exprs.append(pl.col("drivespeed")>=0.5) + exprs.append(pl.col("drivespeed")<=4) + if "driveenergy" in datadf.columns: + exprs.append(pl.col("driveenergy")<=2000) + exprs.append(pl.col("driveenergy")>=100) + if "catch" in datadf.columns: + exprs.append(pl.col("catch")<=-30) + + if exprs: + datadf2 = datadf.filter(exprs) + + if not datadf2.is_empty(): + return datadf2 + + exprs = [] + if workstrokesonly: + workoutstatesrest = [3] + exprs.append(~pl.col("workoutstate").is_in(workoutstatesrest)) + + if exprs: + datadf2 = datadf.filter(exprs) + if not datadf2.is_empty(): + return datadf2 + + return datadf def getsmallrowdata_db(columns, ids=[], doclean=True, workstrokesonly=True, compute=True, debug=False, for_chart=False): @@ -2105,7 +2240,6 @@ def dataplep(rowdatadf, id=0, inboard=0.88, forceunit='lbs', bands=True, barchar df = df.with_columns((pl.col(" AverageDriveForce (lbs)") * lbstoN).alias(" AverageDriveForce (lbs)")) df = df.with_columns((pl.col(" PeakDriveForce (lbs)") * lbstoN).alias(" PeakDriveForce (lbs)")) - if df["driveenergy"].mean() == 0 and df["driveenergy"].std() == 0: df = df.with_columns((0.0*pl.col("driveenergy")+100).alias("driveenergy")) @@ -2113,7 +2247,7 @@ def dataplep(rowdatadf, id=0, inboard=0.88, forceunit='lbs', bands=True, barchar t2 = df["TimeStamp (sec)"].map_elements(lambda x: timedeltaconv(x), return_dtype=pl.Datetime) p2 = df[" Stroke500mPace"].map_elements(lambda x: timedeltaconv(x), return_dtype=pl.Datetime) - + data = pl.DataFrame( dict( time=df["TimeStamp (sec)"] * 1e3, diff --git a/rowers/interactiveplots.py b/rowers/interactiveplots.py index b71a7260..f5aa990b 100644 --- a/rowers/interactiveplots.py +++ b/rowers/interactiveplots.py @@ -459,13 +459,13 @@ def interactive_forcecurve(theworkouts): columns = columns + [name for name, d in metrics.rowingmetrics] - rowdata = dataprep.getsmallrowdata_pl(columns, ids=ids, + rowdata = dataprep.read_data(columns, ids=ids, workstrokesonly=False) if rowdata.is_empty(): return "", "No Valid Data Available" - rowdata = rowdata.fill_nan(None).drop_nulls() + rowdata = dataprep.remove_nulls_pl(rowdata) data_dict = rowdata.to_dicts() @@ -796,10 +796,10 @@ def interactive_histoall(theworkouts, histoparam, includereststrokes, ids = [int(w.id) for w in theworkouts] - columns = [name for name, d in metrics.rowingmetrics]+['spm', 'driveenergy', 'distance', 'workoutstate', 'workoutid'] + columns = [histoparam, 'spm', 'driveenergy', 'distance', 'workoutstate', 'workoutid'] workstrokesonly = not includereststrokes - rowdata = dataprep.getsmallrowdata_pl( + rowdata = dataprep.read_data( columns, ids=ids, doclean=True, workstrokesonly=workstrokesonly) rowdata = rowdata.fill_nan(None).drop_nulls() @@ -1190,7 +1190,7 @@ def forcecurve_multi_interactive_chart(selected): # pragma: no cover 'workoutstate', 'workoutid', 'driveenergy', 'cumdist'] columns = columns + [name for name, d in metrics.rowingmetrics] - rowdata = dataprep.getsmallrowdata_pl(columns, ids=workoutids, + rowdata = dataprep.read_data(columns, ids=workoutids, workstrokesonly=False) rowdata = rowdata.fill_nan(None).drop_nulls() @@ -1379,7 +1379,7 @@ def interactive_chart(id=0, promember=0, intervaldata={}): TOOLS = 'pan,box_zoom,wheel_zoom,reset,tap,hover,crosshair' columns = ['time', 'pace', 'hr', 'fpace', 'ftime', 'spm'] - datadf = dataprep.getsmallrowdata_pl(columns, ids=[id]) + datadf = dataprep.read_data(columns, ids=[id]) if datadf.is_empty(): return "", "No Valid Data Available" @@ -1572,10 +1572,10 @@ def interactive_cum_flex_chart2(theworkouts, promember=0, datadf = pd.DataFrame() if promember: - datadf = dataprep.getsmallrowdata_pl(columns, ids=ids, doclean=True, + datadf = dataprep.read_data(columns, ids=ids, doclean=True, workstrokesonly=workstrokesonly, for_chart=True) else: - datadf = dataprep.getsmallrowdata_pl(columns_basic, ids=ids, doclean=True, + datadf = dataprep.read_data(columns_basic, ids=ids, doclean=True, workstrokesonly=workstrokesonly, for_chart=True) try: @@ -1803,15 +1803,15 @@ def interactive_flex_chart2(id, r, promember=0, columns = columns + ['spm', 'driveenergy', 'distance','workoutstate'] columns_basic = columns_basic + ['spm', 'driveenergy', 'distance','workoutstate'] - datadf = pd.DataFrame() if promember: - rowdata = dataprep.getsmallrowdata_pl(columns, ids=[id], doclean=True, - workstrokesonly=workstrokesonly, for_chart=True) + rowdata = dataprep.read_data(columns, ids=[id], doclean=True, + workstrokesonly=workstrokesonly, for_chart=True) else: - rowdata = dataprep.getsmallrowdata_pl(columns_basic, ids=[id], doclean=True, - workstrokesonly=workstrokesonly, for_chart=True) + rowdata = dataprep.read_data(columns_basic, ids=[id], doclean=True, + workstrokesonly=workstrokesonly, for_chart=True) + if r.usersmooth > 1: # pragma: no cover for column in columns: try: @@ -1825,18 +1825,16 @@ def interactive_flex_chart2(id, r, promember=0, except KeyError: pass - try: - if len(rowdata) < 2: - if promember: - rowdata = dataprep.getsmallrowdata_pl(columns, ids=[id], - doclean=False, - workstrokesonly=False, for_chart=True) - else: - rowdata = dataprep.getsmallrowdata_pl(columns_basic, ids=[id], doclean=False, - workstrokesonly=False, for_chart=True) - workstrokesonly = False - except (KeyError, TypeError): # pragma: no cover - workstrokesonly = False + if len(rowdata) < 2: + if promember: + rowdata = dataprep.read_data(columns, ids=[id], + doclean=False, + workstrokesonly=False, for_chart=True) + else: + rowdata = dataprep.read_data(columns_basic, ids=[id], doclean=False, + workstrokesonly=False, for_chart=True) + workstrokesonly = False + try: _ = rowdata[yparam2] except (KeyError, TypeError, ColumnNotFoundError): # pragma: no cover @@ -1867,14 +1865,6 @@ def interactive_flex_chart2(id, r, promember=0, if rowdata.is_empty(): return "", "No valid data", workstrokesonly - workoutstatesrest = [3] - - if workstrokesonly: # pragma: no cover - try: - rowdata = rowdata.filter(~pl.col("workoutstate").is_in(workoutstatesrest)) - except (KeyError, ColumnNotFoundError): - pass - try: tseconds = rowdata['time'] except (KeyError, ColumnNotFoundError): # pragma: no cover @@ -2205,9 +2195,11 @@ def get_zones_report_pl(rower, startdate, enddate, trainingzones='hr', date_agg= columns = ['workoutid', 'hr', 'power', 'time'] - df = dataprep.getsmallrowdata_pl(columns, ids=ids) + df = dataprep.read_data(columns, ids=ids, workstrokesonly=False, doclean=False) + df = dataprep.remove_nulls_pl(df) df = df.with_columns((pl.col("time").diff().clip(0, 20*1.e3)).alias("deltat")).lazy() + hrzones = rower.hrzones powerzones = rower.powerzones diff --git a/rowers/tests/testdata/testdata.tcx.gz b/rowers/tests/testdata/testdata.tcx.gz index 6170c42cfcd250386d352fa73fcb8a039af814e5..76c3e0acba2bb4cbbe0ea247f1e10a75215adc2a 100644 GIT binary patch literal 4000 zcmV;R4`1*fiwFoaXCY<+|8!+@bYx+4VJ>uIcmVC4NpBoC7J%>m6@m`QVPN)K__!#F z0>{o^4A`C^F*3OgDpE)88L1_x%iF&nvfGhmTZQCFQv|C(9}IPMZCzhI^6fnQ?%ny> z-kassW_7Xtb{`Gw?|pal;PB;Dw_2~(r;nHGetGqz>$~5&&2lg8yxx4f|EljV4;G8J zZ{M2J?&^HCHpds|i}dB@_|@{f+bnKBeDQ|%cc(b$-u;EMJUZ@IZ&v;LAKvxLb-uwH zCwRA6eR77kY%Z=&7X*NdU!Onwd3%b7tM#Vu*2l|q)Mj5_~8Mh zpOHU)-0*(y=gW(${@eY3_OF&#`@1jG*H2H52+%zQc}VaQ`M}}|L$-6h;*7?y0h2IBedxvyFWi(VR6ePZhZZ`d$&4&eg1#C z0dcq7z%7@!@%7X7_EqiQA8p@?Wg>gM<>rGtS*8bkM3^7r?(Z8EuU4n0>DHG0tBaGP z=iT~sca|uAvi0WP-SK}vx%=&w%Zt^zzqvVok=}v7@Gfq5`}D}qUI{UF_%+!ca8Cj^ zDUkny!j8o42tIkhaLe=E?r=SNZN9kd%O8LJ%QL~PWBBK1Oy~Sp*Y`Fz3^KQSz{36_jS4*A2Bxx6`ppI;C8rbuLu82Fiinu3;2rekNQ_`CXau0XsPIbw< zm5R7KLyW3XJ`E9Ekb9gBcV{c&9+SLV@fz;2q`WiChP$tbyQS?o?-bl^FjSB`-h(?< z#XUep+$C{OLZ#8>y>S!aUYEQFu86xg!GqK2^Ujqt)f;mWy$kB$WvY9 zshTOC4j`FwqBG?XrKInM?R;dOr6Tgy#H7@#WS+!?g5+gJkIWc-)(-1+&;^8yu-XJZ+KHRKEX?&6HdQ$^&RL8o)#pz|f?c?)f{ zABB7{k;H1kbLc}!|J}XFx6yu(JfU#}UNzc>vgGFzp+l+2sv~Oz5LFspj3w>)Igz)a zBJxR3XF^r-6~vNJI!iMe?Xe>A(F6o_d0x;YkyqLC3{W&; zGo#VICYcv&(z6B9kVl-1Zr`G`!{nXuf{jjd{)XGRkzLLS>Z ze>-{3xwi~R8c$AoQA%cATi5*dqmUN^3_&IHLPaAsjr=&|y$MO;QA6HwlleQ6=X%H^ z7zq$mO)@+em3iNE`%%cVVT+=g=K~kj@O+c`QOG-k5k)2Q(H8aIHkluX{5JPi)dJ0y zt*E>gd2BwsA@Vu*R;)_5w`D6TZD=w-3VCbNP#-nAy(L#v=0mgUMj`JFMRXePtz}m< z-#)a_ehl&q+m=`w-QHSPRJRW;K0gY1GL9jtn&3%U@^h-FbQO^oV})5Y&%2_vw6h{F z6_F1naaUzv+Oidu0!`+}AaBtmjnk2;C3!W6OwMPl6hM; z>pHet+EK`RBdm+qh?FJYe0U>}w+XgEZ&g{>;u`t8lD8F+ClkD`x=K7wCO+R}ewe%% zOQ@=@l2Ep;uJynihkP)EsOp>`v22ZvZ@T>$B~H_kwB!0(mI?D zd2XJZLGmtVU8_1FI>Ms5{hY|xG}=c4h(QydcT7dAsC|pik3l}?-UbYszPHX&(cIfP zk*~;}Cu_(dtLAy3$>iR0tDhev?+scc)oAais69U;@{TGZpF1ZIYx>(e&qe3?IgzhO zw@=y{M^vp^<^s1U?F`6^0T5M}a~HWqX{SHhvo|8TPKX{F`MZ*@NVm6O*o(@OT5>LE{%dt25_xMV{Zo;rqU+e?`qP$~F#>s?^YB#t@FEvo$42aox}VlWKAPY#s${-! zHkNIDZ^s~S$vDwfSH;5FSaVvJwxTF)l3!t^Ztq-Ct8Px@EAr&nU?>u3^6jJdMMrzr z)}kJR{4V!aH_wMq)T)~k`HDO_k}NbPP?3+JXiko6^ZY2}oe7p$6QxCfvSHdep|8l5 zBQYo3Lr|e7Xz1_Q|FoL%-sR4>F6zd61jeE%_xC{$m7$l4&`WNZ0=fa;(qzzM>vuhf zJ{UnlQlU?(Ay34pZ&E*qKDRR^WF2~Ei`sV5=K4|4TVn-PK3oJZP3rH6o+?7`bBipG zD)hlkM(Qg=&#oTy$xgF@z^c$kUsTPT(}0f^q31N$6R75TgfJQOwGH?+Nxd@+NUFIW zv23}`oX~Sc=wsIKqfY9}u3Fg8kApsULUpXl%)&HdH4T3T=p_&Ms*o+BP|;}poYcIE zj5`j7JO-79kFI(1?rPkHdeEn_8(i*>r$O(F_QAbZMkX}lZVu^D%R z=zT)3Yod{}xr(CahJFbZx^~31Bp&&sw;Cm=qSdUH38v{Ki0ag~=AqAXFM_g=anS^e8pYVO8BpEE24b@lZuW!tGxQ}aea&qe^FPU>Bl z%z%$AL2D4bGZuAg86x{wv~e0WyM7e((NILI)9^#t@C`R_-WcdD<^EYZZyv|k)V%2r z_*M*vvhd0xOeS7mTQYV)*QLR=(jBsy-LlqVw+26`i|Uf%S80Lx!%cS60}-c<^cMzZJ7wVC2!DrX~x~u z?fQ^gCR!0xIdW}i)i-xVA1Xp0atAG4R}CA>hHtnH_+y|)G|b{uQZG<;xq6dmulz#juWfk9$q6?zMkS=+FA^M=uHi>yf% zuTK_1E3%$})VthOBkIsQDr(%dt{UT@4_WHHZov0E8L4jr{vdir6NGi~dI`3uhTnp= z0QhD6)EdDC~~vNxW|s)Du-++CN*{ zb~$5cy>v!F?=t$RvgbYzQ2DG5< zDCi|OR3V+zTPYg3gQn(imBi>w(`;Ts=(v+HMoJ{n2ZolfebpNw5kP0brZ&t$|}UHB$IXmBK{m(HxvOGW5oUegwIQtz?ofbUxg#u(^@jI$V2 zSvnyUt+c%t`c`u{j9zl}O;BazhFG-H_Fe;i>!mY{KIAoRb&6TjMzG3&j7tOX{J{#fIPUOk6!n$F0O98 zCXX&J&sN7>pMLuFi*bGKFHV-e`~LOn?Bx3A{*vCubm)zT^I^;L?my|0*;#P&^Ow7= z)@ieskNa&mTtC{+&VTqa-Sp8<%d=&-Stbto(}zcQH)dC%=76M-R5BSL^f~9w(H`t0(Dc{FW|#^x%KNZ#CD% GfB^twggg%b literal 4001 zcmV;S4_@#eiwFo8Wg%t)|8!+@bYx+4VJ>uIcmVC4NpBoC7J%>m6@m}RVG#CP__!#t zf^35^U}u8Z$mBMtNFBLnq?XWKUh?mU?6zguRv~`U6u~Oc2SZ(5Th~{Qd^?Z6d2?}o z@M^i)t}fS)9-@JV2j3h&IC{S6R_oRJ?EB@qUv7Tr`tJ8`yF5rcueXmLzUcd_!^PtD z>(}P2+gz;H=H&8Xk-prXyjWgz+r{mNFJ95Z{V5K+H=l8qZ%_Kwt5yH!k8k?rI^W=p z6TI23J~+c`HkX^T1p(mVx2I2@?N0I4YQ62d^~o|FwOt+V&->)^r0aLT_D(l&pg-yr znK#?hzdrEqUH^S|ez{pKw~JT!*Zcp`^(XuNm)DoV&%4whpt~Nr|BVli9~>MUy?p@b zXXKCXH~e|5_haQ|ic`tj*80lKe1z9M*zd}#5IoQaM&!ruT60ge`F zr+525>8=i*_1&g_zPeaGdiYnzUSB93`uwu*&huf6ct@ ze!Tu}+4b9>mz$@n^~-+w@BYS*NT>O=JAb)6Mw>3O|MUG77PnmD#@A1~H>-=67yqXl z5O>QB+;WK&Ue*53F_9Sjc@WBIyTb}Q3hwIU6`^jaW|NPr$&jh!Q;UAwdo%3H^-zTw&M~nSmA768k zF8cA$yZsjPi*OrWEp-Nyq+P6nIPg?r}EUovny_O!98UYq-af^3E_D?!F@KmUiR3Q*gJzP(ki^5AIkM z_W%`fm&82@l}4NQ#!ZBKUGg5dBJSP<4^E@cJ6F<-W10_lUk`Ty;{#fayn8V+5$=_Z z^A*`~!X|=^8ty?#+V69c_fQdcXB-e{^tvcT_4)b8J3>W@okYW|bJHS}CC~RDPj!)} zYNmKPgk;Kz&XhxxlD-?Z^O1R$ipX0NlTxpec@h%}l9w5gj}?(m=A0l0Rd~h~wd37< z$Xizrc?9Dul1k>2@n6z&=jSue3sgj&jb#+pkT2}Ji!&lm6_Ixaoz97a&X=6$Ews^o z6!O7D5~~T%p${efclRRSM*BhXgvJqg)o35elAlk64y7Wij;s+tRB3oImbB;RMBavq z$R|CW3028g5KBhsEX`=N$BM{D6A-*rrQJ31Ga%1oT*Rm%pObS7lAn>^zNXQh#jpdf zndixtC2wa$zM_Z?2Lk|JG~@+JM(6E){Px_ex*_sOeyu}Q8Wu?`>BE~9`HFmamJK6X zO_Y{=(G2@WekAhFfK`$YFd3PzZJw_vN^67h=~YzO^AU`18s60Bc|nszUS-cSK+%ZJ zj7IyKWL~UE&lX5S9&s|deT&i#lXu1o22GTf35(kEb0S}nJ)bQ!wyJra8BLrBd2I9i z?c_P<-ZCI*JUQt_DVcR`UGv+ILS76o1eMGS6^+<5^5c;ACM1bR4SB~+=I=miR| zBtTF#$?#lM=6%!cMgB=55)m z>)2{(MGorgM=%kIHQrkX*m7^DFZ0Mo0+A+4>u@^c zxp{I1$-A6&t?Gp62#f0Wb0S~UXdewA22Fh4F%_+%_ANd?2Kk(O8!%}4-a1P~b8qKF zz9M^`tRaW2n&*WklY7gpetwX=H)xSmqrIb|_WX>2L2m7oF$lM7|>3 zK51(lQMGEB3*4f#GaxSpKvZ4MUE~&}o&IRg-iYWrA$n-!?@GQR-QI#>FDg%tM_W?E zyH_4w)9{ALbMAhs^4|Jn3C5OzIeGG8ymhFF)%m0w6pYmg&d8pxC_Zn4vB;>&$?@rJ zu@fO*n-4G4LtZe~(yG?9@KUsXeoiu9lZGd3ybx-9c;2~^e*0O?^A*|iNq#w!s>a59 z*D5OSNWP*dtwaM%pc?Iin~ZMXytf17g15oSSD3w zMnrDpryj9Mmsax;c@r$dhA(p-7;~w~yWz9qnCP zi+T+5``laIJRd?)t8Px@EAr$>ve1}7MLve2IXSM)^P`Y=CRk!klokQXhH2-7z9LtS z#GG&sL4}^6p}%AQ(`v?hmpk9Os2lGQ7>lOd-v>QZhF&T{FS%g~=mvaClR=NI-}NB+ zU<3(Cg+8fqkLvjTKb+a1p#TslOw7stCQ$EwVhS z&<8gesjmz@yL!+kJIw|Ht3n@rQ8jN)13p%Sp3_`UpqlFu!er3bHsIGJ_0BLLspfjb zvgI~&LeCYUk6FWyI;k(aYGFe^4*J{))v+ow3)76%H2fK$mptIBLbixPMWgj|Qu8V@ z?l>6o7*rZQy5`Nht8o|VL7&EMaJfI82E8xZ2lrkXnb3^8LG;!D5bLBq*vX*BX50;; z_X)kOiAK)mDvF*P`ccqF6RpyvL*R-UchtOjW1vsUdPH5>4N`XeH>saIdNIyx(No!0 zX>*c#s0e*9BBW}7jg%}x-_TFnrQ>sFOzT8c#vLUEK_=2Z2ivmrgXpbcuWO%^7si-O z=cH!W52N>nf&|@w@1ba6JvKFO4D^@;SXpp~6mT*faj~Uk4x(p+PS-UHDHdJ0bXEgC zRfOIdu%ye#W#GxIant7dVf4X7QC8nD+M+I)cSQxf=t0&af2J)z`C>ZKpy_%^L+h8v%?usdr&A z13tC{twHq8Sk$d$i0ot0#%a{-`ccqFLlLb`!w+S{H{85=W1zQ``)BFAc^qR?^QJ%G zTQMNY!Yhj~nRtC|$=DSIZEe!nkYZG!Cu(#3^wA5Ls9d-hEmP6RUGw3OgPu)vNILYE zC$qd+bL0-A-?dEiDmBlFZ7!YZJ92j|6VdDDdMA@f&}wO!1L(u9Wg_U7yg}=w8Fy2+ z>qBmtXhl%v$hDzW-`o{_s0e+?9kg^^HEb*!zTr0DkAWW1FpE=3y+GOJ>UlUs=mSosp-NNp2GGacG7%7?&aRKt(lTd&o=pIyH_t=a z<&2^A(isK4%jl!ZuJ_1AYwn;KccY+>#-lP1-&5HV3}`{yG0-FBhAK*_=dvXj(1Nz3 zpqJcGg>+JHrD)_1nwmEX`e0Zv=%ij;(NvA5=8b}$l5t1M)(YOa=E$A8ODE=*iJ)E6 z)_Xshxt?0kb`X6ovSv_)Z+!60uAc$=Xe3#8I;oF-GIl*RHE#$#lM!om;hO-V+4VC( zFS+mvbm*}t`n*(78}LU!9}F<+`r(J1ygw22)Vz5^=<`a}fI)|z<7Cjc_3MYxvk5}F z$X$9-N^1DI6GcNm0(xgW+o*Es1cxPEI4Z?U()M2HTg}}tddbx{L6wmkV$n+5dky%lm(DQykk_=;MeZULUF>XD8h$NVIs&_u ztO32sksDDa1o=8b}$cZ-;-jJqhYt!_DSZ(c_4a^aN{y^BSqzBNi62Yt+iS30{s258OG zChpBk=yxkwD$?p{iznUH^@pD>`{nld zi}U{LAJ3N;%XNS78D@67}w|iH*IC%=14#}9U>SL^f~zE3EZn;+8C_&r_v_`&}GkWP6b H#DD<+{U|`a