Abstract:Randomized methods are becoming increasingly popular in numerical linear algebra. However, few attempts have been made to use them in developing preconditioners. Our interest lies in solving large-scale sparse symmetric positive definite linear systems of equations where the system matrix is preordered to doubly bordered block diagonal form (for example, using a nested dissection ordering). We investigate the use of randomized methods to construct high quality preconditioners. In particular, we propose a new and efficient approach that employs Nystr\"om's method for computing low rank approximations to develop robust algebraic two-level preconditioners. Construction of the new preconditioners involves iteratively solving a smaller but denser symmetric positive definite Schur complement system with multiple right-hand sides. Numerical experiments on problems coming from a range of application areas demonstrate that this inner system can be solved cheaply using block conjugate gradients and that using a large convergence tolerance to limit the cost does not adversely affect the quality of the resulting Nystr\"om--Schur two-level preconditioner.

Ideal Preconditioners for Saddle Point Systems with a Rank-Deficient Leading Block

A Note on Constraint Preconditioners for Nonsymmetric Saddle Point Problems

Preconditioned iterative method for nonsymmetric saddle point linear systems

A Block Positive-Semidefinite Splitting Preconditioner for Generalized Saddle Point Linear Systems.

A preconditioned iteration method for solving saddle point problems

Some preconditioning techniques for a class of double saddle point problems

A New Block Preconditioner for Double Saddle Point Systems Arising from Liquid Crystal Directors Modeling

Schur complement based preconditioners for twofold and block tridiagonal saddle point problems

Spectral analysis of block preconditioners for double saddle-point linear systems with application to PDE-constrained optimization

New preconditioner strategy for solving block four-by-four linear systems: An application to the saddle-point problem from 3D Stokes equation

A null-space approach for large-scale symmetric saddle point systems with a small and non zero (2, 2) block

Shift-Splitting Preconditioners for Saddle Point Problems

On the improvement of shift-splitting preconditioners for double saddle point problems

On the invertibility of matrices with a double saddle-point structure

Permuted preconditioning for extended saddle point problem arising from Neumann boundary control

A class of new extended generalized shift-splitting preconditioners for block saddle point problems

A generalized relaxed block positive-semidefinite splitting preconditioner for generalized saddle point linear system

Alternating positive semidefinite splitting preconditioners for double saddle point problems

On Preconditioned Generalized Shift-Splitting Iteration Methods for Saddle Point Problems.

Two-level Nyström-Schur preconditioner for sparse symmetric positive definite matrices