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: January 31, 2017
Published on: March 24, 2017
Abstract
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.
DOI: https://doi.org/10.22323/1.256.0273
How to cite

Metadata are provided both in "article" format (very similar to INSPIRE) as this helps creating very compact bibliographies which can be beneficial to authors and readers, and in "proceeding" format which is more detailed and complete.

Open Access
Copyright owned by the author(s) under the term of the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.