Multigrid deflation for lattice QCD
From MaRDI portal
Publication:778340
DOI10.1016/J.JCP.2020.109356zbMATH Open1435.65058arXiv1909.12234OpenAlexW3006815883MaRDI QIDQ778340FDOQ778340
Eloy Romero, Kostas Orginos, Andreas Stathopoulos
Publication date: 2 July 2020
Published in: Journal of Computational Physics (Search for Journal in Brave)
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.
Full work available at URL: https://arxiv.org/abs/1909.12234
Recommendations
- Deflation as a method of variance reduction for estimating the trace of a matrix inverse
- Computing and deflating eigenvalues while solving multiple right-hand side linear systems with an application to quantum chromodynamics
- A Multilevel Approach to Variance Reduction in the Stochastic Estimation of the Trace of a Matrix
- scientific article
- Hierarchical probing for estimating the trace of the matrix inverse on toroidal lattices
Monte Carlo methods (65C05) Quantum field theory on lattices (81T25) Numerical solutions to inverse eigenvalue problems (65F18)
Cites Work
- PRIMME
- Templates for the Solution of Algebraic Eigenvalue Problems
- A stochastic estimator of the trace of the influence matrix for laplacian smoothing splines
- A probing method for computing the diagonal of a matrix inverse
- Hierarchical probing for estimating the trace of the matrix inverse on toroidal lattices
- Preconditioning eigensolvers -- an Oxymoron?
- Computing and Deflating Eigenvalues While Solving Multiple Right-Hand Side Linear Systems with an Application to Quantum Chromodynamics
- Deflated and Restarted Symmetric Lanczos Methods for Eigenvalues and Linear Equations with Multiple Right-Hand Sides
- Nearly Optimal Preconditioned Methods for Hermitian Eigenproblems under Limited Memory. Part I: Seeking One Eigenvalue
- Monte Carlo methods for estimating linear combinations of inverse matrix entries in lattice QCD
- Effective noise reduction techniques for disconnected loops in lattice QCD
- Estimating the trace of the matrix inverse by interpolating from the diagonal of an approximate inverse
- Extending the eigCG algorithm to nonsymmetric Lanczos for linear systems with multiple right-hand sides
- Deflation as a Method of Variance Reduction for Estimating the Trace of a Matrix Inverse
Cited In (3)
Uses Software
This page was built for publication: Multigrid deflation for lattice QCD
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q778340)
