diff --git a/triangle_spectrum_plots.m b/triangle_spectrum_plots.m index 2fbbd2afd660a1d7fb4b82a8b5c7240057f93cf4..e525f1aa7a49ff7c4cda6e15d7c26d53a00a96b5 100644 --- a/triangle_spectrum_plots.m +++ b/triangle_spectrum_plots.m @@ -1,6 +1,6 @@ (* ::Package:: *) -gsraw=Import[NotebookDirectory[]<>"cpp/graphdata_spectrum.m"]; +gsraw=Import[NotebookDirectory[]<>"data/graphdata_properties.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. *) *) @@ -9,19 +9,20 @@ gdata=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}]; (* gdata[[ tau index, n index, run index , datatype index ]] *) (* datatype index: 1: {n,tau} -2: etmt +2: avgTriangles +3: edges +4: dstn +5: { HH A, HH L, average A, average L } where for each there is (average of) {lambda1 , lambda1 - lambda2, lambda1/lambda2} +6: switching successrate after mixing +7: initial HH triangles *) tauvalues=gdata[[All,1,1,1,2]]; nlabels=Map["n = "<>ToString[#]&,gdata[[1,All,1,1,1]]]; taulabels=Map["\[Tau] = "<>ToString[#]&,gdata[[All,1,1,1,2]]]; -ListPlot[gdata[[1,1,1,{4,6}]]] -ListPlot[gdata[[1,1,1,{5,7}]]] - - -Histogram[gdata[[1,1,1,{4,6}]],50] -Histogram[gdata[[1,1,1,{5,7}]],50] +Histogram[gdata[[1,1,All,5,3,1]]] +Histogram[gdata[[1,1,All,5,3,2]]]