Abstract:This work proposes a new class of preconditioners for the low rank Generalized Minimal Residual Method (GMRES) for multiterm matrix equations arising from implicit timestepping of linear matrix differential equations. We are interested in computing low rank solutions to matrix equations, e.g. arising from spatial discretization of stiff partial differential equations (PDEs). The low rank GMRES method is a particular class of Krylov subspace method where the iteration is performed on the low rank factors of the solution. Such methods can exploit the low rank property of the solution to save on computational and storage cost. Of critical importance for the efficiency and applicability of the low rank GMRES method is the availability of an effective low rank preconditioner that operates directly on the low rank factors of the solution and that can limit the iteration count and the maximal Krylov rank. The preconditioner we propose here is based on the basis update and Galerkin (BUG) method, resulting from the dynamic low rank approximation. It is a nonlinear preconditioner for the low rank GMRES scheme that naturally operates on the low rank factors. Extensive numerical tests show that this new preconditioner is highly efficient in limiting iteration count and maximal Krylov rank. We show that the preconditioner performs well for general diffusion equations including highly challenging problems, e.g. high contrast, anisotropic equations. Further, it compares favorably with the state of the art exponential sum preconditioner. We also propose a hybrid BUG - exponential sum preconditioner based on alternating between the two preconditioners.

Convergence analysis of inner-iteration preconditioned GMRES

The Convergence of Generalized Lanczos Methods for Large Unsymmetric Eigenproblems

Hermitian Preconditioning for a class of Non-Hermitian Linear Systems

New Convergence Analysis of GMRES with Weighted Norms, Preconditioning and Deflation, Leading to a New Deflation Space

A block quaternion GMRES method and its convergence analysis

The Convergence of Krylov Subspace Methods for Large Unsymmetric Linear Systems

Convergence on Gauss-Seidel Iterative Methods for Linear Systems with General H-matrices

Utilizing Cuda For Preconditioned Gmres Solvers

Convergence Properties of Nonlinear Conjugate Gradient Methods.

Convergence analysis of inexact proximal point algorithms on Hadamard manifolds

New Convergence Analysis of GMRES with Weighted Norms, Preconditioning, and Deflation, Leading to a New Deflation Space

Convergence Analysis of Gradient Algorithms on Riemannian Manifolds Without Curvature Constraints and Application to Riemannian Mass

Preconditioning Low Rank Generalized Minimal Residual Method (GMRES) for Implicit Discretizations of Matrix Differential Equations

The Convergence Discussion of the GSOR Iterative Methods with Two Preconditioners

Convergence analysis of iterative solution with inexact block preconditioning for weak Galerkin finite element approximation of Stokes flow

Preprocessed GMRES for fast solution of linear equations

Field of values analysis of preconditioners for the Helmholtz equation in lossy media

Convergence Analysis of Kernel Conjugate Gradient for Functional Linear Regression

Analysis of a multigrid preconditioner for Crouzeix-Raviart discretization of elliptic PDE with jump coefficient

Some Investigations on Convergence of GMRES in Solving BEM Equations for Slender Beam Structures

Improving performance of contour integral-based nonlinear eigensolvers with infinite GMRES