Polynomial preconditioned Arnoldi with stability control
From MaRDI portal
Publication:5856683
DOI10.1137/19M1302430zbMATH Open1464.65035arXiv1806.08020OpenAlexW3119186146MaRDI QIDQ5856683FDOQ5856683
Authors: Mark Embree, Jennifer Loe, Ronald B. Morgan
Publication date: 29 March 2021
Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)
Abstract: Polynomial preconditioning can improve the convergence of the Arnoldi method for computing eigenvalues. Such preconditioning significantly reduces the cost of orthogonalization; for difficult problems, it can also reduce the number of matrix-vector products. Parallel computations can particularly benefit from the reduction of communication-intensive operations. The GMRES algorithm provides a simple and effective way of generating the preconditioning polynomial. For some problems high degree polynomials are especially effective, but they can lead to stability problems that must be mitigated. A two-level "double polynomial preconditioning" strategy provides an effective way to generate high-degree preconditioners.
Full work available at URL: https://arxiv.org/abs/1806.08020
Recommendations
Eigenvalues, singular values, and eigenvectors (15A18) Numerical computation of eigenvalues and eigenvectors of matrices (65F15)
Cites Work
- The University of Florida sparse matrix collection
- Implicit Application of Polynomial Filters in a k-Step Arnoldi Method
- GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
- Numerical methods for large eigenvalue problems
- Title not available (Why is that?)
- Computing interior eigenvalues of large matrices
- Title not available (Why is that?)
- Title not available (Why is that?)
- Solution of Sparse Indefinite Systems of Linear Equations
- Ritz Value Localization for Non-Hermitian Matrices
- Variations on Arnoldi's method for computing eigenelements of large unsymmetric matrices
- Any Ritz value behavior is possible for Arnoldi and for GMRES
- A harmonic restarted Arnoldi algorithm for calculating eigenvalues and determining multiplicity
- Thick-restart Lanczos method for large symmetric eigenvalue problems
- Title not available (Why is that?)
- Approximate solutions and eigenvalue bounds from Krylov subspaces
- A Krylov--Schur algorithm for large eigenproblems
- Chebyshev Acceleration Techniques for Solving Nonsymmetric Eigenvalue Problems
- On restarting the Arnoldi method for large nonsymmetric eigenvalue problems
- Least Squares Polynomials in the Complex Plane and Their Use for Solving Nonsymmetric Linear Systems
- A Hybrid GMRES Algorithm for Nonsymmetric Linear Systems
- A Newton basis GMRES implementation
- The Arnoldi eigenvalue iteration with exact shifts can fail
- Title not available (Why is that?)
- Convergence of Polynomial Restart Krylov Methods for Eigenvalue Computations
- Convergence of Restarted Krylov Subspaces to Invariant Subspaces
- Dynamic Thick Restarting of the Davidson, and the Implicitly Restarted Arnoldi Methods
- Newton interpolation at Leja points
- Stability of polynomial preconditioning
- Harmonic Ritz and Lehmann bounds
- Shift-invert and Cayley transforms for detection of rightmost eigenvalues of nonsymmetric matrices
- Polynomial Preconditioned GMRES and GMRES-DR
- Quasi-kernel polynomials and their use in non-Hermitian matrix iterations
- A thick-restart Lanczos algorithm with polynomial filtering for Hermitian eigenvalue problems
Cited In (7)
- Stability of polynomial preconditioning
- Polynomial preconditioning for the action of the matrix square root and inverse square root
- Proxy-GMRES: Preconditioning via GMRES in Polynomial Space
- Parallel matrix-free polynomial preconditioners with application to flow simulations in discrete fracture networks
- Coarsest-level improvements in multigrid for lattice QCD on large-scale computers
- GMRES algorithms over 35 years
- Parallelization of the rational Arnoldi algorithm
Uses Software
This page was built for publication: Polynomial preconditioned Arnoldi with stability control
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5856683)