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.

A single-sided all-at-once preconditioning for linear system from a non-local evolutionary equation with weakly singular kernels

A single-sided all-at-once preconditioning for linear system from a non-local evolutionary equation with weakly singular kernels

A bilateral preconditioning for an L2-type all-at-once system from time-space non-local evolution equations with a weakly singular kernel

Optimal preconditioners for nonsymmetric multilevel Toeplitz systems with application to solving non-local evolutionary partial differential equations

An efficient preconditioner for evolutionary partial differential equations with $θ$-method in time discretization

A preconditioned MINRES method for block lower triangular Toeplitz systems

A parallel preconditioner for the all-at-once linear system from evolutionary PDEs with Crank-Nicolson discretization

A Note on Constraint Preconditioners for Nonsymmetric Saddle Point Problems

A fast normal splitting preconditioner for attractive coupled nonlinear Schrödinger equations with fractional Laplacian

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

A preconditioning technique for all-at-once system from the nonlinear tempered fractional diffusion equation

On Hybrid Preconditioning Methods for Large Sparse Saddle-Point Problems

Constraint Preconditioners for Symmetric Indefinite Matrices

Block-splitting preconditioners for indefinite least squares problem

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

Nonlinearly preconditioned semismooth Newton algorithms for nonlinear nonsmooth systems

Preconditioned iterative method for nonsymmetric saddle point linear systems

Preconditioning Technique for Symmetric M-matrices

A Class of Preconditioners Based on the (I+s(alpha))-Type Preconditioning Matrices for Solving Linear Systems

Preconditioning correction equations in Jacobi--Davidson type methods for computing partial singular value decompositions of large matrices