Abstract:Riemannian Gaussian distributions were initially introduced as basic building blocks for learning models which aim to capture the intrinsic structure of statistical populations of positive-definite matrices (here called covariance matrices). While the potential applications of such models have attracted significant attention, a major obstacle still stands in the way of these applications: there seems to exist no practical method of computing the normalising factors associated with Riemannian Gaussian distributions on spaces of high-dimensional covariance matrices. The present paper shows that this missing method comes from an unexpected new connection with random matrix theory. Its main contribution is to prove that Riemannian Gaussian distributions of real, complex, or quaternion covariance matrices are equivalent to orthogonal, unitary, or symplectic log-normal matrix ensembles. This equivalence yields a highly efficient approximation of the normalising factors, in terms of a rather simple analytic expression. The error due to this approximation decreases like the inverse square of dimension. Numerical experiments are conducted which demonstrate how this new approximation can unlock the difficulties which have impeded applications to real-world datasets of high-dimensional covariance matrices. The paper then turns to Riemannian Gaussian distributions of block-Toeplitz covariance matrices. These are equivalent to yet another kind of random matrix ensembles, here called "acosh-normal" ensembles. Orthogonal and unitary "acosh-normal" ensembles correspond to the cases of block-Toeplitz with Toeplitz blocks, and block-Toeplitz (with general blocks) covariance matrices, respectively.

Some results about the Equiangular Algorithm

Equiangular Subspaces in Euclidean Spaces

A Refined Shift-And-Invert Arnoldi Algorithm for Large Unsymmetric Generalized Eigenproblems

Solution of Linear Systems of Equations Ax=b and Ax=0 using Unifying Approach with Geometric Algebra: Outer Product Application and Angular Conditionality

Equiangular lines via improved eigenvalue multiplicity

The eigenvalue decomposition of normal matrices by the decomposition of the skew-symmetric part with applications to orthogonal matrices

The Jordan algebraic structure of the rotated quadratic cone

A Refined Subspace Iteration Algorithm For Large Sparse Eigenproblems

Equiangular lines via matrix projection

Computing Eigenvectors of Normal Matrices with Simple Inverse Iteration

The Refined Harmonic Arnoldi Method and an Implicitly Restarted Refined Algorithm for Computing Interior Eigenpairs of Large Matrices

Improving eigenvectors in Arnoldi's method

A Refined Iterative Algorithm Based on the Block Arnoldi Process for Large Unsymmetric Eigenproblems

A Variation on the Block Arnoldi Method for Large Unsymmetric Matrix Eigenproblems

A linear algebraic approach to orthogonal arrays and Latin squares

A simple, randomized algorithm for diagonalizing normal matrices

On matrix balancing and eigenvector computation

Research on Least‐Square η ‐Hermitian Solutions of Split Quaternion Matrix Equations

Riemannian statistics meets random matrix theory: towards learning from high-dimensional covariance matrices

Arnoldi Reduction Algorithm for Large Scale Gyroscopic Eigenvalue Problem

Exact Algorithms for Computing Generalized Eigenspaces of Matrices via Annihilating Polynomials