SU(2) Simulations Summary

Author

Laurence Sebastian Bowes

Published

November 29, 2022

Code
proj_path = "/home/laurence/projects/RipRep"
include(proj_path*"/src/spectroscopy.jl");
include(proj_path*"/src/utilities.jl");
include(proj_path*"/src/parser.jl");
include(proj_path*"/src/wilson.jl");
theme(:dark)
#plotlyjs()
default(size = (800, 500))

\beta = 2.1

m = -0.08 L16T32

Code
e = load_ensemble("/home/laurence/projects/RipRep/data/Prepared/b2.1_csw3.0000_m-0.08_L16T32_expclv/out_0")
w = load_wilson_flow("/home/laurence/projects/RipRep/data/Prepared/b2.1_csw3.0000_m-0.08_L16T32_expclv/out_wf")
default_thermalise_bin!(e)
default_thermalise_bin!(w)
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Plaquette

Code
display(plot_plaquette(e))

t^2E(t)

Code
display(plot_t2e(w, title=_ensemble_to_latex_string(e)))

W(t)

Code
display(plot_w(w, title=_ensemble_to_latex_string(e)))

Topological Charge

Code
display(plot_tc(w, title=_ensemble_to_latex_string(e)))
display(plot_tc_hist(w, title=_ensemble_to_latex_string(e)))

dS History

Code
try
default_thermalise_bin!(e)
display(plot(1:nrow(e.analysis), e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = "\$dS\$", ylabel = "\$dS\$", xlabel = "Trajectory #"))
display(histogram(e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = 
"\$dS\$", xlabel = "\$dS\$"))
catch error
end
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Polyakov

0

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = "Polyakov 0 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = 
"Polyakov 0 Re"))

1

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = "Polyakov 1 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = 
"Polyakov 1 Re"))

2

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = "Polyakov 2 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = 
"Polyakov 2 Re"))

3

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = "Polyakov 3 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = 
"Polyakov 3 Re"))

Eigenvalues

\lambda_{min}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "Minimum Eigenvalue", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = 
"Minimum Eigenvalue"))
catch error
end

\frac{\lambda_{max}}{\lambda_{min}}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$"))
catch error
end

\beta = 2.2

m = -0.12 L16T32

Code
e = load_ensemble("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.12_L16T32_expclv/out_0")
w = load_wilson_flow("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.12_L16T32_expclv/out_wf")
default_thermalise_bin!(e)
default_thermalise_bin!(w)
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Plaquette

Code
display(plot_plaquette(e))

t^2E(t)

Code
display(plot_t2e(w, title=_ensemble_to_latex_string(e)))

W(t)

Code
display(plot_w(w, title=_ensemble_to_latex_string(e)))

Topological Charge

Code
display(plot_tc(w, title=_ensemble_to_latex_string(e)))
display(plot_tc_hist(w, title=_ensemble_to_latex_string(e)))

dS History

Code
try
default_thermalise_bin!(e)
display(plot(1:nrow(e.analysis), e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = "\$dS\$", ylabel = "\$dS\$", xlabel = "Trajectory #"))
display(histogram(e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = 
"\$dS\$", xlabel = "\$dS\$"))
catch error
end
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Polyakov

0

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = "Polyakov 0 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = 
"Polyakov 0 Re"))

1

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = "Polyakov 1 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = 
"Polyakov 1 Re"))

2

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = "Polyakov 2 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = 
"Polyakov 2 Re"))

3

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = "Polyakov 3 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = 
"Polyakov 3 Re"))

Eigenvalues

\lambda_{min}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "Minimum Eigenvalue", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = 
"Minimum Eigenvalue"))
catch error
end

\frac{\lambda_{max}}{\lambda_{min}}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$"))
catch error
end

m = -0.18 L16T32

Code
e = load_ensemble("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.18_L16T32_expclv/out_0")
w = load_wilson_flow("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.18_L16T32_expclv/out_wf")
default_thermalise_bin!(e)
default_thermalise_bin!(w)
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Plaquette

Code
display(plot_plaquette(e))

t^2E(t)

Code
display(plot_t2e(w, title=_ensemble_to_latex_string(e)))

W(t)

Code
display(plot_w(w, title=_ensemble_to_latex_string(e)))

Topological Charge

Code
display(plot_tc(w, title=_ensemble_to_latex_string(e)))
display(plot_tc_hist(w, title=_ensemble_to_latex_string(e)))

dS History

Code
try
default_thermalise_bin!(e)
display(plot(1:nrow(e.analysis), e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = "\$dS\$", ylabel = "\$dS\$", xlabel = "Trajectory #"))
display(histogram(e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = 
"\$dS\$", xlabel = "\$dS\$"))
catch error
end
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Polyakov

0

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = "Polyakov 0 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = 
"Polyakov 0 Re"))

1

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = "Polyakov 1 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = 
"Polyakov 1 Re"))

2

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = "Polyakov 2 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = 
"Polyakov 2 Re"))

3

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = "Polyakov 3 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = 
"Polyakov 3 Re"))

Eigenvalues

\lambda_{min}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "Minimum Eigenvalue", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = 
"Minimum Eigenvalue"))
catch error
end

\frac{\lambda_{max}}{\lambda_{min}}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$"))
catch error
end

m = -0.22 L16T32

Code
e = load_ensemble("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.22_L16T32_expclv/out_0")
w = load_wilson_flow("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.22_L16T32_expclv/out_wf")
default_thermalise_bin!(e)
default_thermalise_bin!(w)
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Plaquette

Code
display(plot_plaquette(e))

t^2E(t)

Code
display(plot_t2e(w, title=_ensemble_to_latex_string(e)))

W(t)

Code
display(plot_w(w, title=_ensemble_to_latex_string(e)))

Topological Charge

Code
display(plot_tc(w, title=_ensemble_to_latex_string(e)))
display(plot_tc_hist(w, title=_ensemble_to_latex_string(e)))

dS History

Code
try
default_thermalise_bin!(e)
display(plot(1:nrow(e.analysis), e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = "\$dS\$", ylabel = "\$dS\$", xlabel = "Trajectory #"))
display(histogram(e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = 
"\$dS\$", xlabel = "\$dS\$"))
catch error
end
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Polyakov

0

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = "Polyakov 0 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = 
"Polyakov 0 Re"))

1

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = "Polyakov 1 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = 
"Polyakov 1 Re"))

2

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = "Polyakov 2 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = 
"Polyakov 2 Re"))

3

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = "Polyakov 3 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = 
"Polyakov 3 Re"))

Eigenvalues

\lambda_{min}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "Minimum Eigenvalue", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = 
"Minimum Eigenvalue"))
catch error
end

\frac{\lambda_{max}}{\lambda_{min}}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$"))
catch error
end

m = -0.231037 L24T24

Code
e = load_ensemble("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.231037_L24T24_expclv/out_0")
w = load_wilson_flow("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.231037_L24T24_expclv/out_wf")
default_thermalise_bin!(e)
default_thermalise_bin!(w)
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Plaquette

Code
display(plot_plaquette(e))

t^2E(t)

Code
display(plot_t2e(w, title=_ensemble_to_latex_string(e)))

W(t)

Code
display(plot_w(w, title=_ensemble_to_latex_string(e)))

Topological Charge

Code
display(plot_tc(w, title=_ensemble_to_latex_string(e)))
display(plot_tc_hist(w, title=_ensemble_to_latex_string(e)))

dS History

Code
try
default_thermalise_bin!(e)
display(plot(1:nrow(e.analysis), e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = "\$dS\$", ylabel = "\$dS\$", xlabel = "Trajectory #"))
display(histogram(e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = 
"\$dS\$", xlabel = "\$dS\$"))
catch error
end
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Polyakov

0

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = "Polyakov 0 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = 
"Polyakov 0 Re"))

1

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = "Polyakov 1 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = 
"Polyakov 1 Re"))

2

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = "Polyakov 2 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = 
"Polyakov 2 Re"))

3

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = "Polyakov 3 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = 
"Polyakov 3 Re"))

Eigenvalues

\lambda_{min}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "Minimum Eigenvalue", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = 
"Minimum Eigenvalue"))
catch error
end

\frac{\lambda_{max}}{\lambda_{min}}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$"))
catch error
end

m = -0.239 L16T32

Code
e = load_ensemble("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.239_L16T32_expclv/out_0")
w = load_wilson_flow("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.239_L16T32_expclv/out_wf")
default_thermalise_bin!(e)
default_thermalise_bin!(w)
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Plaquette

Code
display(plot_plaquette(e))

t^2E(t)

Code
display(plot_t2e(w, title=_ensemble_to_latex_string(e)))

W(t)

Code
display(plot_w(w, title=_ensemble_to_latex_string(e)))

Topological Charge

Code
display(plot_tc(w, title=_ensemble_to_latex_string(e)))
display(plot_tc_hist(w, title=_ensemble_to_latex_string(e)))

dS History

Code
try
default_thermalise_bin!(e)
display(plot(1:nrow(e.analysis), e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = "\$dS\$", ylabel = "\$dS\$", xlabel = "Trajectory #"))
display(histogram(e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = 
"\$dS\$", xlabel = "\$dS\$"))
catch error
end
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Polyakov

0

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = "Polyakov 0 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = 
"Polyakov 0 Re"))

1

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = "Polyakov 1 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = 
"Polyakov 1 Re"))

2

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = "Polyakov 2 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = 
"Polyakov 2 Re"))

3

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = "Polyakov 3 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = 
"Polyakov 3 Re"))

Eigenvalues

\lambda_{min}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "Minimum Eigenvalue", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = 
"Minimum Eigenvalue"))
catch error
end

\frac{\lambda_{max}}{\lambda_{min}}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$"))
catch error
end

m = -0.239000 L32T32

Code
e = load_ensemble("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.239000_L32T32_expclv/out_0")
w = load_wilson_flow("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.239000_L32T32_expclv/out_wf")
default_thermalise_bin!(e)
default_thermalise_bin!(w)
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Plaquette

Code
display(plot_plaquette(e))

t^2E(t)

Code
display(plot_t2e(w, title=_ensemble_to_latex_string(e)))

W(t)

Code
display(plot_w(w, title=_ensemble_to_latex_string(e)))

Topological Charge

Code
display(plot_tc(w, title=_ensemble_to_latex_string(e)))
display(plot_tc_hist(w, title=_ensemble_to_latex_string(e)))

dS History

Code
try
default_thermalise_bin!(e)
display(plot(1:nrow(e.analysis), e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = "\$dS\$", ylabel = "\$dS\$", xlabel = "Trajectory #"))
display(histogram(e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = 
"\$dS\$", xlabel = "\$dS\$"))
catch error
end
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Polyakov

0

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = "Polyakov 0 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = 
"Polyakov 0 Re"))

1

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = "Polyakov 1 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = 
"Polyakov 1 Re"))

2

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = "Polyakov 2 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = 
"Polyakov 2 Re"))

3

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = "Polyakov 3 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = 
"Polyakov 3 Re"))

Eigenvalues

\lambda_{min}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "Minimum Eigenvalue", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = 
"Minimum Eigenvalue"))
catch error
end

\frac{\lambda_{max}}{\lambda_{min}}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$"))
catch error
end

m = -0.246570 L32T32

Code
e = load_ensemble("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.246570_L32T32_expclv/out_0")
w = load_wilson_flow("/home/laurence/projects/RipRep/data/Prepared/b2.2_csw2.299498_m-0.246570_L32T32_expclv/out_wf")
default_thermalise_bin!(e)
default_thermalise_bin!(w)
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Plaquette

Code
display(plot_plaquette(e))

t^2E(t)

Code
display(plot_t2e(w, title=_ensemble_to_latex_string(e)))

W(t)

Code
display(plot_w(w, title=_ensemble_to_latex_string(e)))

Topological Charge

Code
display(plot_tc(w, title=_ensemble_to_latex_string(e)))
display(plot_tc_hist(w, title=_ensemble_to_latex_string(e)))

dS History

Code
try
default_thermalise_bin!(e)
display(plot(1:nrow(e.analysis), e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = "\$dS\$", ylabel = "\$dS\$", xlabel = "Trajectory #"))
display(histogram(e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = 
"\$dS\$", xlabel = "\$dS\$"))
catch error
end
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Polyakov

0

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = "Polyakov 0 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = 
"Polyakov 0 Re"))

1

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = "Polyakov 1 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = 
"Polyakov 1 Re"))

2

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = "Polyakov 2 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = 
"Polyakov 2 Re"))

3

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = "Polyakov 3 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = 
"Polyakov 3 Re"))

Eigenvalues

\lambda_{min}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "Minimum Eigenvalue", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = 
"Minimum Eigenvalue"))
catch error
end

\frac{\lambda_{max}}{\lambda_{min}}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$"))
catch error
end

\beta = 2.3

m = -0.282 L32T32

Code
e = load_ensemble("/home/laurence/projects/RipRep/data/Prepared/b2.3_csw1.9059_m-0.282_L32T32_expclv/out_0")
w = load_wilson_flow("/home/laurence/projects/RipRep/data/Prepared/b2.3_csw1.9059_m-0.282_L32T32_expclv/out_wf")
default_thermalise_bin!(e)
default_thermalise_bin!(w)
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Plaquette

Code
display(plot_plaquette(e))

t^2E(t)

Code
display(plot_t2e(w, title=_ensemble_to_latex_string(e)))

W(t)

Code
display(plot_w(w, title=_ensemble_to_latex_string(e)))

Topological Charge

Code
display(plot_tc(w, title=_ensemble_to_latex_string(e)))
display(plot_tc_hist(w, title=_ensemble_to_latex_string(e)))

dS History

Code
try
default_thermalise_bin!(e)
display(plot(1:nrow(e.analysis), e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = "\$dS\$", ylabel = "\$dS\$", xlabel = "Trajectory #"))
display(histogram(e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = 
"\$dS\$", xlabel = "\$dS\$"))
catch error
end
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Polyakov

0

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = "Polyakov 0 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = 
"Polyakov 0 Re"))

1

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = "Polyakov 1 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = 
"Polyakov 1 Re"))

2

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = "Polyakov 2 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = 
"Polyakov 2 Re"))

3

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = "Polyakov 3 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = 
"Polyakov 3 Re"))

Eigenvalues

\lambda_{min}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "Minimum Eigenvalue", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = 
"Minimum Eigenvalue"))
catch error
end

\frac{\lambda_{max}}{\lambda_{min}}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$"))
catch error
end

m = -0.2987 L42T42

Code
e = load_ensemble("/home/laurence/projects/RipRep/data/Prepared/b2.3_csw1.9059_m-0.2987_L42T42_expclv/out_0")
w = load_wilson_flow("/home/laurence/projects/RipRep/data/Prepared/b2.3_csw1.9059_m-0.2987_L42T42_expclv/out_wf")
default_thermalise_bin!(e)
default_thermalise_bin!(w)
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Plaquette

Code
display(plot_plaquette(e))

t^2E(t)

Code
display(plot_t2e(w, title=_ensemble_to_latex_string(e)))

W(t)

Code
display(plot_w(w, title=_ensemble_to_latex_string(e)))

Topological Charge

Code
display(plot_tc(w, title=_ensemble_to_latex_string(e)))
display(plot_tc_hist(w, title=_ensemble_to_latex_string(e)))

dS History

Code
try
default_thermalise_bin!(e)
display(plot(1:nrow(e.analysis), e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = "\$dS\$", ylabel = "\$dS\$", xlabel = "Trajectory #"))
display(histogram(e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = 
"\$dS\$", xlabel = "\$dS\$"))
catch error
end
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Polyakov

0

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = "Polyakov 0 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = 
"Polyakov 0 Re"))

1

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = "Polyakov 1 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = 
"Polyakov 1 Re"))

2

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = "Polyakov 2 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = 
"Polyakov 2 Re"))

3

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = "Polyakov 3 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = 
"Polyakov 3 Re"))

Eigenvalues

\lambda_{min}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "Minimum Eigenvalue", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = 
"Minimum Eigenvalue"))
catch error
end

\frac{\lambda_{max}}{\lambda_{min}}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$"))
catch error
end

m = -0.29 L36T36

Code
e = load_ensemble("/home/laurence/projects/RipRep/data/Prepared/b2.3_csw1.9059_m-0.29_L36T36_expclv/out_0")
w = load_wilson_flow("/home/laurence/projects/RipRep/data/Prepared/b2.3_csw1.9059_m-0.29_L36T36_expclv/out_wf")
default_thermalise_bin!(e)
default_thermalise_bin!(w)
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Plaquette

Code
display(plot_plaquette(e))

t^2E(t)

Code
display(plot_t2e(w, title=_ensemble_to_latex_string(e)))

W(t)

Code
display(plot_w(w, title=_ensemble_to_latex_string(e)))

Topological Charge

Code
display(plot_tc(w, title=_ensemble_to_latex_string(e)))
display(plot_tc_hist(w, title=_ensemble_to_latex_string(e)))

dS History

Code
try
default_thermalise_bin!(e)
display(plot(1:nrow(e.analysis), e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = "\$dS\$", ylabel = "\$dS\$", xlabel = "Trajectory #"))
display(histogram(e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = 
"\$dS\$", xlabel = "\$dS\$"))
catch error
end
┌ Error: Thermalisation time does not get rid of all missing configurations!
└ @ Main /home/laurence/projects/RipRep/src/utilities.jl:24

Polyakov

0

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = "Polyakov 0 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = 
"Polyakov 0 Re"))

1

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = "Polyakov 1 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = 
"Polyakov 1 Re"))

2

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = "Polyakov 2 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = 
"Polyakov 2 Re"))

3

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = "Polyakov 3 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = 
"Polyakov 3 Re"))

Eigenvalues

\lambda_{min}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "Minimum Eigenvalue", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = 
"Minimum Eigenvalue"))
catch error
end

\frac{\lambda_{max}}{\lambda_{min}}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$"))
catch error
end

\beta = 2.5

m = -0.25 L16T32

Code
e = load_ensemble("/home/laurence/projects/RipRep/data/Prepared/b2.5_csw1.672_m-0.25_L16T32_expclv/out_0")
w = load_wilson_flow("/home/laurence/projects/RipRep/data/Prepared/b2.5_csw1.672_m-0.25_L16T32_expclv/out_wf")
default_thermalise_bin!(e)
default_thermalise_bin!(w)

Plaquette

Code
display(plot_plaquette(e))

t^2E(t)

Code
display(plot_t2e(w, title=_ensemble_to_latex_string(e)))

W(t)

Code
display(plot_w(w, title=_ensemble_to_latex_string(e)))

Topological Charge

Code
display(plot_tc(w, title=_ensemble_to_latex_string(e)))
display(plot_tc_hist(w, title=_ensemble_to_latex_string(e)))

dS History

Code
try
default_thermalise_bin!(e)
display(plot(1:nrow(e.analysis), e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = "\$dS\$", ylabel = "\$dS\$", xlabel = "Trajectory #"))
display(histogram(e.analysis[:, :dS], title = _ensemble_to_latex_string(e), label = 
"\$dS\$", xlabel = "\$dS\$"))
catch error
end

Polyakov

0

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = "Polyakov 0 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_0], title = _ensemble_to_latex_string(e), label = 
"Polyakov 0 Re"))

1

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = "Polyakov 1 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_1], title = _ensemble_to_latex_string(e), label = 
"Polyakov 1 Re"))

2

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = "Polyakov 2 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_2], title = _ensemble_to_latex_string(e), label = 
"Polyakov 2 Re"))

3

Code
display(plot(1:nrow(e.data), e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = "Polyakov 3 Re", xlabel = "Trajectory #"))
display(histogram(e.data[:, :Polyakov_3], title = _ensemble_to_latex_string(e), label = 
"Polyakov 3 Re"))

Eigenvalues

\lambda_{min}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "Minimum Eigenvalue", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = 
"Minimum Eigenvalue"))
catch error
end

\frac{\lambda_{max}}{\lambda_{min}}

Code
try
display(plot(1:nrow(e.data), e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$", xlabel = "Trajectory #"))
display(histogram(e.data[:, :eig_max]./e.data[:, :eig_min], title = _ensemble_to_latex_string(e), label = "\$\\frac{\\lambda_{max}}{\\lambda_{min}}\$"))
catch error
end