Parallelizing the QUDA Library for Multi-GPU Calculations in Lattice Quantum Chromodynamics

Graphics Processing Units (GPUs) are having a transformational effect on numerical lattice quantum chromodynamics (LQCD) calculations of importance in nuclear and particle physics. The QUDA library provides a package of mixed precision sparse matrix linear solvers for LQCD applications, supporting single GPUs based on NVIDIA's Compute Unified Device Architecture (CUDA). This library, interfaced to the QDP++/Chroma framework for LQCD calculations, is currently in production use on the "9g" cluster at the Jefferson Laboratory, enabling unprecedented price/performance for a range of problems in LQCD. Nevertheless, memory constraints on current GPU devices limit the problem sizes that can be tackled. In this contribution we describe the parallelization of the QUDA library onto multiple GPUs using MPI, including strategies for the overlapping of communication and computation. We report on both weak and strong scaling for up to 32 GPUs interconnected by InfiniBand, on which we sustain in excess of 4 Tflops.Comment: 11 pages, 7 figures, to appear in the Proceedings of Supercomputing 2010 (submitted April 12, 2010

Diquark correlations in baryons on the lattice with overlap quarks

We evaluate baryon wave functions in both the Coulomb and Landau gauge in lattice QCD. These are constructed from quark propagators calculated with the overlap Dirac operator on quenched gauge configurations at beta = 6. By comparing baryon states that differ in their diquark content, we find evidence for enhanced correlation in the scalar diquark channel, as favored by quark models. We also summarize earlier results for diquark masses in the Landau gauge, casting them in a form more easily compared with subsequent studies.Comment: 9 pages, 12 figures, v2 as appears in PR

Light hadron spectroscopy in quenched QCD with overlap fermions

A simulation of quenched QCD with the overlap Dirac operator has been carried out using 100 Wilson gauge configurations at beta = 6 on an 18^3 x 64 lattice and at beta = 5.85 on a 14^3 x 48 lattice. Here we present results for meson masses, meson final state "wave functions," decay constants, and other observables, as well as details on our algorithmic and data analysis techniques. We also summarize results for baryon masses and quark and diquark propagators in the Landau gauge.Comment: 11 pages, 17 figures. Combined contribution by C.R. and R.B. to Lattice 2005 (Hadron spectrum and quark masses), Trinity College, Dublin, 25-30 July 200

Exploring strange nucleon form factors on the lattice

We discuss techniques for evaluating sea quark contributions to hadronic form factors on the lattice and apply these to an exploratory calculation of the strange electromagnetic, axial, and scalar form factors of the nucleon. We employ the Wilson gauge and fermion actions on an anisotropic 24^3 x 64 lattice, probing a range of momentum transfer with Q^2 < 1 GeV^2. The strange electric and magnetic form factors, G_E^s(Q^2) and G_M^s(Q^2), are found to be small and consistent with zero within the statistics of our calculation. The lattice data favor a small negative value for the strange axial form factor G_A^s(Q^2) and exhibit a strong signal for the bare strange scalar matrix element _0. We discuss the unique systematic uncertainties affecting the latter quantity relative to the continuum, as well as prospects for improving future determinations with Wilson-like fermions.Comment: 19 pages, 11 figures; v2 includes additional references; v3 as appears in PR

Light hadron and diquark spectroscopy in quenched QCD with overlap quarks on a large lattice

A simulation of quenched QCD with the overlap Dirac operator has been completed using 100 Wilson gauge configurations at beta = 6 on an 18^3 x 64 lattice and at beta = 5.85 on a 14^3 x 48 lattice, both in Landau gauge. We present results for light meson and baryon masses, meson final state "wave functions," and other observables, as well as some details on the numerical techniques that were used. Our results indicate that scaling violations, if any, are small. We also present an analysis of diquark correlations using the quark propagators generated in our simulation.Comment: 28 LaTeX pages, 41 figures, v2: minor updates, version published in JHE