(* ::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]