$ \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}} $

Adjoint Jacobi field

Evaluates the adjoint of the Jacobi field along the geodesic at with initial conditions .

For the geodesic from to a Jacobi field is a smooth vector field along that can be evaluated with the function JacobiField(x,y,t,eta), where denotes one of the two initial conditions, i.e. .

We interpret in the following the Jacobi field for fixed as a function of and obtain . We introduce the short notation The adjoint is then meant in the sense that for any it holds

and hence . The adjoint can be computed employing the same weights as for the Jacobi field. If the weights of the Jacobi field represent a differential, this yields the adjoint differential.

Matlab Documentation

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
% AdjJacobiField(x,y,t,eta) - evaluate a adjoint of a Jacobi field
%    along the geodesic geo(x,y) at point t, where the
%    'weight'-function f(k,t,d) determines the boundary
%    conditions of the field and hence the meaning of eta,
%
% INPUT
%    x   : a point on the manifold Sn (or a vector of points)
%    y   : a point on the manifold Sn (or a vector of points)
%    t   : a value from [0,1] indicating the point geo(x,y,t)
%           (or a vector of values)
%    eta : an initial condition of the Jacobi field, where the
%           following weights determine the type of initial
%           condition.
%
% OPTIONAL
%    'weights' : [@(k,t,d) = (k==0)*t
%       + (k>0)*sin(sqrt(k)*t*d)/sin(sqrt(k)*d)
%       + (k<0)*sinh(sqrt(-k)*d*t)/sinh(sqrt(-k)*d)
%       provides the weight depending on the eigenvalue (k) of
%       the curvature tensor coresponding to the ONB basis
%       vector, the position t along the Jacobi field and d the
%       length of the geodesic.
%       For the standard value, eta is the Jacobi field at 0
%       and the second boundary condition is J(1)=0, i.e. the
%       Jacobifield corresponds to D_x\gamma_{xy}(t)[\eta]
%
% ---
% R. Bergmann | MVIRT | 2017-12-01

See also

used in