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
Catalogue identifier: AABW_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 36(1985)409
Programming language: Fortran.
Computer: CRAY-1.
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.

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.