diff options
Diffstat (limited to 'figs/fig-susmag.gplot')
| -rw-r--r-- | figs/fig-susmag.gplot | 133 |
1 files changed, 0 insertions, 133 deletions
diff --git a/figs/fig-susmag.gplot b/figs/fig-susmag.gplot deleted file mode 100644 index 901576f..0000000 --- a/figs/fig-susmag.gplot +++ /dev/null @@ -1,133 +0,0 @@ - -set terminal pslatex rotate size 3.417,4.222 -set multiplot - -cc1 = "#5e81b5" -cc2 = "#e19c24" -cc3 = "#8fb032" -cc4 = "#eb6235" - -Tc = 2 / log(1 + sqrt(2)) -Delta = 15. / 8 -gamma = 7. / 4 -beta = 1. / 8 - -Ch = 0.838677624411 - -t(T) = abs((T - Tc) / T) -h(T, H) = H / T -X(T, H) = h(T, H) * t(T)**(-Delta) - -poly(A, l, m, X) = A * X**m / (1 + (l * X)**(m + 1)) -polyint(A, l, m, X) = A * l**(-(m + 1)) * log(1 + (l * X)**(m + 1)) / (m + 1) - -G(i) = i == 1 ? -1.3578383417066 : \ - i == 2 ? -0.048953289720 : \ - i == 3 ? 0.038863932 : \ - i == 4 ? -0.068362119 : \ - i == 5 ? 0.18388370 : \ - i == 6 ? -0.6591714 : \ - i == 7 ? 2.937665 : \ - i == 8 ? -15.61 : 0 -GC(i) = G(i) * (2 * asinh(1))**2 * (Ch * (2 * asinh(1))**(-Delta))**i - -M0 = (2**2.5 * asinh(1))**0.125 -B = Tc**2 * M0 / (16 * pi) -C0 = 0.0255369719 -A = pi / 2 * C0 / (B**2) - -n = 1 -c(i) = i == 1 ? 0.00856277 : 0 -lamb = 10.487 -A2 = 1.70034 -B2 = B - -n2 = 5 -c2(i) = i == 1 ? 0.00639522 : \ - i == 2 ? 0.00143461 : \ - i == 3 ? -0.206879 : \ - i == 4 ? 4.39278 : \ - i == 5 ? -56.1389 : 0 -lamb2 = 13.2489 -A3 = 1.91747 -B3 = B - -susfunc = "figs/fig-sus_scaling-func.dat" -susfuncm = "figs/fig-sus_scaling-func-mod.dat" -magfunc = "figs/fig-mag_scaling-func.dat" -magfuncm = "figs/fig-mag_scaling-func-mod.dat" - -num = "data/data_square-4096.dat" - -set samples 20000 -set key off - -set size 1,1 - 0.52 -set origin 0,0.52 -set xrange [-40:1200] -set yrange [-1:33] -set ylabel offset 1,0 '$T\chi|t|^{\gamma}\times 10^3$' -set lmargin 6 -set xtics format '' -set mxtics 5 -set mytics 5 -set bmargin 0.2 - -plot \ - num using (X($2, $3)):(10**3 * Tc * $10 * t($2)**gamma):(10**3 * Tc * $11 * t($2)**gamma) with yerrorbars pt 0 lc black, \ - susfunc using (10**$1 / B):(10**(3+$2) * A * B**2) with linespoints pt 0 lw 3 lc rgb cc1 dt 1, \ - susfuncm using 1:(10**3 * $2) with lines lw 3 lc rgb cc2 dt 2, \ - -10**3 * (sum[i=1:7] GC(i + 1) * i * (i + 1) * (x)**(i-1)) with lines lw 3 lc rgb cc3 dt 3 - -set bmargin -1 -set tmargin 0.2 -set size 1,0.52 -set origin 0,0 -set yrange [1.18:1.74] -set ylabel offset 1,0 '$M|t|^{-\beta}$' -set xlabel '$h|t|^{-\beta\delta}$' -set xtics format '%g' - -plot \ - num using (X($2, $3)):($6 * t($2)**(-beta)):($7 * t($2)**(-beta)) with yerrorbars pt 0 lc black, \ - magfunc using (10**$1 / B):(M0 + 10**($2) * A * B) with linespoints pt 0 lw 3 lc rgb cc1 dt 1, \ - magfuncm using 1:2 with lines lw 3 lc rgb cc2 dt 2, \ - -sum[i=1:8] GC(i) * i * x**(i-1) with lines lw 3 lc rgb cc3 dt 3 - -set logscale xy -set tmargin -1 -set lmargin -1 -set size 0.7,0.325 -set origin 0.26,0.5 + 0.29 / 2 -set xrange [0.0015:1900] -set yrange [0.00002:0.08] -set ylabel offset 4,0 '\footnotesize$T\chi|t|^\gamma$' -set xlabel offset 0,0.5 '\footnotesize$h|t|^{-\beta\delta}$' -set mxtics 5 -#set xtics format '' -2,10,1000 -set xtics add ('$\footnotesize10^{-2}$' 10**(-2), "" 0.1, '$\footnotesize10^0$' 1, "" 10, '$\footnotesize10^2$' 100, "" 1000) -set mytics 5 -#set ytics format '\footnotesize$10^{%T}$' 0.00001,10,0.01 -set ytics add ('$\footnotesize10^{-4}$' 10**(-4), '$\footnotesize10^{-3}$' 10**(-3), '$\footnotesize10^{-2}$' 10**(-2)) - -plot \ - num using (X($2, $3)):(Tc * $10 * t($2)**gamma):(Tc * $11 * t($2)**gamma) with yerrorbars pt 0 lc black, \ - susfunc using (10**$1 / B):(10**$2 * A * B**2) with linespoints pt 0 lw 3 lc rgb cc1 dt 1, \ - susfuncm using 1:2 with lines lw 3 lc rgb cc2 dt 2, \ - -sum[i=2:8] GC(i) * i * (i-1) * x**(i-2) with lines dt 3 lw 3 lc rgb cc3 - -unset logscale xy -set logscale x -set size 0.65,0.325 -set origin 0.31,0.23 / 2 -set yrange [1.18:1.74] -set ylabel offset 4,0 '\footnotesize$M|t|^{-\beta}$' -set ytics format '\footnotesize {%g}' 1.2,0.2,1.8 -set mytics 5 - -plot \ - num using (X($2, $3)):($6 * t($2)**(-beta)):($7 * t($2)**(-beta)) with yerrorbars pt 0 lc black, \ - magfunc using (10**$1 / B):(M0 + A * B * 10**$2) with linespoints pt 0 lw 3 lc rgb cc1 dt 1, \ - magfuncm using 1:2 with lines lw 3 lc rgb cc2 dt 2, \ - -sum[i=1:8] GC(i) * i * x**(i-1) with lines lw 3 lc rgb cc3 dt 3 - |