From 5e8f0e15e05e6a4e06a3aac3398b9ce5e4f56602 2017-07-13 16:28:06 From: Tom Bannink Date: 2017-07-13 16:28:06 Subject: [PATCH] Update exponent plot --- diff --git a/plots/triangle_exponent.pdf b/plots/triangle_exponent.pdf index ebdb0f4fbb796b478a361cdf2ef07cb0f725df4a..880b9b773fdef1f3dc2ff8a7fa79e11a3f6ab346 100644 GIT binary patch delta 7016 zcmZu!WmuKnwgm(vMY>^wl;mzeIyQ}TNOw0#D6plwq+v@pNN$4?NkO_xQW}vi39tH{ zd+#~lci!iDpEch(#~O1zYmOh|?T^RoS4dPv8Y(9+d?NlZT(ec^HRIjNESV_t0`MY; z0fpBqbzp8@gF~D3w2CVsH#dP-m-4xd_nY}UYFLCq z9P12+=}r#SftUTcX|CS~k5-uPHqULPA(_7%?X~ktFfSrWQl12gUQBkc4=cOTFi`DY zBMqDk;+5+w8{F%tehPSBqSu-%QpJjJj|N7gce=*iBBZ zSmTiFXyOuLtqkD=aMQ&gznStYv)ZmKB&_G-V)=~O955{N_~hnf24>)q?#`Rf5^Ws! zluomDTC1)dZB78$+~!=j0vW>{*NyKM9I25fO=oDP!wTXuBZVVSKQ*!TkKJhrd{yv+ zA3EacT9P%VbHOv-v03scc;_CHJ5#)c zCgGc7ygT1R$5?rydtp7jwNazTCod$magSbCM_RwBCMelU?GW;|(y)sq%BB%Z=*oJ3 z;f|Y|{Z?-{z6zMe_=JM)K_`dYpyiY%(QLG&Jdy7zH$AF-IFC!#1`A)Kc5Z-lJ(9U< zFI)ccy?B;3HMV3(2erASH62R=@*(Cx^3-ys{FaepQQ_ODD9-B1G5hyjvY|Km@?&Y% zCEPSY^1(hzM-fBR&N}XK?c}CjIN$0!9tNqYo7b)ZDd|<{OMsp^F2Iar<3N`e+jwob zhEQzI`=u5k_MH?qLy))NVDoL+h>o~2u6TTZi^!}aZ71O40!OnAQ4bY+z>O4TZ_DJe@Mhh;4;!i&Zrb8YlI`X4~ z1~k=!{zDw`R7$qP0}jVfqOAL3F+xuI!n8g~nypoA979gzq=^AdEmCMHf&!au_YwVl zp*wdfG}mTPf0% zSXFK8;^=@%Sr{yZEX+Zr07VkLDjJTnfik2yWT9>d-L`>XB4r9CGU|d zO*C_*e5y6c>@}KA8 zhy_j+ma3&1g|X?IHa(KLt`ZGegLCC_JXXf-v6$ODKvK+m(%cG<3z84#%BA&L?%`I2k2O^fSH#g0> zBuM}BRF3^Lv9ZS?b=DkdJ6vWOsf~#Z3$6!z^lx6_+G3lWEgD$l(@0C#&pYq$Wp%%r z3TBY$!+)-5KbqkcR`+gO&xg*zS@yAe_!eFAdP)+vcQ*1_e~DjB`(ojl3X7*O9mB@^ z;BuGQ#(kI&*=!~mgwpBla#fA+3cU&QH<(Z%72YV$2iwR+IjtJtMYs>mb7EZLyegGN z&9CYnWz%Z~Ka*opY%7VaM5&{#JZkk2EYI7`1c?Gi0|tp*svvDe5r-q2OQsFOt5r*nZJKR3@y2d9X9 z4#J~Zu4$zu5jtRFBw2Bn5sPvCB#EJ^uj*CyQrLGWjSiQEL{p4!xZ=0cc@r;wIJC9DYoT1U2Oc#ahiBPgt4j4b38WYC3;5lq%3&mRW&Y+(0{0!ada&Qu*;|4RX{AiT&vFeYWM+k>2i50S%X;? zCMUF%e{C!-!HYCxA_;f)->7_pqoH*nu_S=p*@aO9MZyFZ(pSn8LljjvFj#S%Y{QOK zNj>RFS#`wuJ4i^X>>9h-0})PQrY&j3-5XOgqmm%1%${;fQ=cgaDI_Pz` zK4e<$&e(1d%%cT}3>Ql&hEI8&8HAD4x=W~-vqv~pgznQ8Baf~($QO#H&VDd$ZW~HG zXjq*vTB|)p$(Aovrddurxy2(z2*r<#@3GxT68tCvxu>Sh-IWX9KV05;2nwbX5kho* z&D`J%FkD**G6*`;z`)faCJ3+gcfot?^l0Pb}G^%P*Zi=M7aR)s;LT;EAW^1U!hgcTtl0(YrX+`!6uS>JC;|h65zujkO zs}~gdai8$&yPbIyC6JsxG$U#cs3^lELsxC0K>GRt4#GkUz>?IV$3ws^knh%BhsD#4 z34l7<$#5f=*9%(4SaDR<2t+85`cU}>Y8;qEJVP+axJI5Cju{5H;$iQsVmtlEXgY=V z-0Y-SU-oBV(gy z;ufzj9k(tKH7g$rn40INL`Kh|{JY->i*|8nccqzbJY(h;#P~Xel0Roa;&QT{I|!>_ z3&H2+iF|M>b%bQL$mc}HIsAWa+?tZ)?n{Wa+}tj_Gj)C0efE)GT*Y;j@|EbMCH4`; zGq);EULO(OQv5;d{sa<sy29(7eju3-tqmttyS#a1;kyiLs$>XlP087lvFNc=oMZ1 zJYmD_8vjsnAD_%!FAYp?Bq<1k_3;bncgB(KX&Sxl)^z*5+3}t*t*>R`(K}Y*qSVt@ z-!Ql9gWnIN&)ylc#z}&EOEz6|pP!t-tnwWWgn7m@Q(2C1PJDgSCw~;}dF|xPyQ=k- z)Vd$D@I<1$?u@PLJ&8h+@%5+-{_xaIgYamq?EvWN4RXEAe%0{Zdz`NWdI>y~Qj7!R zQt!QzL3A$X9P1;8V5Pk22F@iJ2fHTdf`}5w0I2GlaOlynf34q9yY>!s(?(Bt5eqkm zd~t3S7E}Xs*?-U2|4w8X#;E6=2~1XQN$EQ4RkX?4be0;k81h2%+vSG`4HD*opN*!^ z-b9s6rQJEVVgwXmQSb=55d#WvwO#0N^OQ(H7$?U1>KifF2FE;YS4H2yPWWG~Z68Ze zuOQXNsHUx7^_ABP60jcOYF@z3LMbjcoyn#5?)_K|0M%2q-iWR85d@kEQ&*1O^ zRPiuAU9?A60}>WF$ugSMcQjK}8JUXMg`QDnHsB=oN+&Qd6x?b-Bl)$6(($PzaR(C@VwMpm+*}BCB&pd( z^C2DeS_0Z9tqK%t1(v;)O~Y{1&P^Pnry8Hh?U{c6sW?$Bw@}_k_A-R8Y<&b0Wa_aZ zXt2NJe|dd)d)QeT@y@4uMv&y^<&{9fMKkjFqS$A0iSd>Y|A}8Mfwha3r%bz+s0ci-LPJN?AGnJP64_@9N zX99z1s93s2MA)0V!Y{=vV5F6PWMg0ELV1p_;B@`hn6Q*p5@Ka1x~zF=8>RWW+poem z1p;o{B(g<$UU7+01*#jS;ufU+K>7XLrK0yAu}Zw^gM?dh+n|oUS zCISfjS86ZXx&PkNiZ1^Dr6p#Ri~vyP*>&NQVs^-PnK0bW)ih>h67AvQQTn2r&^ zx{v9HkBhy+9hsKSx;j5aU}SO(BaDha0Dan2^0+DC=n8d|sT5 zHr^zA^>&=uCwKPz+ncHL*Qa~P{S2sx$eob|_fvGZMENtD{^5)B^KFs7t%iQ_TwNZ4?#x?rON!5h#EW{HmALn+? z0C7R5%P$5aJYj=s*%derAcI=?6T%VcH+X}lF%unc0<%Ik%C_e)ob!-p(Nm_omE8{R z1}#*@tu~!y-_QXDh4!Hj2Tw?{DYCWTKF+V4kDckA_h)eCRU0bnlUeny8>V+mNS@7~ zKi+C>J%38aPBLTYz?ZAN%S=(-O{We|`z}c<1t$>2tI&AVOKnbQQ#y9tSXoa{I`-~< zcH1MPQ7ng^I)*NCD!c*YkSUd>n4Wz0>gSE= zAa4u9!`s24@Lm9;FpjdvvEAd1*kb%V^&&<|zLk^wqv!Kw0lU?4>+^UnO|90iM2U6xM(Fg%kD!{#AGh8?@?!~F-#ukI+`k1} z!WC}Ge6r4ZO`ue)0zi-$4MoIuv2ZH-;OiuFjRc}w^~Z+%Cci&nQyihx0V_p;TL4Qr zxe;Qq)oHdc#LJiF@LF9Y^6GJM{|( zF@c+$vYIxnjnLj1_Q7l>x%aoIXHt4;;_-L1%QKAWwU|Dhn11524$k`6hc%a?cn`}M zF?cJ$Plna~H0O{kCQ}iW`4RhKpr(KW7?qUaPdJvNI3u{&v?L?F@!5Q!Lt$9<`$0&| zq1}o__r(AP^1u|DfPdY3``Nf~@{>)lc?I%GC+`rJbN;D}0_?T#jouE7nm%jo>dbxW;*TinHmgNXu=Ma-vs zvF8Nbwka2fHz3KrfV7aNXY_?edDp^m$mWk3qEN_ly77HJojMoN6n&D&8oYwGig&onu2jeL&i=I}N* za^cJz$1CN{`iV?>s7*gD0kF*33GsvIi>{xRr&{^xwJ-%oetUwSLXopM_n!?fYl)UKIPpT^xu@FRPeQZ{9A&wj#(r>4v+GCc_V;yH|1;^2Hs+KoR;i$C= zYPp4JzX4RbXS9W1>^Yn6eaS>v*g3ihncXTPXA%qYPwsEqIg~tw$ckIL%xQDKB;Q3| zjJ}6_(t39{$wS-Zhfi5tD37_e?3V9HGoEWlfUcH|<-F$%6{j-uje+&Hl>Z>DvzVba z(sU~2f!uSMSu$y*@tVh_hvkq5?e_|PNSq}O#kSKXhyZBk2QxPp#iXOxs9I@t`ta5g zL)W76q7=%Yns6OGlTCvPIahfZq8`#7xh~%qE%)7VeeitfcIHjdZaGe)?<5(x$r99) zV?jjzN(?bhLS&WV*F16IGc?9vM_XJ*BvBLiiu5Z@aD;PW4e2AB>Dcl-W$WCjW%{q0 z@ogojj8A44&uX+fiAgCLm-L&~Nzhjc*4y@-m5YuPulw7bZrwVs?ij8}6>Ed{Dy<(v?{gI0 zt>RuBs6PO&out)fO&NCNGE1$}cUsJ_t=_XS#@_4szIwDoS{<3Xt6ODCp!BJNjKKm= z>C1Yus;=(kNINs1V-6kWW(J0uIv#RB-aY7b)4Nl|xjN=>3I7rD z4C9?FZu??&vvaC5+R|&c^7&Ctky0puUNeHz756h6%iS_IraW~PhR)($?YAzc&0q~$ zL$S8Lik0VM*nK%XH?2i+&&byjwqhsTbSVs9B`SBE?2iIJeDW~O5?y5cvReky-St5{ zLOI_&f&9d6@TSNAr|Sd3{@?b2goKmdN+=WMBHSkxMcA${mK#7NLNI~2eOoT^RXYPz`{|+_s>cwDXWMdO;>0$y zIB@>tHanlNKw-l>(}sa=Urs|u-M^&UKAx~DXV15X6BDERIhy7;M^Z z=LfBl#h4WTRq~$er)jeLRynR03N~n=upzT=*Ek5BG?@0# zBv$LqI}RTwVG+;0g0i!0^c&uB0C(#!D!gc4e<_LoNqUjv6s*qrZI|H&-(k|FLhOJ- zRR(39rFU2OESoM$S%bgNMA1a2-b9_ulONGYk^5vcb28{sa-fl7UDGJ0z<^@9(xa6@ zUW;MyxmR`iSG(v!VOP2>@cuiwX{8(FN&goMmDisW*v=!H#-IDQ6asNQ*7vE|`paUI zA4xI)vQ!3A!W8Pdx>Aqj|IM9&_<(=h*)MYj1_=TL;b4Fu__sS_VFwHU@n5pmo|e|X z?AP!5-|gBjZQZ)2mye}5;C;O>N42L`Z7ijU_c-M27~{3 z|B8SBFc|zx8-J}r5J>W<3^kt6UnpSmcbQBa5I>L~lpHHdPo`z<=4I_{VeJ6`gLna8 zAPAn^Aj``Q0l~pA7>ox9<^%#C{yCeftCg&|moKT!86V+53lxIF1p&a{ z55NUL@FxI5U_jyISOhx~Dg^&~4F-b$BPRp||C03IIN*?f@(RH~f5$-~;9q}1|G5_m zf&S@%{P#ma;qZTMh70|3voQGA4gMz~3<4JZdrlAnhWs-J{r5v7{?G^t!vE1kP*CXK hfq!s9zq?j(3{cer`?KKI=7?epC+ zx@XON!T&kEh3_$KxOp6hOLil|*gmTe(KFa=tFUw*kTvM)OfUOb zh{3(fqueN?)6-hTN$R{0G9+_IhL2_E%eE=%cZ-F0CxfSgDaK=#_tvUY=t2$UZO>=S zptY_S;})flpDO|mA8oB?FVU+3ZLg2SxQ%b_SMAa3xj)kHbvG^5HaqFp-$0s|uy$Mo zXE|zXf=){d0*{l=4KK>lJG(iT<3H<$ab=QD4U@XfT7Yn>J>`Oa!312VE85wf=PU4LdoJ{)CA*b{Pq2mJkBH^esJXIy5aUIdDBi$D-J#}b^kwdgIw zZ!8nka*U!_H)7lN9nV7au8Mlm5jM#DLo}^m7ph@!#3y#;*s`p#b4~e$&#{TaRFaco z=MD{{JJ}(-Hv^iS+?byevH02g3L-&Kr8lIobbHLi&aALt4Ud6kvF_``k{kg_m<) zacpzntbQY5da(}B8{ImFpA%@ClkWjAxaH)2-R2~aJl2dcfbgxrNV z4Zd%BL)R~PN3z8iyZh>G`so(wEu>>p?7MlE&W8V&LkaFttuj!n^p|A*a%G;CI~79Jc=K^QsnE8%dr#!oLsl2n_4nx+N#vXK+wc>0Mh8;s%&C(7Rd z57#+tq*?ea*_)&mQ?6+Cke~DJYN(|~dr=5iMgpla7LK{)T3)@HMUF<)c~eyx8xx^9 z#J)XcTXBlc^f1#{H`J*Z8Ll~un;iI5&$xu()TB4TfxgK5%9|YU%ig4_B12%+PD&#b zDX~hmn4&lxF;7vD)gS=B-8%CMd3i19OqPs9Bi4+DmBdSBhyhi zlZjPz<0$;5@GVr2cOJS7Z8G`)>_F z|0<&C%o1+o@6En(Ln)t-EZa~oMrB-YTj5)VO*ICq!-swk;qEEadpbm@UKx1y5&;fU z+tl+wZoBE@EI3wIT1=S3bZl^wu-RUHq{HPrL0nQG@dc)xl$??Y5FbvD}0Rez4x= z{55A8i(4Y0v|lA`(GmcvQ@i@KUdJY}N|ikE7140`@{9vn5m~*5?%$@WSc09ef)$*8B_YWL7Js<9U~n8~UIwUoL}Eo4${q?L35{gS%Rw8$3#A@p zakg@U#2efew%mL6LEV96(>CFR8rT6C=t zL>hBy%e9%O+zahUwl!Qd-XNd;ld`SLc~4(i6NT0BuJCpu9NW9;C?WR~`&7+*YibGJ2_WudoI`|VRLvDlKJsAJvX}l&;Dm3)#lH>L8R9-E z-&?+A!?lY#AnwWJ0>t z)e?Uji%5V2bHXEb@Tr`G0fHI--XfAYwKg1F|WpC`6 zuHo^APwNsn$!!U`ba7k}m!-$}LBVl5Qu9oouZNK@>yfV=BF7zU5%h*7a?Dz!e6$^U z`b@FsMk2tsp)kivCt}80XMZc9OmEW#2`r))xbJjq9YH&?)BqJl9{iAeqQJqax z!ikl!IV%3@M<^JZ3>6Lfy4X8@%VmqFoIP4)>rfSN8x95Rk#M0&NP6pU(#d{-RAg~Z z%TRUb+HD*nf=}~WP06%KE1Q>wvo_d!2gW4x0CP-k5O~a?5K#^z+_=?UuG+z)N{@C9 zwTYwg7+cUAE=WE9#MGw`2SAIW2q7z919p{Vz4bHkh^q&bHI#w_8Y&#MW*1L-* zUhI~I{+VCN+I-a$y)Fl zn*6WIK^ej{lhTpUBMxGG9&pCfEt;&4Eq8Gvuwd0?nmT_TYsEHbP-6vyu%#N{lA9tY z#P$+vWR^}32i^~E^2R)Cl$z9<^G{E@`}1B^DoLQp{4!2!qZPFWH~H8J!1+uO|EkUm zqT`Y*Nt*@!RN$%;EsN@|j5NZ*f})H8?$tXH1QodCSmdM>l)RIY(4{_ti3#iBN2}6>ag^rxD-U*LrlNN0ez~Oc0Vfc^1 zhrRCeq--?Hx^BC6$Pd}_bYc)BFwNjTFCoNdaX9lS45F5am=#(ml(tSK*Rcd>KlKu8M~woS5xskpmbBj zJ+&khzT2d^!DW&ooMgVa>07DC4eIDFgc`bcmHT}}AOI>%TtQ#TLrJ*r2DD>vF#1)c zT9~$}bYOT?8J&f&z*R}_6RReiY;;AY=nz8UdhAHx$XnQ+u$H*)S7&ng76&$0Oo`Da4(Tf@wD zQ#wv`@A(F+)YF;FJ9?egAoL?_QU$TNVivL)Tnc9ZMd;t6n2$G)4nzo>SxPK6rCrfC zNvU8d6fykCe^w%T9l?B#%9-LE`i!{EMHT^d4RO3XckwZ`-bpk!j6jy}C;Y4};=}wy z>{+{7g6#6!ZGEN%yA+=#9J{pWC6q@X%$?)zCyzi~71S)u1!$iV&{h{g9@7i}KKt(@ z{sMt1+i44~-dz(giniL=13@U>Mg>(F zBKuMKt(VSxKb%$TW?s+~oDT+FSzrG7Pz*K^Fi0Z2%z73`uJ0Zbg+s*lS#r26;_svR zEp}zzEHD`K=Y>D&WY`g@;Q}yu54m_8@QrF4bpr4Tp&BEk9DM!K$8n6hkoO0e_WMsD zzPW&MR5h&a%$+c87jc=$3S=>^@0w!iX}yocT6}tzJV*JxifJjaNG!_!PjWM6CV`zh zCz5xC4#ms%v2CSt)aW2EJ|RrjK>X|D!TNGH=vW!ASj67 z=I&}?;(*|lRjsw@M5u+`&BoVg2;#ev(p8c}>bw(&zFZ{5oXE$C1#y8U7No1GWf+06Kljb_2Pc^Bm^*UR2G5pv9W%42Vk|Uyom}G%F#IHbHBL(?uULy<1)> zN#%8ui{FI0;dH5aTRaUQ4Fv>!pnf@s-oO6=+MpD-L%i*U#`VdyJ_#3u>2&(_dtA<7 zd;a}X_M3;kSG_vc@__e|+{|~TFKJ&ms`nX>YK7B4(P;-bLK*3=M7IxlAHqiwas#Qm z*y76ojEiApg`?GpUZPQ=^1Zk{SNRQtrNPVyt@o#;1|9TF?#_4mn~tl*70Tb0e7m}w zg&xg#_gh&3k24W2mpzV~{JGydUx)76DnTWzGbhUjhHr=(w^wZ^PJX_eTpx`}wK;9? z&lUjCGnP#>!CF!MaBJk#Zskv7Q2O%9%>~Jz;T)A&bPx~Q35b20qee0-9Y)b?!tbHx z>+(FP;-cGEL8#F8_pe-|-3NuM5^p)Uo47;_Q-}PUF(a$V_~i>>p`QU(4sjkBvPs4*@Tl+#8{R+&EPUu$Oc5tGPS3^q~b z%A-j&<~nm{pYh(TH1BN`9J`WpGj5VR=PG=X5TKc53cDz6B4!gV#_J^MQf)}&Jrw3X z#iO#^c=m$f3Uref)Gt44Q^a!s@y*$BwF56+zS!Vx|B%7z@ZO$Qm6v%0s!Bka|AdYQ z$*~D+0A)+xWrZXwkRD!?ji0DaJM2a;N5$rc-lO0>p7p3fEO}MAxC+~gwzB#8LmX&fE zC<&QM+pCRbELzRZj>vUU1pXfCsQ!Qh>MB2w7zaS%tirGwHAi>AsicGqrcl_Bj>OW? zrZz;IHPdF2Tyfw1qBYr+b2t%Wq?5HSZlg_oNJ`DVogZaru25dSQD=&?C|3xip=@rQ zzS;fMva3dK$2OyzZ&6O<$Xk#kt)3U+T|S9y&XKSXR$x&srZDM~uR?tDQs%J*G8~&j z;ovw^$gaYxrZ~+`^e;t_X!jjMTwkyEk_ca8+sL!lqqhGQuj!Yt6r2)Kj9@r2DAuBg z9X&o@0_G%M@6mL~A`#hc9F(HMi~xOJdP)0(j);peYg4Gp6q$&vFkjI8N1Y7ICt#`2 zb3M`9R#J@9TFNbHxJhl+J!ON!eM>Bt@0oKE^O)eCe z8ecwR{M3q|klINX)7YfG(sK=uO{E>67#hp)m0qMeV;XIQ7^;B)h2pbGhB$!XfebB` z^RvlnxUhB!8r*mQRNwUtT*$uQUnxTwUKk)tFTaj1ohIg7m}wTzKdJyqQAwUs#$T7n zSEJq|A9KJ7w-+pZr-yt?Dn32Hw-zXxEAl`UAKh=!+<(t6P39*YZ5nQj3$ z7~H<{^F|L#p5|ZDRug2LN^QZpzi`hs59TuvDYqH)>+n8ymmSq%9EI0m6;M;Uxdwfu zlJ)A)&90-aOA*Lg=Dq;c8KBzge#JQuiXkt8-Hh*!5RG-T?8SsC)uI z3ic*pU0M=n-)nl{ftM<_NKWpE1a~mFY0e@X;?3;ZP|v~NP&3wOYuZItOVafi%(TBo ze{t@Cod}-BLdRgj@FldfG7EFrgGQhH3|f~rU5;agY{zNlnw16}h#CG$WRsN;*=bbf z>pPzer37W#0Zf1|CmeO*Dw+8s^6>Wfw8LlcIWbx%8D-{K!kq=COXH&^`aXm%n>0MB zu}Q;fGjUFxns$iI28T;pLkIT)s_{LW+L_-jHeuB^w8yMc7~S-nU9x`m2QUq^BXDkA zmbBF^`rQWHw`76PCBoU(c?C+WS=0sxHW(uAO;A!PrSO!cxcKY~e{Se=>I#G_sJcB6 zq978^RdO*hR@<;kXs_+p|cf=WsC!Z5l!z!WyY7 zPqS&{TH@QtHRS>KC+hV6ydFZB`3V!2{w(xw;J=7yW4aeHf-Ay+liD7+7c3^^(F2DML< z99597a1Y|qB)b<4Nq@V81KxEMSaFtb8&eh-Mxd!NOGcD$tGkLO5)4?;nV?75aaQV9ZT(+ zkL$KS7N0hc^uUIov|+cau{9CkszWvQ`<=@VY>@Zbv)PW^V#AZ3%i5#z2FnN2Bj19V zDOQ83uOY6k{;j0G{+9wHKg1EIE;OCX-ybh{9k@0^0~RWL;4fk(y!7qtiI#P1$uk%n z3)u0>8*RngVY`w$Fo$1CqMYSpLjBiLZ0w8Z4RPt?L*J@4Wf>4BaXVG}=!s~&OdhoK zFxYkFo=aNZ{aFqBp~IZqKOh`n_FoVz9L@i$5BNKCGPOnE0k8o6kpKXIm4&tO3N0Ff zOGj5+N?eLw!NlF#!okGd#>_;`#L?{^2(zSGALz7LWj1V^#%>~IUY zKoIF$B~17w;0?NWyE$>eU2t`8zjRI?zTG{84i?Jy)i*1gvUxqGpaGxqg3(w23g4Go zcXnmEpKyD0V1=ob^ctks8g}C4KIx->61=E^W;@&j-_4zkJzbSO>}hYoQHzC&mG@Ck2~#Bvhv*I=8>367Dc^1Pa>+U=9QitWw6;w76IzT$xey#=GElOo|E zz67MEPvu&m)C3NAenz8#G={Fh_=)2gD>F;J1%*t zdY}|A7chOByR)!O8+TU!(8W-|-ZafTLVwyG{3#q^@Cn(~p-Rrv+r1#{Q0b?^sN)a0 zT~A<$t-j;X5jy#8n|~qW00!yQo1_veYz~|N!*_U&BV9r(7r};UBWaX*4M3puGiVfM z{N0r|bgGn|Rt|L*=`u;dcbQ(Uy4zKz1D^r2R$nzuH=w-zyS%hDTB&FwDeCr>@rX~f521R)$Xse{IK=8HAY z4@jk?V3`k}kJH6`=B_9CY^#U!C{Tk`QT9d-AWCU+sk-mW711HjGwDbEyIJG-zDb*T{OpnS zEM5Ac4ZO7k?qtYrMz)-(Q-!foWZK~dzyd%%(|Gy%?^vJs*_O#2XeS0lUF{L#A2%cNgK zkX%7{CVnv%Rb?i|J9vTf$9ND(FojnX{Zwp71NOEwK-7JM*bV3c|R(xy{X0&>OkTQ%_&987KA z=KxB&(rUu^%|%wZ$N~^0Wp!sTLq+bYy<{@X;vylG*=eJvT_h0WMR*A(ndP`s7MGK; znjxZl&8JKXk7pHd0+K_@2U5cDD4R-uCBM+XP^Y}S{F3-tE^yMxgdxQ7*iZpUCyXhv}Ye(@X3E6`dxgp>bM>c8tfIq4G5T1#8`h^sueI}9eiIT+9Tn;tfp3tnQ74)j^>WM(-6v8>GI|Uvq z_DzbV2nwHvd12EN*YFkpTtK7h;ud#~b_`3bN_H|vq(*o{L{C&jF3NTfjmhgyn=hvI zTem1@Eos9?B3oTTlqCXfMAahO8H%DQ~d;}$lVp$>GR ztLA2HEPTmf162ye8@^78ZJ)k3wR_zTsUYV19=yL8S`L4=JdB+nxO}6Y$1nJJHCr#n zF3MyeXYQ;&namZc2DNcsFDwc8qVe*#?V`714NO4mae zB<|N`F__`lWZoenN5?+!EQt}U_S7s^&DC;NaR{Dzfu>!_Ml<68E$x(PeH64>CY?p! zf~rbO%|FFsB9-TC@%e(cG;_lbeY~n9mBv*VXC;s1b&y)pE#OT!Fj=B4ZQ~N8Q1qn9%cnb&XU3sMT=H-@#n|AXq7X2b zlNmolgJ#N;bo#o2-58DTL8>VO(AXz-IG7LX5?J)0i?G*Uzef{Ds38Iw z7|I@8!TU6c);_f^QFL|qq#20pG+2J<8< zgl`A^T&8;>hb>E&kmEUu0Z~i_f@gt?mBnvPllHlS`=eVq%Uqt4phpA~<#SAuhtFFl z+oLv%HNNv}b1f8|>FyJk`M-A6V>w4{Kx3Z3HGU=~5ZnbeOfMGLozqTrP6o`gnw)IU zN`S7yGg@@BVm_~k$&oS95W8aUNN>1Wo1x>XDMG;Px;pK%Avuoa($%WELD4`)cA}o> zf#Dpb_-BJbCq99U1w8Cx^xXXY{(~>1^9R{0;*hlgHe0^6t79j$Ob3Ob{o*YO$kp0Z zr1RHqNnzW&Wf}T5MegJy<~Qr-op zY)3>Ta)n*oi_?A-b72*73eC}t-VFYi{AE?)S`YQC?8^b1w&;}?ScaP}-n~66ivYZO z(2Tt!!CLi<$jl(`!xLnv?tUnCI(#u20H!dG#Tchqq%rPCJ6@-be7L(QsL2vvT7E{| zrV>8{T2y+4p$;J3r%1}+j$j34$SrOa48Gm;P5a3+sl^nPhIM%T8{5JP>u_Avl(8NQ zE8>Et5gCS!2w^3VEV{d_N4088w*o%Dz)pO1p^c|f>telIh<2s_<`%`JC&uci30XXy zz~4*kp;q#sHTlKP>y%&#G>P(!SXq8)i8@}xwPZET zIW^`uV#ii<*0%X>w=5^Qj!4V=z@6aQI((HpFYs>t(iC1Yr`Ul~%a#r{(C$m!Loy3L zke^$CxtN~seqD;#F2DuNiRF~)e**K=Ao{p{c2B8jQL}{N)wg37EiNu^JL)J9zS@Zs|Q^!LSEdB9n3WU!GAo9$wm zIe0cDvAn%~*`axI>fmnCb@y^huDzzZt(xSiDO~Rq$odsz-A`2L)Gj5!j5^e1!c+(0 z+;;-o8;Wb~5GP7Uxr#V>v^*wIJsR3ek(4jfKp>rDFdaMV+QeJla!B7jVFIRLrn zuM*hU!4LpKz8@uUuRS%D0MFVNb>#ix$O`)ZB`J36Xb6CdXc@Ld zo-LA-p@=NeV};+)W{Hc0n8R!RhQ3L2lCT>Iy~{iJOpoYNzrUY0gVT8)7xbuDY4hYF zVfQ7;XE0vfkvKfB0>P=V%GYDGIkj2^fgAXiVzjJ($U(9TYOzWh*@)qbyy! z97fC_tm9kBwogAZXRh<3ULRC7ro2Z=)%Hw2cwS;;P-6s0di;6CtU>f)kii84;og77 z#i88SmMR%F)A*zcNNd#)1NSgbeJg#`kP;pW`y~iX8*ABiiH2fsk!rKBr^ucnFz14eJ{St7|iA57m!@=FrUU-yCx&>GF4g(N_-k`*E?Ek zKh)rVHhe(e_z^f;Kk+pp5_b0+Q=3TbVUpWz#whu%#%Kp8BPF90iS*P}+R+OY>N2Wg z4#mRJ_gkUkfVvEy7qR-*OmW8cb0|Ainm_K`PtDgMT`hjG70$2Lm2H(q(|b$P%$K_1KuSuWWzam9189V8o1kr|OzfTWzU ze14E)o2?pb>r;6ebyRGK;$SlM4PG~h>#v7k~Q)9@d0ul<&=$~Ak znz0i{ikQg73}9uB;}hy?0&41kqs{3v9hCT~NVQTo`+~H^Tn^?ib{~<=UxUM>b@W7a{b66>2)+WBJ%Y$w#I+533>DI zNL#F129?n_Qr?w+U8Ey~6lRY6S?5+tWzDTDlE2GSc^=UvcP57O5Om=+9J!aXJfUf- zs6_jVyncEQl9m!Yskg$I=}>ZZpjZRGz}n)ln{4PD7^JXqkWAASm%XbTk?f>e_Tkf` zU-)Auuar>zv-)}2df>9}4o^kl90%ewr;vvXYX&c3&GtfQ!x0P|}o> zKtf>4f1%6@sUU7x8l*y<^9@s%&N7G-Pk3B;7k+?ezK4PUDMqo<0(UH0vat8l&sj}V zB1?S}>E#Rp?=u7)9MbaI(N|o-0XDx7lq3a=#-Nh;p-_}Y2ELrfqXXIFnhQ@J=18L` zEhc}Bcu9o;>YNtW&fvT1@%SJ2o-ikvQ{IuQ8x=fTwixllY9R9<83y|Ku}j_$(D8!T%Fr*7?rUG=O^73y>l zHshYnlXvin9H>K~>dUGtErUw!!u*cX_>cWkWO+S+Ln?c5+}y&ixL@VXnS4kdaX%MU zGkL;ve)A7!Ehrw_UzbiP8CVJ6fSG-mq=BR7R>7t-nggj(!K=T1n_tPj8O;|IY?=1W z=Nn#l{*XZ}NseLTCHrG;5X@9$<`FHe3c*-s`Pm(#Ap>+uVKOJ7%}h&OsDiB4Ar_S_ zE&%H#7Khb7iz!7OF)Zxj$fTpbsU_&b$sBS}G{9IZo;r5eDx4bHHoDg=sP$|;|C6l| zC$}4PU9y#mlgi+~KPbWITiDpTygn4SiMIZ1f4laKUc9-SeyB(UdQiH^d1d>p5!phH zc*QD;&;cS{)`K>PGJz9PXN#;|oY=FnVGlJeZ(cCmtHP+PUc9HXBGM`)qF&NFIor*; zY4u(^{>^{F|KffqFi^73vr9Aq2$B_^l*6T;;24jk2H7C$9U^_ zSui325S#j4Y(8v#ytPwn$o0^me8T-!T$tc1Jb;f`6@rVWgg>g<5X3L5G^14XEtxt+ zwi|R;JZw7l1GFMDln7&Qlw+;bcb`OX_Hh@#!zW0(J__EW_4B04!F=pc*m`&JrSqxv zJDmHC``M5Pf1d>Z+lSEYQ}@R)=0+3+S;j&V^k~di)HSLFvi)syBwPQvkDG@N(3O}G z8T9C2CmgLU19W|DVwZ}SyXGlDvB9rfJ}gii1)MR?DOIB}7hx!u49@);By2w977VO< zAGi*QSj;ISd{p}eEv7%b|9MziwhGeqm(48Pk+^!h<%*uz!Wd|b(%&AM+PrX2OL^-^Zjy*JJY-h4~3K084+k6BUv({CWD zQ_Q1K0pf56=N$fuO;pw{%8;b)z|gwvaw{2o%TwyCK-BYIo;nr`#Y9E?fmy4bEfe z_Y$X1*iY03fzz@FP<-5+Xdb@E`~`W)-Qm-b(R&4LIpV%^Pq(y|q&anjr|YmK{R{)X z^l|Ta@JWyh!b)t@4YhU4k4;CR>F{4EK6Qs}_Ye$3(5FZvdVlkieiwzErM(3$-y$9W z?n!mc554|OR#`RCi89GFI zBX~=6VL6shpfw#!*R>R!sz_OlnGQH`TN=!dQ(`~gk$%V2|k_YWa8y8~%S*JZOfpPFDnkPbgx_{JL=38dv z?$&JfKOD*KSn4M<&KIG5_Td*=s@-UQy^Q>N>ere~kVxPdn3<)N6LW8IE8L< z98dRK3iqSbP2p{_fxUDm5#>|vjTykFww;1oh>o--ZP{(aK7)$JXGffVJhL8r|3vf( zg&&xa6`8jPVZ0q17d%>zW2|RSwL#>I<&!U99n08{{w$T9l@T)5jAeTs?M zftjD**wJQQlN2i0Jx$sAFTR(=}*Sz}Fm_mF@LdaP&mRt@NjA7RoImxa!NSJuCh#Z?)tOuptl1 zwc<6HpvXCobOPp8Ec{+|>|L_Tx}&^+l)+?iN6(^9HJFZSla97<&*0;G5n)=glQf13 zU2uaTk7|k$LT<&a96GDFJ@Nrb&F!-6FFD(QqRj0o&_Qha!&6OvbfnCKTzr_1o;segs2;!_he?34wz3iVq^c(5NpUAn zwgUkfh)muA4!`83lxt}*Bx6W-pC}d3^@H`a^`hV=jSR(+NPXVN&5mwfb3coccCe*I zleWXnesZ!}kg~dZ5Lm}NNw_h?iEGpeAyP%m{A>Lvb1joaC%85_IfbS&s@5s2TejRL zoiVieJ@CDxWfi%Vqo#!=rhR$&+;tSz45vo}MA@uRcJS*sZcXs@q66DMP`R?FsOYKCo;bqM zhvJK+PYwtfYWz7y7YWuaUqZ8S*0#sRSF-aE-2SLjqgX1wRb6G?#F$4)gK=Vf=q~mf zzeRFDQp9EM{R0nPUvvM_lwUDu#LdvFylWyFtAvjg7U{*xc7;kl7J#r=5@rzsvQHrU z_KZ_2i{*qaudru#RcdWjo7HQ5w6+HGQrzo+>=2u%ZCzo6@#X+~{QY*a1zH1n8Jik| zuSt-D-%E_qbWvQ^x7=5D`(yogi~9qX+bU|NdaREa|4mwzedtAZ_1jtQz56tl9LLxP zk(wzdMz`pv5Qte77U_`OdvU}8(46=ivk+grhCO2PueSiWCvgg#IcbVpXcN_fbN&`9 z07={}z#nfZ;fFENfnFkv>AT3@B7>wuiJ?T!#dDTzON!;jtreW-(A9qHr+l z+VG4RO-y)T*e>pAb}yPU8%nUpZ)GH?`r2j`$ZmtO1qJTS>amY%H7%22P!SaLB8^#*p#6y_0iFDlYbB%*3UBu2nm?&&%T?&?P4Y^Kfb6HKEUA z=yeK#&n%ox8taR(?2!u#2r~#}L~#hlgktP65n4?f_tn^7TueprKSP8h z(sZ3d623wvBn$yK0F1)Ffe)kyf+1*gmPOgpCKSXkei8nH1A^3vTVl^>*9qoDKxurR z5+Wna{L*zu$W9Qp=WWe)pY*bv((Dh7cVdD3>#xTb8t*8ssMp|I@%CXvY#)2rQ*tt# z#6wd)sqK|Cle-~0hFhAx-9dJoazY#Szga0r^9#t%ub28(bB~yCN5XE(t^d+X>Hair zEu3?f0vW74OhE~sycNbH?1(9WK&>KkaEF$OG{Jb_ZN|fqLaA}GY(>aR%Jw~ zWvsd2n-i43g|N!d1cy+lGDwxx`fCS{Rj9Yrit=cUI1wZ`4e>mqsgP7OqsG6f$SaW# zE+oZ>s`}`24;EhD`mgR^xn9h@K4FU;15dq>e*D5+Yd1SOJ2YMU=8hjJGAR%iOwwP? z3zIFKmS`%?FI&1hz*Sz3H?Lm=axJ3$d#3*$mOt)HA#6ROxq%5;x%P&aRp`F zPYVTk6x^dlTxc3c6()#T*D!BpOzYNXpLH1V z6y5Ufyv_R0zEtyx0$*ZZ>K>+5LnzO-giI#sVLM1#Nv^xV1hdAIiNYZOK8o@J*v+?& zcBP7f$~qI@_OGFRC|+sp;a-QWi(zXUCdGe;&#p6&d0ZsCcAS$QIxVfyq-iu28u<*mEK3q&jAc z5%}PtyjdG|e8Ya*G8p#JF%UW;U77_eGhK{TO0vccnT9HUw;b)Pa?FoD8F`LsVxa`} zaK2KJJ9Sz=j}9SG_akTL_dp8!>aaiXko|8#i|C}BO!|aXD2QgD^RiXBks;kGUa{GnHalzA_Zu##=Ve7urP}Wj7|XGRZBnJ38y@9$%s%32)xmQLCHg$=4{gy&mr*xsgwk0%dFRS-rFfs0kSkyXR$=kuGvlBv;D58ZfZ!}bL`WVs4sZpb zENI-ZAAlzD;5W!HFJXR4F(aK?Ots80<*#?|upnx8^df8gI_@GD=1$NVbqE=kliv4hGyNas zXSKiG?%*^8(-=l4zTj*VAa%Cql@;6P5r;|UHGa4*R`2ttUkFC{Xceo{&QtwCU@K)D zq#UNa%JEjJO2@t;zahVh^Oo{67DP*%urQr;W1mzZ)Z&@C%lyM}GFv1?)OH{}v1Wy3 zqn4?INnj~9;ijI`zs7De?}wxCRRwanj)x0pgvw{S{o?>%v_VkF=RRwLm}u(*QONht z*f#TTgPe{3OKrLSKdKxXkozCNKg2m8H<%d<8*GMwry-}Qp{65AFXC!rVy~e7zxZq)50ZTs?I->g}=xc;RB9|JNNSeRIuSeRG+zQy$jR_G10V;uiIjovuViE8V5etcrDkFIcm4`a=He#q766*RATAaz z;D15D2xPW+KrVJRR={8O-rs*2903Ubny|3_^#G0~qXe^ZaLaKrT)mAoz}4 T4hhJ`#m"data/graphdata_exponent_hightau.m"]; gsraw=SortBy[gsraw,#[[1,1]]&]; (* Sort by n *) averagesGrouped=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}]; @@ -69,6 +70,7 @@ mediansFitsExtra=Map[LinearModelFit[#,logn,logn]&,mediansLoglogdata]; canonicalloglogdata=Log[canonicalDatapoints[[All,nRange]]]; canonicalFitsExtra=Map[LinearModelFit[#,logn,logn]&,canonicalloglogdata]; +canonicalNLMs=Map[NonlinearModelFit[#,factor*n^exponent,{factor,exponent},n]&,canonicalDatapoints]; averagesFitsExtra[[1]]["ParameterConfidenceIntervalTable"] @@ -94,21 +96,28 @@ mediansExponents=Transpose[{tauValues,mediansFits[[All,2,1]]}]; Show[ListPlot[{averagesExponents,mediansExponents},Joined->True,PlotMarkers->Automatic,AxesLabel->{"tau","exponent"},PlotRange->{{2,3},{0,1.6}}],Plot[3/2(3-tau),{tau,2,3}]] +canonicalNLMs[[1]]["ParameterConfidenceIntervals"] + + (* ::Subsection:: *) (*T(\[Tau]) including error bars*) (* For visual, shift the tau values slightly left or right to distinguish the two datasets *) tauValues=averagesGrouped[[All,1,1,1,2]]; +canonicalTauValues=averagesGrouped2[[All,1,1,1,2]]; averagesExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&, Transpose[{tauValues-0.003,averagesFitsExtra}]]; mediansExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&, Transpose[{tauValues+0.003,mediansFitsExtra}]]; canonicalExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&, -Transpose[{tauValues+0.000,canonicalFitsExtra}]]; +Transpose[{canonicalTauValues+0.000,canonicalFitsExtra}]]; +canonicalExponentsErrorBars2=Map[{{#[[1]],exponent/.#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-exponent/.#[[2]]["BestFitParameters"]]}&, +Transpose[{canonicalTauValues+0.000,canonicalNLMs}]]; plot2=Show[ ErrorListPlot[{averagesExponentsErrorBars,mediansExponentsErrorBars,canonicalExponentsErrorBars}, Joined->True,PlotMarkers->Automatic, +PlotLegends->Placed[{"average","median","canonical"},{Left,Bottom}], Frame->True,FrameLabel->{"tau","triangle exponent"}, PlotRange->{{2,3},{0,1.6}}, ImageSize->300],