From 9f4b087c6db1b4a8014020a477af9f6a81cf4952 Mon Sep 17 00:00:00 2001 From: Joshua Coles Date: Fri, 17 Mar 2023 15:06:36 +0000 Subject: [PATCH] Separate out figure gen into different python files --- notebooks/fig1.png | Bin 29915 -> 0 bytes notebooks/fig1.svg | 23275 ------------------------------------------ notebooks/rmax-n.py | 30 + notebooks/sp-fd.py | 30 + 4 files changed, 60 insertions(+), 23275 deletions(-) delete mode 100644 notebooks/fig1.png delete mode 100644 notebooks/fig1.svg create mode 100644 notebooks/rmax-n.py create mode 100644 notebooks/sp-fd.py diff --git a/notebooks/fig1.png b/notebooks/fig1.png deleted file mode 100644 index 40aea211ce99d60be1e93b935ca1130b3fcf8006..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 29915 zcmeFZWn7eP*FK6OHejHHlma41h=6nnk|NSMpd!*BEsYo`DIq9IcMVcQgGEcnAVVl2 z3@I=}!~EC9{XF}9pWplL{rl|?``iA2_rO)>b)M@y)^Qwb&HeiB9VN=+49CgH$S75m z6|~66$Y;pN4!k*f82%!qu6iH-Bj$R`z*XDvk*mjjXA3g*`>v1e99`|KAF#SxIJ;Oo zItcQK@Ll6!wQ_ZR>>|$3Z~uRtz~|^}$uCKZJ_R>9_E_1_g^Y|EZb*IqK($k7Jf6pHM?CPPLR7am5qUp;2eeI0FTQ#x3E3AAw zr(8?8Ll(J%8B;I43VcfOkU~X)Vwi`U=hUlXhg1S(m2aMjdjHWy!%VkxDIy+!n^Umm z_Ph6s=`u@LTUOMr4pCdGUy4{Toikk&Dxe6%A|qqXo43>7|3r4G>j3hh`Pd=X{XaO9 z9biH}1Tf0V9zi~^9{Kkv2mb$4{=b};r?4c?+Eyj)he}l%8XAguuHyQUwRywQNGB{J zV!QG8cd`ADgUf*fQk`~iPFY#mPc#y;0gOfWza0)^6uXY=4>&sGnN?9?5JkGk!^PF& zY>zDALp5#fXd7GGUit$+m1H{w6Gn#AR8@~c(36o>5QT#tr|ri8nOHyo74kP(D%NCj zloLV)fF1hX0m!fpPhr zG{vwUk-B9{l+SVletNnysj!W&n#Vlhv1q;l&y7>n@R$JV6v2hVG)&3c5DZTN0{v{CBbbcr~-Fs1AuGNtK?9L5!*Oy9UT1tTGBE0xXUUGA`Ty!pZq$oyXaBh`;A(Lr*l=vyz!*{a79tm z!l-i_FE=l*I;n49fNFqmfAzNl0utHkhy{-LP{EF_*PMyj_pW7^bV?h?nK-{@Z^7)E z^%C&CZp#g{wkwk@;!`o&gvn@)+~Q&_8PaB__wZtk=W5^J;Jdw_pYDH4o`=~a>oWTT zYgt@eoL{q6=wFgBLoEb%Q?$?rEaf$?;hBJndh^~|&ksd9F1PFGtyZF{j=142JH zKY!jm0Upn#ooPj!?U!;8VWhgsNJXtv21D&qa&)j`EX9gHdT5>0YUYK{@cR4Z$hv{F z?EsI|_I$~%E-#rBYx49vxIrYFl4=yQ_q>4bN_5ig^5mrsH&pHBpAX`zKLn4I1j86b zUnTf$;kWt={oWSYIZ{fuJkX-RCwe4I>?C-McEcp+dyM(&tHttkZ%!u_p1Z0y<2|-*TJ^)DZk1`b zle7~7OO}Q;^;S>X9-b!bZ7aAg!ZH?io!e@S9SNR~~cqhoM#gVdpcMT_O-NU}yE3+H@py-{c-=;+vz2Ifk zqv?)>s`22fLVi0>5})nyrdA~P-^P7@Fl9Q)-!$YrZ%@3iZBOrUk0IaM#L}`deD&Ban%QgS zVlAOZ&Bn^A3sxn`uMXDHGSM?L(R)#>c5iD=e3#HG?7cQ0>f)eM1q1H=vOSwuIH0SK z$|xW6lc){T@LP!%pUqAywflJ;S%a$IFE4E2QRHi35EYYs+8S>KGMf<8N?Lxhv$4Eg zBknPs;4urRFds2Hb<+H#HnMjo-^lFzrsn46ZZ%kQB&|7amHTg(6D9y4@5Z6r4He%j zTuKRI_PQeF*OPiHgaB`EvoKVKfBzhL0VZEU_ZDt!FB8_R4oCEgx_z10o3y9Xo2een zF17LN8H0>pO)QUL>e)M@a?zL7d(u^+d!|SeGJEnT&hvC>Ye;3-qAM@g?)(t!?CaMz z<5k((!q%fDmTI*9A)167hDD{eX3}6r*T?+pFGaozm)^1hsL<3$x#=d2%HQ)py-yR; z6WoWatELmZugB`@=sc866nC3to(p4^5@&l@PaH4?*lK_GXFWSDT+7h#{ixSmfvfBo zEXJXOPYL2tu9J~UndNR@k8x(h8o;`i>}@p3*sT6(SIy9}adaF@qwpR7{u=LZ9Gd|- z#QtY72WtBDSh$dS?#9<{`2q8d{Gb%^)=n-Hp(Cz@=(y`-ckuPiVs zQ>CP&+!lVR#r?Ji;<`6t%x~LXSV+hU)-x(S&eTf zq}`?9h6i|y!kQJ;q`Bks{Qnvs;CfY6RcMku4ny>m5VY**6KG;g3*oC$F05UPM^-9? zzs?}*qqAx`YI2=}?&R!Mv&~>rU!pFrG3GY%usf2=3JH(D9*ar`T#=OO8ouCIlLFE5 zsgrkt*7QjVPg~D9eUlX3JUu)fmLGHqX=Jegk06jB*WDSHCp z5*RI1l4R*t7j^Sd>V%8~MtC15Ab%6(@k@uv$t&_K279 zuLHs&V+COuf9?KFOkRx0=#vRb)Z;G{5Z<~%4E-C+6L7W1g61MjT3)!+n!$^2_mOaf zYxjQXpC`ZkX5lz%eRtsg^GCgxr26~%tQte2Ir5zb+M$QCl$&psQk3-_TGl5?+uMa&vOJ7s|)#I_>|RJrGV1#vrT;2dlPb z^QNlM4z)x*U44Dq!B6+MfGm6{GL@HqR`>S>e>UWH4i1j3!6xZUm;<)m9q*;uEF|CP z>gocL^(;5B6q0twNbQ4-|9b{0)L3fR*<95HQYf>F*FaDzI5^mK3!h~R)9X8Q$LvRf zgucg!TihM7DH^c8?VzDbQtCyZAbFANW|5pyJLf38C zpxpX(y|x$1w+Vhos>0wLdiM$C3tH2zE#jOT0NT7;LTYMi?VgF})LY+Nd2{nLWYQyO z`AzfR3tB#kc`s}md!FBPIBVgOM|ZXonJ?`+RXM+t$Q=^4YX!`Kn)OX456>m1L%|a;sot1u*NRh`i)x%35PVK zTHolu)8udUC6K~vHYclRd1@(!jzMn%bjMTtNS)SKE7{s6m2CmG`GMef*?a_5;<9b zPj{lLTc5f=0|AR7CKz%bT^O?_Aqx))7v3c!4l)^mdnO!q4p2{g6U4SNqY}p6k_T68 z5yI~OrZhJ8JSF{gc5{h^g$1vr+I99#+LI1-n~!#vBmDb7UA$GzfCtwpo!I*R`fM)= z!=nbn!ej9TK>1xo0;)z3rd0($ZjdXLSLAn*G+D;X(ipTbUgtce>N{`F z(#sKNQrtx56(!`m5Hpw0aI0a&;Dh68>c!kw|Hsh3cG zn^RTGjTeL;HlG49c6)52QRI46&Uzf|1Dh=ZK3>o=9Dy+=r|+kSO0A=j;F$gB($U?Q zqwTiwJ7^BFU?1o`3&U7M;Seh=d{}+o@hx7}%kOHEUK0Rr}blaQo6;)9Es%mQegM))vn(RkZ&5gAv?kmX!+1^63@ag?5 zY!t>Mp%=Z7u2!__TSBVRb{>El|D{({^i8SiZRGirW z1+7Uk=nE4NFjaqkyiXHz`C|q9ki6^yc?Vq!@gIuu`O4z6&oD59PqHq5e&CCKgwhMT z0vo~y$?V+s`l(;PrCSRb(+dSlKS}tGZ4XugB z1bgn26gsAKt?s%o=q>f>=AY(`p_)Pr?e1=Gnuc(=GKUdyMJE~AuE0dL@#PI>mUmV{ zxg2V-SKi-NjRIKz4~&fe%iqshglYmJ^Hols837Scys*18^SF(-d@(HpK;Nh>~Z!BxcI;O0~$%vs7|T%kLc>fD%Z_hA?K|DU8EfN|FIve5@hnz zr!A@x=a{_FiSg!{0MmlJpljER5}aGP5cTOcU}%X4TG0zK#BI!XopNiB zNNfI`%n2Sz#;WMVa!oCcuRn%M8DYlMjmj+ciEq z(|o-c)|HHG13n%*dkS`mkjJ8l{aBsU?hf#43q!p|4fZqipKIw%h7r|4 z#+n2cjk2w69ufo9OLf)yQ(%;FAO>sI&PrPu{OQvNgTC zKF0ED){+8|7i6;1CmvQGJd!8aAu4bA(Mhqejo-A!0C23key6QoIjS+7DUaV5ueR@b zYIYR!fR0$VMY&L{{`1b3rODH;AhRPvd!I{&YWjj5@MQcL;`-1}85tP~c@eDo2L@*6 zOF9wsF439TQ`fmnfW**zKqjEMa#m;pg-eCXEp&@ z3_e>xIBRR`j)&xsGz0SLKRS*V)_zK>Ol}d=> z2YTx^;)caPKYOR^@9%FMDr1`+9}Q#6s&YlQo(w6$saJD z$PFnmVS@v#GLuF($kvnDg08hAJ8geUcl=Y@XtqSQnD02^Bz*{|JvP&=q_2mOUZ0m75|*+84NgbbKlSi@VNybviuLf zSfN3#zzZlNJ1byE0gk9t`JJe`ueJXG&zEpn)a^~)8G`G6{5t|jZ~;6Cy!;pRF#*m7 zTfVbF{Lz1&NJn4U*xa1Wk=W|KkF@u?07|qaHLf!i%q*w3mD<9YWGjc^xed+Df{HCM znmxk78N+t^$B+YukymJ-BXK4=zjc-3(Eb-tQ6q(@!34|Dp)-DMo~^6JCWt=`#Rr%a z8gq6wuBlDbi1zS5J4>Uhfk}*Wih2Ty=mm62$UesI$jT~ps8l41yNDs2AJ&ouRpKohku@`o?CxK`lpxanH+Crr{z<%$b=2|W$915(Q0X#lv%R)JZ^7)p zq^5sPW^+L-+@P zPojYc&jJU>f(CYNIH#Z9Rs)e)^(*j%E7&=6kd(2&zl1D17z+%F9wfV)%NC6GxIESM?uu+s}HYSMuRH#FZcUT3F2{Y!tFKQ! zXopD)cxJ|DLB3%n*Lk`(O9Kmz6fz5KbZJ{#UQ$5dkOE1s#lLB2z6kJ2UelV2riI1{ zCvc)Hb|#Ov6?^H1pWlz&O{8{oM;ycEBrKz0JqcSJUFsGQ}zqjUd$iu6ik#EcAS zRR$?Gfg=Ls(lTG)%&ZWEA9MKys-SaJK}Sj#xARJZw(E?#q+Y5Z0r%zzGEMgO_QeZ7 z%z`FFqtZK^c94VtMmS&E4@7?XZV|rZJe0wm>sn|2EV1et9J@a+M+A%@4Xdez=oOkK z_|>>C81b6bnVdU!jzQQa2B}L@OK-?hOA+3%Td&RkL@FY#z&@>@OoSL%Ktjb43z$LL z2fP|&K9J0WbF4=p1|d?_;zd>I05dkzTU-8C+*f_s6)UKkF4>@s4rW4gA^S;GG}8U> z|8uEkABd;;YYAEiE8PM=As4vd3}SU5k_-~zb<`tBtzLv~C4`Y$Lvg_i5>F+abwolk}2bCquzX( zor+NM9|P61iZN7L3vg<{Gvj|^QE4rGT5fDzIQ5T#&AqvOm^B0zHXwxAL=A`te_t~k zpd&S{&iqAQ?bk0r(1Cp1DEK0zjIvtOATf|eWMqu_M`uKn*^lj`2bI<#SH`<0|3JB@ z^64%o*^6Kb!NeFKpHr0M>9#b-UHyHnjnzzaId-V_M`N%eoty-$ z9?`{UWbHsBCW3kYJBTHbYCb^*m{FU%F11a2k#BjZ+MdZw6ic7?@fr06m4DJG`JK@t z1V_Ks4|som!(s(crrpnvorO{=Qd$M;k)@bN4J<2gJ+#l;2H)Z8%odK+vX;AMUMZlc zzQENGmH0qKKQcKwg8`XO*|8~K*`?XsJao^9$5??&hZtw=JW@YETsxF=3x{DoeU)le z;?Nni7MimO1J7aSN?g$yV~n2kaADCruP%e+CbG(Kx$&wea|Pu>o-31~7DzqKZ>=N= z58&8wCu2py%-Qw21&vyjjgNYz|SF=VIuS70$W2AI^LX%1C=(iD!NY*KR;!MTLN(^ z^VxJ(=|LP`Yw+0-Qy)qoM!SDRpiDq*L9LkX(?S!;H1pdU_e+82-FsKT!yQ9xE195C7d8s)>@ zI?{Cjc}-yV-RUpeN&Rsxo$K;Z`zr1iS=2ge*xNSw8~YXuRtLF`Q5iP6m0S3DqN|iA zv-Lg?b6`uJNVW)_^io_<(LZ(*DJlgpD#>CDc!3=DSWPQq87^h_Cm|k_wMbhp=Y=>jN{G8i+z!^H=k+_N9qPqUR5RG zh@l7Xkf<85gRwf^J#PibMRL96y|s)#Oss@zs$09tnu7IU@+cHc`0I_$f_vs{T=$Hx zD;D!u?8yA~f(-^YPK>-6Uht$<_jAz{ecHeiRk9{TELsE=69u@BE>9oa z%#ixc9l64;7R-UVB_Dneld<+Z5A!^4Tv?VQ@>hVbhsPTC+@hmOE0IQH^MJ_c@;c`k zR~<&Gc};1&mrYn!-x-aJCX2}Qaph}LH$u#f;T;c#^Wvfc!+C)-vZ=+!KTzq9QGD2< z^|`3;jM^q`{Zzfk!ZDHE+U}7B(}XWIt7E=JwH%Rq<^7r>Si$SpSKO-?9Srw6yHY=0 zl2f{_rOD3PxlvyI6GKBKdc6@gZ)Kv2*|npCv7tK)9!#%$7bxSbA=ugv{>+5=cpCM}-?{5VQMki~eDGJ_ zF}Yvz@U}M=?evW=ao~z$oIQN4suJ(UDUdTrrC&22fv9UQN^ecQ_=0>*&Mgv}bqsnO z-_H4bRO1{(&AEI+{hG7W4Lg)zTK1B`ypKkW3no;fLe*2Tk?vNYU-BKVMv^Ip7m$1? z8XbtqXdVA#bQa)wEU>zS_jmc0yU(EqEw$@_wQ9MnAMfg_C;(!H-d1|G`Y@E{*pN+F z^-3tq!4rW8bk7G((MMB#Q(tHyj0GF?yd7DVP_Ie^aHU|n%A`EORmdh8PZ&!W*_0d6 z6a6)RJX~tWex#dYW3AEK*l2yr2fGx$NN7UYH8$5<6~|47MPh`Al7DVr`{m#1rq9wA z)p|m{n7u(eQPO-Ts4CRjxOZzL_vClkT}E4lk)szyAIv)t5$LbcXMw zd*Ux;Kga0Q$4hl}s7@=$>5a-@pC?+W#TEgAt9MiqJ2!??VMwZO|cECKUF@#>V5YRH`fTS}N?>Uq7Ys zwmaPl42QWH5;xNOc7QZlRV+F@K~dSEH~94v%MA(szitloC1vQcvqIIGumRaz&~u^7>wa1%@*w$dmQycobj393=`M^!JJY* zG5$SMo;f!{0fxr@8o_}%9m8XnjIQ&Fxb>j)Ht4DK1&%CiW zwh5`;+Y<{1=z`FXf527Ns6FyB9gV~Dg*sS=`_FU@p7$tow_X2Iwz}G#Bk)$o+`Cu< zr%|c6uu>&Cx^2L{H{CNQ@cCQ4!Pr^zE%)E4Luk7?$u8Y0|%Jj7!AmlIM6( z7Uof>)V9V#xUkt=vA6l7*LF<}V|D1Lt_&@s;*~eYI8y=&V&9QQCJkcq!Nf-HIPxxDz? zf{E=x%x=Boi*GTV zT6B7nQ?08z(Ip)x%#7pRU#P2Stp?>ytTEp)30r=q0BM(~zE>Pzo|7v&UalvOX_ep$ zB~BSE*wHWAm077(1RvfV;51K9Hg6GDyx64omDtCy?wuBTB_tCYy%T$_X>>JE?$T6w zg@gTTu2_P`FZZSExWT%sI%O4H5cl+RYFaKl8TUEvLoD37-%O2t9W(o39RFtECvl*# z-2dw^!+JFZ?F(`nwM;K6s#?*kvEfg)Z&aHwYogX;10Ua5hF+4BNzXDV!NxYDf8Rib zJe{cYDuL?Byrm8$KqTAW@LHySuIk3_#bVr`>hbro0P?prNxD~{7IWX zEU&ljw~>uO#|B|SHo)w@*^|WXfjJiG>LI>0bhM_6<(xdw4gL-lHb{aO?^4?X(_~w-v>FOPN^S)t$aeqmCPCje) zfP(52#YJ73m!)aKDLq!7u;qLNMdM2u>2n_PBR^wjL`U7F#!$*DbGfco*5N7=MrY~2 z6w7GDZ{@BEM{}ixp^QRiKRh^2G4HCTBsbXImHN0!&B=v_!_aU2M}U0}2e3Nr+j31y z?GKFoJYe|*eV4MMe+t*9@V=KmJ;!sFa&yE_!uw(cQwhJTZ^YXAuCr7$RkN`~rpvIe zG0L9>Sb%~Jwu)sB*u-!I5 zd-{jF>{d&syNoBe>}~AEcv%-b&kYGH;*8GJioZ82Dbux}kmEpUz>F=}skJd;FM&c`LQ;t=|Eu7ddxF-7%t^T8Y4~Bydrx zRn-TLEy+fQ7_n=ysNa=UeJXKt@tgK-T0N^YvlXnp+HKyJ`m>n^S`bF-YOh9j*|}&n z;zRMGr{}FykPz@zyEL-kxZX8lt>32HL{ZH5+r^PrC1FX<^8J_DHMP8>P&>K?*8eQh zivz`UVK2f==+si(Hu@zs3v^4TLOo#)45V)GuDyouNHXh7*3AuBX>0ce^jFI4%oi@< z$2u5gk_&%(DHa^lzjdEKt=isO1XXdPe%QdDiqhNotp zT$r8Z+@w9Kh4yzezCYOgCiT;}iatkoGGm*TeA6pZfs_94rkI&Z=<*jI6~zm+{pY?yp)XV2~)iD~HolyVNgm z53Hqxvs@am*Mbixd}OFBOX2MnknwqRYEJz45}eEZ`minym6f-eNyBMUosGZH6zW+H z9b3Q1V(nQ|?DnEpkJzW5d(0LI%93H}Q5ZqNq5)$Zt*2BpgyoIO%Zg%%CpEvJ9J%Wx!pboV}MrqP^A zzM9qW5aSX?*zww2w6nT+ABVQ0Ee6AV&v+$tjZV|gL%k^d3{HxYKWxora^oo`Lk!a? z%e%Q;q9 z^YqsZE+fTY>BoPtquoR0o`y1^9owK6S_q%-RI*RS`-W!@wXE*o{fx$Q{Ep2etotlQ zR5cphsY`;J%c3Q&z4Yr~>y#nYfF(F$9`Zz=mzNi5DguLA@Tl+w#y#I(&w7XGG`wtr zIrSj;@^v?Gn?XKH`Wdx?dSLV(U=$`f3+nh!^681|h(;r9&&w03iqjKrxFNd96EF{h8m>&l=ylp|_-AJ16os=?+SkXj zMd^7l)14<+D8(5}&-cfM53E1_EIl)-0b{&(4e4^v$x7rdiOc*BLt84!hvpZsfu+!^ z(4<3D>hEg!GB|{bQs%;Y<~jz$s&{vOJD?-fZY>qMoC5q=kwSHGh9&!h(&) zM!_tXNQN8E=-eD#q|dcK?-0|8>d_(BU#ul1!bMFlf~egWNd&h9m-71hT4AORiMXbB zd?(?g8Ij>g|4KX5&4_`vgn;T$RqD9!AdLa-yv6Vk|5w zV>iJM07zphI0k&(sJvDe+>+~2;dFZW$JkVLP7T#-^J4{zCHx8?6WD}vSQy9N5XM|R z<0_PXE&cIl&Ci%g`m`9!(h}a37~4_@?>s~k`)c3*aKYT>brxt*ZvLY!bpD6K(``YE zfS^^3RXRtx7+)+n?DX45%l^GdXErO%U2SizaC>Vtm)jCX(Fm!KYt}UW`gWTB#X$@R zU&gkmNSQY?fLV~ZZbb$dJ>as(KT_@Z`i9L=%hpFfq5$%Ju_%w)KgCnc7mSF6FT)v^ zRbI?|&8DVnN~dKgl&!&tDp(18W764-eF)rxrqTM&mja4_zQ?f9=EldxFj|S6 zZN)uH5%z{KejUbR%@`hf^M?yfmfghAfhDW{@PYK!FXYLvg>$a|f z7Xn)FcDQ{GZ={xrP&Qj$$T=#X4ZPmEQV)^Fb!Yw3l-JNe)Loz?IHPf!#dkgK6CLi9 zT*q^+FwTf(@=`DvPrZ=SeXF3uK-n7gc25deJi^gA;P_pk$!z#nuK8!-I))!6j3coJ z)Oo3wZ3w(Y!U9udasfCu;Df+4F2`h9b!j4EEH^1y8sT+3vvvo zCvC;I)_WG5CugL8eSS#KH!=AZ3er&S{#QZzCpsJgxqEi&JZiJKICQGM{?9{jbQM+o zf0E0huQjG?^sTKw_xDi2_GZ3xx3{s<_jNHTME^ve`w-8r!pW?a^2WRk!^?{*7aS=b zOnCfx;M4f0Q1+%`nv;V@wHdcao@<9d`0M#~Wye@TkmT}&TaB#g6>KDHsu#uFn78VU z>~8>9%HCFGrqAgu(=*j7FAud4X(yex9f^`RVU~-ea_SxJ?zRdBA7;As43W|{SC7*i zAKZ3}WQ}zt-tjHLlwa04PnlaRnt$3t6Iv=@hCpY)tT4FJJkBTB3N`I~8ua{XMmKsc zaDq#}Q91HnXG#N9!!j#ioqL;CvOz`{4EANOmzUe@hPxNwMssJQUIC>De3J&X42W=X z2iMBu^5L;{_1cgzZ8WO3=ZkbS>9jttkt?f7X1p=agP_C;a$h}avw*n)4K+Qv?dsPR z!6Ar^>>jc&m<{T4Wt9S%neWUZUfv)HmchZ7nCCBPUh3$b(&{0`C~wkU4h$!~caci7 zFF0mD+>MxJaH9<}Y6f^2-8cnAV$*yQl=jf&I4u^N@uXGV?8h;2Oz6&fI_qRM$D$dl zjaqt)bJ*;M2^%}LW?{txIO?(22NKuT2Q@9cR%W9Q6NJ6ku*C=9TaKJq>F0UnrP)%y zXGyz{Dj`Hw_|2XWq2Tvgq_mbJCCPuUqBhZhX=ZQk?K<&8#sxx`1!s26wU~uj)(ATD zQXiyu0lq55sDY1e>YGX6ZTHHmyLDD6W|qaRok6fv*0#LNX6^fLakt@XkaVrJ&2MLP zC&ObjL6+QaIIVucI7LAIgvD(V#lWa?Ria&6l>H{>5L}f>{=m|pgqNh?zy5iKzJ4dZ zfOjkTVkU=ex$U49l<(J(5ztV;;gCvU`?5%`kNCll|LshjdYAGyz`z7--FTl;+e~cMg27 zJ$A&-v+)X^YA5q%e#Yo%Gide@v-<7S=f#8{-EZ!$m^~AAB_uobT7Kk| z=3Cn5rNEHa(myyjQ6Jd};fZ{cMkz6*Ge#3_VF(trmqojK87B(1S7x!*UHOo2kTF!F zUY=9ZF+HWD&65FtoO+SsY9-j%tjryDxIgxBKVYqbVITEWY{tb~7a>yL*sDN^DRc*~ zp`Dl|e>J49FHuKLlA}|3$-#O6_Sx|xj1Ag;Z9YPiRu4U&j&qR$dwSMAtR_0vo;^^E zgzXGj0fYJOJ7m<%Y>p8I}7YDupXI6|UY2-Cq+hPA8_h{RnWQPbf(Y7pK zI|AzND}6xw5y}WFLHtRXog<`NC?S7m@Lf&|dXyKu3lO8H`}Vo+sTB(rv;DXo1Eni7TZxF3~2POVY-S$7ZtQ&j>ewEtpf2ynRBVU^yny&!uF~ zQnwZ()YLP=_Q9Oo`=Gw=nSs-~4fp{{tRRD#cIwwRX#f)5^RYAZGJnyqX2QnCr4L7_ zfGX;op|sh2c^F8ymA3i3eE-6Rz!bzdsK1--MSTI~ZwR^&LN!9Uh?_mm1=yaQGqUtY zJYO_lo;yS+19Lmtrj2H%;#Aa1=;-iZHXs}Kae1TKvRhCY-6bXzfJ1In_S5kNKRAvY zv0C7I+rLKYgt~Exx|*kS%&W(2ZdmA4+NiP(^>v{Ewz%E~yHylc21hcr)DpDQ7H!Mx)ZpT>jZOm&UiYRM<7wqjJ?QF=#UPX9pjA z!AKH34TaaR?U7dF@tp7hYkiK$WXPi0Ux_O+zWhef6)Su1b~Z*3EPam9T{tSCPNT<`z2+>P@$rJZ6n22epJ~?taFU zW{9CZ%7v~nKA!YSl&G&EkG zeuPmiF>>td#}(*|*6AwJ)jQR5QT_zDkvV|>*BF&Oh)WuY%iaW^t}b+f{V|}3K+uZW zJ69dK28Wi{JLuA?E|bP$cS1@NIv%*!P}7?F2^=V!*VVy??9M4Y10F!Z zWLxj3_Va^PbNBsk!CSu7+y5W_hNotDhtu5i5f!a*w_bq(` zI6t_CmPyp8!)Vrfd7>hwps%dH*13_jpGr!g(aPKQ=D8Oa%NkIgjoIuBF_yOrje{hX zyoWykmbqeGmVj&#iM`T+`0_5?xzN*&^0`>l<(;$h1QD6>?gS-((t5 z-i;ms-T^hJR11oWOfs1!Fva+2!X|%V@3zF+$4e!mj9nY?XKn>zY+6%~i(_8KDe6GF z0BQ!F9Fg$hgYYA(6lt%GyXthBxhc1PAE@>B@n3qq&=;$W@F|c0y60aA=d-)~cw^7& zRK{UaZ&4+N1{$pk@%hAB=fZ`dvR}locoXAHwmdZ*Lqs-YPuXCMlACHa*o%-gt^gca>3U7JM}bGK!3HLtP+&cgUFQtO6YjB2wQSLcYW z@eA80%gr_i?{AHjV`d7%^;+(H3;s5OS@YfO`&L>Yj3zBQo1Dx9^tl{kDdZv=H&h?4?TCCgh8xONr(78NRHE+h>pd5gyZ)*; zJN2UfuRM?PQY$6fn^R2BI1m*0Q!$pX2EeQ9`g{tK_=K%$U*cc#cn&pEwrvmSSW}}~ zbwJ`f(luJ=DCFo>%BCn7^ijk7#wi065MXLGt-K=V=jBu z6v}{(>0!i1kao~|4mPux+c#Wmlh(e=a|<~syoVqLu(YmHTSZ>QumlyPa~?l-z9vWh z<@sC)v6Nf)!}9nBuG(U3|7>upGaBZm*%Y&!fhfo@smED=E+ULq<4S$ZpD_$Y;jbrs zG)Utrn=0nD{?YZ-^#rtKAGz=!FMGk9kW=kVRYmgJiWJ9xN_nPRX8M`~<@;wk#2d9& z?Q(f0N90CM9ap9?$ulbTV*QlEnughrD<19EK{+W>IYDO2)#B7wscdK89hchu^T+l5 zwzU7py@?1%YgPOA=#A!!yPrfV^LT$pxbCzLF>QS*A*j#}$pJ(6F>`K0wa5)(MLUF5{7|%b zX|eK)-vOq^keRNRgy$l;Ol@o|r}H*bPEOj@Qpoh@3xA@4JXnvv()x*cORN2m-n_IZ z1!eD~5Bq~#w(~TZY$9v2IdqJAUIto>BCjg7ygM?oORF(5k=SG>zEsy>Gxw-(I3zQy zXIvOp_9)L^V*w9C%GIV`5O@M=ZhCli#Sc4if-|E3XSCA%`WATsM{Ey8y^lr2lzGtk z*^bUTi)(%h20i_bu?}5La& z)oVWZtl*!$_H)!VrFth!N%v$izOARWE^v!f8C>&)ujFlpR~*JRZHi_b$( zowy_+y{0fe_&)r$Qc^;GQ=EO@*|r~4>qg6zs;0iZ3k1* z`}A*l!QQ%%EY5(Fzs`eVAn)%&?s?BvNEPhs-y1yx$aV>=(1Gd%VaijmjYF`S% zn3+qlU+u$MW>m3XJ%$9Ud#b$1E{48sG z#L=pq$s(irG$RAO!;Ew56e z<+a;B%P@)~f(TCVEjD{j-=$XYO-ucINyL<7`WhKHXA^!cm+aGTZ)k&!?zYdl-kHy44`vs2#q!s8TCk-$GOK#uRx1aDn0JSCg*7Nb?+g^JOx>4SIvv z0eATlr8BVycBAm-1PGWY_AKgx;(FhOXVR!cPgZj2#m@1yosiGHa?%>oMiimnU;HMt zP0kjyt5*><9Mx8`F3vFu_^x|SWrcEMV@CKy8D%nWf&*mGGpCKJ^8+LlYs&nNG(0*p zLH{-0l9nlNXyowa4{kbd_evuMa9MReQP7R}r_p`4GuNw(WBud70*S_VuGwQBUg#&K zuZ7$H(Gp!kUU=?9y>Uu!cy?dT#7mVZY3C;y@H6G-xVp^W-W6SF-{)wC%ho(cVN{f( z5^dWC(rF;IeAew7XJZRda%8i9$;QX6$YHOq8`_QfkzPc?Im_;5E}81T22g4Ce(s`D z^~XFYi+qb##`)FacD=B748A6!l68`rdmvGHL9>yLjm^TfcCX!mciN&m+0U~@N2)sY z4mgouHvu?1*&mtqayA;9nKH=#KJo|8z zg-Sx2yBLN8`mLx2p=s;Ny)19e!*<38a9XW`#XK37F-}rx$>v;>g8*k-g^p5ZBh8#e ztXcW{_2zxW3t9G76`(YSwVrrsy_WtPhQJwX-}S8!jV&Na3ZT#~zbdmv2h-IoL;EF& z?OU@Gjx((h%&A%WeOyys;!7Kv;`SyPl5A2tt=o9FYIky(>bUhn^3V5*co2^d98Zv} zxZ)mK>a>eW3&-?%mt#O)KA@_2pfK$BKi5%{nso!VVR{fVw` zZ@-snHHz&+G0$6XN70-#V8Y!ve#9{C@ES9JE={U9BErCgP2i(1*sVM{qznU&mncmda}v{!e9R z0#0SS?r~D6q*37`6^(`rDH%co7EOdgq|BK!MW!#61}#Ii49PsCgb0}uGK9(y5|R{4 z!e=hy`9G^~@3Z&0&OX<5uB$<pq_6UhUr7S|nYG?Xis^VJ8%n^rHEqO@$| zbvAADX4=8=%Vmj-_wkVc%H+0wn*mhyprW99)Q{o%-%WOSd9<^~EMVR=R9MBoP^b_iN7pcW4xOt_bZ9f+zyX_Yh! z(}4LGayf}PS1u`?4*Fim{^f(N><4{X=Lj#y2e>+mYz%OAF(xV7X|@h3%Q{hSGd)D+7qT|H$rv-;2eOLpD3 zD`45Z{aBm97aR*Q^bz^uX`5=P9h_`md0tCxne`5%>Z&54A?ih^6|qoy`>cESZ8bcG zkR==yr`qKYysNSs*Rh7Ap#soRkLDv$=F=$>PMwrSEd9dmbK6e_BWdt-&5!3XE0*m= zx&Yze7+be1=jZd@mgu#uPP3`eDs|zav-{6oT`lLeF?Mu6XOJ-Io)jCDQClX=IqCi& z{SEIPCTZpjMQm)ODM~Zu=P#-{qZ5(Xct;4ehhPbLb*Z+BO@ zxz(^-j(YD@!x1;bWo~0jEN@J&V_%gPTHn!t1A!Sh>9M#$R ziT=Z^PDh85a6g18!w1E{+YTMLRGixjBUa~tf!49Yp5>fw<0<~$4};%MpP!_Ky~t_b zmwfD=9E;_4tF{j*6x^h;zST^uZoSM>;VeP58 zqpz3eNA#j8zGb+9`2oCIR|CWov<`q^KFpnkRb>5N^K*`g9>w{r#Y!VfK` zn$_27hB66toK928L4ly$GLAf=4;vKhsnLln9X+J0n+;__q6SdnH%A_i>+ zt2tiYJs?t3{GppCJ}yIXiLDUbZ8%$QP?*Zfj~FLTcal9L{@NPTXJ1`%s$B=X^;RuD zuW{n^o8g3;Wy(CiSPr{KHHbEceC6WK+p?)Lvtd${jo*OMVqY@&W5N2t=K|RUwAQ!H z#SNC?@!3za)@!{JjXWxl<=8^2$WSzWvS}jC>|XiVW!pJ^O*Z!E2}{C6KSoN`73LEA zJWNddsDJ20HF=z(DSDl5xx`Vhuf^&V>k?ud`f1ud71Z z0K0Bidq+~08k8M>Du6i*qOaim$R9 zl})|9Eih@`P=~jkdRyzoSS|hsnSRHm|D3&TuA9!=yKzD_b6#5aH^p$BbnUri?p zB}mIgw0ER>B7uQ;^+-iic1{B_V{hm)JD!hrLTTHNyJqRCiVdlHS%i5>fKM!$0`vkH z?bBk<*3FfND}=6F$d4+0Sbn4);%)O4Se+d?Nuer!kDmOcx&py=^avn>6MH~-X(LN8 zvKx=o-6bV6MKxfRT`JWDlm;*roM(nB`w6Q>==G>P=_V)rEpy{%=Ce<> z*)p1@Q+WUQ15MEl=sdCioizl^bhh-VD6%Bd36{qCk!}_Nq|V&=_IKVkXuk_K8n6y4 zIJO@SRQi!0xKH> z??*(irZu?R)z)2J-jrRGE>fd~!ezdo3QkdV|MTypG7GFE1_#)`c4Y5%J1S^yK!e6b z!PC{V(1|fxQHtL^H)UIldLY>PovJ$@<-8~Rl$G+3N)hoqIky1Qf!no30o|RXpRcio zutHm0bO7tn1Q$|$9^f30_Ga`q#gTe47TuFSHwlcq7+wC@`r)8h+`7DNwjIt=zr5;pw zQsd+?exA^Tcxn#B29cg}&|{qX%%Xbx0E{JoQ?`N#_jD)TESO7>;Nu<%4DtqjamqPg z!e#}Hw6D^-C9P@XTaAvUW~^qsiq)q=_o>f*%g5>tEI0f2(yF#X047~tZw|Hm{ys&v z>&oh}c!629Cr_S`CN)-1pCBts%~M^zPb{nYT!W-=p-gmPKxXywol}6|dN$VGpI&Sf z;!F;JNplHiB<(xYT3&j z9CD~QfL+KUyMeTE`_7cglOSP_Z0)#vom{0%|IFQmEHvycZJCT-&rf&@de#`bmdyjn zAmloQ2lF46{mlSw z*poUM*e0k3S!kR8l*B%&b}Wz`K_A9LRD4x%mol3EQ_)sFbc%);*rN?}XZbCWj5qF+ zYvW>KHIOkzu|{N+3hA_Gx`57g_pRtJSS+c}=vr^zucZ3v#=^Zn$Q6HGAF^ab9TK`Q zsP8>_``7hQxED*5hK4rpl!?E#~ejjq}YP~VFz{Lu4m$@ArCtvjC5Dj(%;(YX}O=x znBT!kzI+zU6PBZT%E#dm-z;N3B-bOglYK9`o}?2skvc+tzptcQL<4S@r2@IbL5Q=fj{&jD`CugL@&9V}=yYnh_uh z7CZB7NORJjjIkAl(XsvZ4<{xDET(F=*3!)LlM|y+yzk{g{TH~(+6?VuKb@>r=ko%N z&Mx}4RbvIkWQCN(ux`oD>+(f4`ugf+hon32zqK1h=%-Np{ro)`0ZhKkfC}gid^)%O zN0P8d-LSzMQH>>RQ?6As6`nlOB={1vHEQH4dXkGg{gk)ba@?X7&&g8qMdAkTu+A;0 zD=1(qaymHSYnwQy<;!Cv*sa<(YVf=2=(FczV?2B{+RRZcrUF;ur&y@kk z`%gfZvwGnj?fhG%u7s|;vSx8z_U}szTsy}16pgQhk^&NgeunX%&_+ohE$X-B+P9>8 zj|)>8o;!9EFrg*5u;*=T-V@LukCeIPx9WSmdoP&g64cAh%I>~OuulCl6^oi$T@<3m zO!Y#t3_@q)9F%s8_F|owMwmqDdrL>}`Q!=&MzD5%CKbOJ#D)iCw- zovrMhr{DY{qXY+5UWt%DlwaNw(OdB`m}%7 zDc$^A1)DXx#8^u7nxiHY$liKca`SJd0^z}<`M+5&fEYA(yOwoj zpl?^w`@+q8bH^%Y=>-c4IwbB9;lcEk1vkQ~y2|rp8*r!-G8vG zI$VNetXqoch^OvS`^|2dDN%qpNUqs+zAtSsC(+kdxN6OBpb9&oj=JRy6-uSD*%)lB zG++1|OmY>?xtY(b$0pK5=_O2EWDS(OH^vA)mej8sc)QlwI&ftN5fPt_D>pb+?Wh@H z^Uq2(WG)h~JxhQr8|LdSQ&T@mU!G3$UBCF_-z*imT{@c{qv&<>_OCx#$|zbwpR1(o zT!k48>wW>?aG`;iD?@Y&!?;jw zT{za}T8AHD@brQ#AoTLz(niyQBbC44ESSC`ubP=2KVE44IPZ$=-bQSEN?io@>Y-wO=OYU?*s9ru*0GeyN66=bMpLf?7$RJ6W~2oI$f&;tD;H6e%O z7%krNs=zmZRrrsWw*DSjL11AQF8_U<|KqTczmda);`RTL*8d3}!y)!JONbl}1|GrA zz|9;qN*`@~=dDYaa)UrV5~%9bc)HIM$L@;6XdlpGTuK(LqmR_=Y=Nnl1`mgfM*teG zT3c;va5g<7epL`W-oZC^Tp3NFvE$92F73?p{JuzA2g>~ z^$_6wl7Slx7{?#@xux5+d-(a6;fth=qw(brawMkD2av#2y+~@U$EpMa>4Y^EZV{D2 z;FtiU7^?_-oFTm|c(fT43nomSkf(@g_UeCIQ2}0eB zj|W`>tL^Wx04)%)n<;-g2HD=T9u)7I2`H3EF`urh>&{mNZ6qfDv-3?miU8SWuuMU4 z15$$00Ezk-SbklC?3)8pmQFyH1poI&x`ktEj&Z3|>!7+tKzlinbhsY$KOqXil@gUg zjQ)^kP$|x*LKA>W71e*H2ADY!s$%){O@!R7h@&=|LAl4O3V0D-Of;dw5vb_Wbz<*9 z6Jd2@Xi277DyV-PPW&W07lh4m{Wsl*fxaTQ%lE%$!Yl-1ha>InMmFp$Kq$$pt?$DZ zbJIlJRv_c=tG_L>_9cOlXus+TF$rc=7% z`8S#Y2~Ru*S{YZ^h<08Qe6hHDqx{?yNS&z|HuCV{LpA8V`BbuFC4q2cFsLtIUO-s4 zAaBQOrD&^!7tf8S%T9bSf$iu5#69W#Vkx@By|8nBZW2u{1X(#IN#$|w1x=-nsAxco z4?=-g@;L>7^%3Fyi94R`7T}w1s>Da^bxQgK?ni(vMoAo4XXqRZD?%HfZo)~0xxr&L zhj>!Ipun%|c(%d4a{=vzvwAMpCok=>y;r{hubgWtu&spEJX`5A_q` z5!rNQ{>wwn^z@%QoUIu3bN7@7ZBQ0)D6Sg^oogj3>0r(bPR<^p@ffdAz%^}Ivm~5t z{l%B5r}G+4JC<`LrKZ-t2b`(kJhw40bdtpw&{2DnOSl$~*&E0PD~zZxBM+<@;iR81 zFt8@$9I~AIw48t@ihzlj5Pd9B%e?f;Vsijrdw}HHLpZ!Vvd#|?hQLe?u+$F&rI-f? z{SC-~wV3j-x?6_oKGABy=HgRd3`YMsP>)PN`>_O|^W)eBD$ftQN-DoO`qFw3Pn1r7 zG4%*C@%AqO$xR91No;g>$fxrzjTV6zRRq?p%=7z+Xi1Ywsx*w_KFQ>0@BDD+1>n3- z))H=dcp0Z$?P-AJ`+8Afbp%0}JbK(M#&v=$Z@v1q${Qe1N0mFycwQv<{x?Az57dJG z#g+wj%qYM)_wm3xp+zdO#-a9WE{^}v6Ap$r*y#zurjW-$!!ySG0$@9ZIPz7tgT!jy z^hl@Ye&lyK@VslSg=}z9$7s?#g!lr-hhB)2OH3oproCZdc&06USNmchvl7D*{kqQ z3arN-26|Tyn8C&VWj%lTiAaqkD64QWB``?#_yyP2Jm{nkD&N3J>jQ z;!dxWGm|S@53O_KhyWlQ@G&dwpFXo3TWkESUwF}Rt@SkaYxdG1WgY3$o zdnT+M?qH+ua`~=JVzz`+OVF;fm|38vSd4YXZ7`d8tD&KR`0cZRZaf)>Pb42b5*m&N z7*ddc>BUCt1HH9lk8d+H2Re*+fZ3Rs^5>@`WAS(r=p?d*Q6GX7kM7%h#=x-qjtM_O zuO$46{;2KtD?SO&g{quo+qOF|?Y?822n8n&=wgl5lb~gYi{d7{ZlpJ&lMGUZE$5{p%Z+qb$U?}c92Djdri$`$lNhcJrn|f#Jm^J! z%4FIM_WJ{(+9-Td9|}t3_j&9n8+LM73~PjGZ z&n$yC1)~ejkL60QJJF6WCX0qBAw-3uVU`XvsR=EJA>Tf3T*sg?5$GsUm-zh*jg)Eq zX$40i(|W_0#p#yHl`%w;@Kv^&T=@oK9N647a0AvcS1KQApyeqS+N@%izOh8sY)AS zfVYR_i3I~0WPa_6Oft7*c4|Dtb@tQbMDD}84*xPIBXBfCqpD-0CK=A@5@Muvel&LG8 zfF|>_(M8H6TtELpL^+Mj@FBC4<_}|v5Oj^^R3{d%1Ka9ZJA^*$sHMxuxq_K>rH;ehaUgBejWoA11K4kv~g;1`%{B9JkAVBXr=`i`wd z)G9fu;JZ;F$6Dhd*3qAzhjgN|y2!kLohvsAm!dqu%;HyfEI8Do zfXPNwS6AC5*AmY5W|RyY|f(}QCUx~eAtJr@T4ep4ahBt4kj@GfjsN^ES2wvJ9X z(GKPSnfKt4lLDaoaJUOo{?`zV+UqL?8yZs;eto@TjwpGL^_gk?e2&V)-qv<@y9or{ z9XQt8LQa+ZZecJaI||{&qwAXkOgmC+sHEY7am#eh-RJhkf zuD|(-6bQoCZm7gF&WJF~C|mBp%76_N21n)o@olD@d%p$+w7mS&cGCB)gx}x{43Ism zw;UGla`-^RYY+UXnv0)CPDBM$~Q z3{fYQaH5JM!AL{r_WhNDDq$F2GtKhCp6-Kf)&s7hG`Xl?wO4!?K+^L^=TDUi=aYO= zzXtFWYvlHdlL2TF<_>{@1peQ1U63k{8$dFNJ!JNW@4%kSGRcatAI~~SFoMbV=x;RY zzfRHn4S@Z(88(dY>%R=z`@i4f5G03zYoPIwDJ2#n^RJJ2$3aO$=pei5M$Ayoa~0fV zoE2f6lhjQ>FHubeFZ|3>p3BrzR4+M(qi zF?S)2oC)76@@c>pc6xR#C8NclXay-4bp#tCm-y)${E3S#a~+nEb%cy#-2T`gG9HPb zLqig9zYnBn-vqc#iF4!c?*(|dVdD(R+X>A7R4wmzBp<|_H_}f{vu0;vk|h(SWUjju zy_y<{n!j-I7*g|8Nv>_8xX5UO2Guv8!+8Vjcpjx zk_=3n!UTZFMC6bSC$v3SjZTJ+5Z)mMVd==Z4Qgv@YVLlFvopew0t9#80|yPYw(D%^ zm9b2pxn7@UGBFE~og}|$(D@j+V>0JU=w4{3FgJ&Jd^MsZ(uB{*mO{q4;j}O;PJQQe zvd`33n8yb<1#{BvU;IisD${kvMhD)Jb4Gx~mbQC1-PX}N(K82$$irB%=qeXOum&H& zjB7(p - - - - - - - 2023-03-15T18:14:16.784620 - image/svg+xml - - - Matplotlib v3.7.0, https://matplotlib.org/ - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - diff --git a/notebooks/rmax-n.py b/notebooks/rmax-n.py new file mode 100644 index 0000000..b7b95ad --- /dev/null +++ b/notebooks/rmax-n.py @@ -0,0 +1,30 @@ +import numpy as np +import matplotlib.pyplot as plt +import scipy + +from notebooks.lib import read_load + + +def linear(x, a, b): + return x * a + b + + +alpha = read_load("../data/alpha") +meaned_by_N = alpha.groupby('N').agg({'fd': ['mean', 'std']}) \ + .reset_index() \ + .replace([np.inf, -np.inf], np.nan) + +without_prefix = alpha[alpha.N > 50] + +p, pcov = scipy.optimize.curve_fit(linear, np.log(without_prefix.cr), np.log(without_prefix.N)) +linear_extent = np.linspace(0, np.max(np.log(alpha.cr))) + +plt.scatter(np.log(alpha.cr), np.log(alpha.N), s=1, marker='.', color="tab:blue") +plt.plot(linear_extent, linear(linear_extent, *p), color="tab:red") + +plt.xlabel("$\\log r_{max}$") +plt.ylabel("$\\log N$") + +plt.savefig('../figures/rmax-n.svg') +plt.savefig('../figures/rmax-n.png') +plt.show() diff --git a/notebooks/sp-fd.py b/notebooks/sp-fd.py new file mode 100644 index 0000000..12f2907 --- /dev/null +++ b/notebooks/sp-fd.py @@ -0,0 +1,30 @@ +import numpy as np +from matplotlib import pyplot as plt + +from notebooks.lib import read_sp + +c_sp = read_sp("../data/stick-probability") +by_run = c_sp.groupby(['probability', 'N']) + +by_probability = by_run.agg( + overall_fd=('fd', lambda fd: np.mean(fd[-100:])), + overall_fd_std=('fd', 'std') +).reset_index().groupby('probability') + +ggg = by_probability.agg( + fd=('overall_fd', 'mean'), + # TODO Check stats + fd_std=('overall_fd_std', lambda std: np.sqrt(np.mean(np.square(std)))) +) + +plt.fill_between(ggg.index, ggg.fd - ggg.fd_std, ggg.fd + ggg.fd_std, alpha=0.2, label=f"Standard error band") +plt.plot(ggg.index, ggg.fd, color='tab:blue', label='fd mean, seeds = 100') +plt.plot(ggg.index, ggg.fd) + +plt.xlabel("$p_{stick}$") +plt.ylabel("$fd$") +plt.legend() + +plt.savefig('../figures/sp-fd.svg') +plt.savefig('../figures/sp-fd.png') +plt.show()