Deflated GMRES for Systems with Multiple Shifts and Multiple Right-Hand Sides
arXiv:0707.0502 · doi:10.1016/j.laa.2008.04.019
Abstract
We consider solution of multiply shifted systems of nonsymmetric linear equations, possibly also with multiple right-hand sides. First, for a single right-hand side, the matrix is shifted by several multiples of the identity. Such problems arise in a number of applications, including lattice quantum chromodynamics where the matrices are complex and non-Hermitian. Some Krylov iterative methods such as GMRES and BiCGStab have been used to solve multiply shifted systems for about the cost of solving just one system. Restarted GMRES can be improved by deflating eigenvalues for matrices that have a few small eigenvalues. We show that a particular deflated method, GMRES-DR, can be applied to multiply shifted systems. In quantum chromodynamics, it is common to have multiple right-hand sides with multiple shifts for each right-hand side. We develop a method that efficiently solves the multiple right-hand sides by using a deflated version of GMRES and yet keeps costs for all of the multiply shifted systems close to those for one shift. An example is given showing this can be extremely effective with a quantum chromodynamics matrix.
19 pages, 9 figures
References in corpus (9)
- Numerical Methods for the QCD Overlap Operator: I. Sign-Function and Error Bounds
- A study of chiral symmetry in quenched QCD using the Overlap-Dirac operator
- On the low fermionic eigenmode dominance in QCD on the lattice
- Chiral Symmetry, Quark Mass, and Scaling of the Overlap Fermions
- Deflated Iterative Methods for Linear Equations with Multiple Right-Hand Sides
- An alternative to domain wall fermions
- Progress on lattice QCD algorithms
- Alternative to Domain Wall Fermions
- Deflation of Eigenvalues for Iterative Methods in Lattice QCD
Cited by in corpus (17)
- Effective noise reduction techniques for disconnected loops in Lattice QCD
- Light Hadron Masses from Lattice QCD
- A framework for deflated and augmented Krylov subspace methods
- Krylov Subspace Recycling for Sequences of Shifted Linear Systems
- Application of block Krylov subspace algorithms to the Wilson-Dirac equation with multiple right-hand sides in lattice QCD
- A Flexible Krylov Solver for Shifted Systems with Application to Oscillatory Hydraulic Tomography
- Uncertainty Quantification in Lattice QCD Calculations for Nuclear Physics
- Block Krylov subspace recycling for shifted systems with unrelated right-hand sides
- Restarted Hessenberg method for solving shifted nonsymmetric linear systems
- GMRES algorithms over 35 years
- Extending the eigCG algorithm to nonsymmetric Lanczos for linear systems with multiple right-hand sides
- A flexible and adaptive Simpler GMRES with deflated restarting for shifted linear systems
- Efficient variants of the CMRH method for solving a sequence of multi-shifted non-Hermitian linear systems simultaneously
- Application of preconditioned block BiCGGR to the Wilson-Dirac equation with multiple right-hand sides in lattice QCD
- Reduced-basis method for the iterative solution of parametrized symmetric positive-definite linear systems
- Preconditioning the Restarted and Shifted Block FOM Algorithm for Matrix Exponential Computation
- Numerical experiments using deflation with the HISQ action