Sobolev Metrics on Spaces of Discrete Regular Curves
Jonathan Cerqueira Emmanuel Hartman Eric Klassen Martin Bauer

A Comprehensive Study of Open Cluster Chemical Homogeneity 
 using APOGEE and Milky Way Mapper Abundances
    [
    
=========================================================================================================

§ ABSTRACT 
Reparametrization invariant Sobolev metrics on spaces of regular curves have been shown to be of importance  in the field of mathematical shape analysis. For practical applications, one usually discretizes the space of smooth curves and considers the induced Riemannian metric  on a finite dimensional approximation space. Surprisingly, the theoretical properties of the corresponding finite dimensional Riemannian manifolds have not yet been studied in detail, which is the content of the present article. Our main theorem concerns metric and geodesic completeness and mirrors the results of the infinite dimensional setting as obtained by Bruveris, Michor and Mumford. 











§ INTRODUCTION



  
Motivation and Background:
Reparametrization invariant Sobolev metrics on spaces of regular curves play a central role in the field of mathematical shape analysis. Due to their reparametrization invariance, these metrics descend to Riemannian metrics on  spaces of unparametrized curves, which are of relevance in mathematical shape analysis and data analysis. Examples include applications where one is interested  in the shape of planar objects (represented by their boundary curves); see <cit.> and the references therein. Motivated by their appearance in these applications there has been a large interest in studying their mathematical properties. Michor and Mumford <cit.> showed a surprising degeneracy of the simplest such metric; namely they proved that the reparametrization invariant L^2-metric, i.e., the Sobolev metric of order zero, induces a degenerate distance function. This purely infinite dimensional phenomenon renders this metric unsuited for mathematical shape analysis, as it assigns a zero distance between any two curves and thus cannot distinguish between different shapes. For higher order metrics this degeneracy disappears and it has been shown in many applications that they lead  to meaningful notions of distance <cit.>. As a result, these metrics can be  used to define a mathematical framework for statistical shape analysis on these spaces of curves <cit.>. A natural question that arises in this context concerns the existence of minimizing geodesics, i.e., whether the space of  regular curves equipped with these Riemannian metrics is a geodesically complete and/or geodesically convex space. For metrics of order two or higher a positive answer to this question has been found by Bruveris, Michor and Mumford <cit.>; more recently, it has been shown by one of the authors and collaborators <cit.> that 3/2 is actually the critical index for this property, i.e., for a Sobolev metric of order greater than 3/2 the resulting space is geodesically complete, whereas there always exists geodesics that leave the space in finite time if the order is smaller than 3/2. The behavior at the critical value 3/2 is still open.



  
Main Contributions:
For practical applications, one usually discretizes the space of smooth curves and considers the induced Riemannian metric  in a corresponding finite dimensional approximation space. Discretizations that  have been considered include approximating curves via piecewise linear functions <cit.>,
B-spline discretizations <cit.> or finite Fourier series approximations <cit.>.  In this paper, we will discretize a curve as a finite sequence of points in Euclidean space, so that the space of curves is the space of these sequences. 
Using methods of discrete differential geometry, we will define a class of metrics on this finite dimensional space that are motivated by and analogous to the class of reparametrization invariant Sobolev metrics mentioned above for the infinite dimensional space of smooth curves.  To our surprise, these rather natural finite dimensional Riemannian manifolds have not yet been studied in much detail, with the only exception being the case of the homogenous Sobolev metric of order one, where the space of PL curves can be viewed as a totally geodesic submanifold of the infinite dimensional setting <cit.>;  note that a similar result for more general Sobolev metrics  is not true. 

Considering these finite dimensional approximation manifolds leads to a natural question, which is the starting point of the present article: 


Which properties of the infinite dimensional geometry are mirrored in these finite dimensional geometries?


As vanishing geodesic distance is a purely infinite dimensional phenomenon — every finite dimensional Riemannian manifold admits a non-degenerate geodesic distance function <cit.> —  one cannot hope to observe the analogue of this result in our setting. The main result of the present article shows, however, that some of these discretizations do indeed capture the aforementioned completeness properties, cf. Theorem <ref>. In addition to these theoretical results, we present in Section <ref> selected numerical examples showcasing the effects of the order of the metric on the resulting geodesics. Finally, for the special case of triangles in the plane, we study the Riemannian curvature of the space of triangles and observe that it explodes near the singularities of the space, i.e., where two points of the triangle come together. 



  
Conclusions and future work: In this article we studied discrete Sobolev type metrics on the space of discrete regular curves in Euclidean space (where we identified this space as a space of sequences of points) and showed that these geometries mirror several properties of their infinite dimensional counterparts. In future work we envision several distinct research directions: first, we have restricted ourselves in the present study to integer order Sobolev metrics. In future work it would be interesting to perform a similar analysis also for the class of fractional (real) order Sobolev metrics, such as those studied in <cit.>.
Secondly, we aim to study stochastic completeness of these geometries.  For extrinsic metrics on the  two-landmark space it has  recently been shown by Habermann, Harms and Sommer <cit.> that the resulting space is stochastically complete, assuming certain conditions on the kernel function. We believe that a similar approach could be  applied successfully to the geometries of the present article, which would be of interest in several applications where stochastic processes on shape spaces play a central role. Finally, we would like to study similar questions in the context of reparametrization invariant metrics on the space of surfaces: in this case geodesic completeness in the smooth category, i.e., in the infinite dimensional setting, is wide open and we hope to get new insights for this extremely difficult open problem by studying its finite dimensional counterpart. 



  
Acknowledgements: M.B was partially
supported by NSF grants DMS–2324962 and DMS-1953244. M.B and J.C. were partially supported by the BSF under grant 2022076. E.H. was partially supported by NSF grant DMS-2402555.


§ REPARAMETRIZATION INVARIANT SOBOLEV METRICS ON THE SPACE OF SMOOTH CURVES

In this section we will recall some basic definitions and results regarding the class of reparametrization invariant Sobolev metrics on the space of smooth, regular (immersed) curves. 
We begin by defining 
the set of smooth immersions of the circle S^1 into the space ^d:

    (S^1,^d)={c∈ C^∞(S^1, ^d) : |c'(θ)| 0, ∀θ∈ S^1}.

Here we identify S^1 with the interval [0,1] with its ends identified. The set (S^1,^d) is an open subset of the Fréchet space C^∞(S^1,^d) and thus can be considered as an infinite dimensional Fréchet manifold using a single chart.
We let h,k denote our tangent vectors, which belong to

    T_cImm(S^1,^d)≅ C^∞(S^1,^d).

Next, we consider the space of orientation-preserving smooth self-diffeomorphisms of the circle 
    Diff(S^1)={φ∈ C^∞(S^1,S^1):φ is bijective and φ'(θ)>0, ∀θ∈ S^1},

which is an infinite dimensional Fréchet Lie group and acts on 
(S^1,^d) from the right 
via the map (c,φ)↦ c∘φ.

The principal goal of the present article is to study Riemannian geometries on (S^1,^d). To introduce our class of Riemannian metrics we will first need to introduce some additional notation: we  denote by D_θ the derivative with respect to θ and let D_s=1/|c'(θ)|D_θ and ds=|c'(θ)|dθ be arc-length differentiation and integration respectively. Furthermore, we let ℓ(c):=∫_S^1|c'(θ)|dθ denote the  length of a curve c. Using these notations, we are ready to define the m-th order, reparametrization invariant Sobolev metric on (S^1,^d):

  Let c∈(S^1,^d), m∈ℤ_≥ 0, and h,k∈ T_c(S^1,^d). The m-th order Sobolev metric on (S^1,^d) is then given by
  
    ^m_c(h,k):=∫_S^1⟨ h,k⟩/ℓ(c)^3+ℓ(c)^-3+2m⟨ D_s^m h,D_s^m k⟩ ds.



In the above definition we used length dependent weights which allowed us to define a scale-invariant version of the Sobolev metric of order m. Alternatively we could have considered the constant coefficient Sobolev metric
  
    ^m_c(h,k):=∫_S^1⟨ h,k⟩+⟨ D_s^m h,D_s^m k⟩ ds.

and its discrete counterpart.
Almost all of the results of the present article hold also for this class of metrics, albeit with minor adaptions in the proof of the main completeness result, cf. Appendix <ref>.


We start by collecting several useful properties of the above defined family of Riemannian metrics:

  Let m≥ 0 and let ^m be the Riemannian metric as defined in (<ref>). Let c∈(S^1,^d) and h,k∈ T_c(S^1,^d). Then we have:
    
  
  * The metric ^m is invariant with respect to reparametrizations, rescalings, rotations, and translations; i.e. for φ∈Diff(S^1), λ∈^+, R∈SO(^d) and 𝐯∈^d we have 
    
    _c^m(h,k)   =_c∘φ^m(h∘φ,k∘φ)=_λ c^m(λ h,λ k)=_c+𝐯^m(h,k)=_Rc^m(Rh,Rk).

    We note in the above the action of translation (+𝐯) only affects the tangent space in which h,k lie, but not the vectors h and k, since the derivative of translation is the identity map.
  
  * The metric ^m is equivalent to any metric of the form
    
    ^m(h,k):=∫_S^1∑_j=0^m a_jℓ(c)^-3+2j⟨ D_s^j h,D_s^j k⟩ ds
 for a_0,a_m>0 and a_j≥ 0 for j=1,…,m-1; i.e., there exists a C>0 such that for all c∈(S^1,^d) and all h∈ T_c(S^1,^d)
    
    1/C^m_c(h,h)≤^m_c(h,h)≤ C ^m_c(h,h).

  



The proof of the invariances follows by direct computation. The second statement is implied by <cit.> in a similar way as in Lemma <ref> below.


For a Riemannian metric g on a smooth manifold ℳ one defines the induced path length, i.e., for γ(t):[0,1]→ℳ we let 

    L_g(γ)=∫_0^1√(g_γ(t)(D_tγ(t),D_tγ(t)))dt.

Using this one can consider the induced geodesic distance, which is defined via

    d_g(p,q)=inf L_g(γ)

where the infimum is calculated over the set of piecewise C^∞ paths [0,1]→ℳ with γ(0)=p and γ(1)=q.
In finite dimensions this always defines a true distance function; for infinite dimensional manifolds, however, one may have a degenerate distance function, i.e., there may be distinct points c_0,c_1∈ℳ such that d_g(c_0,c_1)=0, see for example <cit.>. For the class of Sobolev metrics on spaces of curves this phenomenon has been studied by Michor and Mumford and collaborators and a full characterization of the degeneracy has been obtained:

 The metric  defines a non-degenerate geodesic distance function on the space of curves (S^1,ℝ^d) if and only if m≥ 1.





The main focus on this article concerns geodesic and metric completeness properties of the corresponding space. Recall that 
a Riemannian manifold (ℳ,g) is metrically complete if it is complete as a metric space under the distance function given by the Riemannian metric. Furthermore it is called geodesically complete if the geodesic equation has solutions defined for all time for any initial conditions γ(0)=p∈ℳ and γ'(0)=v∈ T_pℳ and it is called  geodesically convex if for any two points p,q∈ℳ there exists a length minimizing path connecting them. In finite dimensions  the theorem of Hopf-Rinow implies that metric and geodesic completeness are equivalent and either implies geodesic convexity, but this result famously does not hold in infinite dimensions <cit.>. 

The following theorem will characterize the metric and geodesic completeness properties of ((S^1,^d),^m).
To state the theorem, we will first need to introduce the space of regular curves of finite (Sobolev) regularity, i.e., for m≥ 2 we let

    ℐ^m(S^1,^d)={c∈ H^m(S^1,^d): |c'(θ)|0,∀θ}.
 

  Let m≥ 0, c∈(S^1,^d), h,k∈ T_c(S^1,^d) and G^m as defined above. Then 
  ((S^1,^d),^m) is a Riemannian manifold and for m≥ 2 the following hold:
  
  
  * The manifold ((S^1,^d),G^m) is geodesically complete.
  
  * The metric completion of ((S^1,^d),G^m) is (ℐ^m,G^m).
  
  * The metric completion of ((S^1,^d),G^m) is geodesically convex.
  


  Proof of the above three statements can be found in <cit.>, specifically via the proofs of Theorems 5.1, 5.2, and 5.3.


Note that ((S^1,^d),G^m) is not metrically complete for any m despite being geodesically complete for m≥2.

In the exposition of the present article, we have only focused on integer order Sobolev metrics. More recently, the equivalent of these results have been also shown for fractional (real) order Sobolev metrics, where the critical order for positivity of the geodesic distance is 1/2 and the critical order for completeness is 3/2, see <cit.> for more details. 




§ DISCRETE SOBOLEV METRICS ON DISCRETE REGULAR CURVES


In this section we will introduce the main concept of the present article: a discrete version of the reparametrization invariant Sobolev metric.
We start this section by first defining a natural discretization of the space of immersed curves and then introducing a discrete analogue of ^m that, under modest assumptions, converges to its smooth counterpart. Let  denote the set of ordered n-tuples of points in ^d and consider the subset

    ={(𝐱_1,…, 𝐱_n)∈: 𝐱_i∈^d  and 𝐱_i≠𝐱_i+1, ∀  i∈ℤ/nℤ}.

As  is clearly an open subset of  it carries the structure of an dn-dimensional manifold.
Furthermore, 
 and  are both acted on from the right by the cyclic group of n elements by cyclically permuting the indices. That is, if j∈ℤ/nℤ then

    ((𝐱_1,…, 𝐱_n),j)↦ (𝐱_1+j,…, 𝐱_n+j)

where addition is modulo n.




In this remark we will observe that we can identify the above introduced space of discrete regular curves  with the set of piecewise linear, regular curves. Therefore we  let

    PL^n(S^1,^d)={c∈ C(S^1,^d): c|_[i/n,i+1/n] is linear}

be the space of  continuous, piece-wise linear curves with n control points from S^1→^d
and 

    PLImm^n(S^1,^d)={c∈PL^n(S^1,^d): c(i/n) c(i+1/n), ∀  i∈ℤ/nℤ}

the open subset of piecewise linear immersions.
To identify the space of piecewise linear immersion with the previously defined space
, thereby allowing us to visualize  as a space of curves, we simply consider the map 
    c↦(c(0/n),c(1/n),…,c(n-1/n))∈

  Likewise we may identify T_cPLImm^n(S^1,^d)≅PL^n(S^1,^d) with T_c≅. 

We now wish to define a Riemannian metric ^m on , which can be interpreted as a discretization of ^m. First, we define some notation.
Given a curve c∈ we denote our vertices as c_i and let e_i=c_i+1-c_i denote the edge beginning at c_i and ending at c_i+1. We also denote the average length of the two edges meeting at the vertex c_i as μ_i=|e_i|+|e_i-1|/2. Let h∈ T_c denote a tangent vector to c with h_i denoting the i-th component of h. To define our discrete metric ^m we will also need to define the j-th discrete derivative at the ith vertex D_s^j h_i 
using some ideas from discrete differential geometry <cit.>. We define these recursively via

    D_s^0 h_i=h_i     D_s^j+1 h_i=D_s^jh_i-D_s^jh_i-1/μ_i   j∉2ℕ
    D_s^jh_i+1-D_s^jh_i/|e_i|   j∈2ℕ.
 
A visualization of this construction can be seen in Figure <ref>. 


This allows us to define some useful notation for something analogous to the m-th order component of ^m.


  For c∈ and h,k∈ we let 
  
    _c^m(h,k):=∑_i=1^nℓ(c)^2m-3⟨ D_s^m h_i,D_s^m k_i⟩μ_i,m where μ_i,m=μ_i   m∈ 2ℤ^+
     |e_i|   m∉ 2ℤ^+

  and 
    
    _c^m(h,k)=_c^0(h,k)+_c^m(h,k).


In the following lemma we justify the above notation and show the above defines a Riemannian metric on .

  Let n≥ 3 and m≥ 0. Then
  ^m
   is a Riemannian metric on .


  We need to check non-degeneracy of the inner product and the smooth dependence on the base point c.
  As g_c^m(h,h)≥ġ_c^0(h,h) it follows that if g_c^m(h,h)=0 then we also have ġ_c^0(h,h)=0. As each edge length |e_i| and μ_i are nonzero for any c, this can only happen if ⟨ D_s^0 h_i, D_s^0 h_i⟩=⟨ h_i, h_i⟩=0 for all i. That implies h_i=0 for all i. Thus for any m, g_c^m(h,h)=0 implies h≡ 0.
  
  We next check that ^m_c varies smoothly with respect to the base point c. The only suspect term in the definition is |e_i|. Varying c_i in the u direction yields the following:
  
    ∂/∂ u|e_i|=-u· e_i/|e_i|    ∂/∂ u|e_i-1|=u· e_i-1/|e_i-1|

  which implies the smoothness of the metric as |e_i|>0 for each i of any curve in .

In the next Lemma we will show that higher order metrics dominate lower order metrics, which will be of importance in the next section, where we will establish the main completeness results of the present article.

  Let m≥ 1,  c∈ and h∈, then
  
    ^m_c(h,h)≤1/4^m+1_c(h,h).



The proof of this result is rather technical and we postpone it to the appendix. 

With these lemmas in hand, we arrive at a result concerning the main properties of the discrete Riemannian metric g^m, which parallels the statements of Lemma <ref>. 

  The metric ^m_c(h,k) has the following properties:
  
  
  *  For all n≥ 3, m≥0, ^m is invariant with respect to the action of ℤ/nℤ on , as well as rescaling, rotation, and translation. If j∈ℤ/nℤ, λ∈^+, R∈SO(^d) and 𝐯∈^d, then
    
    _c^m(h,k)   =_c∘ j^m(h∘ j,k∘ j)=_λ c^m(λ h,λ k)=_𝐯+c^m(h,k)=_Rc^m(Rh,Rk).

  
  *  The metric ^m(h,k) is equivalent to any metric of the form
    
    ^m(h,k):=∑_j=0^m ∑_i=1^na_jℓ(c)^2j-3⟨ D_s^j h_i,D_s^j k_i⟩μ_i,m
 for a_0,a_m>0 and a_j≥ 0 for j=1,…,m-1 in the following sense: for some C∈, all c∈(S^1,^d) and all h∈ T_c≅
    
    1/C^m_c(h,h)≤^m_c(h,h)≤ C ^m_c(h,h).

  


  The proof of statement <ref> is similar to the smooth case and we will not present the proof for translation, rotation, or scale invariance. The key difference is the action by ℤ/nℤ. Note _c^m(h,k) and _c∘ j^m(h∘ j,k∘ j) sum over the exact same terms. Therefore, since |e_i∘ j|=|e_i+j| and μ_i∘ j=μ_i+j, we see
  
    _c∘ j^m(h∘ j,k∘ j)=∑_i=1+j^n+jℓ(c)^2m-3⟨ D_s^m h_i,D_s^mk_i⟩μ_i,m=∑_i=1^nℓ(c)^2m-3⟨ D_s^m h_i,D_s^mk_i⟩μ_i,m=_c^m(h,k)

  where the center equality is due to the fact that i belongs to ℤ/nℤ so that the i=1 term of the sum is identical to the i=n+1 term.
  
  Statement <ref> largely follows via Lemma <ref>. Given our ^m, it allows us to collect all our coefficients into ^0 and ^m terms alone to form a ^m with
  
    a_0^0_c(h,h)+a_m^m_c(h,h)≤^m_c(h,h)≤^m_c(h,h).

  As ^m_c(h,h)=â_0^0(h,h)+â_m^m(h,h) for some a_0≤â_0 and a_m≤â_m. Note also that
  
    min(â_0,â_m)^m_c(h,h)≤^m_c(h,h)≤max(â_0,â_m)^m_c(h,h)

  and similar for a_0^0_c(h,h)+a_m^m_c(h,h). Thus we may write
  
    1/C^m_c(h,h)≤min(a_0,a_m)^m_c(h,h)≤^m_c(h,h)≤max(â_0,â_m)^m_c(h,h)≤ C^m_c(h,h)

  for
  
    C=max(1/min(a_0,a_m), max(â_0,â_m)).

  Via multiplication by C we find ^m_c(h,h)≤ C^m_c(h,h) and via division we find 1/C^m_c(h,h)≤^m_c(h,h).
  Together these give the desired statement

    1/C^m_c(h,h)≤^m_c(h,h)≤ C^m_c(h,h).
  

We have seen that this discrete metric shares several properties with its infinite dimensional counterpart. In the following proposition we show that we are indeed justified in referring to it as a discretization of ^m:
  Let n≥ 3 and for i=0,…, n-1 let θ_i=i/n. Let c∈(S^1,^d) and h,k∈ T_c(S^1,^d). Define c̃_n=(c(θ_1),c(θ_2),…,c(θ_n))∈. Similarly define h̃_n and k̃_n using h,k. Then
  
    lim_n→∞^m_c̃_n(h̃_n,k̃_n)=^m_c(h,k)



  To prove this statement we will make use of a technical result  regarding the approximation of derivatives, which we postponed to the appendix, cf. Lemma <ref>.
  We first define several relevant operators and notational items.
  Given c∈(S^1,^d) we let μ_c,n,m=|I_n^1(c)| if m is even and μ_c,n,m=|I_n^1(c)|+|I_n^1(c)|/2 if m is odd. Now, for a given c∈(S^1,^d) and n∈ℤ_≥ 3 define 𝒟_c,n,m via the following recursive formula
  
    𝒟_c,n,0(f)=I^0_n(f)    𝒟_c,n,m(f)=I^1_n(D_c,n,m-1(f))/μ_c,n,m   m∈ 2ℤ^+
    I^1_n(D_c,n,m-1(f))/μ_c,n,m   m∉ 2ℤ^+.

  Note that all the n cancel in these definitions once expanded as at each step we introduce a factor of n in both the numerator and denominator.
  Next define ℓ_n(c) as ∑_i=1^n|c(θ_i+1)-c(θ_i)|. With these in hand we define:
  
    F_m(h,k,c,n,θ)=ℓ_n(c)^2m-3⟨𝒟_c,n,m(h),𝒟_c,n,m(k)⟩μ_c,n,m.

  We prove the proposition in two steps
  
    ∫_S^1F_m(h,k,c,n,θ)dθ=^m_c_n(h_n,k_n)
  
  We expand to
  
    ∫_S^1F_m(h,k,c,n,θ)dθ=∫_S^1ℓ_n(c)^2m-3⟨𝒟_c,n,m(h),𝒟_c,n,m(k)⟩μ_c,n,mdθ

  and note that the integrand on the right is constant on each interval [θ_i,θ_i+1) as ℓ_n is constant with respect to θ and each of 𝒟_c,n,m(h),𝒟_c,n,m(k), and μ_c,n,m are constant on such intervals. Thus for some 𝒟_c,n,m,i(h),𝒟_c,n,m,i(k)∈^d we can write either  
    ∫_S^1F_m(h,k,c,n,θ)dθ=∑_i=1^nℓ_n(c)^2m-31/n⟨𝒟_c,n,m,i(h),𝒟_c,n,m,i(k)⟩ |n(c(θ_i+1)-c(θ_i))|

  or
  
    =∑_i=1^nℓ_n(c)^2m-31/n⟨𝒟_c,n,m,i(h),𝒟_c,n,m,i(k)⟩|n(c(θ_i+1)-c(θ_i))|+|n(c(θ_i)-c(θ_i-1))|/2

  depending on the parity of m. Note all n cancel and we can rewrite further to either
  
    ∑_i=1^nℓ_n(c)^2m-3⟨𝒟_c,n,m,i(h),𝒟_c,n,m,i(k)⟩|e_i|

  for odd m or
  
    ∑_i=1^nℓ_n(c)^2m-3⟨𝒟_c,n,m,i(h),𝒟_c,n,m,i(k)⟩μ_i

  for even m. Thus the only thing that remains to be shown for this first claim is that these 𝒟 terms are our discrete D_s derivatives of h_n and k_n. This is clear for m=0 as I^0_n(h) is simply the step function taking on the value h_n,i=h(i/n) on [θ_i,θ_i+1). For m>0 we see either
  
    𝒟_c,n,m,i(h)=𝒟_c,n,m-1,i+1(h)-𝒟_c,n,m-1,i(h)/|e_i|

  or
  
    𝒟_c,n,m,i(h)dθ=𝒟_c,n,m-1,i(h)-𝒟_c,n,m-1,i-1(h)/μ_i

  so that by induction it follows that these really do correspond to the discrete derivatives as 𝒟_c,n,m,i(h) relates to 𝒟_c,n,m-1,j(h) identically to how D_s^mh_i relates to D_s^m-1 h_j.

  The second step is easier by comparison. We have shown that ∫_S^1F_m(…)=_c̃_n^m(…). We next need to show how the limit of F_m is related to ^m_c(h,k)=ℓ^2m-3⟨ D_s^m h,D_s^m k⟩|c'| which will allow us to relate ^m and ^m.
  
    lim_n→∞ F_m(h,k,c,n,θ)L^∞→^m_c(h,k)   
  By Lemma <ref> it follows that
  
    lim_n→∞𝒟_c,n,m(f)-D_s^m(f)_L^∞→ 0  and  lim_n→∞μ̂_c,n,m-|c'|_L^∞→0.

  It is also clear that lim_n→∞ℓ_n(c)→ℓ(c) as our curves are rectifiable.
  Now 
    lim_n→∞F_m(h,k,c,n,θ)-^m_c(h,k)_L^∞→ 0.

  This implies that lim_n→∞F_m(…) is equal to ^m_c(h,k) almost everywhere, and we may join our claims together via
  
    ∫_S^1^m_c(h,k)dθ=∫_S^1lim_n→∞ F_m(…) dθ*=lim_n→∞∫_S^1F_m(…)dθ=lim_n→∞^m_c̃_n(h̃_n,k̃_n),

  where * follows by the Lebesgue dominated convergence theorem as we can bound the integrand almost everywhere by the constant function on S^1 equal to ^m_c(h,k)_L^∞. By linearity it follows
  
    ^m_c(h,k)=lim_n→∞^m_c̃_n(h̃_n,k̃_n)

  which completes the proof.




§ COMPLETENESS RESULTS

In this section we will prove the main theoretical result of the present article, namely we will show that the finite dimensional manifolds (,^m) indeed capture the completeness properties of the space of smooth, regular curves equipped with the class of reparametrization invariant Sobolev metrics ^m, i.e., we will prove the  following theorem:

  Let m≥ 2, n≥ 3 and d≥ 2. Then the space (,^m) is metrically and geodesically complete. Furthermore,  for any two points in  there exists a minimizing geodesic, i.e.,  (,^m)  is geodesically convex.

This theorem can readily be seen as a discrete equivalent to the above Theorem <ref> from the smooth case. As  is finite dimensional, Hopf-Rinow implies we need only show metric completeness and that both geodesic completeness and geodesic convexity follow automatically. The requirement that d≥ 2 is due to the fact that ^1× n_* is not connected and thus Hopf-Rinow does not apply.  We begin by proving a useful lemma

  Let (ℳ,g) be a possibly infinite dimensional manifold. (ℳ,g) is metrically incomplete if and only if there exists a path γ:[0,1)→ℳ, such that 
  
  
  * γ(t)∈ℳ for 0≤ t<1;
  
  * lim_t→ 1γ(t) does not exist in (ℳ,g);
  
  * the length of γ (w.r.t. g) is finite, i.e., L_g(γ)<∞.
  


  We first assume incompleteness. As ℳ is metrically incomplete there exists some Cauchy sequence x_n that fails to converge in the space. Without loss of generality we can assume ∑_n^∞ d_g(x_n,x_n+1)<∞. This is since we can always find a Cauchy sub-sequence of our Cauchy sequence with this property. To construct such a sub-sequence we perform the following procedure: if, under our Cauchy conditions, for all n,m≥ N(i) we have d_g(x_n,x_m)<1/10^i then we define our sub-sequence as {x_N(i)}_i∈ℕ⊂{x_n}. The length ∑_i=1^∞ d_g( x_N(i),x_N(i+1)) of this sub-sequence is bounded above by 1/9 by construction as ∑_i=1^∞ d_g( x_N(i),x_N(i+1))<∑_i=1^∞1/10^i=1/9.
  
  While we cannot guarantee paths lying in the space between the elements of our sequence that realize their distances, we can have a collection of paths γ_n going from x_n to x_n+1 with a length <d_g(x_n,x_n+1)+1/2^n lying within the space so that the overall distance of γ=γ_1*γ_2*γ_3*… will have a total length L(γ)< 1+∑_n=1^∞ d_g(x_n,x_n+1). Here * denotes usual path composition. γ is then of finite length with lim_t→ 1γ(t)∉ℳ and γ(t)∈ℳ for 0≤ t< 1.
  
  On the other hand, assuming such a path, we wish to show we have a Cauchy sequence that does not converge with respect to our metric. Without loss of generality assume γ is traversed with constant speed. 
  If a and b are points on our curve, then let d(a,b) be the length of the subsegment of gamma with endpoints γ(a),γ(b). As γ is traversed with constant speed we see d(t_i,t_j)= L_g(γ)|t_i-t_j|. The distance of γ(t_i) and γ(t_j) along γ is not necessarily equal with the distance of these points in the space, but if d_g is the metric distance, then d_g(γ(t_i),γ(t_j))≤ d(t_i,t_j)= L_g(γ)|t_i-t_j|. Set t_n=1-1/n noting that this sequence (t_n) is a Cauchy sequence in the reals. The distances between a pair of these points is d_g(γ(t_i),γ(t_j))≤ L_g(γ)|t_i-t_j|. Now given any ε>0 we have some N so that for all i,j>N we have |t_i-t_j|<ε/L_g(γ) so that d_g(γ(t_i),γ(t_j))<ε. This makes γ(t_n)  a Cauchy sequence in ℳ, but (t_n) converges to 1 so this Cauchy sequence does not converge in our space:
  
    lim_n→∞γ(t_n)=lim_t→ 1γ(t)∉ℳ.

  Thus our space is not metrically complete.


  Let (ℳ,g) and (ℳ,g^*) be two Riemannian manifolds such that for all p∈ℳ, r∈^+, q∈ B_g^*(p,r) and h∈ T_qℳ we have g_q(h,h)≤ Cg^*_q(h,h) for some constant C potentially dependent on our choice of p and r. Then if (ℳ,g) is metrically complete so is (ℳ,g^*).


  Suppose (ℳ,g^*) were incomplete while (ℳ,g) were complete. By Lemma <ref> there exists a path γ(t) such that γ(t)∈ℳ for t∈[0,1) and lim_t→ 1γ(t)∉ℳ with finite length under g^* and no similar finite length path exists under g. Let r_* be such that L_g^*(γ)=r_0<r_*<∞. Note g_γ(t)(h,h)≤ Cg^*_γ(t)(h,h)  since  γ(t)∈ B(γ_0,r_*) for all t∈[0,1). We thus have
  
    L_g(γ(t))   =∫_0^1√(g_γ(t)(D_t γ(t),D_t γ(t)))dt
       <√(C)∫_0^1√(g^*_γ(t)( D_t γ(t),D_t γ(t)))dt=√(C)L_g^*(γ(t))=r_0√(C)<∞


  This is a contradiction as it implies γ is of finite length in (ℳ,g).

With this in hand, we are now ready to prove Theorem <ref>.

  Corollary <ref> and Lemma <ref> together imply we need only consider the metric completeness of  equipped with ^2. Via Lemma <ref> it is sufficient to show any path :[0,1)→ with L_^2()<∞ has lim_t→ 1(t)∈ to show metric completeness of (,^2). 
  There are three cases to consider, as there are three ways for lim_t→ 1(t) to not exist in :
  
  
  * Some point escapes to infinity, edge lengths stay greater than 0, and curve length does not become infinite:
 ∃ j such that lim_t→ 1|_j|=∞, ℓ(γ(t))<∞ and |e_i|>0, ∀ i,∀ t∈[0,1].
  
  * Some edges but not all edges shrink to length 0:
 ∃ i: lim_t→ 1|e_i|=0 and ∃ j: |e_j|>0, ∀ t∈[0,1].
  
  * All edges shrink to length 0 or curve length is unbounded:
 lim_t→ 1ℓ((t))=0 or lim_t→ 1ℓ((t))=∞   
  
  
  Below we make use of the following estimates
  
    |D_t|(t)_i||≤ |D_t (t)_i|    |D_t|e_i||≤ |D_t e_i|    |D_tℓ((t))|≤∑_i=1^n|D_t e_i|.

  The first two of these are derived from the same general pattern:
  
    ± D_t |v|=±⟨ D_t v,v⟩/|v|=⟨ D_t v,± v⟩/|v|≤ |D_t v|

  as ± v/|v| is a unit vector. The third inequality comes from noting |D_tℓ(γ(t))|=|∑_i=1^nD_t|e_i|| and using the second inequality.


  Now we proceed by cases:

  (Case <ref>): We appeal to ^0. Let (t)_j be the vertex escaping
  
    Ł_^2((t))   ≥∫_0^1√(∑_i=1^n1/ℓ((t))^3|D_t(t)_i|^2μ_i)dt≥∫_0^11/ℓ((t))^3/2|D_t (t)_j|√(μ_j)dt≥ M∫_0^1|D_t (t)_j|dt
       ≥ M∫_0^1|D_t(|(t)_j|)|dt≥ M|∫_0^1D_t(|(t)_j|)dt|=M|∞lim_t→∞|(t)_j|-<∞|(0)_j||=∞

  where M=min_i,t√(μ_i/ℓ(γ(t))^3).

  (Case <ref>): We appeal to ^2. Let e_k be an edge shrinking to length 0 such that e_k-1 is not shrinking to length 0. Such a pair of consecutive edges must exist, as we have edges shrinking to length zero and others not doing so, these subsets of our edges must meet at some vertex γ(t)_k.
  
    Ł_^2(c)≥   ∫_0^1√(ℓ((t)))|D_t e_k/|e_k|-D_t e_k-1/|e_k-1||1/√(μ_i)dt≥ M|∫_0^1|D_t e_k/|e_k||-|D_t e_k-1/|e_k-1||dt|
    ≥    M|∫_0^1|D_t e_k/|e_k||dt-∫_0^1|D_t e_k-1/|e_k-1||dt|≥ M|∫_0^1|D_t|e_k|/|e_k||dt-∫_0^1|D_t e_k-1/|e_k-1||dt|
       ≥ M||∫_0^1D_t|e_k|/|e_k|dt|-∫_0^1|D_t e_k-1/|e_k-1||dt|≥ M|∞|[ln(|e_k|)]_0^1|-<∞∫_0^1|D_t e_k-1/|e_k-1||dt|=∞

  where M=min_i,t√(ℓ((t))/μ_i).

  (Case <ref>): We appeal to ^1.
  
    Ł_^2((t))   ≥ CŁ_^1((t))≥ C∫_0^1√(∑_i=1^n|D_t e_i|^2/ℓ((t))|e_i|)dt≥ C∫_0^1|D_te_i|/ℓ((t))dt   for each i
    
    Ł_^2((t))   *≥C/n∑_i=1^n∫_0^1|D_t e_i|/ℓ((t))dt=C/n∫_0^1∑_i=1^n|D_t e_i|/ℓ((t))dt≥C/n∫_0^1|D_tℓ((t))|/ℓ((t))dt
       ≥C/n|∫_0^1D_tℓ((t))/ℓ((t))dt|=C/n|[ln(ℓ((t)))]_0^1|=∞

  The ≥ marked * comes from the fact that the average of several underestimates is itself an underestimate.





§ GEODESICS, CURVATURE AND A COMPARISON TO KENDALL'S SHAPE SPACE

In this section we will further demonstrate the behavior of these geometries by  presenting selected examples of geodesics for different choices of metrics. Finally, we will consider the special case of planar triangles, where we will compare the discrete Sobolev metrics defined in this paper with Kendall's shape metric; it is well known that under Kendall's metric the space of planar triangles reduces to a round sphere <cit.>.

All numerical examples were obtained using the programming language , where we implemented both the Riemannian exponential map (and the Riemannian log map) by solving the the geodesic initial value problem (and boundary value problem, resp.).

[12]r0.32  
  

    
    < g r a p h i c s >

    

  Riemannian exponential map w.r.t. to the metric ^2.
  







  
The geodesic initial value problem:To approximate the Riemannian exponential map we calculate the Christoffel symbols using the automatic differentiation capabilities of . This in turn allows us to approximate the exponential map using a simple one-step Euler method. In Figure <ref> one can see an approximated geodesic computed with initial conditions consisting of a right triangle (0,0),(0,1),(1,0) and an initial velocity of (0,1) on the third vertex and zero elsewhere. We note that an initial velocity which is non-zero only on a single vertex immediately puts multiple vertices in motion. 



  
The geodesic boundary value problem:To approximate the Riemannian log map we employ a path-straigthening algorithm; i.e.,  we minimize the Riemannian energy over paths of curves starting at c_0 and ending at c_1, approximated using a fixed (finite) number of intermediate curves. Thereby, we reduce the approximation of the Riemannian log map to a finite dimensional, unconstrained minimization problem, which we tackle  using the  implementation of L-BFGS-B, where we use a simple linear interpolation as initialization. In particular, for higher order metrics the algorithm can fail if the initialization leaves the space of PL immersions, as the metric is undefined for curves that are not in . To fix this issue it is possible to add a small amount of  noise to the vertices of the initialization. We present two different examples of solutions to the geodesic boundary value problem in Figure <ref>.







  
Gaussian Curvature of the Space of Triangles and Kendall's shape metric:
A prudent comparison of the discrete metrics described in this work is with the Kendall metrics of discrete shapes, which may be defined on the same space. Comparing the geodesics with respect to these metrics also  highlights the nature of our completeness result. 

We start by restricting ourselves to the space of triangles ℝ^2× 3_* modulo the actions of rotation, translation, and scale. This space can be identified with the surface of a sphere with three punctures (the punctures correspond to degenerate triangles in which two vertices coincide) and when equipped with the Kendall metric <cit.>, this space is famously isometric to the punctured sphere and has constant Gaussian curvature. In Figure <ref>, we display boundary value geodesics with respect to the Kendall metric as well as for the discrete Sobolev metrics for m=0,1,2. While the geodesic with respect to the Kendall metric passes directly through a discrete curve where adjacent points coincide, each of the geodesics with respect to the metrics proposed in this paper does not pass through this point in the space. Moreover, geodesics with respect to the higher-order metrics pass further from this point.

Finally, we calculated the Gaussian curvature for the space of triangles w.r.t to the same four Riemannian metrics: while the Gaussian curvature for the Kendall metric is constant, this is clearly not the case for discrete Sobolev metrics of the present article, cf. Figure <ref>.
Indeed all of the ^m metrics exhibit negative curvature near the points corresponding to triangles with double-points (which do not correspond to elements of ), but only for higher order metrics is this negative curvature strong enough to prevent all geodesics from leaving the space in finite time. More generally, an increase in order leads for a more curved space both for negatively curved regions but also for positively curved regions. Finally, in Figure <ref>, we further visualize the  curvature along selected paths to further demonstrate the behavior at key points of interest. One can see again the increase in negative curvature near the punctures of the sphere. 











§ CONSTANT COEFFICIENT SOBOLEV METRICS

Our use of the ℓ(c) terms in ^m metrics allowed for scale invariance. These terms may be removed to arrive at the constant coefficient Sobolev metrics

    _c^m(h,k):=∫_S^1⟨ h,k⟩+⟨ D_s^m h,D_s^m k⟩ ds.

These lack scale invariance, but the manifolds ((S^1,^d),^m) are otherwise similar to the scale invariant formulation with identical completeness properties <cit.>. We can likewise drop the ℓ(c) terms from our discretized metrics to arrive at constant coefficient discrete metrics ^m. These also converge to the corresponding smooth ^m like the scale invariant formulation, using exactly the same proof as used for Proposition <ref>, but without the ℓ(c) terms in F. Similarly, the parallel of the completeness statement, Theorem <ref> still holds. The removal of the length terms  does not impact the proofs of cases <ref> or <ref> beyond changing the relevant  M, but the proof for case <ref> no longer works. We make use of the following lemma for an alternative proof, stated here in a much more narrow and simple context than in its original form.

  Let (ℳ,g) be a Riemannian manifold with a weakly continuous metric and f:ℳ→ a C^1-function. Assume that for each metric ball B(y,r) in ℳ there exists a constant C, such that
  
    |df/dh|≤ C(1+|f(x)|)h_g

  holds for all x∈ B(y,r) and all h∈ T_xℳ. Then the function
  
    f:(M,d)→ (,|·|)

  is continuous and Lipschitz continuous on every metric ball. In particular, f is bounded on every metric ball.


With this we can re-prove case <ref> for the constant coefficient case.

  Note that it is sufficient to show that 1/ℓ(c) and ℓ(c) are bounded on every metric ball to handle case <ref>. By the above lemma then, it is sufficient to show that
  
    |d(1/ℓ(c))/dh|≤ C(1+|1/ℓ(c)|)h_g^2

  We do this directly. 
  
    |d/dh1/ℓ(c)|=1/ℓ(c)^2d/dh∑_i=1^n |e_i|=1/ℓ(c)^2∑_i=1^n ⟨ h_i+1-h_i,e_i⟩/|e_i|≤1/ℓ(c)^2∑_i=1^n  h_i+1-h_i/√(|e_i|)e_i/√(|e_i|)

  
    ≤1/ℓ(c)^2(∑_i=1^n  h_i+1-h_i^2/|e_i|)^1/2(∑_i=1^n|e_i|)^1/2= 1/ℓ(c)^3/2(∑_i=1^n  h_i+1-h_i^2/|e_i|)^1/2≤1/ℓ(c)^3/2√(^2_c(h,h))

  So that at any particular c the constant C=√(ℓ(c)) yields the desired inequality. We can show via similar work that √(ℓ(c)) is globally Lipschitz:
  
    |d/dh√(ℓ(c))|≤1/√(ℓ(c))(∑_i=1^n  h_i+1-h_i^2/|e_i|)^1/2(∑_i=1^n|e_i|)^1/2=(∑_i=1^n  h_i+1-h_i^2/|e_i|)^1/2≤√(^2_c(h,h)).

  Thus we can indeed choose a C for any metric ball as needed for the lemma for 1/ℓ(c). Thus 1/ℓ(c) and ℓ(c) are bounded on all metric balls and therefore no path γ:[0,1)→ such that lim_t→ 1γ(t)=0 or ∞ can be of finite length.

The item most complicated by the removal of the length terms is the relation between ^m_c(h,h) and ^m+1_c(h,h). As ^m_c=^m_c·ℓ(c)^-3+2m the statement of Lemma <ref> here becomes

    ^m_c(h,h)≤ℓ(c)^2/4^m+1_c(h,h)

which is a parallel to the equivalent statement for the smooth constant coefficient case as seen in the proof of Lemma 2.13 of <cit.>.






§ APPROXIMATING DERIVATIVES IN L^∞


  Let f be a bounded function in C^∞(S^1). Let θ_i=i/n and define 
  
    I_n^0:C^∞(S^1)→ L^∞(S^1)

  so that I_n^0(f)=g_0 where g_0 is the piecewise constant function where g_0(θ)=f(θ_i) for θ∈[θ_i,θ_i+1). Further define I_n^1(f)=g_1 where g_1(θ)=n(f(θ_i+1)-f(θ_i)) for θ∈[θ_i,θ_i+1). Then
  
    lim_n→∞I_n^0(f)-f_L^∞=0    and    lim_n→∞I_n^1(f)-f'_L^∞=0.



  Note that S^1 is compact so that our smooth, bounded function f is Lipschitz continuous. Now if K is our Lipschitz constant then 0≤sup_θ∈[θ_i,θ_i+1)|g_0(θ)-f(θ)|≤K/n but this does not depend on i. This means that as n→∞ this supremum goes to 0 and thus also
  
    lim_n→∞g(θ)-f(θ)_L^∞=lim_n→∞max_i=1,…,nesssup_θ∈[θ_i,θ_i+1)(g_0(θ)-f(θ))= 0.

  For the second operator we choose n Lipschitz constants, one for each i with K_i=sup_θ∈[θ_i,θ_i+1)f'(θ) being the constant we will use for [θ_i,θ_i+1). This means sup f(θ_i+1)-f(θ_i)≤K_i/n. We note immediately as n→∞ K_i→ f'(θ_i) since f is smooth. Leaning on this fact makes the work swift:
  
    I_n^1(f)-f'_L^∞   =g(θ)-f'(θ)_L^∞=max_i=1,…,nn·(f(θ_i+1)-f(θ_i))-f'_L^∞
       ≤max_i=1,…,nn·(K_i/n)-f'_L^∞=max_i=1,…,nK_i-f'_L^∞→ 0 as n→∞


Define I_n^1(f) analogously to I_n^1(f), but with g_1(θ)=n· (f(θ_i)-f(θ_i-1)) for θ∈[θ_i,θ_i+1) and note that this is also such that lim_n→∞I_n^1(f)-f'_L^∞=0.




§ PROOF OF LEMMA <REF>

Before we are able to present the proof of Lemma <ref> we will need an additional technical Lemma:

  Let m≥ 1, let c∈ and h∈ T_c≅, we then have
  
    max_k∈ℤ/nℤ|D_s^m-1(h_k+1)-D_s^m-1(h_k)/|e_k||≤1/2∑_i=1^n|D_s^m-1(h_i+1)-D_s^m-1(h_i)/|e_i|-D_s^m-1(h_i)-D_s^m-1(h_i-1)/|e_i-1||

 for odd m and
  
    max_k∈ℤ/nℤ|D_s^m-1(h_k)-D_s^m-1(h_k-1)/μ_k|≤1/2∑_i=1^n|D_s^m-1(h_i+1)-D_s^m-1(h_i)/μ_i+1-D_s^m-1(h_i)-D_s^m-1(h_i-1)/μ_i|

  for even m.


  This Lemma follows by direct estimation. Therefore, let k∈ℤ/nℤ and note that ∑_j=1^n(D_s(h_j+1)-D_s(h_j))=0 as each term appears in the sum both in positive form and negative. Then, for odd m,
  
    
    |D_s^m-1(h_k+1)-D_s^m-1(h_k)/|e_k||
       = |1/ℓ(c)∑_j=1^n(D_s^m-1(h_j+1)-D_s^m-1(h_j))-D_s^m-1(h_k+1)-D_s^m-1(h_k)/|e_k||
       =|1/ℓ(c)∑_j=1^n(D_s^m-1(h_j+1)-D_s^m-1(h_j))-1/ℓ(c)∑_j=1^nD_s^m-1(h_k+1)-D_s^m-1(h_k)/|e_k||e_j||
       ≤|1/ℓ(c)∑_j=1^n(D_s^m-1(h_j+1)-D_s^m-1(h_j)/|e_j|-D_s^m-1(h_k+1)-D_s^m-1(h_k)/|e_k|)|e_j||
       ≤|1/ℓ(c)∑_j=1^n1/2(∑_i=k^j-1(D_s^m-1(h_i+i)-D_s^m-1(h_i)/|e_i|-D_s^m-1(h_i)-D_s^m-1(h_i-1)/|e_i-1|)..
       ..            -∑_i=j^k-1(D_s^m-1(h_i+1)-D_s^m-1(h_i)/|e_i|-D_s^m-1(h_i)-D_s^m-1(h_i-1)/|e_i-1|))|e_j||
       ≤1/2ℓ(c)∑_j=1^n∑_i=1^n|D_s^m-1(h_i+1)-D_s^m-1(h_i)/|e_i|-D_s^m-1(h_i)-D_s^m-1(h_i-1)/|e_i-1|||e_j|
       ≤1/2∑_i=1^n|D_s^m-1(h_i+1)-D_s^m-1(h_i)/|e_i|-D_s^m-1(h_i)-D_s^m-1(h_i-1)/|e_i-1||

  
  The case for even m is similar, swapping |e_i| for μ_i and rotating some indices.
  
    
    |D_s^m-1(h_k)-D_s^m-1(h_k-1)/μ_k|
       = |1/ℓ(c)∑_j=1^n(D_s^m-1(h_j)-D_s^m-1(h_j-1))-D_s^m-1(h_k)-D_s^m-1(h_k-1)/μ_k|
       =|1/ℓ(c)∑_j=1^n(D_s^m-1(h_j)-D_s^m-1(h_j-1))-1/ℓ(c)∑_j=1^nD_s^m-1(h_k)-D_s^m-1(h_k-1)/μ_kμ_j|
       ≤|1/ℓ(c)∑_j=1^n(D_s^m-1(h_j)-D_s^m-1(h_j-1)/μ_j-D_s^m-1(h_k)-D_s^m-1(h_k-1)/μ_k)μ_j|
       ≤|1/ℓ(c)∑_j=1^n1/2(∑_i=k^j-1(D_s^m-1(h_i+i)-D_s^m-1(h_i)/μ_i+1-D_s^m-1(h_i)-D_s^m-1(h_i-1)/μ_i)..
       ..            -∑_i=j^k-1(D_s^m-1(h_i+1)-D_s^m-1(h_i)/μ_i+1-D_s^m-1(h_i)-D_s^m-1(h_i-1)/μ_i))μ_j|
       ≤1/2ℓ(c)∑_j=1^n∑_i=1^n|D_s^m-1(h_i+1)-D_s^m-1(h_i)/μ_i-D_s^m-1(h_i)-D_s^m-1(h_i-1)/μ_i-1|μ_j
       ≤1/2∑_i=1^n|D_s^m-1(h_i+1)-D_s^m-1(h_i)/μ_i+1-D_s^m-1(h_i)-D_s^m-1(h_i-1)/μ_i|

  
  This proves the lemma.



  We consider the odd case which is not meaningfully different from the even one. Following from Lemma <ref>, we have
  
    max_i∈ℤ/nℤ|D^m-1(h_k+1)-D^m-1(h_k)/|e_k||^2
       ≤1/4(∑_i=1^n|D_s^m-1h_i+1-D_s^m-1(h_i)/|e_i|-D_s^m-1(h_i)-D_s^m-1(h_i-1)/|e_i-1||)^2
       ≤1/4(∑_i=1^n|D_s^m-1(h_i+1)-D_s^m-1(h_i)/|e_i|-D_s^m-1(h_i)-D_s^m-1(h_i-1)/|e_i-1||^2 1/μ_i)(∑_i=1^nμ_i )
       =ℓ(c)/4(∑_i=1^n|D_s^m-1(h_i+1)-D_s^m-1(h_i)/|e_i|-D_s^m-1(h_i)-D_s^m-1(h_i-1)/|e_i-1||^2 1/μ_i).

  
  Then, due to Hölder's inequality and the above, we conclude
  
    
    _c^m(h,h)   =∑_i=1^nℓ(c)^-3+2m|D_s^m(h_i)|^2|e_i|
       =∑_i=1^nℓ^-3+2m|D_s^m-1(h_i+1)-D_s^m-1(h_i)|^2/|e_i|
       ≤ℓ(c)^-3+2m(max_k∈𝒞_n|D_s^m-1(h_k+1)-D_s^m-1(h_k)/|e_k||^2)(∑_k=1^n|e_k|)
       ≤ℓ(c)^-3+2mℓ(c)^2/4∑_i=1^n|D_s^m-1(h_i+1)-D_s^m-1(h_i)/|e_i|-D_s^m-1(h_i)-D_s^m-1(h_i-1)/|e_i-1||^2 1/μ_i
       =ℓ(c)^-3+2(m+1)/4∑_i=1^n|D_s^m+1(h_i)|^2μ_i=1/4_c^m+1(h,h).

  
  This is symmetric with the argument for the odd case, which we omit.


abbrv

