Accelerating the solution of families of shifted linear systems with CUDA

We describe the GPU implementation of shifted or multimass iterative solvers for sparse linear systems of the sort encountered in lattice gauge theory. We provide a generic tool that can be used by those without GPU programming experience to accelerate the simulation of a wide array of theories. We stress genericity, which is important to allow the simulation of candidate theories for new physics at LHC, and for the study of various supersymmetric theories. We find significant speed ups, which we conservatively bound below at at least twelve times, that promise to put a variety of research questions within practical reach.

[1]  Joshua A. Anderson,et al.  General purpose molecular dynamics simulations fully implemented on graphics processing units , 2008, J. Comput. Phys..

[2]  Kipton Barros,et al.  Solving lattice QCD systems of equations using mixed precision solvers on GPUs , 2009, Comput. Phys. Commun..

[3]  J. Maldacena.  The Large-N Limit of Superconformal Field Theories and Supergravity , 1997, hep-th/9711200.

[4]  Exact 2+1 flavour RHMC simulations , 2005, hep-lat/0409133.

[5]  Greg van Anders,et al.  First results from lattice simulation of the PWMM , 2010, 1003.4952.

[6]  Róbert Mészáros,et al.  Air pollution modelling using a Graphics Processing Unit with CUDA , 2010, Comput. Phys. Commun..

[7]  Yousef Saad,et al.  Efficient Parallel Solution of Parabolic Equations: Implicit Methods on the Cedar Multicluster , 1989, PPSC.

[8]  Y. Saad,et al.  Arnoldi methods for large Sylvester-like observer matrix equations, and an associated algorithm for partial spectrum assignment , 1991 .

[9]  Michael Garland,et al.  Implementing sparse matrix-vector multiplication on throughput-oriented processors , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.

[10]  Zoltán Fodor,et al.  Lattice QCD as a video game , 2007, Comput. Phys. Commun..

[11]  Towards lattice simulation of the gauge theory duals to black holes and hot strings , 2007, 0706.3518.

[12]  James C. Osborn.  Initial guesses for multi-shift solvers , 2008 .

[13]  B. Jegerlehner.  Krylov space solvers for shifted linear systems , 1996, hep-lat/9612014.