Abstract:In this paper, we propose an adaptive fast solver for a general class of symmetric positive definite (SPD) matrices which include the well-known graph Laplacian. We achieve this by developing an adaptive operator compression scheme and a multiresolution matrix factorization algorithm which achieve nearly optimal performance on both complexity and well-posedness. To develop our adaptive operator compression and multiresolution matrix factorization methods, we first introduce a novel notion of energy decomposition for SPD matrix $A$ using the representation of energy elements. The interaction between these energy elements depicts the underlying topological structure of the operator. This concept of decomposition naturally reflects the hidden geometric structure of the operator which inherits the localities of the structure. By utilizing the intrinsic geometric information under this Energy framework, we propose a systematic operator compression scheme for the inverse operator $A^{-1}$. In particular, with an appropriate partition of the underlying geometric structure, we can construct localized basis by using the concept of interior and closed energy. Meanwhile, two important localized quantities are introduced, namely the error factor and the condition factor. Our error analysis results show that these two factors will be the guidelines for finding the appropriate partition of the basis functions such that prescribed compression error and acceptable condition number can be achieved. By virtue of this insight, we propose the Patch Pairing algorithm to realize our energy partition framework for operator compression with controllable compression error and condition number.

Weighted Non-Linear Criterion-Based Adaptive Generalised Eigendecomposition

Adaptive Algorithms for Generalized Eigenvalue Decomposition with A Nonquadratic Criterion

Weighted Rule Based Adaptive Algorithm for Simultaneously Extracting Generalized Eigenvectors

RLS-based Adaptive Algorithms for Generalized Eigen-Decomposition

An Adaptive Quasi-Newton Algorithm for Generalized Eigen-Decomposition Problem

Robust Adaptive Modified Newton Algorithm for Generalized Eigendecomposition and Its Application

Online Generalized Eigenvectors Extraction Via A Fixed-Point Approach

Low Complexity Adaptive Algorithm for Generalized Eigenvalue Decomposition

Online Dominant Generalized Eigenvectors Extraction Via a Randomized Algorithm

Adaptive quasi-newton algorithm for extracting generalized eigen subspace

Fast Adaptive Extraction Algorithm for Multiple Principal Generalized Eigenvectors

A Power-Based Adaptive Method for Eigenanalysis Without Square-Root Operations

Penalized Orthogonal Iteration for Sparse Estimation of Generalized Eigenvalue Problem

An Adaptive Gradient Method for Computing Generalized Tensor Eigenpairs

A Decomposition Algorithm for the Sparse Generalized Eigenvalue Problem.

A Weighted GMRES Method Based on Polynomial Preconditioning Generalized Minimal Residual Method

WEIGHTED GOLUB-KAHAN-LANCZOS BIDIAGONALIZATION ALGORITHMS

An adaptive fast solver for a general class of positive definite matrices via energy decomposition

Fast adaptive algorithm to extract multiple principal generalized eigenvectors

AN ADAPTIVE TRUST-REGION METHOD FOR GENERALIZED EIGENVALUES OF SYMMETRIC TENSORS

A Riemannian Optimization Approach for Solving the Generalized Eigenvalue Problem for Nonsquare Matrix Pencils.