Abstract:A sparse matrix is called double irregular sparse if it has at least one relatively dense column and row, and it is double regular sparse if all the columns and rows of it are sparse. The sparse approximate inverse preconditioning procedures SPAI, PSAI( t o l ) and RSAI( t o l ) are costly and even impractical to construct preconditioners for a large sparse nonsymmetric linear system with the coefficient matrix being double irregular sparse, but they are efficient for double regular sparse problems. Double irregular sparse linear systems have a wide range of applications, and 24.4% of the nonsymmetric matrices in the Florida University collection are double irregular sparse. For this class of problems, we propose a transformation approach, which consists of four steps: (i) transform a given double irregular sparse problem into a small number of double regular sparse ones with the same coefficient matrix A ˆ , (ii) use SPAI, PSAI( t o l ) and RSAI( t o l ) to construct sparse approximate inverses M of A ˆ , (iii) solve the preconditioned double regular sparse linear systems by Krylov solvers, and (iv) recover an approximate solution of the original problem with a prescribed accuracy from those of the double regular sparse ones. A number of theoretical and practical issues are considered on the transformation approach. Numerical experiments on a number of real-world problems confirm the very sharp superiority of the transformation approach to the standard approach that preconditions the original double irregular sparse problem by SPAI, PSAI( t o l ) or RSAI( t o l ) and solves the resulting preconditioned system by Krylov solvers.

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

Prestructuring sparse matrices with dense rows and columns via null space methods

A Note on Constraint Preconditioners for Nonsymmetric Saddle Point Problems

Deflation for the Off-Diagonal Block in Symmetric Saddle Point Systems

Preconditioned iterative method for nonsymmetric saddle point linear systems

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

A randomized small-block Lanczos method for large-scale null space computations

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

Extended Newton-type Method for Inverse Singular Value Problems with Multiple And/or Zero Singular Values

Generalized Golub-Kahan bidiagonalization for nonsymmetric saddle point systems

A Non-Krylov subspace Method for Solving Large and Sparse Linear System of Equations

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

A transformation approach that makes SPAI, PSAI and RSAI procedures efficient for large double irregular nonsymmetric sparse linear systems

An Efficient Gauss--Newton Algorithm for Symmetric Low-Rank Product Matrix Approximations

On the Unreasonable Effectiveness of Single Vector Krylov Methods for Low-Rank Approximation

Large-Scale Optimization with Linear Equality Constraints Using Reduced Compact Representation

Solving Dense Linear Systems Faster Than via Preconditioning

Iterative Methods for Sparse Symmetric Multilinear Systems

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

A Computational Study of Using Black-box QR Solvers for Large-scale Sparse-dense Linear Least Squares Problems