A power Schur complement Low-Rank correction preconditioner for general sparse linear systems

02/03/2020
by   Qingqing Zheng, et al.
0

An effective power based parallel preconditioner is proposed for general large sparse linear systems. The preconditioner combines a power series expansion method with some low-rank correction techniques, where the Sherman-Morrison-Woodbury formula is utilized. A matrix splitting of the Schur complement is proposed to expand the power series. The number of terms used in the power series expansion can control the approximation accuracy of the preconditioner to the inverse of the Schur complement. To construct the preconditioner, graph partitioning is invoked to reorder the original coefficient matrix, leading to a special block two-by-two matrix whose two off-diagonal submatrices are block diagonal. Variables corresponding to interface variables are obtained by solving a linear system with the coeffcient matrix being the Schur complement. For the variables related to the interior variables, one only needs to solve a block diagonal linear system. This can be performed efficiently in parallel. Various numerical examples are provided to illustrate that the efficiency of the proposed preconditioner.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
08/12/2022

Solving Linear Systems on a GPU with Hierarchically Off-Diagonal Low-Rank Approximations

We are interested in solving linear systems arising from three applicati...
research
04/19/2022

The Binary Rank of Circulant Block Matrices

The binary rank of a 0,1 matrix is the smallest size of a partition of i...
research
05/04/2022

parGeMSLR: A Parallel Multilevel Schur Complement Low-Rank Preconditioning and Solution Package for General Sparse Matrices

This paper discusses parGeMSLR, a C++/MPI software library for the solut...
research
01/03/2020

A note on 2× 2 block-diagonal preconditioning

For 2x2 block matrices, it is well-known that block-triangular or block-...
research
10/14/2020

Hierarchical Orthogonal Factorization: Sparse Square matrices

In this work, we develop a new fast algorithm, spaQR – sparsified QR, fo...
research
11/26/2021

Symmetric Near-Field Schur Complement Preconditioner for Hierarchal Electric Field Integral Equation Solve

In this paper, a robust and effective preconditioner for the fast Method...
research
07/12/2017

Discriminative Block-Diagonal Representation Learning for Image Recognition

Existing block-diagonal representation researches mainly focuses on cast...

Please sign up or login with your details

Forgot password? Click here to reset