Main Image
Volume 256 - 34th annual International Symposium on Lattice Field Theory (LATTICE2016) - Algorithms and Machines
Progress Report on Staggered Multigrid
E.S. Weinberg,* R.C. Brower, K. Clark, A. Strelchenko
*corresponding author
Full text: pdf
Pre-published on: 2017 January 31
Published on: 2017 March 24
Multigrid methods have become the optimal solvers for Lattice QCD, first for the Wilson-clover discrete representation and more recently for the domain wall formulation. However, at present, the ensembles with the largest lattices use a third discretization, staggered. At the physical pion mass, propagator inversions with the staggered operator require $\mathcal{O}(10,000)$ iterations using the Conjugate Gradient algorithm. The staggered discretization raises new challenges for a multigrid implementation. It is a first order discretization which in four-dimensions corresponds to four copies of a Dirac fermion in the continuum. Here we report on our investigation into a new geometric adaptive multigrid for staggered fermions both for the normal equations and for a first order projection based on a novel blocking scheme stabilized by a second order gauged Laplacian. Current numerical tests, applied to the two-dimensional staggered representation of the two-flavor Schwinger model, are promising. Further improvements are under way as well as generalizations and tests for four-dimensional QCD using multigrid preconditioned staggered solvers in the QUDA library for multi-GPU architectures.
Open Access
Creative Commons LicenseCopyright owned by the author(s) under the term of the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.