Abstract:Motivated by the Cayley-Hamilton theorem, a novel adaptive procedure, called a Power Sparse Approximate Inverse (PSAI) procedure, is proposed that uses a different adaptive sparsity pattern selection approach to constructing a right preconditioner M for the large sparse linear system Ax = b. It determines the sparsity pattern of M dynamically and computes the n independent columns of M that is optimal in the Frobenius norm minimization, subject to the sparsity pattern of M. The PSAI procedure needs a matrix-vector product at each step and updates the solution of a small least squares problem cheaply. To control the sparsity of M and develop a practical PSAI algorithm, two dropping strategies are proposed. The PSAI algorithm can capture an effective approximate sparsity pattern of A(-1) and compute a good sparse approximate inverse M efficiently. Numerical experiments are reported to verify the effectiveness of the PSAI algorithm. Numerical comparisons are made for the PSAI algorithm and the adaptive SPAI algorithm proposed by Grote and Huckle as well as for the PSAI algorithm and three static Sparse Approximate Inverse (SAI) algorithms. The results indicate that the PSAI algorithm is at least comparable to and can be much more effective than the adaptive SPAI algorithm and it often outperforms the static SAI algorithms very considerably and is more robust and practical than the static ones for general problems. Copyright (C) 2008 John Wiley & Sons, Ltd.

A filtering technique for the matrix power series being near-sparse

Efficient Computational Method for Matrix Function in Dynamic Problems

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

High-Performance Computation of the Exponential of a Large Sparse Matrix

The Improved Signal Processing Method Based On Matrix Pencil Algorithm

Mixed-Precision Paterson--Stockmeyer Method for Evaluating Polynomials of Matrices

A New Stable and Avoiding Inversion Iteration for Computing Matrix Square Root

Two Iterative algorithms for the matrix sign function based on the adaptive filtering technology

Efficient And Memory Saving Method Based On Pseudoskeleton Approximation For Analysis Of Finite Periodic Structures

Low-Rank Matrix Iteration Using Polynomial-Filtered Subspace Extraction.

Projected Kernel Least Mean $p$ -Power Algorithm: Convergence Analyses and Modifications.

A Unified Precision Matrix Estimation Framework Via Sparse Column-Wise Inverse Operator under Weak Sparsity

A sparse-sparse iteration for computing a sparse incomplete factorization of the inverse of an SPD matrix

State space least p-power filter.

Sparsification of the Impedance Matrix in the Solution of the Integral Equation by Using the Maximally Orthogonalized Basis Functions

Near-Field Sparse Inverse Preconditioning of Multilevel Fast Multipole Algorithm for Electric Field Integral Equations

A preconditioned iteration method for solving saddle point problems

Approximating Sparse Matrices and their Functions using Matrix-vector products

A Fast Random Sampling Algorithm for Sparsifying Matrices

Near Optimal Compressed Sensing of a Class of Sparse Low-Rank Matrices via Sparse Power Factorization

An Efficient Two-Step Iterative Method for Solving a Class of Complex Symmetric Linear Systems