Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aeet_v1_1.tar.gz(210 Kbytes)
Manuscript Title: Program package for multicanonical simulations of U(1) lattice gauge theory
Authors: Alexei Bazavov, Bernd A. Berg
Program title: STMC_U1MUCA_v1_02
Catalogue identifier: AEET_v1_1
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 184(2013)1075
Programming language: Fortran 77 compatible with Fortran 90 and 95.
Computer: Any capable of compiling and executing Fortran code.
Operating system: Any capable of compiling and executing Fortran code.
RAM: 10 MB and up depending on the lattice size used.
Keywords: Markov Chain Monte Carlo, Multicanonical, Wang-Landau Recursion, Fortran, Lattice gauge theory, U(1) gauge group, Phase transition of continuous systems.
Classification: 11.5.

Does the new version supersede the previous version?: Yes

Nature of problem:
Efficient Markov chain Monte Carlo simulation of U(1) lattice gauge theory (or other continuous systems) close to its phase transition. Measurements and analysis of the action per plaquette, the specific heat, Polyakov loops and their structure factors.

Solution method:
Multicanonical simulations with an initial Wang-Landau recursion to determine suitable weight factors. Reweighting to physical values using logarithmic coding and calculating jackknife error bars.

Reasons for new version:
The previous version was developed for the g77 compiler Fortran 77 version. Compiler errors were encountered with Fortran 90 and Fortran 95 compilers (specified below).

Summary of revisions:
epsilon=one/10**10 is replaced by epsilon/10.0D10 in the parameter statements of the subroutines u1_bmha.f, u1_mucabmha.f, u1wl_backup.f, u1wlread_backup.f of the folder Libs/U1_par. For the tested compilers script files are added in the folder ExampleRuns and readme.txt files are now provided in all subfolders of ExampleRuns. The gnuplot driver files produced by the routine hist_gnu.f of Libs/Fortran are adapted to the syntax required by gnuplot version 4.0 and higher.

Due to the use of explicit real*8 initialization the conversion into real*4 will require extra changes besides replacing the implicit.sta file by its real*4 version.

Unusual features:
The programs have to be compiled with the script files like those contained in the folder ExampleRuns as explained in the original paper.

Running time:
The prepared tests runs took up to 74 minutes to execute on a 2GHz PC.