$ \DeclareMathOperator{\arccosh}{arccosh} \DeclareMathOperator*{\argmin}{arg\,min} \DeclareMathOperator{\Exp}{Exp} \newcommand{\geo}[2]{\gamma_{\overset{\frown}{#1,#2}}} \newcommand{\geoS}{\gamma} \newcommand{\geoD}[2]{\gamma_} \newcommand{\geoL}[2]{\gamma(#2; #1)} \newcommand{\gradM}{\nabla_{\M}} \newcommand{\gradMComp}[1]{\nabla_{\M,#1}} \newcommand{\Grid}{\mathcal G} \DeclareMathOperator{\Log}{Log} \newcommand{\M}{\mathcal M} \newcommand{\N}{\mathcal N} \newcommand{\mat}[1]{\mathbf{#1}} \DeclareMathOperator{\prox}{prox} \newcommand{\PT}[3]{\mathrm{PT}_{#1\to#2}#3} \newcommand{\R}{\mathbb R} \newcommand{\SPD}[1]{\mathcal{P}(#1)} \DeclareMathOperator{\Tr}{Tr} \newcommand{\tT}{\mathrm{T}} \newcommand{\vect}[1]{\mathbf{#1}} $

The orthonormal basis in a tangent space

This function computes an ONB and correspoinding curvature coefficients, belonging to the parallel transportet orthonormal frame and diagonalizes the curvature tensor along . The set of vectors in a tangent space added as last dimension of the array and consists of vectors.

The following derivation of the orthonormal frame follows [1]. Let . We denote its eigenvalues by and an orthonormal basis of eigenvectors by . Then we can write

For ease of notation we introduce the indexing system

Then the matrices

form an orthonormal basis of $T_x {\mathcal P}(r)$.

The corresponding coefficients characterizing curvature (see Jacobi fields ) are

Matlab Documentation

1
2
3
4
5
6
7
8
9
10
11
12
% [Xi,k] = TpMONB(x,y) Compute an ONB in TpM and curvature
% coefficients corresponding to the transported frame along g(.,x,y)
%
% INPUT
%     p : base point( sets)
% OPTIONAL:
%     q : directional indicator( sets) for the first vector(s).
% OUTPUT
%    W : orthonormal bases ( n x n x SetDim x Dimension )
%    k : (optional) curvature coefficients (Dimension x SetDim)
% ---
% MVIRT R. Bergmann, 2017-12-03

See also

References

  1. Bačák, M, Bergmann, R, Steidl, G and Weinmann, A (2016). A second order non-smooth variational model for restoring manifold-valued images. SIAM Journal on Scientific Computing. 38 A567–A597