The Reverse Telescoping Coordinate System for Positive Definite Matrices: Geometry, Computation, and Generative Modeling

We design a new unconstrained coordinate system where a p× p symmetric positive definite (SPD) matrix Theta is represented by a reverse telescoping map Theta(x)=rmRT(x), with x=(v,d,r)∈ℝ×ℝ(p-1)×ℝp(p-1)/2, representing respectively the log volume or log determinant; and the shape, as encoded by log relative diagonal scales and partial covariances among the nodes. This construction results in important properties not available in other charts, e.g., matrix logarithm, such as Jacobian depending on only the log-determinant. A useful feature of our construction is x contains a lossless symbolic representation of both the matrix and its inverse. Many important computations involving a matrix and its inverse can be performed in O(p2) in the transformed domain, while it is the rendering of results in matrix forms (on demand) that must incur an O(p3) cost. Moreover, two unit-determinant matrices in the transformed domain can be joined by a straight line with pathwise unit determinant. For generative modeling, this allows designing a split volume-shape flow model trained by conditional flow matching for transporting the shape over the unit-determinant path, with a separate one-dimensional flow for transporting the volume or the determinant. The forbidding SPD constraint, tamed thus into a powerful guiding force, leads to the surprising insight that it is in some sense easier to design a volume-normalized shape flow for SPD compared to the unconstrained ℝp× p, with no intrinsic notion of volume to aid normalization, unlike the determinant of SPD matrices. We apply our construction for up to p=200 in generative modeling of SPD matrices on a difficult synthetic bimodal target, and in generating brain connectivity networks by models trained on fMRI data; as well as in intrinsic diffusion on the SPD manifold.