Solving lattice QCD systems of equations using mixed precision solvers on GPUs

By M. Clark, R. Babich, Kipton Barros, R. Brower, C. Rebbi. Computer Physics Communications, 2009

506 citations recorded by the source corpus.

Abstract Modern graphics hardware is designed for highly parallel numerical tasks and promises significant cost and performance benefits for many scientific applications. One such application is lattice quantum chromodynamics (lattice QCD), where the main computational challenge is to efficiently solve the discretized Dirac equation in the presence of an SU ( 3 ) gauge field. Using NVIDIA's CUDA platform we have implemented a Wilson–Dirac sparse matrix–vector product that performs at up to 40, 135 and 212 Gflops for double, single and half precision respectively on NVIDIA's GeForce GTX 280 GPU. We have developed a new mixed precision approach for Krylov solvers using reliable updates which allows for full double precision accuracy while using only single or half precision arithmetic for the bulk of the computation. The resulting BiCGstab and CG solvers run in excess of 100 Gflops and, in terms of iterations until convergence, perform better than the usual defect-correction approach for mixed precision.

Subfields: Computational Theory and Mathematics, Nuclear and High Energy Physics

View Paper: https://www.semanticscholar.org/paper/06721d4acf2544a77fdb4e8972d1e424f55cb00d

How researchers see this paper

The researchers whose trusted networks rank it highest, strongest advocate first. Every researcher gets a different order.

In Nuclear and High Energy Physics

  1. J. Dudek
  2. Christopher E. Thomas
  3. Zoltán Fodor
  4. Christian Hoelbling
  5. R. Young
  6. R. Briceño
  7. William Detmold
  8. S. Meinel

Citations in the corpus