Changeset - eba8261885e8
[Not reviewed]
master
0 1 0
Tom Bannink - 6 years ago 2019-02-27 16:38:05
tom.bannink@cwi.nl
Change trimeevol plot for thesis
1 file changed with 8 insertions and 3 deletions:
0 comments (0 inline, 0 general)
triangle_analysis.m
Show inline comments
 
@@ -4,6 +4,7 @@ Quit[]
 

	
 

	
 
Needs["ErrorBarPlots`"]
 
Needs["MaTeX`"]
 

	
 

	
 
(* ::Section:: *)
 
@@ -179,8 +180,9 @@ plot4=ListPlot[coarseData[[16;;20]],Joined->True,PlotRange->{0*minCount,maxCount
 

	
 

	
 
(* For export *)
 
numPlots=20;
 
selectedData=gdata[[2,1]][[-numPlots;;-1]];
 
numPlots=25;
 
selectedPlots={6,7,8,11,12,13,16,17,18,21,22,23};
 
selectedData=gdata[[2,1]][[selectedPlots]];
 
measureSkip=1;
 
minCount=Min[Map[Min[#[[2]]]&,selectedData]];
 
maxCount=Max[Map[Max[#[[2]]]&,selectedData]];
 
@@ -189,7 +191,10 @@ maxTime=Max[Map[Length[#[[2]]]&,selectedData]];
 
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]
 
plotTimeEvol=ListPlot[coarseData,Joined->True,PlotRange->{0*minCount,maxCount},DataRange->{0,measureSkip*maxTime},
 
Frame->True,FrameLabel->{MaTeX["\\text{timesteps}"],MaTeX["\\text{number of triangles}"]},
 
PlotLabel->MaTeX["n=1000,\\; \\tau = 2.2"],
 
ImageSize->250]
 
(* Map[ListPlot[#,Joined->True,PlotRange\[Rule]{minCount,maxCount},DataRange\[Rule]{0,maxTime}]&,coarseData] *)
 

	
 

	
0 comments (0 inline, 0 general)