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.

On IGMRES: an Incomplete Generalized Minimal Residual Method for Large Unsymmetric Linear Systems

Theoretical and numerical comparisons of GMRES and WZ-GMRES

Obtaining the pseudoinverse solution of singular range-symmetric linear systems with GMRES-type methods

Generalized minimal residual method for systems with multiple right-hand sides

Inexact Global FOM and GMRES Methods for Solving Linear Systems with Multiple Right-Hand Sides

Restarted GMRES Method Augmented with the Combination of Harmonic Ritz Vectors and Error Approximations

Global GMRES with Deflated Restarting for Families of Shifted Linear Systems

Optimal Solutions of Well-Posed Linear Systems via Low-Precision Right-Preconditioned GMRES with Forward and Backward Stabilization

Regularization Properties of Krylov Iterative Solvers CGME and LSMR for Linear Discrete Ill-Posed Problems with an Application to Truncated Randomized SVDs

GMRES with randomized sketching and deflated restarting

Randomized flexible GMRES with deflated restarting

On the backward stability of s-step GMRES

Global FOM and GMRES algorithms for a class of complex matrix equations.

A spectrally preconditioned and initially deflated variant of the restarted block GMRES method for solving multiple right-hand sides linear systems

Some Recursions on Arnoldi's Method and IOM for Large Non-Hermitian Linear Systems

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

Solving Shifted Linear Systems with Restarted GMRES Augmented with Error Approximations.

Approximation Accuracy of the Krylov Subspaces for Linear Discrete Ill-Posed Problems

On IOM(q): the Incomplete Orthogonalization Method for Large Unsymmetric Linear Systems

Improving GMRES(m) Using an Adaptive Switching Controller

Block GMRES Method with Inexact Breakdowns and Deflated Restarting