Abstract:Suppose that $A \in \mathbb{R}^{N \times N}$ is symmetric positive semidefinite with rank $K \le N$. Our goal is to decompose $A$ into $K$ rank-one matrices $\sum_{k=1}^K g_k g_k^T$ where the modes $\{g_{k}\}_{k=1}^K$ are required to be as sparse as possible. In contrast to eigen decomposition, these sparse modes are not required to be orthogonal. Such a problem arises in random field parametrization where $A$ is the covariance function and is intractable to solve in general. In this paper, we partition the indices from 1 to $N$ into several patches and propose to quantify the sparseness of a vector by the number of patches on which it is nonzero, which is called patch-wise sparseness. Our aim is to find the decomposition which minimizes the total patch-wise sparseness of the decomposed modes. We propose a domain-decomposition type method, called intrinsic sparse mode decomposition (ISMD), which follows the "local-modes-construction + patching-up" procedure. The key step in the ISMD is to construct local pieces of the intrinsic sparse modes by a joint diagonalization problem. Thereafter a pivoted Cholesky decomposition is utilized to glue these local pieces together. Optimal sparse decomposition, consistency with different domain decomposition and robustness to small perturbation are proved under the so called regular-sparse assumption (see Definition 1.2). We provide simulation results to show the efficiency and robustness of the ISMD. We also compare the ISMD to other existing methods, e.g., eigen decomposition, pivoted Cholesky decomposition and convex relaxation of sparse principal component analysis [25] and [40].

Symmetry Breaking in Symmetric Tensor Decomposition

Symmetry & Critical Points for Symmetric Tensor Decomposition Problems

Symmetric Tensor Decompositions On Varieties

Rank Decomposition and Symmetric Rank Decomposition over Arbitrary Fields

Symmetry: a General Structure in Nonparametric Regression

Iterative Methods for Symmetric Outer Product Tensor Decompositions

Desingularization of bounded-rank tensor sets

A sparse decomposition of low rank symmetric positive semi-definite matrices

SySTeC: A Symmetric Sparse Tensor Compiler

Stochastic Mirror Descent for Low-Rank Tensor Decomposition Under Non-Euclidean Losses

Computing symmetric rank for symmetric tensors

Symmetry Breaking and Equivariant Neural Networks

Decomposition of linear tensor transformations

Approximate Real Symmetric Tensor Rank

Symmetry Breaking in Neural Network Optimization: Insights from Input Dimension Expansion

On the ranks and border ranks of symmetric tensors

A chiseling algorithm for low-rank Grassmann decomposition of skew-symmetric tensors

General Low-rank Matrix Optimization: Geometric Analysis and Sharper Bounds

Computational complexity of decomposing a symmetric matrix as a sum of positive semidefinite and diagonal matrices

Scalable Symmetric Tucker Tensor Decomposition

High-order sum-of-squares structured tensors: theory and applications