AENC/switchchain Changeset - 843d893407af · Centrum Wiskunde & Informatica (CWI)

Changeset - 843d893407af

Parent rev.

Child rev.

[Not reviewed]

0 1 1

Tom Bannink - 8 years ago 2017-06-02 23:05:16
tombannink@gmail.com

Add time evolution plot

2 files changed with 38 insertions and 5 deletions:

plots/timeevol.pdf

bin+new file

triangle_analysis.m

0 comments (0 inline, 0 general)

plots/timeevol.pdf

➞

Show inline comments

		new file 100644
		binary diff not shown

triangle_analysis.m

➞

Show inline comments

@@ @@ -71,7 +71,7 @@ Needs["ErrorBarPlots`"] @@
 (*Data import*)
-gsraw=Import[NotebookDirectory[]<>"data/graphdata_temp2.m"];
+gsraw=Import[NotebookDirectory[]<>"data/graphdata_timeevol.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. *) *)
 (*Plot triangle count over "time" in Markov chain instances*)
 numPlots=10;
 selectedData=gdata[[1,-1]][[-numPlots;;-1]];
 numPlots=20;
 selectedData=gdata[[1,1]][[-numPlots;;-1]];
 measureSkip=1;
 minCount=Min[Map[Min[#[[2]]]&,selectedData]];
 maxCount=Max[Map[Max[#[[2]]]&,selectedData]];
 maxTime=Max[Map[Length[#[[2]]]&,selectedData]];
 maxTime=30000;
 skipPts=Max[1,Round[maxTime/400]]; (* Plotting every point is slow. Plot only once per `skipPts` timesteps *)
 (* maxTime=30000; *)
 skipPts=Max[1,Round[maxTime/500]]; (* Plotting every point is slow. Plot only once per `skipPts` timesteps *)
 coarseData=Map[#[[2,1;;maxTime;;skipPts]]&,selectedData];
 labels=Map["{n,tau} = "<>ToString[#[[1]]]&,selectedData];
 ListPlot[coarseData,Joined->True,PlotRange->{0*minCount,maxCount},DataRange->{0,measureSkip*maxTime},PlotLegends->labels]
 (* Map[ListPlot[#,Joined->True,PlotRange\[Rule]{minCount,maxCount},DataRange\[Rule]{0,maxTime}]&,coarseData] *)
 selectedData=gdata[[1,1]];
 measureSkip=1;
 minCount=Min[Map[Min[#[[2]]]&,selectedData]];
 maxCount=Max[Map[Max[#[[2]]]&,selectedData]];
 maxTime=Max[Map[Length[#[[2]]]&,selectedData]];
 maxTime=30000;
 skipPts=Max[1,Round[maxTime/100]]; (* Plotting every point is slow. Plot only once per `skipPts` timesteps *)
 coarseData=Map[#[[2,1;;maxTime;;skipPts]]&,selectedData];
 labels=Map["{n,tau} = "<>ToString[#[[1]]]&,selectedData];
 plot1=ListPlot[coarseData[[1;;5]],Joined->True,PlotRange->{0*minCount,maxCount},DataRange->{0,measureSkip*maxTime}]
 plot2=ListPlot[coarseData[[6;;10]],Joined->True,PlotRange->{0*minCount,maxCount},DataRange->{0,measureSkip*maxTime}]
 plot3=ListPlot[coarseData[[11;;15]],Joined->True,PlotRange->{0*minCount,maxCount},DataRange->{0,measureSkip*maxTime}]
 plot4=ListPlot[coarseData[[16;;20]],Joined->True,PlotRange->{0*minCount,maxCount},DataRange->{0,measureSkip*maxTime}]
 (* For export *)
 numPlots=20;
 selectedData=gdata[[2,1]][[-numPlots;;-1]];
 measureSkip=1;
 minCount=Min[Map[Min[#[[2]]]&,selectedData]];
 maxCount=Max[Map[Max[#[[2]]]&,selectedData]];
 maxTime=Max[Map[Length[#[[2]]]&,selectedData]];
 (* maxTime=30000; *)
 skipPts=Max[1,Round[maxTime/5000]]; (* Plotting every point is slow. Plot only once per `skipPts` timesteps *)
 coarseData=Map[#[[2,1;;maxTime;;skipPts]]&,selectedData];
 labels=Map["{n,tau} = "<>ToString[#[[1]]]&,selectedData];
 plotTimeEvol=ListPlot[coarseData,Joined->True,PlotRange->{0*minCount,maxCount},DataRange->{0,measureSkip*maxTime},Frame->True,FrameLabel->{"timesteps","number of triangles"},ImageSize->300]
 (* Map[ListPlot[#,Joined->True,PlotRange\[Rule]{minCount,maxCount},DataRange\[Rule]{0,maxTime}]&,coarseData] *)
 Export[NotebookDirectory[]<>"plots/timeevol.pdf",plotTimeEvol]
 (* ::Subsection:: *)
 (*Plot success rate over "time"*)

0 comments (0 inline, 0 general)