(* ::Package:: *)

gsraw=Import[NotebookDirectory[]<>"data/graphdata_canonical_timeevol.m"];
gdata=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}];
(* Data format: *)
(* gdata[[ tau index, n index, run index, datatype index ]] *)
(* datatype index:
1: {n,tau}
2: { triangles }
*)


tauChoice=3;
runs=1;;5;
Map[
ListPlot[#[[1,runs,2,1;;3000]],
Joined->True,
PlotRange->Automatic,
Frame->True,
FrameLabel->{"timesteps","number of triangles"},
ImageSize->300
]&,gdata]
