Tree quasi-separable matrices: a simultaneous generalization of sequentially and hierarchically semi-separable representations
Nithin Govindarajan, Shivkumar Chandrasekaran, Patrick Dewilde
TL;DR
This work introduces tree quasi-separable (TQS) matrices as a unifying generalization of sequentially and hierarchically semi-separable representations for matrices whose graph is a tree. By modeling a graph-partitioned matrix with an input–output dynamical system along the tree, the authors derive a minimal, universal construction where generator sizes are dictated by Hankel-block ranks, ensuring TQS inherits the fast matvec and direct-solve capabilities of SSS and HSS. They prove a graph-induced rank structure (GIRS) property for TQS, show equivalence with minimal Hankel-rank-based representations, and demonstrate that TQS operations (addition, multiplication, inversion) preserve low-rank structure with predictable rank growth. The paper also provides an explicit construction algorithm, numerical validation, and a framework for efficient linear-system solving on trees, suggesting TQS as a flexible drop-in replacement with potential advantages for graph-based dynamical systems and network problems.
Abstract
We present a unification and generalization of what is known in the literature as sequentially and hierarchically semi-separable (SSS and HSS) representations for matrices. Describing rank-structured representations of (inverses of) sparse matrices whose adjacency graph is a tree, it is shown that these so-called tree quasi-separable (TQS) matrices inherit all the favorable algebraic properties of SSS and HSS under addition, products, and inversion. To arrive at these properties, we prove a key result that characterizes the conversion of any dense matrix into a TQS representation. Here, we specifically show through an explicit construction procedure that the generator sizes are dictated by the ranks of certain Hankel blocks of the matrix. Analogous to SSS and HSS, TQS matrices admit fast matrix-vector products and direct solvers provided the generator sizes are small. A sketch of the associated algorithms is provided.
