Computing platforms equipped with accelerators like GPUs have proven to provide great computational power. However, exploiting such platforms for existing scientific applications is not a trivial task. Current GPU programming frameworks such as CUDA
C/C++ require low-level programming from the developer in order to achieve high performance code. As a result porting of applications to GPUs is typically limited to time-dominant algorithms and routines, leaving the remainder not accelerated which can open a serious Amdahls law issue. The lattice QCD application Chroma allows to explore a different porting strategy. The layered structure of the software architecture logically separates the data-parallel from the application layer. The QCD Data-Parallel software layer provides data types and expressions with stencil-like operations suitable for lattice field theory and Chroma implements algorithms in terms of this high-level interface. Thus by porting the low-level layer one can effectively move the whole application in one swing to a different platform. The QDP-JIT/PTX library, the reimplementation of the low-level layer, provides a framework for lattice QCD calculations for the CUDA architecture. The complete software interface is supported and thus applications can be run unaltered on GPU-based parallel computers. This reimplementation was possible due to the availability of a JIT compiler (part of the NVIDIA Linux kernel driver) which translates an assembly-like language (PTX) to GPU code. The expression template technique is used to build PTX code generators and a software cache manages the GPU memory. This reimplementation allows us to deploy an efficient implementation of the full gauge-generation program with dynamical fermions on large-scale GPU-based machines such as Titan and Blue Waters which accelerates the algorithm by more than an order of magnitude.
The energies of the excited states of the Nucleon, $Delta$ and $Omega$ are computed in lattice QCD, using two light quarks and one strange quark on anisotropic lattices. The calculation is performed at three values of the light quark mass, correspond
ing to pion masses $m_{pi}$ = 392(4), 438(3) and 521(3) MeV. We employ the variational method with a large basis of interpolating operators enabling six energies in each irreducible representation of the lattice to be distinguished clearly. We compare our calculation with the low-lying experimental spectrum, with which we find reasonable agreement in the pattern of states. The need to include operators that couple to the expected multi-hadron states in the spectrum is clearly identified.
In order to advance lattice calculations of moments of unpolarized, helicity, and transversity distributions, electromagnetic form factors, and generalized form factors of the nucleon to a new level of precision, this work investigates several key as
pects of precision lattice calculations. We calculate the number of configurations required for constant statistical errors as a function of pion mass, describe the coherent sink method to help achieve these statistics, examine the statistical correlations between separate measurements, study correlations in the behavior of form factors at different momentum transfer, examine volume dependence, and compare mixed action results with those using comparable dynamical domain wall configurations. We also show selected form factor results and comment on the QCD evolution of our calculations of the flavor non-singlet nucleon angular momentum.
Moments of the generalized parton distributions of the nucleon, calculated with a mixed action of domain wall valence quarks and asqtad staggered sea quarks, are presented for pion masses extending down to 359 MeV. Results for the moments of the unpo
larized, helicity, and transversity distributions are given and compared to the available experimental measurements. Additionally, a selection of the generalized form factors are shown and the implications for the spin decomposition and transverse structure of the nucleon are discussed. Particular emphasis is placed on understanding systematic errors in the lattice calculation and exploring a variety of chiral extrapolations.
Energies for excited isospin I=1/2 and I=3/2 states that include the nucleon and Delta families of baryons are computed using quenched, anisotropic lattices. Baryon interpolating field operators that are used include nonlocal operators that provide G
_2 irreducible representations of the octahedral group. The decomposition of spin 5/2 or higher spin states is realized for the first time in a lattice QCD calculation. We observe patterns of degenerate energies in the irreducible representations of the octahedral group that correspond to the subduction of the continuum spin 5/2 or higher. The overall pattern of low-lying excited states corresponds well to the pattern of physical states subduced to the irreducible representations of the octahedral group.
