Conditioning and interpolation error bounds for second-order Stiefel retractions with closed-form inverses
Rasmus Jensen, Ralf Zimmermann
Abstract
Retractions provide a computationally efficient alternative to the Riemannian exponential and logarithm maps for practical data-processing tasks on manifolds. In particular, second-order retractions with closed-form inverse are well-suited for interpolation problems on manifolds. On the Stiefel manifold of orthogonal frames, there are only two retractions of this type: the Cayley retraction, which is second-order accurate under the canonical metric, and the recently proposed polar-light retraction, which is second-order accurate under the Euclidean metric. In this paper, we study the properties of these maps in the context of interpolation on the Stiefel manifold. To obtain explicit interpolation error bounds, we examine the conditioning of the retraction maps and their inverses. We show that the retractions are well-conditioned, and we derive interpolation error bounds similar to those of classical Euclidean interpolation. The inverse retractions are not well-conditioned in general, and we discuss how data can be mapped via an isometric group action to ensure stable computations. As with all retractions on compact manifolds, the inverse canonical Cayley retraction and the inverse polar-light retraction exist only locally, and we construct normal neighborhoods around any point in which either the inverse Cayley retraction or the invese polar-light retraction are guaranteed to be computable. As an application of the retraction maps, we consider Hermite interpolation, where the objective is to reproduce both sampled function values and derivative information. A numerical example demonstrates that retraction-based interpolation is competitive with classical methods based on Riemannian normal coordinates.
Create a lesson
Related papers
Graph lattice sums and graph zeta functions for long-range interacting quantum lattice models
Andreas Alexander Buchheit, Andreas Rupp
A numerical benchmark for fluid--structure--contact interaction
Daniele Corti, Jakub Fara, Miguel Angel Fernández et al.
Largest-dihedral-angle bisection algorithm does not preserve mesh regularity for tetrahedral partitions
Sergey Korotov, Jérôme Michaud
A Highly Scalable Quantized Tensor-Train FDTD Framework for the Simulation of Three-Dimensional Electromagnetic Scattering Problems
Daan Vanhaecke, Emile Vanderstraeten, Dries Vande Ginste
Pressure-robustness by commuting interpolation operators for Stokes discretizations with continuous pressures
Philip L. Lederer, Theresa Vock
A Reynolds-Semi-Robust, Globally Divergence-Free HDG Method for the Smagorinsky Model
Shuaijun Liu, Xiaoping Xie