From 1676b1bfc1544ba7049f198cb4ae75bea0513014 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Mon, 24 Sep 2018 15:27:05 -0400 Subject: many figure changes --- figs/fig_correlation-times.gplot | 122 +++++++++++++++++++++++++++---- figs/fig_generator-times.gplot | 21 ++++++ figs/fig_harmonic-susceptibilities.gplot | 37 ++++++++++ 3 files changed, 165 insertions(+), 15 deletions(-) create mode 100644 figs/fig_generator-times.gplot create mode 100644 figs/fig_harmonic-susceptibilities.gplot (limited to 'figs') 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 + -- cgit v1.2.3-54-g00ecf