Abstract:We investigate the SPAI and PSAI preconditioning procedures and shed light on two important features of them: (i) For the large linear system Ax=b with A irregular sparse, i.e., with A having s relatively dense columns, SPAI may be very costly to implement, and the resulting sparse approximate inverses may be ineffective for preconditioning. PSAI can be effective for preconditioning but may require excessive storage and be unacceptably time consuming; (ii) the situation is improved drastically when A is regular sparse, that is, all of its columns are sparse. In this case, both SPAI and PSAI are efficient. Moreover, SPAI and, especially, PSAI are more likely to construct effective preconditioners. Motivated by these features, we propose an approach to making SPAI and PSAI more practical for Ax=b with A irregular sparse. We first split A into a regular sparse Ã and a matrix of low rank s. Then exploiting the Sherman–Morrison–Woodbury formula, we transform Ax=b into s+1 new linear systems with the same coefficient matrix Ã, use SPAI and PSAI to compute sparse approximate inverses of Ã efficiently and apply Krylov iterative methods to solve the preconditioned linear systems. Theoretically, we consider the non-singularity and conditioning of Ã obtained from some important classes of matrices. We show how to recover an approximate solution of Ax=b from those of the s+1 new systems and how to design reliable stopping criteria for the s+1 systems to guarantee that the approximate solution of Ax=b satisfies a desired accuracy. Given the fact that irregular sparse linear systems are common in applications, this approach widely extends the practicability of SPAI and PSAI. Numerical results demonstrate the considerable superiority of our approach to the direct application of SPAI and PSAI to Ax=b.

The Application of Sparse Approximate Inverse Preconditioner in the Fast Solution of Volume-Surface Integral Equation

Improved SAI Preconditioning of Hybridized MLFMA for Electrically Large Multiscale Problems

Improved SAI Preconditioning of Hybridized MLFMA for Electrically Large Multi-scale Problems

A power sparse approximate inverse preconditioning procedure for large sparse linear systems

A Residual Based Sparse Approximate Inverse Preconditioning Procedure for Large Sparse Linear Systems

Accuracy Improvement of the Algebraic Fast Methods for the Volume-Surface Integral Equation

On the Preconditioners for Fast Multipole Boundary Element Methods for 2D Multi-Domain Elastostatics

Projection Improved SPAI Preconditioner for FGMRES

A -Preconditioner for Space Fractional Diffusion Equation with Non-separable Variable Coefficients

A New Preconditioning Scheme for MLFMA Based on Null-Field Generation Technique

Hierarchical Matrices Method and Its Application in Electromagnetic Integral Equations

A Simple Way of Applying the Calderón Preconditioner to a Finite-Element Boundary-Integral Method

Mixed Precision Iterative Refinement with Adaptive Precision Sparse Approximate Inverse Preconditioning

N A ] 2 7 A pr 2 01 5 A Residual Based Sparse Approximate Inverse Preconditioning Procedure for Large Sparse Linear Systems ∗

Multilevel Field-Split Preconditioners with Domain Decomposition for Steady and Unsteady Flow Problems

An Approach to Making SPAI and PSAI Preconditioning Effective for Large Irregular Sparse Linear Systems.

Generating Approximate Inverse Preconditioners for Sparse Matrices Using CUDA and GPGPU

A preconditioned iteration method for solving saddle point problems

On the efficient preconditioning of the Stokes equations in tight geometries

Generative modeling of Sparse Approximate Inverse Preconditioners

Approximate inner solvers for block preconditioning of the incompressible Navier–Stokes problems discretized by isogeometric analysis