Abstract:Deflation techniques are typically used to shift isolated clusters of small eigenvalues in order to obtain a tighter distribution and a smaller condition number. Such changes induce a positive effect in the convergence behavior of Krylov subspace methods, which are among the most popular iterative solvers for large sparse linear systems. We develop a deflation strategy for symmetric saddle point matrices by taking advantage of their underlying block structure. The vectors used for deflation come from an elliptic singular value decomposition relying on the generalized Golub-Kahan bidiagonalization process. The block targeted by deflation is the off-diagonal one since it features a problematic singular value distribution for certain applications. One example is the Stokes flow in elongated channels, where the off-diagonal block has several small, isolated singular values, depending on the length of the channel. Applying deflation to specific parts of the saddle point system is important when using solvers such as CRAIG, which operates on individual blocks rather than the whole system. The theory is developed by extending the existing framework for deflating square matrices before applying a Krylov subspace method like MINRES. Numerical experiments confirm the merits of our strategy and lead to interesting questions about using approximate vectors for deflation.

Review of Cyclic Reduction for Parallel Solution of Hermitian Positive Definite Block-Tridiagonal Linear Systems

Generalized Golub-Kahan bidiagonalization for nonsymmetric saddle point systems

Positive matrices partitioned into a small number of Hermitian blocks

On the global convergence of the block Jacobi method for the positive definite generalized eigenvalue problem

Convergence of the Cyclic and Quasi-cyclic Block Jacobi Methods

Computing the Hermitian Positive Definite Solutions of a Nonlinear Matrix Equation

Hermitian Preconditioning for a class of Non-Hermitian Linear Systems

An efficient numerical algorithm for solving linear systems with cyclic tridiagonal coefficient matrices

Solving system of nonlinear matrix equations over Hermitian positive definite matrices

A Parallel Iterative Method For Solving Symmetrical Tridiagonal Extreme Eigenpair Problems

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

An Inverse Extreme Eigenpair Problem and Its Parallel Iterative Solution for Symmetric Tridiagonal Matrices

A Hybrid Parallel Tridiagonal Solver on Multi-core Architectures

Parallel incomplete factorization preconditioning of block tridiagonal linear systems with 2-D domain decomposition

Deflation for the off-diagonal block in symmetric saddle point systems

Solutions of a class of nonlinear matrix equations

On the Cholesky method

Generalized Block Lanczos Methods for Large Unsymmetric Eigenproblems

A parallel iterative method for solving symmetric tridiagonal extreme eigenpair problems

Practical Restrictively Preconditioned Conjugate Gradient Methods for a Class of Block Two-by-Two Linear Systems