Abstract:The deflated block conjugate gradient (D-BCG) method is an attractive approach for the solution of symmetric positive definite linear systems with multiple right-hand sides. However, the orthogonality between the block residual vectors and the deflation subspace is gradually lost along with the process of the underlying algorithm implementation, which usually causes the algorithm to be unstable or possibly have delayed convergence. In order to maintain such orthogonality to keep certain level, full reorthogonalization could be employed as a remedy, but the expense required is quite costly. In this paper, we present a new projected variant of the deflated block conjugate gradient (PD-BCG) method to mitigate the loss of this orthogonality, which is helpful to deal with the delay of convergence and thus further achieve the theoretically faster convergence rate of D-BCG. Meanwhile, the proposed PD-BCG method is shown to scarcely have any extra computational cost, while having the same theoretical properties as D-BCG in exact arithmetic. Additionally, an automated reorthogonalization strategy is introduced as an alternative choice for the PD-BCG method. Numerical experiments demonstrate that PD-BCG is more efficient than its counterparts especially when solving ill-conditioned linear systems or linear systems suffering from rank deficiency.

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

Restrictive Preconditioners for Conjugate Gradient Methods for Symmetric Positive Definite Linear Systems

Additive block diagonal preconditioning for block two-by-two linear systems of skew-Hamiltonian coefficient matrices

On a type of matrix splitting preconditioners for a class of block two-by-two linear systems.

RESTRICTIVELY PRECONDITIONED CONJUGATE GRADIENT METHOD FOR A SERIES OF CONSTANTLY AUGMENTED LEAST SQUARES PROBLEMS

Restrictively Preconditioned Chebyshev Method for Solving Systems of Linear Equations

Modified restrictive preconditioners for double saddle point problems arising from liquid crystal director modeling

Preconditioned MHSS Iteration Methods for a Class of Block Two-by-two Linear Systems with Applications to Distributed Control Problems

Block [math]-Circulant Preconditioners for Parabolic Optimal Control Problems

On a Chebyshev Accelerated Splitting Iteration Method with Application to Two‐by‐two Block Linear Systems

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

BCCB preconditioners for solving linear systems from delay differential equations

A preconditioner for block two-by-two symmetric indefinite matrices

On Hybrid Preconditioning Methods for Large Sparse Saddle-Point Problems

A Preconditioned MINRES Method for Block Lower Triangular Toeplitz Systems

A Note on Constraint Preconditioners for Nonsymmetric Saddle Point Problems

Constraint Preconditioners for Symmetric Indefinite Matrices

Block Variants of the Cocg and Cocr Methods for Solving Complex Symmetric Linear Systems with Multiple Right-Hand Sides

A New Projected Variant of the Deflated Block Conjugate Gradient Method

A Generalization of T. Chan’s Preconditioner

Block-circulant with circulant-block preconditioners for two-dimensional spatial fractional diffusion equations