From 9ef2b1560059b73a9b88628953c6390060d829d2 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Tue, 25 Sep 2018 00:04:10 -0400 Subject: many text and figure changes --- figs/fig_clusters_ising2d.gplot | 177 ++++++++++++++----------------- figs/fig_correlation-times.gplot | 2 +- figs/fig_harmonic-susceptibilities.gplot | 14 +-- 3 files changed, 90 insertions(+), 103 deletions(-) (limited to 'figs') diff --git a/figs/fig_clusters_ising2d.gplot b/figs/fig_clusters_ising2d.gplot index 61172db..6f1a75d 100644 --- a/figs/fig_clusters_ising2d.gplot +++ b/figs/fig_clusters_ising2d.gplot @@ -1,140 +1,127 @@ -set terminal epslatex size 6.834, 2.81576 -set output "fig_clusters_ising2d.tex" - -cc1 = "#5e81b5" -cc2 = "#e19c24" -cc3 = "#8fb032" -cc4 = "#eb6235" -cc5 = "#8778b3" -cc6 = "#c56e1a" - -set palette model RGB defined ( 8 cc1, 16 cc2, 32 cc3, 64 cc4, 128 cc5, 256 cc6 ) - -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}$' - -set xrange [0.000005:12000000] - - -set multiplot layout 2,3 +# data and constants, sorted by column -set size 0.38, 0.5 - -dataIsing2D = "data/cluster-size/cluster-size_ising2d.dat" - -set yrange [0.8:8] - -unset xlabel - -set label "2D Ising" at graph 0.1, 0.75 +isingRange = "[0.0000005:5000000]" +dataIsing2D = "data/correlation-times/ising-2d.dat" betaIsing2D = 1.0 / 8.0 nuIsing2D = 1.0 deltaIsing2D = 15.0 gammaIsing2D = 7.0 / 4.0 -plot dataIsing2D using ($3/$2 * exp(log($1) * betaIsing2D * deltaIsing2D / nuIsing2D)):($4 * exp(-log($1) * gammaIsing2D / nuIsing2D)):($5 * exp(-log($1) * gammaIsing2D/nuIsing2D)):1 with yerrorbars pt 0 lw 2 palette, dataIsing2D using ($3/$2 * exp(log($1) * betaIsing2D * deltaIsing2D / nuIsing2D)):($4 * exp(-log($1) * gammaIsing2D / nuIsing2D)):1 with lines palette, 1.5 * exp(log(x) * 2 / deltaIsing2D) with lines lc black - -unset label -set origin 0.35, 0.5 -set size 0.35, 0.5 - -data3Potts2D = "data/cluster-size/cluster-size_3potts2d.dat" - -set xrange [0.000005:9000000] -set yrange [0.4:7] - -set label "2D 3-State Potts" at graph 0.1, 0.75 +dataIsing3D = "data/correlation-times/ising-3d.dat" +betaIsing3D = 0.326419 +nuIsing3D = 0.629971 +deltaIsing3D = 4.78984 +gammaIsing3D = 1.237075 -unset xlabel -unset ylabel +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 -plot data3Potts2D using ($3/$2 * exp(log($1) * beta3Potts2D * delta3Potts2D / nu3Potts2D)):($4 * exp(-log($1) * gamma3Potts2D / nu3Potts2D)):($5 * exp(-log($1) * gamma3Potts2D/nu3Potts2D)):1 with yerrorbars pt 0 lw 2 palette, data3Potts2D using ($3/$2 * exp(log($1) * beta3Potts2D * delta3Potts2D / nu3Potts2D)):($4 * exp(-log($1) * gamma3Potts2D / nu3Potts2D)):1 with lines palette, 1.2 * exp(log(x) * 2 / delta3Potts2D) with lines lc black - -unset label -set size 0.38, 0.5 - -data2Vector3D = "data/cluster-size/cluster-size_2vector3d.dat" - -set xrange [0.000005:9000000] -set yrange [0.6:200] +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 -set label "3D O(2)" at graph 0.1, 0.75 +vectorRange="[0.000001:10000000]" +data2Vector3D = "data/correlation-times/planar-pert.dat" beta2Vector3D = 0.3470 nu2Vector3D = 0.6703 delta2Vector3D = 4.79539 gamma2Vector3D = 1.3169 -plot data2Vector3D using ($3/$2 * exp(log($1) * beta2Vector3D * delta2Vector3D / nu2Vector3D)):($4 * exp(-log($1) * gamma2Vector3D / nu2Vector3D)):($5 * exp(-log($1) * gamma2Vector3D/nu2Vector3D)):1 with yerrorbars pt 0 lw 2 palette, data2Vector3D using ($3/$2 * exp(log($1) * beta2Vector3D * delta2Vector3D / nu2Vector3D)):($4 * exp(-log($1) * gamma2Vector3D / nu2Vector3D)):1 with lines palette, 3 * exp(log(x) * 2 / delta2Vector3D) with lines lc black +data3Vector3D = "data/correlation-times/heisenberg-pert.dat" +beta3Vector3D = 0.3662 +nu3Vector3D = 0.7073 +delta3Vector3D = 4.79465 +gamma3Vector3D = 1.3895 -unset label +# define colors -set origin -0.05,0.0 -set size 0.43, 0.55 +cc1 = "#5e81b5" +cc2 = "#e19c24" +cc3 = "#8fb032" +cc4 = "#eb6235" +cc5 = "#8778b3" +cc6 = "#c56e1a" +cc7 = "#5d9ec7" -dataIsing3D = "data/cluster-size/cluster-size_ising3d.dat" +set palette model RGB defined ( 8 cc1, 16 cc2, 32 cc3, 64 cc4, 128 cc5, 256 cc6 , 512 cc7 ) -set xrange [0.000005:9000000] -set yrange [0.9:400] +# global plotting settings -set label "3D Ising" at graph 0.1, 0.75 +set terminal epslatex size 6.834, 2.81576 +set output "fig_clusters_ising2d.tex" -set ylabel '$\avg{s_{\text{\sc 1c}}}L^{-\gamma/\nu}$' -set xlabel '$hL^{\beta\delta/\nu}$' +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 logscale xy -betaIsing3D = 0.326419 -nuIsing3D = 0.629971 -deltaIsing3D = 4.78984 -gammaIsing3D = 1.237075 +scaledylabel = '"$\\avg{s_{\\text{\\sc 1c}}}L^{-\\gamma/\\nu}$" offset 2,0' -plot dataIsing3D using ($3/$2 * exp(log($1) * betaIsing3D * deltaIsing3D / nuIsing3D)):($4 * exp(-log($1) * gammaIsing3D / nuIsing3D)):($5 * exp(-log($1) * gammaIsing3D/nuIsing3D)):1 with yerrorbars pt 0 lw 2 palette, dataIsing3D using ($3/$2 * exp(log($1) * betaIsing3D * deltaIsing3D / nuIsing3D)):($4 * exp(-log($1) * gammaIsing3D / nuIsing3D)):1 with lines palette, 5 * exp(log(x) * 2 / deltaIsing3D) with lines lc black +# plotting the first row +set format x "" +unset xlabel +set ylabel @scaledylabel +# ploting Ising 2D +set xrange @isingRange +#set yrange [0.8:10] + +plot dataIsing2D using (column("H")/column("T") * exp(log(column("L")) * betaIsing2D * deltaIsing2D / nuIsing2D)):(column("S") * exp(-log(column("L")) * (gammaIsing2D / nuIsing2D))):(column("dS") * exp(-log(column("L")) * (gammaIsing2D / nuIsing2D))):1 with yerrorbars pt 0 lw 2 palette,\ + dataIsing2D using (column("H")/column("T") * exp(log(column("L")) * betaIsing2D * deltaIsing2D / nuIsing2D)):(column("S") * exp(-log(column("L")) * (gammaIsing2D / nuIsing2D))):1 with lines palette -unset label unset ylabel -set origin 0.35, 0 -set size 0.35, 0.55 -data4Potts2D = "data/cluster-size/cluster-size_4potts2d.dat" +# plotting 3 Potts 2D +set xrange @pottsRange +set yrange [0.5:10] -set xrange [0.000005:9000000] -set yrange [0.6:4] +plot data3Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta3Potts2D * delta3Potts2D / nu3Potts2D)):(column("S") * exp(-log(column("L")) * (gamma3Potts2D / nu3Potts2D))):(column("dS") * exp(-log(column("L")) * (gamma3Potts2D / nu3Potts2D))):1 with yerrorbars pt 0 lw 2 palette,\ + data3Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta3Potts2D * delta3Potts2D / nu3Potts2D)):(column("S") * exp(-log(column("L")) * (gamma3Potts2D / nu3Potts2D))):1 with lines palette -set label "2D 4-State Potts" at graph 0.1, 0.75 +# plotting 2 Vector 3D +set xrange @vectorRange +set yrange [0.4:250] -beta4Potts2D = 1.0 / 12.0 -nu4Potts2D = 2.0 / 3.0 -delta4Potts2D = 15.0 -gamma4Potts2D = 7.0 / 6.0 +plot data2Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta2Vector3D * delta2Vector3D / nu2Vector3D)):(column("S") * exp(-log(column("L")) * (gamma2Vector3D / nu2Vector3D))):(column("dS") * exp(-log(column("L")) * (gamma2Vector3D / nu2Vector3D))):1 with yerrorbars pt 0 lw 2 palette,\ + data2Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta2Vector3D * delta2Vector3D / nu2Vector3D)):(column("S") * exp(-log(column("L")) * (gamma2Vector3D / nu2Vector3D))):1 with lines palette -plot data4Potts2D using ($3/$2 * exp(log($1) * beta4Potts2D * delta4Potts2D / nu4Potts2D)):($4 * exp(-log($1) * gamma4Potts2D / nu4Potts2D)):($5 * exp(-log($1) * gamma4Potts2D/nu4Potts2D)):1 with yerrorbars pt 0 lw 2 palette, data4Potts2D using ($3/$2 * exp(log($1) * beta4Potts2D * delta4Potts2D / nu4Potts2D)):($4 * exp(-log($1) * gamma4Potts2D / nu4Potts2D)):1 with lines palette, 1.2 * exp(log(x) * 2 / delta4Potts2D) with lines lc black +# plotting the second row +set format x '$10^{%T}$' +set xlabel '$hL^{\beta\delta/\nu}$' +set ylabel @scaledylabel -unset label -set size 0.38, 0.55 +# plotting Ising 3D +set xrange @isingRange +set yrange [0.5:500] -data3Vector3D = "data/cluster-size/cluster-size_3vector3d.dat" +plot dataIsing3D using (column("H")/column("T") * exp(log(column("L")) * betaIsing3D * deltaIsing3D / nuIsing3D)):(column("S") * exp(-log(column("L")) * (gammaIsing3D / nuIsing3D))):(column("dS") * exp(-log(column("L")) * (gammaIsing3D / nuIsing3D))):1 with yerrorbars pt 0 lw 2 palette,\ + dataIsing3D using (column("H")/column("T") * exp(log(column("L")) * betaIsing3D * deltaIsing3D / nuIsing3D)):(column("S") * exp(-log(column("L")) * (gammaIsing3D / nuIsing3D))):1 with lines palette -set xrange [0.000005:9000000] -set yrange [0.6:200] +unset ylabel -set label "3D O(3)" at graph 0.1, 0.75 +# plotting 4 Potts 2D +set xrange @pottsRange +set yrange [0.8:4] -beta3Vector3D = 0.3662 -nu3Vector3D = 0.7073 -delta3Vector3D = 4.79465 -gamma3Vector3D = 1.3895 +plot data4Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta4Potts2D * delta4Potts2D / nu4Potts2D)):(column("S") * exp(-log(column("L")) * (gamma4Potts2D / nu4Potts2D))):(column("dS") * exp(-log(column("L")) * (gamma4Potts2D / nu4Potts2D))):1 with yerrorbars pt 0 lw 2 palette,\ + data4Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta4Potts2D * delta4Potts2D / nu4Potts2D)):(column("S") * exp(-log(column("L")) * (gamma4Potts2D / nu4Potts2D))):1 with lines palette + +# plotting 3 Vector 3D +set xrange @vectorRange +set yrange [0.1:100] -plot data3Vector3D using ($3/$2 * exp(log($1) * beta3Vector3D * delta3Vector3D / nu3Vector3D)):($4 * exp(-log($1) * gamma3Vector3D / nu3Vector3D)):($5 * exp(-log($1) * gamma3Vector3D/nu3Vector3D)):1 with yerrorbars pt 0 lw 2 palette, data3Vector3D using ($3/$2 * exp(log($1) * beta3Vector3D * delta3Vector3D / nu3Vector3D)):($4 * exp(-log($1) * gamma3Vector3D / nu3Vector3D)):1 with lines palette, 2 * exp(log(x) * 2 / delta3Vector3D) with lines lc black +plot data3Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta3Vector3D * delta3Vector3D / nu3Vector3D)):(column("S") * exp(-log(column("L")) * (gamma3Vector3D / nu3Vector3D))):(column("dS") * exp(-log(column("L")) * (gamma3Vector3D / nu3Vector3D))):1 with yerrorbars pt 0 lw 2 palette,\ + data3Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta3Vector3D * delta3Vector3D / nu3Vector3D)):(column("S") * exp(-log(column("L")) * (gamma3Vector3D / nu3Vector3D))):1 with lines palette -unset multiplot diff --git a/figs/fig_correlation-times.gplot b/figs/fig_correlation-times.gplot index 7dac658..db62a2d 100644 --- a/figs/fig_correlation-times.gplot +++ b/figs/fig_correlation-times.gplot @@ -73,7 +73,7 @@ unset colorbox set nokey set logscale xy -scaledylabel = '"$\\avg{s_{\\text{\\sc 1c}}}L^{-\\gamma/\\nu}$" offset 2,0' +scaledylabel = '"$\\tau\\avg{s_{\\text{\\sc 1c}}}L^{-z}$" offset 2,0' # plotting the first row set format x "" diff --git a/figs/fig_harmonic-susceptibilities.gplot b/figs/fig_harmonic-susceptibilities.gplot index 85d474f..6129723 100644 --- a/figs/fig_harmonic-susceptibilities.gplot +++ b/figs/fig_harmonic-susceptibilities.gplot @@ -1,5 +1,5 @@ -set terminal epslatex size 3.417, 5.3 +set terminal epslatex size 3.417, 4.8 set output "fig_harmonic-susceptibilities.tex" data4 = "data/harmonic-susceptibilities/order-4.dat" @@ -11,12 +11,12 @@ set cbrange[STATS_min:STATS_max] set logscale xycb set nokey -set xrange [1:5000] -set yrange [0.01:10000] +set xrange [3:5000] +set yrange [0.005:20000] unset colorbox -set multiplot layout 2,1 margins 0.15, 0.99, 0.225, 0.99 spacing 0,0.015 +set multiplot layout 2,1 margins 0.15, 0.99, 0.1, 0.88 spacing 0,0 set size ratio 1 / 1.61803398875 unset xlabel @@ -27,9 +27,9 @@ 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 colorbox horizontal user origin 0.15, 0.89 size 0.84, 0.03 +set cblabel '(Top) $h_4$ or (bottom) $h_6$' offset 0,6.5 +set cbtics format '$10^{%T}$' offset 0,2.8 set format x plot data6 using "L":"X":"dX":"H" with yerrorbars pt 0 lw 2 palette,\ -- cgit v1.2.3-54-g00ecf