A multigrid accelerated eigensolver for the Hermitian Wilson-Dirac operator in lattice QCD

Abstract

Eigenvalues of the Hermitian Wilson-Dirac operator are of special interest in several lattice QCD simulations, e.g., for noise reduction when evaluating all-to-all propagators. In this paper we present a Davidson-type eigensolver that utilizes the structural properties of the Hermitian Wilson-Dirac operator $Q$ to compute eigenpairs of this operator corresponding to small eigenvalues. The main idea is to exploit a synergy between the (outer) eigensolver and its (inner) iterative scheme which solves shifted linear systems. This is achieved by adapting the multigrid DD-$\alpha$AMG algorithm to a solver for shifted systems involving the Hermitian Wilson-Dirac operator. We demonstrate that updating the coarse grid operator using eigenvector information obtained in the course of the generalized Davidson method is crucial to achieve good performance when calculating many eigenpairs, as our study of the local coherence shows. We compare our method with the commonly used software-packages PARPACK and PRIMME in numerical tests, where we are able to achieve significant improvements, with speed-ups of up to one order of magnitude and a near-linear scaling with respect to the number of eigenvalues. For illustration we compare the distribution of the small eigenvalues of $Q$ on a $64\times 32^3$ lattice with what is predicted by the Banks-Casher relation in the infinite volume limit.

References in corpus (18)

• Simulation of QCD with N_f=2+1 flavors of non-perturbatively improved Wilson fermions
• Practical all-to-all propagators for lattice QCD
• The strange quark mass and Lambda parameter of two flavor QCD
• Adaptive multigrid algorithm for the lattice Wilson-Dirac operator
• Chiral symmetry breaking and the Banks--Casher relation in lattice QCD with Wilson quarks
• Local coherence and deflation of the low quark modes in lattice QCD
• Stability of lattice QCD simulations and the thermodynamic limit
• On the low fermionic eigenmode dominance in QCD on the lattice
• Low-energy couplings of QCD from current correlators near the chiral limit
• Improving meson two-point functions in lattice QCD
• Solution of the Dirac equation in lattice QCD using a domain decomposition method
• Adaptive Multigrid Algorithm for Lattice QCD
• High precision renormalization of the flavour non-singlet Noether currents in lattice QCD with Wilson quarks
• Remarks on strange-quark simulations with Wilson fermions
• The moment $⟨x{⟩}_{u-d}$ of the nucleon from $N_f=2$ lattice QCD down to nearly physical quark masses
• HQET at order $1/m$: II. Spectroscopy in the quenched approximation
• Variance reduction with practical all-to-all lattice propagators
• A Multi-level Correction Scheme for Eigenvalue Problems