3 files changed, 165 insertions, 15 deletions
diff --git a/figs/fig_correlation-times.gplot b/figs/fig_correlation-times.gplot
index fec61d4..7dac658 100644
--- a/figs/fig_correlation-times.gplot
+++ b/figs/fig_correlation-times.gplot
@@ -1,11 +1,55 @@
 
-#set terminal epslatex size 6.834, 2.81576
-#set output "fig_clusters_ising2d.tex"
+# data and constants, sorted by column
 
-set term pdf
-set output "test.pdf"
+isingRange = "[0.0000005:5000000]"
 
-set key autotitle columnhead
+dataIsing2D = "data/correlation-times/ising-2d.dat"
+betaIsing2D = 1.0 / 8.0
+nuIsing2D = 1.0
+deltaIsing2D = 15.0
+gammaIsing2D = 7.0 / 4.0
+zIsing2D = 0.231
+
+dataIsing3D = "data/correlation-times/ising-3d.dat"
+betaIsing3D = 0.326419
+nuIsing3D = 0.629971
+deltaIsing3D = 4.78984
+gammaIsing3D = 1.237075
+zIsing3D = 0.316
+
+pottsRange="[0.000001:5000000]"
+
+data3Potts2D = "data/correlation-times/potts-3.dat"
+beta3Potts2D = 1.0 / 9.0
+nu3Potts2D = 5.0 / 6.0
+delta3Potts2D = 14.0
+gamma3Potts2D = 13.0 / 9.0
+z3Potts2D = 0.568
+
+data4Potts2D = "data/correlation-times/potts-4.dat"
+beta4Potts2D = 1.0 / 12.0
+nu4Potts2D = 2.0 / 3.0
+delta4Potts2D = 15.0
+gamma4Potts2D = 7.0 / 6.0
+z4Potts2D = 0.953
+
+vectorRange="[0.000001:10000000]"
+
+data2Vector3D = "data/correlation-times/planar-pert.dat"
+beta2Vector3D = 0.3470
+nu2Vector3D = 0.6703
+delta2Vector3D = 4.79539
+gamma2Vector3D = 1.3169
+z2Vector3D = 0.39
+
+data3Vector3D = "data/correlation-times/heisenberg-pert.dat"
+beta3Vector3D = 0.3662
+nu3Vector3D = 0.7073
+delta3Vector3D = 4.79465
+gamma3Vector3D = 1.3895
+z3Vector3D = 0.295
+
+# define colors
 
 cc1 = "#5e81b5"
 cc2 = "#e19c24"
@@ -17,25 +61,73 @@ cc7 = "#5d9ec7"
 
 set palette model RGB defined ( 8 cc1, 16 cc2, 32 cc3, 64 cc4, 128 cc5, 256 cc6 , 512 cc7 )
 
+# global plotting settings
+
+set terminal epslatex size 6.834, 2.81576
+set output "fig_correlation-times.tex"
+
+set multiplot layout 2,3 margins 0.05, 0.995, 0.15, 0.99 spacing 0.05, 0.015
+set size ratio 1 / 1.61803398875
+set key autotitle columnhead
 unset colorbox
 set nokey
-set format x '$10^{%T}$'
 set logscale xy
-set ylabel '$\avg{s_{\text{\sc 1c}}}L^{-\gamma/\nu}$'
-set xlabel '$hL^{\beta\delta/\nu}$'
 
+scaledylabel = '"$\\avg{s_{\\text{\\sc 1c}}}L^{-\\gamma/\\nu}$" offset 2,0'
 
-data3Potts2D = "data/correlation-times/potts-3.dat"
-
+# plotting the first row
+set format x ""
 unset xlabel
+set ylabel @scaledylabel
+
+# ploting Ising 2D
+set xrange @isingRange
+set yrange [0.2:1.1]
+
+plot dataIsing2D using (column("H")/column("T") * exp(log(column("L")) * betaIsing2D * deltaIsing2D / nuIsing2D)):(column("t") * column("S") * exp(-log(column("L")) * (zIsing2D + 2))):(sqrt((column("dt")/column("t"))**2 + (column("dS")/column("S"))**2) * column("t") * column("S") * exp(-log(column("L")) * (zIsing2D + 2))):1 with yerrorbars pt 0 lw 2 palette,\
+  dataIsing2D using (column("H")/column("T") * exp(log(column("L")) * betaIsing2D * deltaIsing2D / nuIsing2D)):(column("t") * column("S") * exp(-log(column("L")) * (zIsing2D + 2))):1 with lines palette
+
 unset ylabel
 
-beta3Potts2D = 1.0 / 9.0
-nu3Potts2D = 5.0 / 6.0
-delta3Potts2D = 14.0
-gamma3Potts2D = 13.0 / 9.0
-z3Potts2D = 0.55
+# plotting 3 Potts 2D
+set xrange @pottsRange
+set yrange [0.05:2]
 
 plot data3Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta3Potts2D * delta3Potts2D / nu3Potts2D)):(column("t") * column("S") * exp(-log(column("L")) * (z3Potts2D + 2))):(sqrt((column("dt")/column("t"))**2 + (column("dS")/column("S"))**2) * column("t") * column("S") * exp(-log(column("L")) * (z3Potts2D + 2))):1 with yerrorbars pt 0 lw 2 palette,\
   data3Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta3Potts2D * delta3Potts2D / nu3Potts2D)):(column("t") * column("S") * exp(-log(column("L")) * (z3Potts2D + 2))):1 with lines palette
 
+# plotting 2 Vector 3D
+set xrange @vectorRange
+set yrange [0.4:2]
+
+plot data2Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta2Vector3D * delta2Vector3D / nu2Vector3D)):(column("t") * column("S") * exp(-log(column("L")) * (z2Vector3D + 3))):(sqrt((column("dt")/column("t"))**2 + (column("dS")/column("S"))**2) * column("t") * column("S") * exp(-log(column("L")) * (z2Vector3D + 3))):1 with yerrorbars pt 0 lw 2 palette,\
+  data2Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta2Vector3D * delta2Vector3D / nu2Vector3D)):(column("t") * column("S") * exp(-log(column("L")) * (z2Vector3D + 3))):1 with lines palette
+
+# plotting the second row
+set format x '$10^{%T}$'
+set xlabel '$hL^{\beta\delta/\nu}$'
+set ylabel @scaledylabel
+
+# plotting Ising 3D
+set xrange @isingRange
+set yrange [0.2:1]
+
+plot dataIsing3D using (column("H")/column("T") * exp(log(column("L")) * betaIsing3D * deltaIsing3D / nuIsing3D)):(column("t") * column("S") * exp(-log(column("L")) * (zIsing3D + 3))):(sqrt((column("dt")/column("t"))**2 + (column("dS")/column("S"))**2) * column("t") * column("S") * exp(-log(column("L")) * (zIsing3D + 3))):1 with yerrorbars pt 0 lw 2 palette,\
+  dataIsing3D using (column("H")/column("T") * exp(log(column("L")) * betaIsing3D * deltaIsing3D / nuIsing3D)):(column("t") * column("S") * exp(-log(column("L")) * (zIsing3D + 3))):1 with lines palette
+
+unset ylabel
+
+# plotting 4 Potts 2D
+set xrange @pottsRange
+set yrange [0.02:4]
+
+plot data4Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta4Potts2D * delta4Potts2D / nu4Potts2D)):(column("t") * column("S") * exp(-log(column("L")) * (z4Potts2D + 2))):(sqrt((column("dt")/column("t"))**2 + (column("dS")/column("S"))**2) * column("t") * column("S") * exp(-log(column("L")) * (z4Potts2D + 2))):1 with yerrorbars pt 0 lw 2 palette,\
+  data4Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta4Potts2D * delta4Potts2D / nu4Potts2D)):(column("t") * column("S") * exp(-log(column("L")) * (z4Potts2D + 2))):1 with lines palette
+
+# plotting 3 Vector 3D
+set xrange @vectorRange
+set yrange [0.9:4]
+
+plot data3Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta3Vector3D * delta3Vector3D / nu3Vector3D)):(column("t") * column("S") * exp(-log(column("L")) * (z3Vector3D + 3))):(sqrt((column("dt")/column("t"))**2 + (column("dS")/column("S"))**2) * column("t") * column("S") * exp(-log(column("L")) * (z3Vector3D + 3))):1 with yerrorbars pt 0 lw 2 palette,\
+  data3Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta3Vector3D * delta3Vector3D / nu3Vector3D)):(column("t") * column("S") * exp(-log(column("L")) * (z3Vector3D + 3))):1 with lines palette
+
diff --git a/figs/fig_generator-times.gplot b/figs/fig_generator-times.gplot
new file mode 100644
index 0000000..34a44fd
--- /dev/null
+++ b/figs/fig_generator-times.gplot
@@ -0,0 +1,21 @@
+
+set term epslatex size 3.417, 2.3
+set output "fig_generator-times.tex"
+
+set logscale xy
+
+dataPlanarUni = "data/correlation-times/planar-uni.dat"
+dataPlanarPert = "data/correlation-times/planar-pert.dat"
+
+set nokey
+set size ratio 1 / 1.61803398875
+
+set xlabel '$|H|$'
+set ylabel "$\\avg{s_{\\text{\\sc 1c}}}\\tau$"
+
+set format x '$10^{%T}$'
+
+plot dataPlanarUni using (column("H")):((column("L") == 32 ? column("t") : 1/0) * column("S") * exp(-log(column("L")) * (2))):(sqrt((column("dt")/column("t"))**2 + (column("dS")/column("S"))**2) * column("S") * exp(-log(column("L")) * (2))) with yerrorbars pt 0 lw 2 lc rgb "red",\
+  dataPlanarUni using (column("H")):((column("L") == 32 ? column("t") : 1/0) * column("S") * exp(-log(column("L")) * (2))) with lines lc rgb "red",\
+  dataPlanarPert using (column("H")):((column("L") == 32 ? column("t") : 1/0) * column("S") * exp(-log(column("L")) * (2))):(sqrt((column("dt")/column("t"))**2 + (column("dS")/column("S"))**2) * column("t") * column("S") * exp(-log(column("L")) * (2))) with yerrorbars pt 0 lw 2 lc rgb "green",\
+  dataPlanarPert using (column("H")):((column("L") == 32 ? column("t") : 1/0) * column("S") * exp(-log(column("L")) * (2))) with lines lc rgb "green"
diff --git a/figs/fig_harmonic-susceptibilities.gplot b/figs/fig_harmonic-susceptibilities.gplot
new file mode 100644
index 0000000..85d474f
--- /dev/null
+++ b/figs/fig_harmonic-susceptibilities.gplot
@@ -0,0 +1,37 @@
+
+set terminal epslatex size 3.417, 5.3
+set output "fig_harmonic-susceptibilities.tex"
+
+data4 = "data/harmonic-susceptibilities/order-4.dat"
+data6 = "data/harmonic-susceptibilities/order-6.dat"
+
+set palette defined (0 'blue', 0.33 'green', 0.66 'yellow', 1 'red')
+stats data6 using "H"
+set cbrange[STATS_min:STATS_max]
+set logscale xycb
+set nokey
+
+set xrange [1:5000]
+set yrange [0.01:10000]
+
+unset colorbox
+
+set multiplot layout 2,1 margins 0.15, 0.99, 0.225, 0.99 spacing 0,0.015
+set size ratio 1 / 1.61803398875
+
+unset xlabel
+set ylabel '$\chi$' offset 3
+set format x ''
+
+plot data4 using "L":"X":"dX":"H" with yerrorbars pt 0 lw 2 palette,\
+     data4 using "L":"X":"H"      with lines palette
+
+set xlabel '$L$' 
+set colorbox horizontal user origin 0.15, 0.1 size 0.84, 0.03
+set cblabel '(Top) $h_4$ or (bottom) $h_6$'
+set cbtics format '$10^{%T}$'
+set format x
+
+plot data6 using "L":"X":"dX":"H" with yerrorbars pt 0 lw 2 palette,\
+     data6 using "L":"X":"H"      with lines palette
+