Cholesky decomposition for symmetric matrices, Riemannian geometry, and random matrices

Abstract

For each n ≥ 1 and sign pattern ε ∈ \ 1 \n, we introduce a cone of real symmetric matrices LPMn(ε): those with leading principal k × k minors of signs εk. These cones are pairwise disjoint and their union LPMn is an open dense cone in all symmetric matrices; they subsume positive and negative definite matrices, and symmetric (P-,) N-, PN-, almost P-, and almost N- matrices. We show that each LPMn matrix A admits an uncountable family of Cholesky-type factorizations - yielding a unique lower triangular matrix L with positive diagonals - with additional attractive properties: (i) each such factorization is algorithmic; and (ii) each such Cholesky map A L is a smooth diffeomorphism from LPMn(ε) onto an open Euclidean ball. We then show that (iii) the (diffeomorphic) balls LPMn(ε) are isometric Riemannian manifolds as well as isomorphic abelian Lie groups, each equipped with a translation-invariant Riemannian metric (and hence Riemannian means/barycentres). Moreover, (iv) this abelian metric group structure on each LPMn(ε) - and hence the log-Cholesky metric on Cholesky space - yields an isometric isomorphism onto a finite-dimensional Euclidean space. The complex version of this also holds. In the latter part, we show that the abelian group PDn of positive definite matrices, with its bi-invariant log-Cholesky metric, is precisely the identity-component of a larger group with an alternate metric: the open dense cone LPMn. This also holds for Hermitian matrices over several subfields F ⊂eq C. As a result, (v) the groups LPMnF and LPM∞F admit a rich probability theory, and the cones LPMn(ε), TPMn(ε) admit Wishart densities with signed Bartlett decompositions.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…