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] aabw_v1_0.gz(13 Kbytes)
Manuscript Title: A vectorized Monte Carlo algorithm for computing Wilson line observables in SU(2) gauge theory on a BCH lattice.
Authors: W. Celmaster, F. Green, R. Gupta, E. Kovacs
Program title: LINE, BESSEL AND TRANSLATE
Catalogue identifier: AABW_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 36(1985)409
Programming language: Fortran.
Computer: CRAY-1.
Operating system: LTSS/CTSS (MODIFIED VERSION).
RAM: 150K words
Word size: 64
Peripherals: disc.
Keywords: Particle physics, Elementary, Lattice gauge theory, Su(2) gauge theory, Body-centred Hypercubic lattice, Wilson lines, Polyakov loops, Finite temperature, Deconfinement, Monte carlo techniques, Quark potential, Vector processors, Quantum chromodynamics.
Classification: 11.5.

Nature of problem:
The usual hypercubic formulation of lattice gauge theory may be improved by an alternative formulation of the theory on a body-centred hypercubic (BCH) lattice. We describe a vectorized implementation of the Monte Carlo algorithm for calculating Wilson lines and average plaquettes for an SU(2) gauge theory on a BCH lattice.

Solution method:
SU(2) gauage theory is simulated by importance-sampling Monte Carlo methods on a BCH lattice. The heat-bath method of Creutz is used to thermalize the lattice. A local thermalization technique of Parisi et al. is used to reduce statistical fluctuations in the measurement of Wilson lines.

Restrictions:
The main memory allocation is required for the arrays AL1,...,AL4 each dimensioned to 12 * (ISPACE**3 * (ITIME+4)), where ISPACE is the spatial size and ITIME is the temporal size of the lattice. The program size must not exceed the physical memory of the computer. This imposes a practical limit of ITIME=6 and ISPACE=12 on the present machine. One further I/O restriction (with the FORTLIB utility library) occurs because unformatted records must be smaller than 262,143 words. In our program that limits the size of AL1,...,AL4.

Running time:
The test run took 30 s on the CRAY-1 with ISPACE = 11 and ITIME = 4 for 10 sweeps of the lattice and 2 sets of measurements.