Abstract:We propose several methods that address the problem of fitting a C1\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$C^1$$\end{document} curve γ\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\gamma $$\end{document} to time-labeled data points on a manifold. The methods have a parameter, λ\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\lambda $$\end{document}, to adjust the relative importance of the two goals that the curve should meet: being “straight enough” while fitting the data “closely enough.” The methods are designed for ease of use: they only require to compute Riemannian exponentials and logarithms, they represent the curve by means of a number of tangent vectors that grows linearly with the number of data points, and, once the representation is computed, evaluating γ(t)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\gamma (t)$$\end{document} at any t requires a small number of exponentials and logarithms (independent of the number of data points). Among the proposed methods, the blended cubic spline technique combines the additional properties of interpolating the data when λ→∞\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\lambda \rightarrow \infty $$\end{document} and reducing to the well-known cubic smoothing spline when the manifold is Euclidean. The methods are illustrated on synthetic and real data.

A modified Casteljau algorithm to solve interpolation problems on Stiefel manifolds

SOLVING INTERPOLATION PROBLEMS ON STIEFEL MANIFOLDS USING QUASI-GEODESICS

Surface Interpolation Scheme by Distance Blending over Convex Sets

Comparison of two numerical methods for Riemannian cubic polynomials on Stiefel manifolds

A Modified de Casteljau Subdivision that Supports Smooth Stitching with Hierarchically Organized Bicubic Bezier Patches

A Riemannian optimization method on the indefinite Stiefel manifold

Subgradient Projection Algorithms for Convex Feasibility on Riemannian Manifolds with Lower Bounded Curvatures

The Convergence of the Geometric Interpolation Algorithm.

Gauss–Newton Method for Convex Composite Optimizations on Riemannian Manifolds

A Construction of Interpolating Space Curves with Any Degree of Geometric Continuity

De Casteljau's Algorithm in Geometric Data Analysis: Theory and Application

Convergence of Inexact Steepest Descent Algorithm for Multiobjective Optimizations on Riemannian Manifolds Without Curvature Constraints

Simple Piecewise Geodesic Interpolation of Simple and Jordan Curves with Applications

Fitting Smoothing Splines to Time-Indexed, Noisy Points on Nonlinear Manifolds

Piecewise-Bézier C1 Interpolation on Riemannian Manifolds with Application to 2D Shape Morphing

An objective FE-formulation for Cosserat rods based on the spherical Bézier interpolation

Randomized Submanifold Subgradient Method for Optimization over Stiefel Manifolds

A gradient flow method for smooth splines versus least-squares fitting on Riemannian manifolds

Data Fitting on Manifolds with Composite Bézier-Like Curves and Blended Cubic Splines

The interpolation problem: When can you pass a curve of a given type through N random points in space?

A Simple A-Priory Estimate for 3D Stationary Navier-Stokes System Via Interpolation