Fourier acceleration in lattice gauge theories. II. Matrix inversion and the quark propagator

Abstract

We present a new algorithm for inverting the quark propagator in lattice QCD that removes the critical slowing down associated with the continuum limit. Both a theoretical analysis and a practical implementation are given, together with results from numerical simulations that illustrate the algorithm. The acceleration speeds convergence by factors of 3–4 on small lattices and should lead to proportionately larger gains on larger lattices.