Multigrid deflation for Lattice QCD
arXiv:1909.12234 · doi:10.1016/j.jcp.2020.109356
Abstract
Computing the trace of the inverse of large matrices is typically addressed through statistical methods. Deflating out the lowest eigenvectors or singular vectors of the matrix reduces the variance of the trace estimator. This work summarizes our efforts to reduce the computational cost of computing the deflation space while achieving the desired variance reduction for Lattice QCD applications. Previous efforts computed the lower part of the singular spectrum of the Dirac operator by using an eigensolver preconditioned with a multigrid linear system solver. Despite the improvement in performance in those applications, as the problem size grows the runtime and storage demands of this approach will eventually dominate the stochastic estimation part of the computation. In this work, we propose to compute the deflation space in one of the following two ways. First, by using an inexact eigensolver on the Hermitian, but maximally indefinite, operator . Second, by exploiting the fact that the multigrid prolongator for this operator is rich in components toward the lower part of the singular spectrum. We show experimentally that the inexact eigensolver can approximate the lower part of the spectrum even for ill-conditioned operators. Also, the deflation based on the multigrid prolongator is more efficient to compute and apply, and, despite its limited ability to approximate the fine level spectrum, it obtains similar variance reduction on the trace estimator as deflating with approximate eigenvectors from the fine level operator.
References in corpus (15)
- Practical all-to-all propagators for lattice QCD
- Effective noise reduction techniques for disconnected loops in Lattice QCD
- Adaptive multigrid algorithm for the lattice Wilson-Dirac operator
- Local coherence and deflation of the low quark modes in lattice QCD
- On the low fermionic eigenmode dominance in QCD on the lattice
- Adaptive Multigrid Algorithm for Lattice QCD
- Parton Distribution Functions from Ioffe time pseudo-distributions
- Computing and deflating eigenvalues while solving multiple right hand side linear systems in Quantum Chromodynamics
- High-precision calculation of the strange nucleon electromagnetic form factors
- Exploring strange nucleon form factors on the lattice
- Three-dimensional effective theories for the two Higgs doublet model at high temperature
- Deflation as a Method of Variance Reduction for Estimating the Trace of a Matrix Inverse
- Estimating the Trace of the Matrix Inverse by Interpolating from the Diagonal of an Approximate Inverse
- Deflated and restarted symmetric Lanczos methods for eigenvalues and linear equations with multiple right-hand sides
- Extending the eigCG algorithm to nonsymmetric Lanczos for linear systems with multiple right-hand sides