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] aabt_v2_0.gz(20 Kbytes)
Manuscript Title: An efficiently microtasked CRAY Y-MP C90 version of the Kuba- Moriarty SU(3) gauge theory simulation program.
Authors: K.J.M. Moriarty, S. Sanielevici, T. Trappenberg, D.W. Kuba
Program title: SU3
Catalogue identifier: AABT_v2_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 76(1993)87
Programming language: Fortran.
Computer: CRAY Y-MP C90/16256.
Operating system: UNICOS 7.C.
Word size: 64
Keywords: Particle physics, Elementary, Qcd, Lattice gauge theory, Su(3) gauge theory, Monte carlo method, Vector processor, Microtasking.
Classification: 11.5.

Nature of problem:
The use of nonperturbative computation techniques to study strongly interacting matter requires the discretization of the continuous space- time on a 4-dimensional lattice. Then the quantum field theory is formally equivalent to a statistical mechanical model that can be simulated on a computer by Monte Carlo methods. This program is a microtasked implementation of the multihit Metropolis algorithm on a CRAY Y-MP multiprocessor system. It calculates the Wilson loops of size up to L/2 x L/2 for L**4 lattices for SU(3) gauge theory with the Wilson action.

Solution method:
SU(3) gauge theory is simulated on a hypercubical lattice by Monte Carlo methods with the multihit Metropolis algorithm on a microtasking vector computer, the CRAY Y-MP C90/16256.

For the lattice of L sites in each of the space-time directions we need storage for 4L**4 SU(3) matrices, or 4L**4 x 18 real variables. This is the only restriction of the program.

Running time:
The test run on a 16**4 lattice for SU(3) with 40 Monte Carlo iterations calculating the Wilson loops up to size 8 x 8 took 8.2 seconds for execution on the CRAY Y-MP C90/16. This is 19 times faster than the old record on the CRAY X-MP/4 [1].

[1] D.W.Kuba and K.J.M. Moriarty, Comp. Phys. Comm. 36(1985) 351-362.