Parallel-Transported Multi-Grid for inverting the Dirac-operator: variants of the method and their efficiency

R. Ben-Av*, M. Harmatz, P. G. Lauwers, S. Solomon

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

6 Scopus citations

Abstract

Using U(1) lattice gauge theory in two dimensions as a test case, we discuss in detail several variants of the Parallel-Transported Multi-Grid (PTMG) with different restriction and interpolation procedures, resulting in radically different efficiencies. In particular, it is optimal to stop the coarsening once the coarse mesh size is of the order of the largest correlation length ξ. A dependence of the efficiency on the topological charge Q of the two-dimensional U(1) configurations was not detected, at least not for the most efficient PTMG variant in the physically relevant range of the parameters. In order to investigate this phenomenon in a reliable way, we developed and applied a Monte Carlo simulation algorithm using global update steps, which ensures ergodic covering of all topological sectors of the theory. For large enough systems and correlation lengths, an optimally tuned PTMG algorithm yields a very high convergence rate (one order of magnitude per cycle) and it beats by far the Conjugate-Gradient method in CPU time. Similar conclusions hold for other relevant lattice gauge theories.

Original languageEnglish
Pages (from-to)623-666
Number of pages44
JournalNuclear Physics B
Volume405
Issue number2-3
DOIs
StatePublished - 20 Sep 1993

Fingerprint

Dive into the research topics of 'Parallel-Transported Multi-Grid for inverting the Dirac-operator: variants of the method and their efficiency'. Together they form a unique fingerprint.

Cite this