From b9db6f3e92be4894df6a158db2bc2e5fa0355aa1 2017-07-17 17:09:09 From: Tom Bannink Date: 2017-07-17 17:09:09 Subject: [PATCH] Add CCM plots --- diff --git a/plots/ccm_construction_successrate.pdf b/plots/ccm_construction_successrate.pdf new file mode 100644 index 0000000000000000000000000000000000000000..96aeb848e3d2a6fa72c3396c76f19eeeed16492b GIT binary patch literal 22986 zcma%iWk4NE5^e&)U4wIQhl9Jj2G;;V4({#*cXxLuSa5fj;10nZf;)MTyLb1#ef#4v zhneY;?&|95f`YHf_JvHK6pJJyPv9g4LXmFIvQY z-4YH52)rXqhaT+u)g@0#B1YO18v$77>pLh_t2A!-eayv&1(|`HaSLhRJZx44@ zXLg+jZ|}#j71RehA)Ocdj||=k-uD})Y1W{H^NWcVmO*%3ZA7T2$EHp*+-`LGY4jjP z@cxjV9~i>$YP07?5C0;j=20yLLHn!prwqb!z=O9qs=>Y zdpt#n($rwrk9%D1bE%ipo;_dl&i0HgUJ%ZvV>%8^e@bBWY*7(k>(22EqvZZI(@+pc zK`0W#c1}3R;IzkmT#j!iCyF^JCt^}J9_WwCzWb_wJYyd-VcovX2(R%{k7rmHo%Tv@ z-8*yE)%9u)vpm-AmyWR(Qu}RtZ^Kl~Mo(DY+WyR_KZY|^KC#K=WjTLN6ZN~i?%XDq zEJlI;h()!DGs>Qdcb}<#Y6=bR)@hVMtD!-d8v_$_!^RWa66^l%sDIhRVPGi4yA1R= zI-?Fl70{S=wLxcG)UW9BU~c0wGF0MSHhVcbQ;#}a8ZX>i#I8-a^|IOI`055nqVz*# z3SH7 zMBX+~1;1)`W=(PSgVI)?svhP<6BbVcAG7>7#l*&^%Cn{GJ*NW+|WUPOSI!Lr@EU$IV0uI^)px^oAzi&t{PTR!tP@Is}*tt0Hrv z(r~yBdVQ?W-CkCL1a9zkS07Conq1O8M^_Gw=Clatc*&=3bhnwS2f~-`((3x@ZUpPcIS$un;q!E81^(Wz2JI_JHsuU^?%wbt216}(z}^arw@KxkJ;7> zWTr_OQ|o1h9moe8Q0|cj+bir2EfJREiG&@k^`L|0Wyg@)p3B9uC~nHs=pE%3zi&~H zA+H*Uw8S#dGeSv@p|joDvj$^cg)+Ndk>MDpgK`Cz;%M^H<7=W)lo(+ZkKYw+-}iqV zzO|^$!Cnl)sHPR0#JI3KWORFg8uuxFoK`3Kla&-+94ubk(=?C06*W+4J8jYh4?{IP zS!;VmnIk5xmBQBIV)`CbmSnd16FnAoK6*FJ4$-&O;_||zF3}EB2=B0!x+2}@XE=ru zVv|t3s*;G{d2%@zF)VCE0U5G_eDD-bFeRmkDT)X>dZ;wq$K326>Y7c+5CN2Q+QGh? z@w4>GF~aAKwLbxaA8;h$5Ph?e3oxKw^)$^8ePyzfYh?h?un~qC%tfz^0I*j^&B@0T z9hsa)s|KZNx>a;05984S18JId9fe<}@`f|nKWDNBOEZ`%#ncVwvVT@#h0$W7^I)Q< zlq>Dl$Tf}m`c4qjGie^)?t&*Q!8Fn~!e?5j7J>(P#n?a}kyD4waILT;7IU-cVL9g`J3=iDK=UAV&l}?;x(n@2!@dGlq^Qv-^#)C_=vH zr}HFs{@+IfNR{oLtGL7`7e9(5RdnQ13`&-_&7-vrm{&R;?{128*SuUV*R7%jbvEYQ zPDA;X9ANbD@y_3v!eXb8x;^JQ{&h$+l8%3P^E+Kvua_K7Msh{@ z=xoB0g#H`1)_}9SG=1NkmR()@egpPJUyD~;`6pFl*M7&<*rT@uljd_)w_RV9YdFlB z5`7qk6D~<>3eP7Ry+hKl>@duVpScCPxA=5p?iQ|)_H06wDgbdxLO$!!0>Lzd8~m1W z2*&w|gTh#125xn(vgfsV%Gs!pan<{9|`ppQ{Gp8R5GT3w$Yb1Z zR3IF7bl()t_%#8{Pes@%ZHQRe{}x@MO?(W$F`#uFWPVEkq|5y?OpMg!$uE7fqqJbm zmWfPUN>;DALUwo|(N24KA?`$+3ahe1O5`;M+zuO=&1$mxqyTTJGq5m{V+x2d%y83e z(zQY_TzbD|S%g^E5h{7xibYC7L`ojsG1|WFZ4zq-YZ#6^FyFXlt6YB20*Jn+Mz{NZ(|`moT1bO4&`HKH%A-u-z?eVY~8B9^Kq z@l~o@l+KpmzI)$jo47IHO4e}0ZUh^8T-|mjnnuLbRummQYUq&Y&t^wLdEM;af`Wa8 z;^%}#1alf|{fYi$2l%(;LJg#(l+%C_263TBNsqj3`Iy?n*DbH0>G+x?iVKck^FpX3 z9E5N74+%}mKY0WFt*es}6KMn3zH5^-kSErSd(mkAH>4hV1XwnNP!AmRa2#~bIwhmo z;ILL4ba+gIm>yr>K40P}$6=W{DrXJ(Dve%IzOWP0(aC;&XNS^N2hl}B!7hMcjIgh- zsOIP6nCveWSYhd7RNvW(R+a_2Rs>>Ne0Sbki8g13=qV6fbX(&bc@5}NgTkh@4Y`Jx2FL@wNm0?0Y^i$41Ns*=6n$ zu89j1kg@3|>{>Q$dt2AQkEW;$*D+g&beY&M%t5EQX#%`Ll>>m_@4fxJLdJkWCu>KZNfc;-a$>uF zjG$U_V(3sP!-UsedDTF9)i9wV4QBKKQ@o{s4n@j;dCh}CwO|0~P~m!m+~+_l z^qCV#Xo};N(dVi){Vp%%37!R{gc6v?6sE>1Mx&myz5PRm>2DdvKYYf2@|geT0qNP;nf}&)efx(#%Rgl} z|1HD$FPXn28N{ru99{+gw#`JZGqD5M7=ZxoKhnal zmc5GG1Kuj&zsg;4J8MUqe^~oh|9_}|PeAE)R69K@dz-hh3|#*vO9CA19RFDHHYMR# zDG`v9sR2k)T8G!5oZ;Pz>k6CfPw*RXNe>L=PyZ%=h{x;5E^*SSf zk%g6sA8{u~_ALQdn zD2)XKv{8Owh$suh5W#SK_?-JQpr&fG<_Ly5vjxS;o3=v8$g9#FGXBijdwt^ex9)ta zn`7%?soVA8totnZhj%kSpyBT*CQ{UY5B)TLc<%!DZvSp54D+sG#=8o-vjbZohWne> zL;54r%?li02C-ALa6*S?EbHM#`z|o)z2e+>GMDXe$)>qPnsyY3*XWCnY2=OR=vkM? zQaC5vX(if>0q{59uJLFlw<~cpM$h|ZeQ%)qc<84o!lGdPO_QF7X4lNq9W>koT+;I>bMfsT-x^ImKH1LM=#EJ5ewz5*!lY zF)9+_4`_W9W#HI>hncuQ$nXO!lrI~BS*9ICb~gA=v`*rBGEZ7_gLQuIsBDT=Hk@HI zO`fEfugzbOJD}M9+;pOp5kb=m%9t$c{I-ZxUnf7b*x=uc(w+TJ}uALgO)yZ zaTSuLeX*s*<&{V!^@=ruw5H5!8sMXY7q7JUJzKS^< zJGL4;1W8j{`p_!PJmKRsA_u35dSd$Fr~^b9^W8q|Qcz2T!p6nZRS; z)UNmGj6b7=?N=YWlZyN883RqqO8fEl{In78i$mP~s5t*IVZ+(NX&qXNA88)_c0p&| z2I5?oN_CKFx>KN2i&N!|HXVC1cKIXzNOCk?2bu=}IJ)<{IhI*#1bML3gz-X;i~lQp#jlu^fq53jJ!y zWo{ya(Q8ccn`{k(k$#z6&_H2>)0(r>&c|x$4&B?;mkF-0KubI>@7oF}!i9mzsrenm zp_9w4qBwH{e$O8>)#1vqWQJ!)?rV+b8e{FH`^;KMcmi`9a)be`_QYq$$o`)j0hky| z5+9cV%h^BJm(J`tRp>QfrH4?ReC#rf+MRMv-#@yU5aTY)4r8PyEX_8HAD`i_Go>lU zgp3h%NDFy+NBTEE2%{ReX+LTH;u1KzaavxCaA1aLg{7Z1pK`^GEDNTVDhd;7>;gQ* ze1imo>=HmwxY*y|5f}C+f}dkAVpX>`l`1nSm>XNQHXSZw#qU;0>0*UJOba!)y0J&&E=5E_kXTuK}-i4*o!@N$`6 zmt&XB$)>f}&u5S2sU9!>{&DdU+d1i*BR7A0mHkqjcFm^>UQ8&VCH6Ksu?hA2U;abp zs%ybWi#*2x{XsX96y!c|l(q9QduArc&*8<)n4|hh)3Od}v7RYp)bcj*3brEBLsWGh zV#%Z_;F|;W^mG)5@DMAqf!lVy6Crpz>8WgrLaG~09LkOU@Gv8GNgTM+WS-=<;01Je z&^QWHunP7&(?{Pxdt&_epmDI9E?3x0Rt*kI|>vI>7ShMO_^AizG*+JrF=#*#fz#gq5qP?6Afla0@HN0GweZ{43 z+%+dVXUQccKDg@kIiCV@utc%II4Z0NARe%+Nm~yF6Q8bN$9}v_e?9C ztgyzZ^xQ>-pW#z2Y4^3X`5v=M9&BNTO&9s9IpO8LORYCNnyK^2;+G$UByjGg;Y&V# zjqe63i;tf|h)z*wt0Q(z7Z|+AzgOE_YW)=*S+)^*UlTtrOoog&Kbqz(xgW0|ulh6; zBzL>w6(-7z2^^v62T1BMbQ)T6m04` z;?d?Viu^{H?050dWFf&4nyfUf`hBpSd$2uB2;?Y#9KwyZ62Rp*$GCKaT0ajnBmaE<6ozAxIX&_wZbI_+vnPk4x& z{u)i?uX~EjCgHXI5t;*zY5dR)PMCg*rNLW6xi86BHhnRYs2sI!J3j(A3kXt}qpu;d z=`sBDZwX z%Qr*zEut=4T(H-%PxMHeTN_1frLVd6fVv4%vQnZxN!L-;QK7FikP8alne;{?sFi2R z=%AI9PR>+P$a4^GLm#4TJUb6or$pXPN{pyVlYKv4jQ!mOD~&k>wavCww!A(X^@Wf@ zJ4y-#V?$89kyO&tedb;uOsZhVPJW+W^|vuRCDk7tAJ zKAUY(s8jBb>6lQ(wY_@>#;3-FAOcV&>-&% z+&!-c`iJUnW{Y=jB_d2+m+ZW zQ(JQ`4Y-7s?T?$u+=tX3?Z5SkW95>iR8JW2!@Yo8G7qAaJ&w5$RsU8c5?&Y5!P(HO zc8*MeriH12u;6z>=Z?_sq*fLykSdu&Y*-N!1E?Yk3qN1lPVgv zJN&Gn?VzF2*Gw=+>|9wcsKMHtTF$jPTE4VGzIsT%Z276_=CGWB;}%5VL@>s4z!P`x zbNyWg)D&@d%qPylU6!3pH%`!tI7(1f4*N-I9G+pSLh$I(cy&V0K*1-oN=Q!AhnL}t zt;{}C1Ki2x&}7mf`+@WRmYY@3%h#8h%*i9(aP3%}v@IA1U69!#h(;PvqMjLcko?U zOTKeeBszFK*VUfWOWKj9C0=svlF7xT4TO+NumTyhRkw=|6lSFb6`yOCr-Ft<*Nyt< zOpyYyAS#PORzpN&hT=&ap&}7Xu+%w@whl-;^iMAOs89yljOV|tvVWFeRElyvj*}bz zHlD_8-(Mj=++*<2QERpznuFOko{(74-f)Dka*202rxuEKcI@;!;i&9%Uy^U+3OYsD z4svdd_c-#=5L+hfcbt*!@()N$9x$^N*)DmX^12b=hIGZ2_TQLTkVVD-I4iqChnU82l;Gu7WZ$ z)M=n;@2b)0#K*OS`p+EcR!OriIgIMb*0iy_$cuRSzUDSV`HPD*YV~+rhKR;)S*t`> zDa~PWYXWzR)#7Xpe0-$KWJ3m=!v>JUd5G!1n$-<9XXQ@b^$F~#ize7L>lZ3&0Iry> zYcRN|{iT!&zI7{%sg#MqvKHbjOAO#yMVf^I2UrH!`fLfSOVQI$t&CXMH9Z@g*q9U2 zN;AKyEc4P@UzxBZvgqi?aMHB_k|;XiIIIKDyt-sp6q}S#No0E@<2Qdcw~&&> zH|iyStPm7`nU{u?6d`SpVGYGCNbQwXtj;eSrhNoNUrBhGlO>($-0o8J-uzsj(+zgS zGu8Rt+yt6UH8LZ|#{JmWzR&aA5EPuWzvg}R%-?puM)lc%nW-p-JIH*6*o}@_-)5J? zJ|p-0+I{?u>8eh~RKduMX!?oV92Elp^bU=L!XehTt~_#8ocd%9$k}m3xa!_9V)IYN zJ~>e+8fUc*XBvPb!ny&qD(pwM@eLlaCBVD+V4M@Z6$-tf9&u^tZ`F}F%7SEMHI=#+ zGI$r5UB&X}uw+(1L6DuFf}Rc9frDV5B|gWJQGx;VC#Wb%>Q+Pj7o;&gXdk-nv&TKa$oqEGDDQCkpH z5miheW{W)?{mglMJxxdJqw(&7_}_8*Gii$9R8%sF#M8+;=Z}#4{hFsWk&fd(jS6^7 zG_9A$EvBt+kBM3i2)44WjjFy%S3%QbX53ArF^7S`*-5^u1K;Vb8dWCRx@#mLGp)c7 z@e`cTXYp=}Bb;8zRmkB<3v&hLet~Z~D8?Zh^9huVaJU!Aw+}`rEX>I;&uOUw7rJZPHupx z<577uKH$_x&+o=^<2lWk!msrn9kK!e={R81E*r8SpgG}~xE8cF$wYkuI01l92guHY zjs_T=zlS;89?kE^;|xTK{- z%p99P%Y{f5IOXf*Z)3<8C>G3_IB>P;yekNuOU03$nGw7F8tGuY`!q~>Kwnx%U8g%k zKaVVnix20P+l7NehNBQLUb-$HB_dj@z7Rs&rm3_k0Gi&g zgWbH3bU3SMdU3PU)8nJ>=_JtZ!}(dce1ax6E~r`G^dgb$Gpp{> zCF@iPE`ZHFiWn6L89^!$t|f`n#67=DONs{V%))jfV1YO9^qduii1&hSzm@cI>^={- z7&iYrwe;=%mv%ucKvo%UA>AFjB>+X=lOkDOgO_bn>pJ_6z3|h z0#WPg6dt@(RJIJwqWQyVmKw%NqSxPR9VMoYV#8erFrdCUA3KQPoxGOgKRn9lOf8aM z(yH}7y?U&1Iu6qliV0In+Az+Kl+G=!yI7HCwiUo)sg0{4zftkT80 z!DQ5U;{!hT5_M*oI*mFd428Nf`r?XW7x{bT>r%1VNW)1{mCpq?Io3ihn9p7i-^0K9 zDpi?5WR1UX@mn1dTX4D|E;BK%tU_1kcgM7+4s!vn2q;5vzVC8uA zUd7L+y_MWgkOoVZoE0b+wwP|Es0B3mffK&u5V4dySv*D_r&^`}Bq(~% z;~TW=66sz{AtjB3J(??7D{^3>pk#(B>)?C9QCd%V;JzTBnA80&u~Y`KLk#_S>GKJY zj`=h0yds@q?RGXr@R^C@V^a2{|5B+vLB-b*A{zQ2o9kGmW}QPDPJcwsfK+${DoH=d zCF!QuziurM>863hZsN4bE9UZOU{B`N*k#vng@zt;|LE-?DXy^6`KYy-m6xbklv^0u z-@$E2$5yBrhL+T~Pni=78;4^r)=kd8#SSX3<2*{32l97@ek*NMbPI=uy|dQObB&0xa-fSXzgvWuVn)mF+QS?iu%VDZ0s-a%B2 z91E=M(Wld;HJ|Fp|(F|ly7oUEf=NNEz-yANPmLit7Wl(MGx z0E4Q=TcQIeK92Co1JxwK<-nYzFPn&0bjy3gp#CZ&TjfFKORJ-vC3?h6Bj?jmfMcvJ zeg|lIY~uZer=+j6=a0GJXn$NruHS5N3NI@ki6J`HGH1>j=u$Y$4t^2j4@`0fStGxE zlf_5k`xT|l;4QQ6?$z}?=AK0BnHY{r(CGZ@qGQe+-@q7UN0$Lv&;7V;f!|=gUsH#3 z9G%DUZhIORRTmsOQ+;YcovzA(QlB9c5q@7`bITAtMdFlVWYFM>?xBQ#d*b@`q61GN zVxYjfxhgK4v`)Q)MZ>!PVN>h+_gT;v!mQYnZPN?*z>J|^n1RlN+`omnTUvJ&-3i2@ zT^_jekp{B`Dvvt&`+v(DJrE9DALW4`NOOIKbS2_V)M7OPez6d>qKXkOH7q&tzF+%Y zNl+V*4o!b5YBTKIhWs0IgBxbRsJ|obkW9u~VABNk-A$(#6WUFWt&891E>)le z3^+=K1p)D{?@0`8JcyKgAa3+1eYgeYWs94y9S%CaFXizF#OmyZemjA>%8gwFqAiX)g)$tP=%!V-1CxFN%wRRfV?n&(- zjb{imjkJF8nFp0WM46!22fGS?NrDJl0Q|F~hr|<$l%iTA;d(Pd>IN{7vK5jLjn0{@ zo!bjyhX)sO!$Wf;DnBE{u_8D@DqM`-(MP{L1DRnKC<33PpA`DUyhG4->?wkY(X87FH4*9L#zW4v;gN;k_d!OG%u%ip zCt-T%%bk3q-*)F{+QTBkGQyv!Zj^W!+Hu6PdCEH{`QK0c9F6S6zJu<8xQn$B5x&BD zA{FS5hGL{8iHg=XL|T;AtWNsAEq%dzPjUaTNlu|e9>Q_@oopX+{^fz^mGhOU8)T2h zgYyHkfw5u=P-!q8s51H)Zx+XhIZ$}?V|f^FPk+pW>k6h7R6~$^h`V1)?vD#3D#81g zh&`_FT>Mdj7M2#$<(zXkzrk38&_KiBe3V_2pL*qkj@lr|PT!@IXl`eO8OHD#0_}{z zN=Po{JNQIYbL%~H2p(y+vFr{QwwWR|#Cb20TPB5;k`Q}rkTRt9CdaNAm!u}9b+Ixh zfsp;fhK%6JI7ZapziNzH=n;=Tjv9|9jk4}S26B`cd0J58>h~PKa{_CL_6&SZd`>Le zUE69`Eyl1%FT+siNvy) zVrIO8pRDyr?1*5a)4yPuqi6$!Yrzs89d{t(w{i>#4+ADiqn^@2!2R}gY0w#UyLy0O zFBn`lgp~mlL+DOi7XBVJAKC*x+Ne2foT?a7m2Ycv#aZSWO~J;mhJB|IgG$|?t%{ui zE@LDeh*u$xSDD?VyJCxb2ex;^>s8v248gCwm)2t)4Xgaxx59A)&V-M}ggLbsR=k=B zbhUVQZZ|_cTRe@pb!sKdE8h>_g$gA#{v|$)v6!gdU?J4aR}>~S#QlkWM&H|Cwa)}6 z;2=FSd#+3PqHRF}i_rUheP2;e+xb8?&7mLbIoVx5?=MyhUtU9%W5MD9=6jNrWxuAE zEBJZ-wFI`wsC~-jottY>t0FY0r>`3r`$(RRc|UqzAjAYsJ%pboC|6U8^%}ZJ8`0=D zXuW$0Aem0032~*6`X8CEIo9ET#4q|Q;xj-JKHR`AMwiUJFT{kHBg4H@e9rM#rz$Jo zi>m}z${B(+Di&dm_x~k#HqP0W|IDEwYFC4SLN8TIBa2#%kAfh$PPMz>Z#?a~D! z)BgU9A?UsB;QE7*Y)d@6=Za*fi`qM1SPm&>I&G&GaoXn{+b=Z^qxeUX{5f5egANlh zB~eSx_S0;!RfYxY(q?luJY{oWNtbdjK{>y~v1f34_D!e`6cwc-M{zKXx(tSf;Pp~IJPfUhpg$A0u0(hjvf!$p@h7jY z)e=#6<}Xz_N!bkd9@?T9=_kwtF@`@UiXS6=wtGtLOG^vNIpoak(D)n zipe19f39Z(^r@)aAzDaJKhCpV4^*%N&}hW(tJ#}N%|82=chf&l$>eL0?Xt6~6D7>d z-IJ7lS^pVc*sS>(M@{L27(F6O%jZp$+AVF2m%_;{UmuAQ?$`KVvHt-kOY*USL$u@a zCbeP_3c|p*8tOXGP^j$DQvdK97bjVqbjff`1BY}z3ItS+fXGwaGe@kgV7k&Ls?;U4 zRPtL9CXZ6+`cRoJj1+8z>oQAJH_08vN9ZK_56`l{!Q#Ju8CbVwNJ*wu$cOt4MnyF= zM5oa6kai%!fj2Z9E8>3^WFBg8OtJ^pV#E+S^`%PZ7m3hKm(BJ`Z>bLF;ERBO=Fe*d z6mmWCDGvS5#AiX0r$)&3-6|-WG_$qq<<*ffW;oV!D`q;lB9c}K+R)}A*v;3>&}{a7 z;W{Y#kA$T~T%^eLgqj}ZndEGfBGVilq}RMK~bCX ztb#qykjM&b#H&tm$1S9~=~5l7B(m<&>;;eE`|9|Z?<^gCY|Vx!)!e}g)0sUm5Ne{# zAo#ci8|lD|_zcFcA6>Y7Apt ze!ztjp0PRr9t6#W-tNLJoa)%$%ETHmJ+hzvxF1zZq9DwzdDY*F4|1Mj*m&1%`ytRP zWCQL2g0|>mOd4dkG>Q`dPn^Tg23t^&LBH|qhW!Jlq9-ILbiwO1-dK&{ozoLQ%>AON zrNQcWdeQ7FrErJ`7`Y>~aD*!Z8BZDh{+tYRdg^g*Vv26gA%6zml&Y3Em^eBHDSnfe zrZ5Jn4zu$hPwo-BgCJzV?1jaky9=bQC?H-AW)6X%PswF67g7Xt8j2IK2r~M@E>AsY zl%*U|Ym^Ov9g2wm2DFr0J$6gT6!Tb4dpi@Z=DKwr#c^KI0XLIZRy*0&?S2}Q*T?FQ z=Tr}>-{+w~IOR|lBmhk0y)V4ty4I4+_nROwH@wBI9uBA+bUm!@5B{>yux{yvY@g0k zevyf9pr2w8S_Ynhp1_V3ebx8SgrDX6pA^v#9ia+=OB?Gbx%ws}Y@<>NCR+x{vAI8n z5#H|LG(a-ZoW-!WJC)-bUK6i!l1KSqAFKg*ghPEY=%Fiu;( z0%?S}5)>`#5zH{eDVI=|Q;G16Tr)kn-i=jVDnC9*{)zdm^M4MdZO5Ixt1`h6&^st8P1Ebhr|qzVx6q!c+VN z0}F$50PEZT(7yelz4NOr)vKeN{70Z|KR_rJZFaH`Rlk3YY(;+<=n)ip63Xq3&tlm? zcFtb#i%z>(KLYZdeHi%ITfFS|jH17w_4EMWhlJsEm_~7Lr--X0oejzJrm)DGd-xAt zp5qXeLDbDLR=((;gU2LWOYZm3b>l585%I4RG=}1Sf)oxV82u%`MnAUO<{|{)ZoH>mae}u!_Jc~!%x)FK#~@@m((g7^+B|kef|XC*jG_um!lSzlzqao5 z_RaRr7d^h-o>R4^HBib*36`jtC72Z`mdMvgm3}A?eMoau=aQj=6%>M14b|YJlc$nf zNdP003Vl)U9aWwhR`8Sr<$E$I%Uy_9M8JxL;2@yVFC>Y`B*=3%Tjv`$G7S%ii-mS4 zpm^P_ZXu%9RKrS>I9SGKv!amMmY5`HF+E4UBkEU)Y9vc4Vjf+S0n^cW@I4as(qS$N&r4Z@f zxKE-a{Xj`aMFpG80L1kx&KrO)HOg#2QuV%aR!mDG-6TQ>U45+)~9qg70`0Dm1fw$*6jWI_PN!=vM+buE{iq?AP7ZW@X!}w#~5W^ zfx|M1w~zT(4uMlRLu^96(`?iS)f{YwEzDAdywoHD)TiJAVrK*suC#%&WR!Fftt1V&bW(l<)~E)SjTW15OBDdzIpmauAB2; zkud4hDH(ZCf^c7Sk>*sg2|urW4Kx6enk$3y(VqkTc%lmGrs9d)tP)LGFihGhU|aO` z7)qk`tWcytWQr+L38YVFM;*YBoyw|+D-$@BHpBNU%8V%1J0_u}UmT{SO@-|xxd?Dd zZF~_(WpJ$8g$dY5rvPKOQm2sj(zNs}Q$zS!4p%LN^okPJ2ZE7l{SN_)q&sgN~ME{)}SeyedQ6L)}j0CzaZ zo_+eNN3LtyYqc|P8=(dv`Xx;GC9G4>ugV_(XtF&|5a=ef&i(pGM7B`xiQHc^a@Gk(wm`8#HypAcQN-d zjbG@A=Ruo|B78tXVwyROgHT6Nx3MO(;R23GDO(O}S=R*6Dyls~s97!LML_ljcBVEC)^=}8{Raar zr)T*JRaTG{(~yy(5wtVavyfH(H&aZ{($vBgK>3GG^%vPf&-fLZ%=B*=p}%mxi8=Cv?;Qfye z{_4(A-{B8NTlp0o{l@!471FZ@{lSI*Pk)%&**gfE=-K@}c^SQbwHcZIGW-AbFo5wt z_-#?E*AvR7R>lAZHB&1=D|^$w34exFGc|NDvHuqq{0&_F-xm}68|3`2J1Z+2;P3BO z!217e{kJ|4z|Qm*C;xqMFu$4j|8u>K^Ew*vP4?~mFI!%1`bWQi+wzwL^K0L4@^6xV zY5&oErTjs;zoONxRjf?k4sHYfWSM_19OK_75dN8V{*BH5pK0eUoc~6hvoQV_==?-Y z+X7Dwquu#t;K;k^MKX?RL$2R^YPChHAw-+P7S*}moPx3TD63Nj%z@mT8xrp@5H=;~ ztFFML)Q%4-`**+Ilx0s!XIrngq^-@F|*hl|!&~2puPla1(V0Rfln_L2D7@ zOg2|;h1O$3{eqgU?D*ijlcb}QRTKUd6Ync={m9#GmHTE7v~~0AGw!>!#0LJ_J#o)% z@0rKnS-(&AlJ45cb);9(i5)3YF!b$99*1@J6wjriK_5!#l#?ZXQE&$lq~GYju3Ce*_@m^8h*v)fQU_1u2KgGm$`~I z49DYfw0H^BHx@{l?xq_rCTd)xCyZo@8nl;~_)i9(l@*J@?j5|7ag?0|G$1iusxB(o z|tOfEEwkw)iR4y3JcpsC9}gUuP8007K?;bT0E@%7(Zsye``T#Y4Du@M@C29 za}KP4g|+G4F+PBU*0i|ZfTXn^E!Ng6#t_0!iCQh2{M$m8YJ(jK)i@1-Pd``|RhP^0 z7qLqAI@pHo?Q*f~wCiUd>yilKOPj7fGjNsCB_eAfJTg5yQppq>yJcx9!xNiUHqj%J zW+i?%n9I|t)03IPEr@Ro^`ugnprypD@31o!Bt4Xq9vnK^hzPlgGrO6)r|;^wl(@~A ze=Ox-6r=Sx^bsip3qV{d2H_T)Wfy}B(#R)ere5Z5KQ*p@I)S?}^Km*TF|{_UH)5e% zXfb#+;6lwf#b!7>?%fSQbll|>jH#EK>E}PWBU%i~kT4w4>Fw$<%b;dUDjL&ZYkqaa z+_-l1t;KnDtmGX(%+a&fbudD>|9(9RT~hD4BNeDsZB4)@lzSaqPsb9}3JPo|^LRF9sYJk{mTi;FKKQQ@$j$w>>C{YDZ2WiLo|50hwmgPSWtK50leBHNlDEQnh9nc_<42Sb2SiDfC}v5zErUqd=1Y)3 z09{yRlVhJ_k_?hXZY+ih7i;BHYd5K~oP`NN8^eQ*E?W<$^!St9vm#qM0)DF~kxxT9 z<=dv^xRIyF6g6ykR6m%T8{KOKp!h>W;mf{yL2D5`)RY#J-y0J97nRORx}%CO3b%>1 z_1N-Kf~=8t358$|gm1PO3gy42=@w&J7KoK$ncrD`r=IK2Chh+qtgg)6N7!_%UYp7?s(p!jW8O=Pq51`S;+Jo#MKm4 z&TK}>oy~wv4d5DL+(MmVV()?Iv64jojBuX5h>ilA6IjQNCQJ+W$NDFYNFQG^R{1C& z-q0W6y+*@7e^9dOeY!Ag)e72*J_kEQ!LZj+h}fJmh&mGMYqtpE$p~;dTr;n$SKtO& zc^X^Dt&=z+t(<&=xm)wJ30c;1RbDO~R7ew;#lH?y(7E~fGlFe&$9-(z)c9H89-Rv5P@SZwK5t|QL^2!3L0Y$wDO@e!^MX~valbyy=ah9c zTpTfFJK?3}wJILuSKa27G;h6l{cdJ{QvHG>aYc(ey|LLAu_f$^WnBLNM_xgUW>-E_dDzx3voi5QCdJ4v! zJ3rlXU}P09U8-d>*Yn}{s^=CgG^r@&W@sHPvJNm?vjjc)Qs-GeIkQ{8c$3o0I0L4Q;wqNwvwdY60((cKOQaMXYlCnTK>?#RO+h5w!(zjw;gT}4;-o2^ zY5JS!H`8*WQ{~4DJ_;}=aQNC!q@$THKin~&+mKJox>#sE1wLF6&Xlf4&^ttQ?ebJ3 zTI97w6BBsGRO^H%-oc%){Z{=#!?(Gz2ZFgG7ly9e@}`KH#>VqU(ixN_818|h&Szb! zLuyV%eZoJn#2I}GeYb!RJXezyQaBq(W0f@U?KjYP4gUey0KlZsU!;ZtW!7-v&(1QC zUp0NdMHib$O;o;96=-*;X9}qncNsJqCayb5x^-5!koCfVGjWV2dAh$)YPZ<+qg2lp zmstZeOc!sYDA3I}x3R>(pm&W>GwwMQfuv<0W-2oKVp`DGEcAzAKh_A^EITj!gl=A| zTJQ@U!J6~baYm{rne#6%3L*VW2)_0Pn1vVIB0=mX^eq0xpE@&TU;;1cQhoHlMipXS zEb*oqDHJ_f@QcJ{R0(p2JJ=}Px;`fF`_~_b=g_@Z+IsO-3MYz5I`?@yOFYGNE9Nd+ z=#!%}o6?l;M$%W-#cQiyxBJLR2I1v@6Gcj7esV)7^j9$`WP7>eX85~;Y$>Kyf}yK$ zN<&{;>@3#Z#?le%Zfk+X0L6;GoYsiOmgE0IIiLV>m4HYCOaqF5Awf_$4DoNALk5Lf zx#>FRfUc5|bPGCkadZ8bq~o7r^ta>#@q6w;S)SF}(KV{^zvLeO z&JLuJu-{<%)ws}kDoLS~i?5FIcd3Mpm>wpLLToA|&DHP28=u;jTI}grw)E8g z;g|ed2~A#SpB*c4A$7Z8j<`(D_{1c(7YxSw#>h`vz0tT39Itt6>cR)YrONwFY7J^3 z9Wn>9-5f$me><^jFd#QmudlV!qj~17@&cDU$E&1^-^eVWUmif)lH@HM>aof>%DQ>46Sd~!u4!}?mK6wt=m z?~U15OWqSq2@mv~#B=@rl>1>cf)w53D6p)8=NyN&-{&d**&haa<7|XG;mO)zlD4YZ4g#;bKe}*-&!li@|Cy4P zJ!~j~7mhX-)5Cmd1!Qh97h1sxc=jY9~ubecwzDl34qb zyNjUCA0LKK1KnIM1ds-Tt}^AsfWxX$4?l=Wd!G8**rG)m_l~@Fdo+z-Ji8VZnd2~X zLzTFXLw~6BRz_^ zB#Frl&&TeE+L@7^;VF)dz9$7uOr~TfEF2bc2p8e>9e)snH_W>mh;aFjzL^hr*T+b`+C#fChl^L*quUvzBWz}XtAKC1#o*mLj#R-a8{^*I01n0;=EWC^(lh@?rV9+8BQ6Ln|0{u zpix1Vi_}Apn4DtWl8Br|cRz=^dcPVd;?l7*lojZEcJy8``iDN;#p2h9=}1 zuR_xZ?b;nrprRaQ>4}wDbq5L~!u)&Gy$>->G|t)=6l)#Rb+WgH-x^`K$dEEJK9)BYU>4h;OzdHZrN?p$}R}CJXmPzEOh~l1kIu5e=#X_`O zo`mwNprd8F*qc2=q;}z%LoBLp&u+Adx29Fdi^MIT>#Y0ab}pgahOc+_*kGZbJCgh6 zbO~>uN+Ek|n1$E}?%*W50~*GCNCJ)W(TK+eKB4{`WqK?ut!({|Ze{GggR&R=pxZ=B zdE5A!yK3n4ogL>O_6nuUV0sX@UrkQ^%g>_nxOq`G-V1NytI{fY>Ao~#^PVTgP+2Eb znGpfF4k5hw*qA`~jmM44s}b7|BfRh?7g2G>SM?uq_&+Hu98EdH+V>I57qckqYWZCoJiFKVtYnypy?6V{p~yuH!5%5>ih-!Uh=ZJ!7&JYH#B^ zr+DZYP1}FWJjfh>)kkCvRLb{swrygV=9frsW!?@xR{+8q_yhy#;umVI!>9x-oRaL6(Vwm8Ze2BSH6D%w)-BVlo;+A`;FO%?Y<9i|!`UeZ6zpjl7RV z3c=#8LzVW>;fg=kz?A>!WJZAMC_#vkTeS13S$(}5Ua0Km6?x8Y9focGrszvQP~}}p z&XWHAk-=^4NzCjw!O?c2_c(Ztbc!<0rj}!Zs%{RD9@Ln~4!F5Ye!1E!H_s5l^!S6w zP2rPI8RS|m3`K5V6&>bOT6RwOO&m9<1W*9*u7&ZxR0&knw9crW`oDm9_%Fo2CL8>J z!0=mz13<)BXNZ3X;uTeuDD;MMcvTEpR}cOnPQY0hZ#O`{P<67xQ-VQBacf&z+XA8) zMQ8y!?O3j`9t!wHyW%WeDYAjEDh?1xY~8(o9m04xJ3Cl9Qe=}AO$6B9lruJSFc5TI zg3t#UfdPqul9GqR;UE)}4NYFwC>~`je>YuHe<0qAPCoDBqG+nP}N>#R~pHCr6zN-L!2Cj(T(dr(^LY?viO$x{H8> z&v*dB7(yJ#EI_11R(YzXD<0=z0kowRUj4sv16#b45~WF!FyJTymVp3YunZU}28JF1 zgIE2|;&K1vvK|8Hc5z6+x(;<(p%>Sk%l0*js+sUeJp?%5RonSK!Pb7 zcJqCaFv`+wwZXTI1&N|a?VJ37_a8<7+G+zzk0oWlP+ln7c)UB1UHQe+DX%Ti%9FH`^l)>>y1K9IFd#Q1jbaoO JR7R^X{tuZ@mks~` literal 0 HcmV?d00001 diff --git a/plots/ccm_initialtris.pdf b/plots/ccm_initialtris.pdf new file mode 100644 index 0000000000000000000000000000000000000000..185b06292544271ab8e0448732ca90402f3e58a2 GIT binary patch literal 24819 zcmZsD1z26pt|+CrLvh_W#ogWA-8b&;TD*91clY8B#ogU0?(X_%&-u?i@4dUfy|ZSL z%w&>TD-8LPMJg{WO2a_Q0z+DVTKNpa1fU1l8d$(^a|7t4jcrVw%mBHf!|?LLI665P1Fd1)(&y5p%69wFg7)v#a0ho-^Nr`ERQ!7b zl9A}$e&n0>2lOA>7`*&Cl@)yRc9E!G&~4clr%ZuEfxDY)(c@{TRJN4T+Re)92ac1< zoUtM((bmdc1xPgqdObh@MpiMO?B(CqEy>iN#R ze#Gbd7{_d1Kg6ED(wX^kD)l(#DyQd5ycX5i$rB_X`f%j3{9EstY_H!QYTGNE0Zu{v zbP$n)J99xJeBG}=XVe`>fSZ?e-`WH6-X}RXE~=jpBb;!ML)OOhVurOnAOzsc#o^Zs zm06L;_`rS&=d)bsdV!gQKr8aM#XLkHS-CufSR(Cm*p!oB$Mj&Hw%%354YqH~P(WDZ zA3jh3ym9~6R3jVK%??WV>5WPUg@mJ?MXz--=bOffRBE!oSOlr+s4JA&?Jl_$CxkQO zW{Q!oxR891PRxkrSJ4B!*3(>d%<19O$_^85tq#MRdl})Ug(!yam-_eZS4!S3qm@~`=iXpS<&gMg>xd2nPhlpckn>p0I zFDvVo(_GBB?1*u-&V+*mbtlGDMIRPagUpOwrP&d(uP!iLSht?xYT~p@#aQSmAXbat zU>g5p&$7>YSTBDRx6^O7@ur8o3uj4hy3_p-@cY-S0>{5;n{Q})Z8H0tKR035-RJXN z&2BfF)Ctf_f$8^c*ZCic|I~*8;178HA*tInx@Oir9M72AzCK|d9vkl9rX5D2J^UtH zs*}=g*_c$&D<~isr3m&mxw)x>QEEK3fizukS{tp@xc*fr({LS-9siMv8 zh@aDSA-;=9L(UqeqN?H;1CRoP1PGzP%=aqbjJF`W9_bO+2e>m0Gh?NeSjWUkW71eOor7=~63?=1&qswH_;Dtt(rM1H3p z0y>B=Hbm9y>NHh~cdeu=>lfjZ4{Lg{3A2Am z?oKVf#G@7sOlP(v zhJm@An`X_;Mw+|&3sKl|vnSg$YeY~{kIIa-garso3qBx&hbz<3kPWX=N2>>n7}bdqop)1$vAq~Q4p z*+4{w{1j60FpUpCQS3k_^N>QhSp=6^>RBC&mQ>*w_anfJBMh(UA?}@urZflmw;rD? ztw$=#iRu{*XKxN?4<<$`>ah+-V47H%*kYO&HCp1`tr8DV$1jXT124N zo>;&CSag5vyYY6tHqc10-u|>VmXzSyU)A{@`4;T)baQ!b>-%!?+Nwv%$k+K=m#gP> zIrr2#0+O6-F(iaJ5Y;Mr$tX5IKTsta|SAHqc7Y`RXczm`<&Ktt}@A zKN^+Jtuo{YEDyK9-%0VxOoBUp99TZGS$1@LdkrvO*yuIAKdN_}$b0e3_if$ZSpI5P z(Ra@4`|;;xjXmn@T5O*@!pAlH;v~r?M6D_=LL(>j( zIx@6AG9RK(V{}<>U!T39S)Ev#A83XOvmCf@G~&Vu9JMi@>W1n)Jb#STo4}~5;_Oac zc8H@qEo+xIh3^lCiPzIyV)!SrtXMl#Z=J=Z4_Ud6j3)`@ zRE_n8JtqY1mJjuLv{uoX99$0%=uHJyl3|pcs-2yEtlILF*W-D-YJM!KC_awGx@f4+ z;QpUw6r8~dWw9=mbUaZs7KE#6tN7{Wtuf??lQB+=78@0)1nXLeUfL=96%eX?m*U?i zRVN_nrY>UFq3Bw6L+Yk)zA)(z#R-l_vHH;0%XwUodT2Da47uZ^q<&^# z&(5yeA0K!{C5Kk&(ndH>GMQ^DD0ZV1JL--y(I%gqz$bT&m!bH@P1b!qzF>vn%UnD6 z!1{+Pdy?Ru^k^e21WJb*e2h35Ij&Zo$B@ksS)y;TR?q}@Ik(4dTAzXonrKL7!~l~5b48gxfPeHQR{f=C0{IT6<|*4y$pUnugw~r_uJyxj;xp>uez?uD zy!{bfaFf{bW2S+Pz^I2~wTCMqAE$+LG^5+sgs~5|S#(EXA^n9_8oz%lMFZTVhCp%W zkB=Uv=al>FHCwtQ(N}K{IyTA%Av6luz|GN$i7Yw3g0C)_ec*EN{McO~c{sS=!zGCpG@>jdpGlv4u zii;9Sk->tBKj6o1vB(gKgQqG88>Rj?R1 z2z$Xqg~LMSs8O&O5=|iw@c$Y>qy_=5)qgUS6sgE0qPJ6+tzVQqC>NRzf~QJT!tT-xpdO-DI&5;Rb+er~{H`sUk+s%0q2APszCDOi-h`r(Kb@a#+3SE#2H zR)M4F&*vV5L!fGKpQ=mUJG`uZz|+JXY3ka!#_qz>tRCpBptx(srLoT#=%iQm-P@WL z9k_;>B*!Q}csDY9@Ah%j#Tfa?;!dkPUO^XnUItNJHq=LCwN}BSG*;l_P{a1&Mat5i zu4{Ono0L&JHMQ3!<<@$){p3DcqXdb;1Ky$szTFzKG*DKB};vsU_I#wT1MIU^5ts*$8+T{!9jKAslY`~UFl&T9-?UhF; z4u)EtQ8(bHv3p1HbIrRB6T!7wR-A3jU%m@zjS*HDTxdB;1F0q8j-H})9v{5E(-*y` zME}e5SPj22y>R33A7N`0uwY!K(VuH!PJQ$tpG~R->nq~gqVUiuEkG&@t1ueIDwp_9 zp>?9nAFY?bPNB4;Qq&1j&90Xmqyt|JWR>B}2c$=ZoKT&(coAD{NPgUM@GjG|sKr|& zCMU~s6&{MGHkF7q8j{VkeHJYr5Ths9k`4N$Buq=Jdrv0-hL)2PRiXfbRSkd$IE74F~Cju zXuBCrrMyV*t`(;0mx2qij=dJeCL9CFh|ms=mt*Q~H^wm@eRkU}SWb28h@Cp+A#JQa zf=jpvzJ|hP2oQ!ZFR%&b2k|!r(idXpN0&ngR73v;&rpyrFH3Yr2J6dCbo~q6KO73# z+!bkVwz?XPArVxV0FCL#!3;2#O0gee5N9!jAV8@uMzcZ)@j=-j{GeM73q`;%=Oh*( zf%S1(q`(Ns!&x-{g-bAeSS5mS9*@Pr@F|Hdi@Z~b)QhYNMz9S^)djE(3doQCr{k8; z0oLO}{n>De5XKayAbAo2<0&-|#+$dxh$-!icAmzAXI;e`0SrBA5K&cd8oZddUKva& zyqU(pxn*R}dIgJl(G;x!f1Ln_dtYxRtGv-_fvy_otCJ;BmOK%)y*hpU`fr1l0$o-D zF)P>!4Jrm#TJ-eVqv^N{#d~oFLlG!&=_p7V1qzQMk;!8y@hzl)S}DU>ItE8ZY0u{= zZ3Lc_F+4T(jC2=r-+nH+*Q}DC<5JLW!-or~>N+t6##e@KzNJ;;ntHF-Grk~j-k&jG zjBSkm9UAy({Rt}ki6s2Fe#8=(85kJ(g)yMbW#$o@DaR7A6?smoiI(eX}@n1v9109S%q%r(6j!wqd$Q&qW>ju#JK+yx3 z7})@<4D|ki@2>7$bivOh*$4CFq*8f+S z>0_b)tm}V^GyReCuYeQ-fPsabiS_@AJneeAd!ec{Tr@Y9ws@Q$dCSlp=6{*~B#15u zh=n1c4)%`ICjxE+AqmU>$Sz0z-;{bG`)+0FoUfx(Q%69}jfT8uj+z;b%wu7L`_KCvb7; zPCyGbs12DtoW^OpoZedN@&<*(0|)3W9Mk1?em-p3zl5#=$*)T0HSXToCRJkvbxMo0 zUMn&TFpl!Una}F+=)d8BMW^x5SqQv!t?n?pTb20o?PZ}8B2rcl0p zO~`Kh7Vqx6rt%Ix&(|l>@Yy&H?V^7vY_cdQ7M7AhmTg28=>ZnkOymvpQ==xScJ?Z8 z8J^jQ=@7KC*r}kd-RZkivy+|Y*La}D6L_4BJE7xXdMNZKd>?3siq1h$?I6HL*ht7! zsI6P=Bn(P^R|mSvB+cVisMHv9@w3hrH#-tFRAH~XJ5iOU%)-HZ@;Va?w!(gH!o0wp z)|z*FUOCGh_C3@&wSA&ZtaF2RV(WNo@_L14aqB>xX1#=$w#S5XVx~gIl()f^+e@4N zH4fZ3L?kW@#-2Hp843KH6eLTQQrGAQ4`Y6v@eqMEco+f&B zUWn#KA4A}$Wa1h4BSJ_Yt%3?W50yazF8Jfxm?B9!ZK>!4XS+4<+Mkz6j_-PED+KXG zyT)jL{5f$|Q%dej=IYh(Lpxq*Egl}8n5(r%to1hT;3D960pzmj`i4@U+>j5CLTq==+5>e< zg}RaF2ruz&J$dc3J2d%**`@fF+@g2!k2ObtBWjsdZAsiu#4QWuvs91y%tLe~K1oer zL@0WZ&1CixcMtWkngyI^*Na}f(aZb4p;SD*(!1#T4jny1QbRdpwV&a9(m>wWc}UWr z<}7j2i|lzL=68^E5CndRBWl>LtBDTGG&^DK9txy>`cwmyi(O2LS)5WlbV%zKG`2pi z#~4L7g$52;7$z0lmpH*uk*-~BE+gE0(^TZ*o52r5bv=*r6$ujeNJyT`I1OJbR1AF> z2=7EyGTp1yqCH#mI)9am`8@CY`uir~6J1Ahdw-{;n)_k?{Bx8yi}QS*9E*mW&&z9N zH65Oh%~+t#FSpr%y7<5ViY|yIYS+eyIk^Sfsq@*+A`A%o0JzTI9I%R@HI0F6`C4`G&kLa{kAoEWs5q^PVH@Yh(#$HY{S4rw)usW*1Q68W>SqsS5T`M2li$=X7$p zV~30w2^;S!uslKFLMPH_NQ2$Bh&9yYy9-=te?QdKO;hO3)ys#o*T~QvOWs8GNZvT) z3um46v^h#Bck%zGAp+k6^;{fHf%b0m27f)xv+Z1)RuJyzn>c;no!{rm?27a%_f8b7 z>Y|u&*=c7U^fTeGOiWY$Uen&KgN$p6g_hnYgHB#%s!0ZkQG9+?p*|}YAHv4(Is57W%|L+7=w3b%qkv+Dg{O&kvrY5 zeJU{h`6XX?K9OU9Yo9ZT!`(F9V;_$W^L!gGim`JD61bx&c#vpAZDpi4*f_Lv@TnRf~2+fCtGYm5`c zA#se0O@79g?jMl*J6lAMhdY=Dn3#{6?>d>;L&Q?RuNvKS@(ag%$LEZz+g?bigPcAI z9?dpEJ_)}G_aXNhD-xYZ3~K2PfEx6R3C0-#$f6p)McMc-w3D8%It|W*I@E-nl7n|z zaBsoH*twxI<~69+^g{w9Wuf0B7!{P#x!d?Q=8(3bLSDzP;O~x8r)gJ8ybT+VQ?DGX z)0{L2?W|+F6JwVvt-p3iBXdX`Jl@L75yj+o$H#{uNPmKkq)}agXYPpEG}>!h(jzc` z^_fqfLv+@DeAtenLX^{Hv3we?S1PJFpT$|UooEz0dR!S|7Z^F{~c%&ia}7c)s}$#_;TS^cHx$9 zV4;*=Np+FN(pd}^fJPA^Ad`b8vPC7RM}5J|O)+`LLCd$!#{1Y`Wv|iErq?)e52>!3 zFS&M_8@?o)t5P@S?gT zz<^A8G@9)m1Y`^+(9d5jrjAeBjse~Xs}!18E}e3luNiQO!Ixdl0@zllI7DGqB2dEl zP?1(ACxRuGE?H#Ky&wx^&vVLq{X@i_mDj&mu%xQ_TZs&<6%`jzWh(9biL7qtxAgAS zsNbB+l%AeN=f6joLvyoQE}fl9!9pxiw*)_wrZCMfX$o$1wj4Z^m&BcKEMA_bu%YaH zTG*Sbeoo>x^h~AI%>PuqqmeRR^EfD_^PG(K{fM&hDV6=H?R9$>9O7Vmi)LtV9KGJ){jS)R(j$1A zq@k+iEOx+LIRi)$8E2LxNt3$FohXY!5oATL@ zZ;v(L^>vKt5#)^*Lz%C}+RI&Za(f5MCCz!yqp1*sMBIWUTq@ZMVmZvnwlnZ-XHf8uUbfGk>t!+lqHp`MXnL<$0-0zyY66$j0%#w)3 zhS`YyMHEk-b1OdG##OJZO<;4I9-xivDKg# zDU>S1OMmaI=M@J6b&QiM5>>fUqk)}fujP7nnS~BMj(Pm1hj+i%(9!a0e>Ay+C-&~w zfe&33ZjIOXF|!LJBa8DBFDd8tvbQcu_wEWEbe)!#t3>>rTC1#ho~Gk(Eh{1G4vpOb zJt!Lv)=eS!?1c9V2NJ&)7M2${ra7lIOe-9DC9-xJFf-#auMg%R8ePaTy2en?rYD$~ zv2`cUXHFK*W{&ZmRG#R*no1Mcrdn#6Y)J>k_~_%8W&!0Pul0$ogra)YL8Qs>T6)p5 zWSt(}GA4wjvHG3p#)*V5;Ik)D1RF-(oFdM(F5c7Yy@7(k1rG$lW2W7T8i==f44cU zl>6MeHqu|=e%Or0v=Emtz91XRoa-Glx!am`IiBi$bCbe$u67U~^8h)9%$k>J!40CR zO#53rC-_pD108#f-L+jnirr2;wOFFe0nnO~?<7tRAyt@3BJ5mHSR1UnU+X%A(m#2@sK$K)|@N}AlTP#k%N zXV0LoRaG00(j6i}Ha(-&pl{v#xqQm%ok^=;U8!M4;`lpz>voRwX*9z4Eq$|204!x+ zQ()hr)b4UlYGkT)YI329qpYI@u|HXMN|E`jl?+S90txdnNPnyF&Qz9y{XO!Hm2#UZ zsR_(4+ug=XhK>Zji89ScDQ?)TfBp)waK&=iO~2ur=ZP%Vp59779Wj(z+%nI>)wjqWnI9ff zIH*XfBvOhMl!+EkMf|X!sgFXgJ2J1Opr2xj?xx8B0gwcy| zX;nEtxlNtS99m%!t%zERngBD!`Qz5l07yp*9u2%{D1VBSxJt$D+{y3=0ecxX1k;eVUiev|u_q*){n(mL zT4kG;T16`hi#GhVCSyvWO&udM{dK`zt-B9$jx}!#cnNLBiwjZFRm^5`V8o0Gjn~Du zhJ>M2QuNakIO;qkvftd$Y*49+alU2dw-&=lpXkd(KOS%YRRp%wOxTHPGST3rA{om- z0$U+uxO}R#J~vTG0=ahhfRwDFPYUFe0$cG;fkacQ3NtI7k~w57NWz|_YTakS)Oi9f zu8xMI7J|8LP%Lv%hr`zHZTb$P}0ZOhSV zP>Q+fM36;_p7PCyZjU5^xIfD{l1~|GDOvwiOCv5*0}4ULvi+a1bW1|a;hK!oMWgp8C*_bUbZVz}$ z^;74WvyJR*ssqy2e4?A6-Zy}*UhQ`cs~W8TJUD_it8iW4at^?f#ox*v6^79!TUEJ? ziA6I_G*2Wcc#C#XUI<2$h7-hsr9ha}f&~{Tr!4beX0C=hZ&Xnx-R*uzkQf3*J(fj) zGR|r3o`BKD`nqYwa>bITvH`YU8Y~4;J{X5AAu%URgT^VeQrE-E83SWPK2kJqaAXRH zwirGJw5(^TeNmHtune!(*Ol;#%-sfp)9R4evYK=nmj$6tPkz#nE!uDu!Q*Vy)XVla6SErCtcB=tp2bM_QEc|_I_Vo85pEsPfNgCwd z#Upl-Z0pCU$SW-Qq(Y88BQ#o#$g#G)8iK%}#BNAhx38nJ@;A={W`+#Qo@9LqB}RDh zdYaE;5`6~!q#QD2$rl}*QyK(XB-140_GCxqM68;Bqv*JIjv1B6P)(b z{$4F#|GlmeP{X=sGlJF{(&9=va$0nZQF_8m$y_=EXRm0cSk$uGzgnxY5$pq{H_FI` zN9mbg^Ayo!T-npwyJXhqH%zSy zh9comP|fu0{j^8X7ia%1L}S7P^B9B446{d69qeX>B-yo!k@HslL+qp&>IiR#!*M=t zr-w|8La&-r*9}`A;9)FY>zPKKVY#WB$=i;7Up1t}p=lbFMq|lkQp!7L^MNt#K4OI! zFZe8- z6a#552a^m8U!^e7L4syDGv7=JQ50UwuP#?h*jB#e#W@@+lFYhLF!(X-B`1oRZ;}Z2 zrI}#78t%1$XYh;rj^U#&?^@!zl4DsRk)iZc{0(X9Z7LcTox0w*cN~*hI(`yU73F<| z(#7;tyyTq|l|%rqHJ_hjcpDcsR_+$k^#9lR;LiA`Les zy&cd7$Ya=6I8}NQtH~PlratF5m3^(%m@qz%Y(ulpXW+5COo~u4H9o3uBgd@1ZlN{D z@8|6J);h{F%d_6T;yqo#*!=md`OMLVcXNHS@Jz^w`3OQEmkc>YJ~W*O^}8?y(5Z$a zIoBrPRn(D?P{(50q-cscOJ~Ph$LSN|rBRb{S=nS{rdMrpOh9tgqonNa4_rB395jLE zsuA018@GB^zTN7%x=o znHSgyv3)&pkGrn}`0#1zb4vuG3K7Yq;YJDUDlLQu;4%#dq=Kf>FRbOk5kR{95O@z~ z@`@xQQ}en~iu0gv@%_^J@#Y&wf{W<;!bx;`<0piUQ9kgkwL}^^KdxIb$UCjH?@dz7 z*k73sRfUIBWh~EiYKkeYryAK1Y0hJ=F9r*$qL~J6;F%Ifo!DDE)NtJmMk5-NGGCHc zx6)XACHrWAY9?>Lc7>kk#u0=eSGmSma!rKo7n4lpOi<(LO0XJl8`NSMo@RU5x&7mW99zd;rZ+y?AmX!BpyQkTpmn2HMQIZCY@ zqv$XXY%vX6$9VgBbdnQXs&)lN&~0k%!)+nmqt58xr_8`{+(^=PeosmJuBYqdTKTBD zgc^Z%dsg`H?B8Bve>G@9BNAAZmHUMSZGd`mCs!(Qw!PUvNWmUZ<(b&XUA}8uR6@ltz2`rm2r{l+<&J|-eS$(9(+>2mD z`1+h9_ACD@|Cfc(3FaN_0f)wGSorXWGI$F-v{_72oRDeCO2PQK0C?|+Wf>=?`=*Dc zkEF%i;u`t5sbq5bG=O~ly8hsoS#{-Ov1C#D&>V|;(XtsVIMF;X3;2p!X3gSVz}hE> z=?V>~((U}T@;S>o0gX<$_O<{gILuH4PqXv`qsUPHNoGWNyYwh`8j9JRhOI#Mwd#e{ zvSl=9VhOs1#1?uafldK~jd1P-Dh!LRKHfxewpIL+iqj3R@1@B$Mk=Gx0xQSaa@yd| zTX2)TDfm1d%5kO1_|G*7-^PDDEa2xQ))&g+@h%fkGVBG%$A!H1k@mkB6f{*aH5AHX z}4yga$ob%=8)?bWUrQ^Hb z-ssfq&FaK9YH;QsfIH(pQh@fl;Bt?P5;c!1!>Gs&zb$3-Olci%pXwZ=L)Cc(y#J zn-z8o7z?nkQZN>}2aAw3{A46K{q;+XbWO?ltl5?f>?oF4Z`>9K1L4(zPAv1g%pIdS zwT97m>G3uzsVib3nr8|d*`T&HCG^#>Nm`joFLtvmjTV z6^4LPPk2I|?@blkei-xbt4S-)TfI}RyqwpvigyJ3T#oqOmKbe5Al#!xf?8q$s$w~z zLD(784c)~FM#U8fJk^{|l#{6pH$k9cg(D?_n6lGHB;$nmyToX)4BScFPMM}@)}b4 zMy?(%nW?VWk9WG`KF8)YBzpB{aALNW6d=%N5m6i#QKNlfhnURdIrBQfXr({d#gbzQ zkXHOPIYW)2T6ory8X;gLwqRu16szU96(=krZLGygs=@hN-m;=~MZ?6;(Vyc#=~QW$ zsaZD+TL(I(CyjJXHB5CZSMt|#xTx_M3F^G7Uiq@!X#`NkC~h`(`1Y4HK{$?W&0ji$u1;#m*-r2#k}n`8zzhrXliI=#fDlK!>w|* zX~g7pdL^=6$zl-_bZL%jlKnX~5RR^Igl<>qy#?iNzgIVwL1RdbXts5cD(NW9f}GwC zMq-N78*j$a>eI@@OKmTl7rbtLEjaKQ4CSZd(pDx!_l5QeMpG`@zZtybFH%jJ2mHIYt3 z+*HTg)ivuVCSqUj=pmgOgq$B%=+2x2Nl&$Jh@L~RNfEBKE)cKSQT>w&V)eefIDUog z=T8lX@s7^Q=_g+h!UHuZbHY~=chp=aZIlow!tS#$U9@~8T+ez#pCE|?*1p)ty%u&= zq_^*9(^8%sV(t9sY2>njmD@Gp+H)%mI=`g8t-9tu*w|V;c8}s(|>sq*p+nF2)%RjbI0BcmP-u^mxMH%n8)bamrH5T z*fgK!(6oOv=Vp2n~ZlK`5BjO-+WV`_roQf^)u3eTD0BQmHIL;sVbz+`5yA z>zhi3+V4*j``_ydwnw6ToxEnM@8^iMf5$yX^VHM2qexC#n(_%xeU=N+KNA}p_d|B2 z5D<-pHbQQW<`%z_~9YSW+Cq{_Jw_J7R;Z-_z8M|e&ue5zQ~b7+6}?SSgoR;mf& z7Rl`Q)H?FPId>gueIH@ZLRqQ4BpR2GJnIa8b%RhoFnxXAeL=LYe0SFE)b#-S?2pE= z9Na8nE^>}DDvk%jcJRxO`Rbd<;pH1a*7N&Z&|3yh7jYYuO&`AC7EYk8#ySJ>CC#xZ z$18r!_E(U_PpyM+i{N6LG7ivUr@@~mLC?2*RDUV`1oQb^Wl9@AgDz8mXD4uFlb0O|VN230B)u(`D>6@YE3;6D15W7m8VM`D3D#m0_Z(Pnpk5w_;SqZ{2` zi`{cqBn&i{OydAB*KJ%kOgM6Js>rRGqF-#iB8hh1mtQBy3_=aR?*?clBWP^`z7oBf zix*44$)@vH9>IlQl%(Zg-EHHAEc~XPf=xo1Af+YT*zM-mmVI88D7^BiI3ey(i%n#6N?`vu~i=8mt=PJM?y=+g{&$CVc{6R0W1$ z)&Hh>fW7=tX+(h~#5L7?V{XGMsX-fc3K0i{nD+d_`|5~7`dpHkxs$o8w^vm<%{fjU z(L{{whT;Qji}bmYSsZ)tbqDS-_O4+2FzyESjO|7CO~@X!@(CgNGh)}sg`5E)&!C+H zSBh-nkGUOSO-lN&ym#px za@8)ooseDj=YdUg?WIrb&;;Aw2ps^>DP-bx0TiOVp?0?oF?^zxA%vQ^x=k(F7`j2% z9fBP(hlE_k9Ll2t!qh0&z>j?*kx{(;n}0j^XoO%z`Nye-kp zD*+Hm`@1g}nW4|G><+6o>@qjvT!=%2?L&b4Q8&i8**u^2ta5(q$2P`$wIN{sbk?9q zdad76JcnVv)bnJH+lD~Vq(@5(ub|ce))N!C1Lt_}QZsHPpW2AbZ*;rx=7DC+`%m{+ z_u#HU&*0}^RlA?PI0=5|dHMQq`(nH^zr(*r$`iNIc@<2T?iH(+Y?XF~ul22O(c3qv zQ`sQhhZB?>|KN)5=5_nXfF|E>-qj{;M4X=iU6=8d^NoJt#wT+O`07BpCNMu)kRN0M z6|UTFM#l(7dx@NXyW&LEF6e0=Ac2Mj&1BxVGo#GDoQcW(gn^zM$ zyS-F5#EUzyEzkpqIz;}yh4L%s)zZ)4Ov(+!WQtD+XZRA^&0l)2sLO)_Ahn6(_NgE@ z&aQ}3%$WRq3d`$DeJpzRi>$T|G=xIBegIm2;!19dIA3t#Hg7d9I;`O%;t=!ZAm-rCM4pRQ_wF^}Xi|=9i9kjJLBTWN$j`xszM<)d zBLtydgeZJ9__<6mx!X=Q7R7dB4*D}sU9spCyY!a`uzCcLcb@?$%p|u@tmK7V8KGj? zOh*B8p`;JJco%3^-|&o#2ioZot*-_sRssTjZt^GTrczj^q7TKN9>$)-XU9osAMnrB zatLnF5d^-xlVd^~hdPJy5xJV?1+T$<`K{XYt?cE~2K_h7Apa;rC%Ie#HpB%_gkW-A zX|V6;h2$nZ^6YRE;h@q;c@FN)+Vjte)4DuT9;g^-=sT_>$IAtdHq1pzE~+F<#Oe4j z5HIGuRRgCEoF}YQPoiXb0#!2>v$>L+9Fl-W#?>lX%<$R12XVJw-^R=FGrz%qu`w$S z2d^&K7-O&2&+tcY!%`9*EqRAmk~J7!U7>v4+RNEr87h2&!C7=9&xM1-d<8u2};l?R&2xp*evKu&54eS|ERdo}}IT~+fn08fk zGl1*%HYhCi2lY-^=vDSM%)Lt8gloWdqi-T_tT)%u_78Pdc-^_!fo3r@Tn)6iLvTdh z&_Z@xckH^gL?f^{OO1zBM0u|o#FFG*SNAB*>qa#0WWp3emQr+f` zH^Xi-<7SeBGYDOE(dVBA9%bspyyPIQY(j1`O7}{Z+dnKC@j!iP9|O2Xst%ap9hdmmne=dyHf8Ll%U)?NYXJHk2xEuVwFW-raTlWCVMz&D-)|x(vU3Dko+3N&RVk>x4I6UXZaCVZuf{3Q0q||BmQ%^}} zCeSU{-hiDVC=Fqf?J!N|n)0}kL6+P8Lcq3sqs1bB&S|=~wx;L?(f9fCsCH~yrhu=s z-J(y>>Yb4)WcS&|+~=3_zZ380*MsxzcJ1FB`#NVDkJ&=Mx zTDZZinCd_<`TFSXA-@G$vdBUFqz3&7#5E^?-TYNNA3EXmj#YF`2WN!LcS@|i=MI-w zO-EAu&ZbbjIgk%WjPNN#5xT;cz^hevCF&BC1~DMiv0Ee%3>5N{bSn{@2xj!pPe3Yv z$u2`0N9068T!E(io5A$j#~4tKsHFhO2G_k4uL4d&jfxVLQcv45*}Z@e^6M}^iMmXQ zI-8)OJzDQrcU2q4Mx#Bp8g0u#TN+n81oj|)Abw#W{^u`*2fwVgPmqMN(<$eO9248I zVqPl5*zuvnIQ&0}0j7pyOaeGMwWYq?5Ct!PdircFbvN12F6#I;f2DPZ@e}7a6sTC| zuV=a&x}Oo5Viu|NG(wAEc3{6rx<$ayc%tt)+W3mHi04fw(Q@E`+})~y#84Q>!cSO* zD6jFSk+C%nSfB-229#5Bb{bjw7P+R?PKlX)B^frfv_pK^sJS{yLHs7An>U3p=wd$U zK8r6(Mf~YL&1{(PR#X#t1Ojk=TO+*PA--S&fj}r{sB0u;b{aX>hn>qYbz=#mS7LPo z38Owv>S_{)35O5J=t-EUj*`OtMT$B!0wjN>BbUE!*466R{ZCx@VO}E zTBLFdPKguS1?J0dc)m+38#&hf9Y9~h;FK(fUw=UUCA~SW0Cr@6+N^?o8^ zwwV2RA{Uc~*_hws(Z9#ZjKv7~IdHjyzHmF>ef5JH>TQE$ zrnn`#DACY`;(cP=dmQUvU>}{u9fHRt5oNXT-6Y=#`sLce-nZ4*A>_vgzRhip@Pe$R zyH4;h*vNXJyU>KKg&?*CRvRDye*3LTOn<>1+?m70${b;Jpyt)DpySO?gSWs?vSvF9 zqI{_E*Z0AuUm#w(LAD0`Hor{I5aS=upyJ7@!xocY{kYQwnKeechFq?`(!oC`Kt|C< z#T`0SJ~Xkb;)dn8btF>-_kK$z^AJf}&P|*Zt=FwkDVw8YsVcOPUs>!V6;7$HFCj@R zaZ9|C)X1fl@sR?ToJHNPPfjSyRUv-+8b^%*rS(I|~f}|h_$Zts3!y#?d6S|M`&MZy1 zup`>WiR*i@K4|ZRB~k@!FW-lql{jeb+~x+WtCOZn4Ij&Yh?iA{@eCV;eEvA`D;oEG zx^hiaiRSrrOZHxj+qK-hjy5-^V$Uer6##hT+9Diy(dklpCLFu~6yK3wHjrF0+H-Fi z|Af0_#vW9@U}V7tk0yt+l@k&*?#qmT10TQ&ga9wS^ThwY0`6Zx>C$3cKlSaPHR{kM%7CosdUZx^)8fqof z*`dlcL5rojYn5FeesKvdi;G+&cR_A*(rOtacTTQ{?K?22kFDME^wt#(;TaJTrypN+ zaDA@xfzJt}_U!FMr|<9gb@y$!mH+o8(D1%EM~o)eeuiOMVN96>G>dzf+eE zysKZCgQ&<=)h9YDO*;%?H^xZT7debO6@zr`sh=R(l!;6=l&Zn@efSE!^PCO&>JygM zoI0c}_e1@(i;HnA4y(?>^Qyzk%&Dkbo#*d5CE6S+R{0@!C`)qa{6Z`k*J_cu#W11; zxyb@`S!dou!hIPxIo0ZC6MMA(!xH(`#n6-%qZea(qOP_pm zg6aR46#lp$lPA2^dICHz{hh-!XB%W)E&c=}0^9&cI@NdnRR8K~_%jR4IT(K&%$$bLJ~?n*_s zzFiKf-t4?CtiESC`)k6^K2yaXr)Xh?_WhO_LY3t5Ck>T!dn1#jmabQRk8vc<22aW( z9Abo{Zkx#GsYloi?+geS_yr^%ebCOUSZ6#o{kXLeOx*!!x0t@?^?b)hn`e^ek%v`K zyO?{l`h3&RQkX2gQ+nU(jh8Oxm!tgH-p6}mNKDH*2L_k1l|~q8%F%vs(0rw(p#6cf z1)+f=wIjAAwiDJQSQ%=$^ZAmnnxsOwg2*@MEw&D&{?9M$h_ak28kLVo%OG5uN9sf) zrc=;Y!r8d8>^oWw*?h8GJ?gy0_q5UAMSXE2yN z*de^=?jMgVBh2d_6?Rj$m@`72TVWmfRo;9)9gx+~|NpfZ4F6PJeC!Sr3)6qL@xNaH z@vnM-KP3&q#*T&#=5|iD4u8`7|5iE30!ntl`?F#d-}@L!b*G_0)502(GnW&k5SBRhbBo{8<>6$=thKr3@Y0UJ{* zV*vf%VFHeZe<~h0=vn@D|E%gypN5fv;V(i6XeVxLZff>X6!9_5UzC!Qv9&6I{a@q% zbj1JZ{L$;8PMS$<)wlD?lUm{YuW#|g#irzQ=uVZ^YM}ma~o3totn9gfQ_U1KZw6# z)y$2Y%>GX^*Bwvg`~D*#tBiz%h^)hL9OIC^viB;Kz2|W3;sRZ1C>5;a)a=ol1h%= zxqzINfBFMbRKD*$*#ek`2Z9kZ^M&NTWpzvM?2&R-u!s9}AaUjl_d z!M~-KENJRl=}27=gl{!{@yy)hGv}6{o^q9c#b1tiR*xXMX1}!YNg99llk~_i#>`M| z$4WAqO*6d9kUSmev6pZE56&Z1ST7jqHLy7BZpMz$?Yn&`_m zr^S38g|Yko=N0Eae{GMkHEEB7nVPVlrL{qxNYPv)|r(CdCi>QxtR(Jd8+_9~+i9{>a)Z;16%eB*)0ddf~6;`X* zE=letNbc6%Mu)pG=KC)N@~n8|?0cMB8W`K8==I!(P$hR9+&;+psJK!k7p06vQcf%( zz?L{$wLWe5@|?z^@aJj=)2~;*#71;G4uYIj6b6%mo6lYra;@ztS8*+$)->8E$q90t zZ8|HXRB9D36O+N97;YHEMq^3&x!pXuaGeEa%?JZ+*x5?XM-TP9D|BWb5Gv$kl6{u04LNcMfc_;H)g z{zq*r{biiUvjK%gUflcUO(CR}eUEOG&0tbq&+W}OfV+HpY@R&8tkCnR6n8ZNM5O^C zIB3$ZbjqQ6^#Dz|*vAYOOTdbKrM7ZTgCta$cWO=n8F^+Xz}u>(z6#_ncVMF#^;Y1{ zfVz+~YeLV(SJ~Fss^@xp?K4(2}OMMmJnYS}1@ zJe(wEHZv(z53}hId`PEQ8HaKMkYR>V=V?g z`6DO^)|gVk)=})2@7A0(45sL>k$$Sq7v4i25h$7wjP&y!V8+*WQ}`x}xk{f24-z7kTWRY1iy8R(Cm^a-VdP#xl%={V}+%u6#r&husNPmg773_-; z@qT@$0>SXoLQ0#YC$O(6^u|jp$p-(56vr71kH8|j`9+<+RZ|!yQKXiq0A?>(l#?VD zfi6Ke8O_xmzGAa*R;&>>ODZC`)CKlfHbfx*68-KgidAQdPoqU2Ev!>}k%K$!T5(nP zzEJVz+!|l2ymvsPS?rMb$|xa~+gH`84~eDy-mCVKO?04T+wPL_3nVOQeR00bi%QKE ze5;{8j&*rcgy0SRi`PlPP>2;*S)2~3ia-c$HvSeR0WM`{NNB+LLEgOK)>j?!H+GJ> z!2jvz?3%&B8XuRIcjR-Z?(gtX$Wk>t2`?I}?+N0YZ~v019O&cM*JRflkYMu-Z@Hwa zoQA-@8^=F~vQD6gBC{^wgEBR`n;z)A02UaHhRjU^H4dZ@`+p4c*MBpVo z1&%goN4uZes&Drgh8N7~S2G>MiN1Vh3Tzis_tucWi@4SPlmf)B=3-3Uh-|C=vhfrp zhI1GU&ut3K$l+JRl%8nqUD_Hjj0L@%qMsB4PE4BjQ&D!!su66Hc_&CrYvM^SF6S;i zD}X_yEaZk~jKPIf)*VQIf7 z*^^_Rg>Ga}IqLhvBTus8C0?U7%=-jAgy50hBBE6m6IeOlu<)!Agg1ycX!JQPj3T(xtn(V1e9C>c zQV}Cp=#1%ti9v$6mI)SIF`ke+sQZ@h9J<>pme1$qzF$NwO{2VOBW+d%}(; zF-MyS6&v;tey%ZxW;K5B+M|_ysA!E1ty;NLfT5uMQZ0-%jS!PE-n0qx@rPOEo?D8{ zkc3f>cw&Y4;s?@a_xO3H#aZUt8SqL>Pra>vM*|b8UEJ!jRs2RYPOUIQjNSX@PUn<8 zHM7U_vkx@4J7Wf&%ms|*)3>$hL{-0xn@D4M!PO|FzY@lfq?|FgMs_|c)_YCp)#D@~ z$DNr=LW{+qm@x$-SxjamZqZEaT^Ehe?QvJ*Z}4B+*3rg$skE_oI?BZU(O-gz3ea~p z3gJ9tIjCP8eGs<*CV$`M7mF3h~dZvyoLt`J3Q%6UIRHIIzM?0s+mKdf5U4IVeqfG_}4M^ zza%#QDMtTDaJc`@Z&0`AC(XI*`)Kj2i}7AXM8JnU44CK zsC)$7COWY`Eqzvjq%wvWqI}vu4#cg$&dh?r> z#yp{4W?;NZ0GqjVjul-?-xquGz>Ne_m756#xi(+$Nc-tF=P52NDb_np8{cYFNW=7N z43?vm-MHbNWEMm%z#R%R`TV(hrGH@aia>1UAD6X?P(Q41ylXZ7^`s2MaM2*dE zQrj@k{#@8;?Js#-^{PVEjZ5WYL)*Gn#?Fnz_mx2GzPt!r385Ce?=d#&6@;4`4i2b? zy_jI(a*|qa?NA<@cL@p3m`$Q)&1r;%=pw{y>ZE$&R&K@+T7?}^as>Mjc` z7+US*`aOVRb>A!!wimOu@xZ4iNvB6s=ckzMQFxjwi@a(X9XD2rUPRtNaJQ_8pTk~y z1rK_8Tkfn0EmpWmUY=?R{k5om8OL&$>noSrw2i3PO8mK%(tYz+8N{`T%-gSKL*oY5 zpUmgxw?m_5+#1KBk2A%0l+6aH=OD_$mhqSGHY3vzjUC1~-PWQ=X5~sg)g-FZDJtEQ zRkb0f&6+Z8dP$QIVk}DR8MA>ZDWMN2defu(-?XtVdgKy;8&heeEflJ_AJ7^0t%(lP1ZT-e2utW{qJmzL2I zGBpcr0lJ|wzU)%PZB2sF0-<6?4p&uEmF*f9v|Ze0LY2qE5>rRsb$tjQ-<~gF z9BAldPtm&$S5fr!Gr|d@lpM5kG!;@V>S;#5%+a8$dWK^6wwjt2{(5R34o8H#g!)wY zBa#GjXiM4S)cNGwsh>wnf~8|9>UYM?u83uM>`9PVIIYPBrD#|@e^AOmg`Z;pCkW72 zC?GlG7M(d=iIDAQuCxj2yXWqANrxy%MaoJbFm+PUMFHIg(#B>}Dxw>{pz($5&7cJZ z?+E7F1gCn5cc6HO!4gy>T0-aid**0cKMBXzPaEu9!*+RbnD)DKqr>Q0Bo#(pcUkF& z^wvgFk5e*gjnU;%&M>8!Kp?6_ev;YRx>F6##d?ob8qo?N%%qhmv9;+H!X`wmHE|eEKQ+NODmCBi&MQAZcBpXV^bPcMup4{y(aGwnR^D&QC_^xs<)gRV*prGhj;+I%V*nSP@I^VsnIQL$v!GJY8RGRjEy4Hr zwU(=|J87i`0YN9>$gAYk^O-P9`U&;6v4l^pAIR$_Z)Bp~=CBm4zyo8^kh?YHHoseC z7B>3I45z&RF}g0PmJ)NVRlneET=*%|n7TUvL~jT1`Ipuh%xwxdnRBcD2mtO~rb#6k?8Cns9o6_&6Axh}fdsjNluwXbWqM3n0PMmRcuxL?#?CuEnbbhpr z!L3&By4Bw61RWsu!Y`bo*p^P=Iiigjwmr)fSU7rJ-!}5o=39os&UuQzbh0yV!(O2V zrF7!`2|o`3!ZG{=1B=qf1&#iqg&eMR{xg_3f#oggzyOx_z_~miX5Y0yhMrnJQGj3_ zIL-Sij)HW6vHgmzPy?{ya>r}^An_~6eO9l0gmdOxaGxsXw3l)0ch_onwvwz^%=Juv_|ydgc5VVkwNptitRA_OjV`;-E~R=%&-yXokTdtZ zxzpSXw`Hh&_DJ~RcDMF>cRj_2-(Xd};&yEA?jhdqyW;SOo`IiOFlUoeTA&0q`3tKgbdU&KNJe3(bYw*zkeSG$0QHc zmw&okl>JEK`B&HZlXo0<9dzTchey=QNFR`ffasyLmBVPO9M#bJDLFt;Q~XzI_3M;B z>LjX8*Aa=7H8(=FbVz%?b3id$1gdn`FINdUUKlv49>y*P6jcZGHX24vTPqU>eQTOu z=qDg=u|e3tY+xD_z(da52H{NO%*#!~$PEP=Y4|wNiW6p z0S(XA^f0uC@%Ub??B_ww4#p-pD4hof{NqK#&BevZMPov9WP@<>@&Z&0^`f!)ZUaLh zKy|YdHZX9pCv8v&7?>`<9Sd9}5AekK%?5aJ1Gn&Z8z;|S+T%hgr@#4up**Nd|IG%3 zk@IAGfDNUEo$%x31;)h*8w{mr{chvoMh&^&Z7^Qc;5=!Ao@@^_22REbu>BEJ2$%;7 z%*fx`gFtwJ>G+3@`;T*Q{w1cI9I%u31HuXA`O^l2{$T?~A@JNg?u(;?KA_Y(9M+3d zHFq^ejTh7dMa9"data/graphdata_initialtris.m"]; +gsraw=Import[NotebookDirectory[]<>"data/graphdata_initialtris3.m"]; (* gsraw=SortBy[gsraw,{#[[1,1]]&,#[[1,2]]&}]; (* Sort by n and then by tau. The {} forces a *stable* sort because otherwise Mathematica sorts also on triangle count and other things. *) *) @@ -21,7 +21,7 @@ gdata=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}]; 4: {GCM2 starting triangle counts summed, GCM2 number of successes} *) nlabels=Map["n = "<>ToString[#]&,gdata[[1,All,1,1,1]]]; -taulabels=Map["tau = "<>ToString[#]&,gdata[[All,1,1,1,2]]]; +taulabels=Map["\[Tau] = "<>ToString[#]&,gdata[[All,1,1,1,2]]]; (* ::Section:: *) @@ -33,9 +33,9 @@ taulabels=Map["tau = "<>ToString[#]&,gdata[[All,1,1,1,2]]]; (* Consider all runs *) -getIt[x_,avg_]:=If[x[[2]]>=5, x[[1]]/x[[2]]/avg,-0.5] +getIt[x_,avg_]:=If[x[[2]]>=5, x[[1]]/x[[2]]/avg,0] getAverage[run_]:=Module[{avg}, - avg=20*run[[2]]; + avg=run[[2]]; If[avg>0, { getIt[run[[3]],avg], getIt[run[[4]],avg] } , {3,3}] @@ -56,11 +56,24 @@ TableForm[histogramsTotal,TableHeadings->{taulabels,nlabels}] (*Exporting plots*) -tauIndices={2,5,8}; -nIndices={2}; -makeHistogram[datasets_,n_,tau_]:=Histogram[datasets,{0.05},PlotRange->{{-0.5,1.5},Automatic},ImageSize->300,AxesOrigin->{0,0}]; -histogramsTotal=Map[makeHistogram[#,1000,tau]&,totalStats[[tauIndices,nIndices]],{2}]; -TableForm[histogramsTotal,TableHeadings->{taulabels[[tauIndices]],nlabels[[nIndices]]}] +makeHistogram[datasets_,n_,tau_]:=Histogram[datasets,{0.05},"Probability", +PlotRange->{{0,1},{0,0.5}}, +ImageSize->300,AxesOrigin->{0,0}, +Frame->True, +FrameLabel->{"fraction of average #triangles at CCM start","frequency"}, +(*LabelingFunction->(Placed[NumberForm[#,{2,3}],Above]&),*) +(*LabelingFunction\[Rule](Placed[If[#2[[2]]\[Equal]1,NumberForm[#1,{2,3}],""],Above]&),*) +PlotLabel->n<>", "<>tau]; +histogramsTotal=MapIndexed[makeHistogram[#1,nlabels[[#2[[2]]]],taulabels[[#2[[1]]]]]&,totalStats,{2}]; + + +tauIndices={1,5,9}; +nIndices={3}; +(* TableForm[histogramsTotal[[tauIndices,nIndices]],TableHeadings->{taulabels[[tauIndices]],nlabels[[nIndices]]}] *) +plotgrid1=GraphicsGrid[histogramsTotal[[tauIndices,nIndices]]] + + +Export[NotebookDirectory[]<>"plots/ccm_initialtris.pdf",plotgrid1] (* ::Subsection:: *) @@ -68,16 +81,33 @@ TableForm[histogramsTotal,TableHeadings->{taulabels[[tauIndices]],nlabels[[nIndi (* gdata[[ tau index, n index, run index , {ntau, #tris, ds, greedyTriangles} ]] *) -successrates=Map[{#[[3,2]],#[[4,2]]}&,gdata,{3}]; +successrates=Map[{#[[3,2]]/100,#[[4,2]]/100}&,gdata,{3}]; successrates=Map[Transpose,successrates,{2}]; successratesDelta=Map[#[[3,2]]-#[[4,2]]&,gdata,{3}]; -rateHistograms=Map[Histogram[#,{10},PlotRange->{{0,100},Automatic}]&,successrates,{2}]; -TableForm[rateHistograms,TableHeadings->{taulabels,nlabels}] +rateHistograms=Map[Histogram[#,{0.1},"Probability",PlotRange->{{0,1},Automatic}]&,successrates,{2}]; +rateDeltaHistograms=Map[Histogram[#,{10},"Probability",PlotRange->{{-100,100},Automatic}]&,successratesDelta,{2}]; + -rateHistograms=Map[Histogram[#,{10},PlotRange->{{-100,100},Automatic}]&,successratesDelta,{2}]; TableForm[rateHistograms,TableHeadings->{taulabels,nlabels}] +TableForm[rateDeltaHistograms,TableHeadings->{taulabels,nlabels}] (*TableForm[Transpose[rateHistograms],TableHeadings->{nlabels,taulabels}]*) +(* For export *) +makeHistogram2[datasets_,label_]:=Histogram[datasets,{0.1},"Probability", +PlotRange->{{0,1},{0,1}}, +ImageSize->300,AxesOrigin->{0,0}, +Frame->True, +FrameLabel->{"successrate of CCM construction","frequency"}, +(*LabelingFunction->(Placed[NumberForm[#,{2,3}],Above]&),*) +(*LabelingFunction\[Rule](Placed[If[#2[[2]]\[Equal]1,NumberForm[#1,{2,3}],""],Above]&),*) +PlotLabel->label]; + +plot1=makeHistogram2[successrates[[1,3]],"n = 1000, \[Tau] = 2.1"]; +plot2=makeHistogram2[successrates[[5,3]],"n = 1000, \[Tau] = 2.5"]; +plot3=makeHistogram2[successrates[[9,3]],"n = 1000, \[Tau] = 2.9"]; +columnplot1=GraphicsColumn[{plot1,plot2,plot3}] + +Export[NotebookDirectory[]<>"plots/ccm_construction_successrate.pdf",columnplot1]