Files
@ 5cec3a409ef3
Branch filter:
Location: AENC/switchchain/triangle_spectrum_plots.m - annotation
5cec3a409ef3
732 B
application/vnd.wolfram.mathematica.package
Add triangle creation frequency plots
5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 5cec3a409ef3 | (* ::Package:: *)
gsraw=Import[NotebookDirectory[]<>"cpp/graphdata_spectrum.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. *) *)
gdata=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}];
(* Data format: *)
(* gdata[[ tau index, n index, run index , datatype index ]] *)
(* datatype index:
1: {n,tau}
2: etmt
*)
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]
|