We present new parallel software, SUSY LATTICE, for lattice studies of four-dimensional N = 4 supersymmetric Yang Mills theory with gauge group SU(N). The lattice action is constructed to exactly preserve a single supersymmetry charge at non-zero lattice spacing, up to additional potential terms included to stabilize numerical simulations. The software evolved from the MILC code for lattice QCD, and retains a similar large-scale framework despite the different target theory. Many routines are adapted from an existing serial code (Catterall and Joseph, 2012), which SUSY LATTICE supersedes. This paper provides an overview of the new parallel software, summarizing the lattice system, describing the applications that are currently provided and explaining their basic workflow for non-experts in lattice gauge theory. We discuss the parallel performance of the code, and highlight some notable aspects of the documentation for those interested in contributing to its future development. Program summary Program title: SUSY LATTICE Catalogue identifier: AELS_v2_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AELS_v2_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 72821 No. of bytes in distributed program, including test data, etc.: 1294754 Distribution format: tar.gz Programming language: C. Operating system: Any, tested on Linux workstations and MPI clusters with InfiniBand. Has the code been vectorized or parallelized?: Code is parallelized Classification: 11.5. Catalogue identifier of previous version: AELS_v1_0 Journal reference of previous version: Comput. Phys. Comm. 183 (2012) 1336. External routines: Lapack, Bias, Primme (optional) Does the new version supersede the previous version?: Yes Nature of problem: To carry out non-perturbative Monte Carlo importance sampling for maximally supersymmetric Yang Mills theories in two and four dimensions, and thereby compute observables including Wilson loops, fermion bilinears, eigenvalues of DtD and the Pfafflan of the sparse fermion operator D. Solution method: The central application is 4 rational hybrid Monte Carlo algorithm with a two-level Omelyan molecular dynamics integrator. Gauge field configurations generated by this application may be saved to disk for subsequent measurements of additional observables. Input parameters for either configuration generation or analysis may be entered manually or read from a file. Reasons for new version: The code is parallelized to improve performance and scalability. In addition, several new features and measurements are provided. Summary of revisions: The program is completely rewritten on the basis of the MILC code for lattice QCD, and several new features and measurements are added. Restrictions: The code is currently restricted to two-dimensional N = (2, 2) and four-dimensional N = 4 supersymmetric Yang-Mills theories. The process of topological twisting on which it is based can also be applied to a few other systems, as discussed in Sections 1 and 5. Additional comments: Further documentation is provided in the distribution file, including a set of test runs with reference output in the testsuite directory. Running time: From seconds to hours depending on the computational task, lattice volume, gauge group, and desired statistics, as well as on the computing platform and number of cores used. For example, rational hybrid Monte Carlo generation of 50 molecular dynamics time units for a 16(3) x 32 lattice volume with gauge group U(2) takes approximately 16 hours on 512 cores of the USQCD bc cluster at Fermilab, while standard measurements on a saved 8(3) x 24 U(2) configuration require only 8 seconds on one eight-core workstation. (C) 2015 Elsevier B.V. All rights reserved.
