Adaptive multigrid algorithm for the lattice Wilson-Dirac operator

R. Babich, J. Brannick, R. C. Brower, M. A. Clark, T. A. Manteuffel, S. F. McCormick, J. C. Osborn, C. Rebbi

Research output: Contribution to journalArticle

66 Scopus citations

Abstract

We present an adaptive multigrid solver for application to the non-Hermitian Wilson-Dirac system of QCD. The key components leading to the success of our proposed algorithm are the use of an adaptive projection onto coarse grids that preserves the near null space of the system matrix together with a simplified form of the correction based on the so-called γ5-Hermitian symmetry of the Dirac operator. We demonstrate that the algorithm nearly eliminates critical slowing down in the chiral limit and that it has weak dependence on the lattice volume.

Original languageEnglish (US)
Article number201602
JournalPhysical Review Letters
Volume105
Issue number20
DOIs
StatePublished - Nov 11 2010

    Fingerprint

All Science Journal Classification (ASJC) codes

  • Physics and Astronomy(all)

Cite this

Babich, R., Brannick, J., Brower, R. C., Clark, M. A., Manteuffel, T. A., McCormick, S. F., Osborn, J. C., & Rebbi, C. (2010). Adaptive multigrid algorithm for the lattice Wilson-Dirac operator. Physical Review Letters, 105(20), [201602]. https://doi.org/10.1103/PhysRevLett.105.201602