Computing the Riemannian logarithm on the Stiefel manifold: metrics, methods and performance

03/22/2021
by   Ralf Zimmermann, et al.
0

We address the problem of computing Riemannian normal coordinates on the real, compact Stiefel manifold of orthogonal frames. The Riemannian normal coordinates are based on the so-called Riemannian exponential and the Riemannian logarithm maps and enable to transfer almost any computational procedure to the realm of the Stiefel manifold. To compute the Riemannian logarithm is to solve the (local) geodesic endpoint problem. Instead of restricting the consideration to geodesics with respect to a single selected metric, we consider a family of Riemannian metrics introduced by Hüper, Markina and Silva-Leite that includes the Euclidean and the canonical metric as prominent examples. As main contributions, we provide (1) a unified, structured, reduced formula for the Stiefel geodesics for the complete family of metrics, (2) a unified method to tackle the geodesic endpoint problem, (3) an improvement of the existing Riemannian log map under the canonical metric. The findings are illustrated by means of numerical examples, where the novel algorithms prove to be the most efficient methods known to this date.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
10/19/2012

Learning Riemannian Metrics

We propose a solution to the problem of estimating a Riemannian metric a...
research
05/06/2021

A Unifying and Canonical Description of Measure-Preserving Diffusions

A complete recipe of measure-preserving diffusions in Euclidean space wa...
research
09/21/2020

Operator-valued formulas for Riemannian Gradient and Hessian and families of tractable metrics in optimization and machine learning

We provide an explicit formula for the Levi-Civita connection and Rieman...
research
08/16/2019

Hermite Interpolation and data processing errors on Riemannian matrix manifolds

The main contribution of this paper is twofold: On the one hand, a gener...
research
09/20/2021

Integrated Construction of Multimodal Atlases with Structural Connectomes in the Space of Riemannian Metrics

The structural network of the brain, or structural connectome, can be re...
research
12/14/2022

Toroidal Coordinates: Decorrelating Circular Coordinates With Lattice Reduction

The circular coordinates algorithm of de Silva, Morozov, and Vejdemo-Joh...
research
02/20/2023

Simplifying Momentum-based Riemannian Submanifold Optimization

Riemannian submanifold optimization with momentum is computationally cha...

Please sign up or login with your details

Forgot password? Click here to reset