---
abstract: |
    Counterexamples to some old-standing optimization problems in the smooth convex coercive setting are provided. We show that block-coordinate, steepest descent with exact search or Bregman descent methods do not generally converge. Other failures of various desirable features are established: directional convergence of Cauchy’s gradient curves, convergence of Newton’s flow, finite length of Tikhonov path, convergence of central paths, or smooth Kurdyka-Łojasiewicz inequality. All examples are planar.\
    These examples are based on general smooth convex interpolation results. Given a decreasing sequence of positively curved $C^k$ convex compact sets in the plane, we provide a level set interpolation of a $C^k$ smooth convex function where $k\geq2$ is arbitrary. If the intersection is reduced to one point our interpolant has positive definite Hessian, otherwise it is positive definite out of the solution set. Furthermore, given a sequence of decreasing polygons we provide an interpolant agreeing with the vertices and whose gradients coincide with prescribed normals.
author:
- 'Jérôme Bolte[^1] and Edouard Pauwels[^2]'
date: Draft of 
title: Curiosities and counterexamples in smooth convex optimization
---

Introduction
============

Questions and method
--------------------

One of the goals of convex optimization is to provide a solution to a problem with stable and fast algorithms. The quality of a method is generally assessed by the convergence of sequences, rate estimates, complexity bounds, finite length of relevant quantities and other quantitative or qualitative ways.

Positive results in this direction are numerous and have been the object of intense research since decades. To name but a few: gradient methods e.g., [@newmirovsky1983problem; @Nesterov; @Boyd], proximal methods e.g., [@PLC], alternating methods e.g., [@Beck; @wright2015coordinate], path following methods e.g., [@Aus99; @NN], Tikhonov regularization e.g. [@Golub], semi-algebraic optimization e.g., [@jon; @BNPS], decomposition methods e.g., [@PLC; @Beck], augmented Lagrangian methods e.g., [@Bertsek] and many others.

Despite this vast literature, some simple questions remain unanswered or just partly tackled, even for [*smooth convex coercive*]{} functions. Does the alternating minimization method, aka Gauss-Seidel method, converge? Does the steepest descent method with exact line search converge? Do mirror descent or Bregman methods converge? Does Newton’s flow converge? Do central paths converge? Is the gradient flow directionally stable? Do smooth convex functions have the Kurdyka-Łojasiewicz property?

In this article we provide a negative answer to all these questions.

Our work draws inspiration from early works of de Finetti [@definetti1949stratificazioni], Fenchel [@fenchel51], on convex interpolation, but also from Torralba’s PhD thesis [@torralba96] and the more recent [@bolte2010characterization], where some counterexamples on the Tikhonov path and Łojasiewicz inequalities are provided. The convex interpolation problem, see [@definetti1949stratificazioni], is as follows: given a monotone sequence of convex sets[^3] may we find a convex function interpolating each of these sets, i.e., having these sets as sublevel sets? Answers to these questions for [*continuous*]{} convex functions were provided by de Finetti, and improved by Fenchel [@fenchel51], Kannai [@kannai77], and then used in [@torralba96; @bolte2010characterization] for building counterexamples.

We improve this work by providing, for $k\geq 2$ arbitrary, a general $C^k$ interpolation theorem for positively curved convex sets, imposing at the same time the positive definiteness of its Hessian out of the solution set. An abridged version could be as follows.

Let $\left( T_i \right)_{i\in {\mathbb{Z}}}$ be a sequence of compact convex subsets of ${\mathbb{R}}^2$, with positively curved $C^k$ boundary, such that $T_i\subset{\mathrm{int}\,}T_{i+1}$ for all $i$ in ${\mathbb{Z}}$. Then there exists a $C^k$ convex function $f$ having the $T_i$ as sublevel sets with positive definite Hessian outside of the set: $$\operatorname*{argmin}f=\bigcap_{i\in{\mathbb{Z}}} T_i.$$

We provide several additional tools (derivatives computations) and variants (status of the solution set, Legendre functions, Lipschitz continuity). Whether our result is generalizable to general smooth convex sequences, i.e., with vanishing curvature, seems to be a very delicate question whose answer might well be negative.

Our central theoretical result is complemented by a discrete approximate interpolation result “of order one" which is particularly well adapted for building counterexamples. Given a nested collection of polygons, one can indeed build a smooth convex function having level sets interpolating its vertices and whose gradients agree with prescribed normals.

Our results are obtained by blending parametrization techniques, Minkovski summation, Bernstein approximations and convex analysis.

As sketched below, our results offer the possibility of building counterexamples in convex optimization by restricting oneself to the construction of countable collections of nested convex sets satisfying some desirable properties. In all cases failures of good properties are caused by some curvature oscillations.

A digest of counterexamples
---------------------------

Counterexamples provided in this article can be classified along three axes: structural counterexamples[^4], counterexamples for convex optimization algorithms and ordinary differential equations.

In the following, the term “nonconverging” sequence or trajectory means, a sequence or a trajectory with at least two distinct accumulation points. Unless otherwise stated, convex functions have full domain.

[*The following results are proved for $C^k$ convex functions on the plane with $k\geq 2$.*]{}

#### Structural counterexamples

-   **Kurdyka-Łojasiewicz:** There exists a $C^k$ convex function whose Hessian is positive definite outside its solution set and which does not satisfy the Kurdyka-Łojasiewicz inequality. This is an improvement on [@bolte2010characterization].

-   **Tikhonov regularization path:** There exists a $C^k$ convex function $f$ such that the regularization path $$\begin{aligned}
            x(r)= \operatorname*{argmin}\left\{ f(y) + r \|y\|^2:y\in {\mathbb{R}}^2\right\}, \,\,r\in(0,1)
        \end{aligned}$$ has infinite length. This strengthens a theorem by Torralba [@torralba96].

-   **Central path:** There exists a continuous Legendre function $h \colon [-1,1]^2 \mapsto {\mathbb{R}}$, $C^k$ on the interior, and $c$ in ${\mathbb{R}}^2$ such that the central path $$\begin{aligned}
            x(r) =
        \operatorname*{argmin}\left\{ \left\langle c, y \right\rangle + r h(y):y\in D\right\}
        \end{aligned}$$ does not have a limit as $r \to 0$.

#### Algorithmic counterexamples:

-   **Gauss-Seidel method (block coordinate descent):** There exists a $C^k$ convex function with positive definite Hessian outside its solution set and an initialization $ (u_0,v_0)$ in ${\mathbb{R}}^2$, such that the alternating minimization algorithm $$\begin{aligned}
            u_{i+1} &= \operatorname*{argmin}_{u \in {\mathbb{R}}} f(u, v_i) \\
            v_{i+1} &= \operatorname*{argmin}_{v \in {\mathbb{R}}}f(u_{i+1}, v) 
        \end{aligned}$$ produces a bounded nonconverging sequence $((u_i,v_i))_{i\in {\mathbb{N}}}$.

-   **Gradient descent with exact line search:** There exists a $C^k$ convex function $f$ with positive definite Hessian outside its solution set and an initialization $x_0$ in ${\mathbb{R}}^2$, such that the gradient descent algorithm with exact line search $$\begin{aligned}
            x_{i+1} &= \operatorname*{argmin}_{t \in {\mathbb{R}}} f(x_i + t \nabla f(x_i))
        \end{aligned}$$ produces a bounded nonconvergent sequence.

-   **Bregman or mirror descent method:** There exists a continuous Legendre function $h \colon [-1,1]^2 \mapsto {\mathbb{R}}$, $C^k$ on the interior, a vector $c$ in ${\mathbb{R}}^2$ and an initialization $x_0$ in $(-1,1)^2$ such that the Bregman recursion $$\begin{aligned}
            x_{i+1} = \nabla h^*(\nabla h(x_i) - c)
        \end{aligned}$$ produces a nonconverging sequence. The couple $(h,\langle c,\cdot\rangle$) satisfies the smoothness properties introduced in [@bauschke2016descent].

#### Continuous time ODE counterexamples:

-   **Continuous time Newton method:** There exists a $C^k$ convex function with positive definite Hessian outside its solution set, and an initialization $x_0$ in ${\mathbb{R}}^2$ such that the continuous time Newton’s system $$\begin{aligned}
            \dot{x}(t) &= - \left[\nabla^2 f(x(t))\right]^{-1} \nabla f(x(t)), \,\,t\geq 0,\\
            x(0) &= x_0
        \end{aligned}$$ has a solution approaching $\operatorname*{argmin}f$ which does not converge.

-   **Directional convergence for gradient curves:** There exists a $C^k$ convex function with $0$ as a unique minimizer and a positive definite Hessian on ${\mathbb{R}}^2\setminus\{0\}$, such that for any non stationary solution to the gradient system $$\begin{aligned}
            \dot{x}(t) &= - \nabla f(x(t)) 
        \end{aligned}$$ the direction $x(t) / \|x(t)\|$ does not converge.

-   **Hessian Riemannian gradient dynamics:** There exists a continuous Legendre function $h \colon [-1,1]^2 \mapsto {\mathbb{R}}$, $C^k$ on the interior, a linear function $f$ and a nonconvergent solution to the following system $$\begin{aligned}
            \dot{x}(t) = - \nabla_H f(x(t)),
        \end{aligned}$$ where $H = \nabla^2 h$ is the Hessian of $h$ and $\nabla_H f = H^{-1}  \nabla f$ is the gradient of $f$ in the Riemannian metric induced by $H$ on $(-1,1)^2$.

#### Pathological sequences and curves

Our counterexamples lead to sequences or paths in ${\mathbb{R}}^2$ which are related to a function $f$ by a certain property (see examples above) and have a certain type of pathology. For illustration purposes, we provide sketches of the pathological behaviors we met in Figure \[fig:pathologicalbehaviors\].

Preliminaries
=============

Let us set ${\mathbb{N}^*}={\mathbb{N}}\setminus\{0\}$. For $p$ in ${\mathbb{N}}^*$, the Euclidean scalar product in ${\mathbb{R}}^p$ is denoted by $\langle\cdot,\cdot\rangle$, otherwise stated the norm is denoted by $\|\cdot\|$. Given subsets $S,T$ in ${\mathbb{R}}^p$, and $x$ in ${\mathbb{R}}^p$, we define $${\mathrm{dist}}(x,S):=\inf \left\{\|x-y\|:y\in S\right\},$$ and the Hausdorff distance between $S$ and $T$, $${\mathrm{dist}}(S,T)=\max\left(\sup_{x\in S}{\mathrm{dist}}(x,T),\sup_{x\in T}{\mathrm{dist}}(x,S)\right).$$

Throughout this note, the assertion “$g$ is $C^k$ on $D$” where $D$ is not an open set is to be understood as “$g$ is $C^k$ on an open neighborhood of $D$”. Given a map $G \colon X \mapsto A \times B$ for some space $X$, $[G]_1 \colon X \mapsto A$ denotes the first component of $G$.

Continuous convex interpolations
--------------------------------

We consider a sequence of compact convex subsets of ${\mathbb{R}}^p$, $\left(T_i \right)_{i \in {\mathbb{Z}}}$ such that $T_{i+1} \subset \mathrm{int}\  T_i$. Finding a continuous convex interpolation of $\left(T_i \right)_{i \in {\mathbb{Z}}}$ is finding a convex continuous function which makes the $T_i$ a sequence of level sets. We call this process [*continuous convex interpolation*]{}. This questioning was present in Fenchel [@fenchel51] and dates back to de Finetti [@definetti1949stratificazioni], let us also mention the work of Crouzeix [@crouzeix80] revolving around this issue.

Such constructions have been shown to be realizable by Torralba [@torralba96], Kannai [@kannai77], using ideas based on Minkowski sum. The validity of this construction can be proved easily using the result of Crouzeix [@crouzeix80] which was already present under different and weaker forms in the works of de Finetti and Fenchel.

Let $f \colon {\mathbb{R}}^p \to {\mathbb{R}}$ be a quasiconvex function. The functions $$\begin{aligned}
                                F_x \colon \lambda \mapsto \sup\left\{  \langle z,x\rangle: f(z)\leq \lambda\right\}
                \end{aligned}$$ are concave for all $x$ in ${\mathbb{R}}^p$, if and only $f$ is convex. \[th:crouzeix\]

Our goal is to build smooth convex interpolation for sequences of smooth convex sets. To make such a construction we shall use [*nonlinear*]{} [*Minkowski*]{} [*interpolation*]{} between level sets.

We shall also rely on Bernstein approximation which we now describe.

Bernstein approximation
-----------------------

We refer to the monograph [@lorentz1954bernstein] by G.G. Lorentz.

The main properties of Bernstein polynomials to be used in this paper are the following:

-   Bernstein approximation is linear in its functional argument $f$ and “monotone” which allows to construct an approximation using only positive combination of a finite number of function values.

-   There are precise formulas for derivatives of Bernstein approximation. They involve repeated finite differences. So approximating piecewise affine function with high enough degree leads to an approximation for which corner values of derivatives are controlled while the remaining derivatives are vanishing (up to a given order).

-   Bernstein approximation is shape preserving, which means in particular that approximating a concave function preserves concavity.

The main idea to produce a smooth interpolation which preserves level sets is depicted in Figure \[fig:ideaBernstein\] where we use Bernstein approximation to interpolate smoothly between three points and controlling the successive derivatives at the end points of the interpolation.

Let us now be specific. Given $f$ defined on the interval $[0,1]$, the [*Bernstein polynomial*]{} of order $d \in {\mathbb{N}^*}$ associated to $f$ is given by $$\begin{aligned}
                B_d(x) = B_{d,f}(x) = \sum_{k=0}^d f\left( \frac{k}{d} \right) {d \choose k} x^k(1-x)^{d-k}, \mbox{ for $x\in [0,1].$}
                \label{eq:bernsteinPolynomial}\end{aligned}$$

#### Derivatives and shape preservation:

For any $h $ in $ (0,1)$ and $x $ in $ [0,1-h]$, we set $\Delta_h^1 f(x) = f(x+h) - f(x)$ and recursively for all $k $ in $ {\mathbb{N}^*}$, $\Delta_h^k f(x) = \Delta\left( \Delta_h^{k-1} f(x) \right)$. We fix $d \neq0$ in ${\mathbb{N}}$ and for $h = \frac{1}{d}$ write $\Delta_h^k = \Delta^k$. Then for any $m \leq d$, we have $$\begin{aligned}
                B_d(x)^{(m)} &= d(d-1)\ldots(d-m+1) \sum_{k=0}^{d - m} \Delta^m f\left( \frac{k}{d} \right) {{d-m} \choose k} x^k(1-x)^{d-k-m},
                \label{eq:bernsteinDerivative}\end{aligned}$$ for any $x$ in $[0,1]$. If $f$ is increasing (resp. strictly increasing), then $\Delta^1f(x) \geq 0$ (resp. $\Delta^1f(x) > 0$) for all $x$ and $B'_d$ is positive (resp. strictly positive) and $B_d$ is increasing (resp. strictly increasing). Similarly, if $f$ is concave, then $\Delta^2 f(x) \leq 0$ for all $x$ so that $B^{(2)}_d \leq 0$ and $B_d$ is concave. From , we infer $$\label{e:born}
  \left|B_d(x)^{(m)}\right|\leq d(d-1)\ldots(d-m+1)\sup_{k\in\{0,\ldots,d-m\}}  \left|\Delta^m f\left( \frac{k}{d} \right)\right|$$ for $x$ in $[0,1]$.

#### Approximation of piecewise affine functions:

The following lemma will be extensively used throughout the proofs.

\[lem:convexInterpol\] Let $q_0,q_1 \in {\mathbb{R}}^p$, $\lambda_-<\lambda_0< \lambda_1 < \lambda_+$ and $0<e_1,e_0<1$. Set $\Theta = \left( q_0,q_1,,\lambda_-,\lambda_0,\lambda_1,\lambda_+, e_0,e_1 \right)$ and define $\gamma_\Theta \colon [0,1] \mapsto {\mathbb{R}}^p$ through $$\begin{aligned}
\small
                                \gamma_\Theta(t)=
                                \begin{cases}
                                                q_0 \left( 1 + \frac{e_0}{\lambda_0 - \lambda_-} \left( t (\lambda_+- \lambda_-) \right) \right) & \text{ if } 0\leq t \leq \frac{\lambda_0 - \lambda_-}{\lambda_+ - \lambda_-} \\
                                                q_0(1 + e_0) \left( \frac{\lambda_1 - \lambda_- - t (\lambda_+ - \lambda_-)}{\lambda_1 - \lambda_0} \right) + q_1(1 - e_1) \left( \frac{ \lambda_- - \lambda_0  + t (\lambda_+ - \lambda_-)}{\lambda_1 - \lambda_0} \right)& \text{ if } \frac{\lambda_0 - \lambda_-}{\lambda_+ - \lambda_-} \leq t \leq \frac{\lambda_1 - \lambda_-}{\lambda_+ - \lambda_-} \\
                                                q_1\left( 1 + \frac{(t-1)(\lambda_+-\lambda_-)e_1}{\lambda_+ - \lambda_1} \right) & \text{ if } \frac{\lambda_1 - \lambda_-}{\lambda_+ - \lambda_-}  \leq t \leq 1.
                                \end{cases}
                \end{aligned}$$ The curve $\gamma_\Theta$ in ${\mathbb{R}}^{p+1}$ is the affine interpolant between the points $q_0$, $ (1 + e_0) q_0$, $(1 - e_1) q_1$ and $ q_1$. For any $m $ in $ {\mathbb{N}}$, we choose $d$ in ${\mathbb{N}}^*$ such that $$\begin{aligned}
    \label{e:deg}   \frac{m}{d} \leq \min\left\{ \frac{\lambda_0 - \lambda_-}{\lambda_+ - \lambda_-}, 1 - \frac{\lambda_1 - \lambda_-}{\lambda_+ - \lambda_-} \right\}.
                \end{aligned}$$ We consider a Bernstein-like reparametrization of $\tilde\gamma_\Theta $ given by $$\begin{aligned}
                                \tilde\gamma_\Theta \colon [\lambda_-,\lambda_+] &\mapsto {\mathbb{R}}^p\\
                                \lambda & \mapsto \sum_{k=0}^d \tilde\gamma_\Theta \left( \frac{k}{d} \right){d \choose k} \left(\frac{\lambda - \lambda_-}{\lambda_+ - \lambda_-}  \right)^k\left(1 - \frac{\lambda - \lambda_-}{\lambda_+ - \lambda_-}  \right)^{d-k}.
                \end{aligned}$$ Then the following holds, for any $2 \leq l \leq m$, $\tilde\gamma_\Theta $ is $C^m$ and $$\begin{aligned}
                                \tilde\gamma_\Theta (\lambda_-)&=q_0&
                                \tilde\gamma_\Theta (\lambda_+)&=q_1\\
                                \tilde\gamma_\Theta '(\lambda_-)&=\frac{e_0}{\lambda_0- \lambda-} q_0&
                                \tilde\gamma_\Theta '(\lambda_+)&=\frac{e_1}{\lambda_+ - \lambda_1} q_1\\
                                \tilde\gamma_\Theta ^{(l)}(\lambda_-) &=0  &\tilde\gamma_\Theta ^{(l)}(\lambda_+) &= 0.
                \end{aligned}$$ Furthermore, if $\gamma_\Theta $ has monotone coordinates (resp. strictly monotone, resp. concave, resp. convex), then so has $\tilde\gamma_\Theta $.

Note that the dependence of $\tilde\gamma_\Theta $ in $(q_0,q_1)$ is linear so that the dependence of $\tilde\gamma_\Theta $ in $(q_0,q_1)$ is also linear. Hence $\tilde\gamma_\Theta $ is of the form $\lambda \mapsto a(\lambda) q_0 + b(\lambda)q_1$. We can restrict ourselves to $p = 1$ since the general case follows from the univariate case applied coordinatewise.

If $p=1$, then $\tilde\gamma_\Theta  = B_{f_{\Theta,d}} \circ A$, where $A \colon \lambda \mapsto \frac{\lambda - \lambda_-}{\lambda_+-\lambda_-}$. We have $\tilde\gamma_\Theta (0) = q_0$, $\tilde\gamma_\Theta (1) = q_1$ and $\Delta^1 f_{\Theta}(0) = \frac{\lambda_+ - \lambda_-}{\lambda_0 - \lambda_-}\frac{e_0}{d} q_0$, $\Delta^1 f_{\Theta}\left( 1 - \frac{1}{d} \right) = \frac{\lambda_+ - \lambda_-}{\lambda_+ - \lambda_1}\frac{e_1}{d} q_1$ and $\Delta^{(l)} f_{\Theta}(0) = \Delta^{(l)} f_{\Theta}\left( 1 - \frac{l}{d} \right) =0$. The results follow from the expressions in and and the chain rule for $\tilde\gamma_\Theta  = B_{f_{\Theta,d}} \circ A$.

The last property of $\tilde\gamma_\Theta $ is due to the shape preserving property of Bernstein approximation and the fact that $\tilde\gamma_\Theta  = B_{f_{\Theta,d}} \circ L$.

\[rem:alignedInterpolation\] [(a) [**\[Affine image\]**]{} Using the notation of Lemma \[lem:convexInterpol\], if $(\lambda_-,q_0)$, $(\lambda_0, (1 + e_0) q_0)$, $(\lambda_1, (1 - e_1) q_1)$ and $(\lambda_+, q_1)$ are aligned, then the interpolation is actually an affine function.\
(b) [**\[Degree of the interpolants\]**]{} Observe that the degree of the Bernstein interpolant is connected to the slopes of the piecewise path $\lambda$ by .]{}

Smooth convex interpolation
===========================

Being given a subset $S$ of ${\mathbb{R}}^p$, we denote by $\mathrm{int}(S)$ its interior, $\bar S$ its closure and ${\mathrm{bd}\,}S=\bar S\setminus \mathrm{int}(S)$ its boundary. Let us recall that the [*support function*]{} of $S$ is defined through $$\sigma_S(x)=\sup\left\{\langle y,x\rangle :y\in S\right\}\in {\mathbb{R}}\cup\{+\infty\}.$$

Smooth parametrization of convex rings
--------------------------------------

A [*convex ring*]{} is a set of the form $C_1\setminus C_2$ where $C_1\subset C_2$ are convex sets. Providing adequate parameterizations for such objects is key for interpolating $C_1$ and $C_2$ by some (regular) convex function.

The following assertion plays a fundamental role.

Let $T_-,T_+ \subset {\mathbb{R}}^2$ be convex, compact with $C^k$ boundary ($k \geq 2$) and positive curvature. Assume that, $T_- \subset \mathrm{int} (T_+)$ and $0 \in \mathrm{int}(T_-)$. \[ass:curvature\]

The positive curvature assumption ensures that the boundaries can be parametrized by their normal, that is, for $i=-,+$, there exists $$\begin{aligned}
                c_i \colon {\mathbb{R}}/ 2 \pi {\mathbb{Z}}\mapsto {\mathrm{bd}\,}(T_i)\end{aligned}$$ such that the normal to $T_i$ at $c_i(\theta)$ is the vector $n(\theta) = (\cos(\theta),\sin(\theta))^T$ and $\dot{c}_i(\theta) = \rho_{i}(\theta) \tau(\theta)$ where $\rho_i > 0$ and $\tau(\theta) = (-\sin(\theta),\cos(\theta))$. In this setting, it holds that $c_i(\theta) = \mbox{argmax}_{y \in T_i} \left\langle n(\theta), y\right\rangle$. The map $c_i$ is the inverse of the Gauss map and is $C^{k-1}$ (see [@schneider1993convex] Section 2.5).

\[l:curv\] Let $T_-,T_+$ be as in Assumption \[ass:curvature\] with normal parametrizations as above. For $a,b\geq0$ with $a+b>0$ set $T=aT_-+bT_+$.

Then $T$ has positive curvature and its boundary is given by $${\mathrm{bd}\,}T=\{a\,c_-(\theta)+b\,c_+(\theta): \theta \in  {\mathbb{R}}/ 2 \pi {\mathbb{Z}}\},$$ with the natural parametrization ${{\mathbb{R}}/2\pi{\mathbb{Z}}\,}\ni \theta \to a\,c_-(\theta)+b\,c_+(\theta).$

We may assume $ab>0$ otherwise the result is obvious. Let $x$ be in ${\mathrm{bd}\,}T$ and denote by $n(\theta)$ the normal vector at $x$ for a well chosen $\theta$, so that $x={\mathrm{argmax}\,}\left\{\langle y,n(\theta)\rangle\right\}: y \in T\}$. Observe that the definition of the Minkowski sum yields $$\begin{aligned}
\max_{y\in T}\langle y,n(\theta)\rangle & = & \max_{(v,w)\in T_-\times T_+}\langle a v+bw,n(\theta)\rangle \\
&=& a \max_{v\in T_-} \langle v,n(\theta)\rangle+b\max_{w\in T_+} \langle w,n(\theta)\rangle 
 \end{aligned}$$ so that $$\begin{aligned}
\langle x,n(\theta)\rangle & =& a \langle c_-(\theta),n(\theta)\rangle+b\langle c_+(\theta),n(\theta)\rangle \\
& = &  \langle a c_-(\theta)+b c_+(\theta),n(\theta)\rangle \end{aligned}$$ which implies by extremality of $x$ that $x=  a c_-(\theta)+b c_+(\theta)$. Conversely, for any such $x$, $n(\theta)$ defines a supporting hyperplane to $T$ and $x$ must be on the boundary of $T$. The other results follow immediately.

In the following fundamental proposition, we provide a smooth parametrization of the convex ring $T_+\setminus {\mathrm{int}\,}T_-$. The major difficulty is to control tightly the derivatives at the boundary so that the parametrizations can be glued afterward to build smooth interpolants.

Let $T_-,T_+$ be as in Assumption \[ass:curvature\] with their normal parametrization as above. Fix $k\geq 2$, $\lambda_-< \lambda_0 < \lambda_1 < \lambda_+$ and $0 < e_0,e_1<1$. Choose $d $ in $ {\mathbb{N}}^*$, such that $$\begin{aligned}
                                \frac{k}{d} \leq \min\left\{ \frac{\lambda_0 - \lambda_-}{\lambda_+ - \lambda_-}, 1 - \frac{\lambda_1 - \lambda_-}{\lambda_+ - \lambda_-} \right\}.
                \end{aligned}$$ Consider the map $$\begin{aligned}
                                G \colon [\lambda_-,\lambda_+] \times {\mathbb{R}}/2\pi{\mathbb{Z}}&\mapsto {\mathbb{R}}^2 \nonumber\\
                                (\lambda, \theta) \mapsto \tilde\gamma_\Theta(\lambda)
                                \label{eq:mapG}
                \end{aligned}$$ with $\Theta = (c_-(\theta), c_+(\theta), \lambda_-, \lambda_0, \lambda_1, \lambda_+, e_0,e_1)$ and $\tilde\gamma$ as given by Lemma \[lem:convexInterpol\]. Assume further that:

${(\mathcal{M})}$ For any $\theta \in {\mathbb{R}}/ 2\pi{\mathbb{Z}}$, $\lambda \mapsto \left\langle G(\lambda, \theta), n(\theta)\right\rangle$ has strictly positive derivative on $[\lambda_-,\lambda_+]$.

Then the image of $G$ is ${\mathcal{R}}:=T_+ \setminus \mathrm{int}(T_-)$, $G$ is $C^k$ and satisfies, for any $2 \leq l \leq k$ and any $m $ in $ {\mathbb{N}}^*$, $$\begin{aligned}
                                \frac{\partial^m G}{\partial \theta^m}(\lambda_-,\theta) &= c_-^{(m)}(\theta)&
                                \frac{\partial^m G}{\partial \theta^m}(\lambda_+,\theta) &= c_+^{(m)}(\theta)\\
                                \frac{\partial^{m+1} G}{\partial \lambda\partial \theta^m} (\lambda_-,\theta) &= c_-^{(m)}(\theta) \frac{e_0}{\lambda_0 - \lambda_-}&
                                \frac{\partial^{m+1} G}{\partial \lambda\partial \theta^m} (\lambda_+,\theta) &= c_+^{(m)}(\theta) \frac{e_1}{\lambda_+-\lambda_1}\\
                                \frac{\partial^{l+m} G}{\partial \lambda^l \partial \theta^m }(\lambda_-,\theta) &= 0&
                                \frac{\partial^{l+m} G}{\partial \lambda^l \partial \theta^m }(\lambda_+,\theta) &= 0.
                \end{aligned}$$ Besides $G$ is a diffeomorphism from its domain on to its image. Set ${\mathcal{R}}\ni x \mapsto (f(x), \theta(x))$ to be the inverse of $G$. Then $f$ is $C^k$ and in addition, for all $x$ in $ {\mathbb{R}}$, $$\begin{aligned}
                    \nabla f(x) =\quad& \frac{1}{\left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle} \, n(\theta(x))\nonumber\\
                    \nabla \theta(x) =\quad& \frac{1}{\left\langle\frac{\partial G}{\partial \theta}(f(x),\theta(x)),\tau(\theta(x))\right\rangle} \tau(\theta(x)) \nonumber\\
                    &-  \frac{\left\langle \frac{\partial G}{\partial \lambda} (f(x),\theta(x)), \tau(\theta(x))\right\rangle}{\left\langle \frac{\partial G}{\partial \lambda} (f(x),\theta(x)), n(\theta(x))\right\rangle \left\langle\frac{\partial G}{\partial \theta}(f(x),\theta(x)),\tau(\theta(x))\right\rangle   } n(\theta(x)) \nonumber\\
                    \nabla^2 f(x)=\quad&\frac{\left\langle\frac{\partial G}{\partial \theta}(f(x),\theta(x)),\tau(\theta(x))\right\rangle}{\left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle}
                             \nabla \theta(x) \nabla \theta (x)^T \nonumber\\
                             &-   \frac{\left\langle  \frac{\partial^2 G}{\partial \lambda^2}(f(x),\theta(x)), n(\theta(x))\right\rangle}{\left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle} \nabla f(x) \nabla f(x)^T,
                    \label{eq:gradient} 
                \end{aligned}$$ where all denominators are positive. \[prop:diffeoMorphism\]

[ Note that $G$ is actually well defined and smooth on an open set containing its domain. As we shall see it is also a diffeomorphism from an open set containing its domain onto its image.]{}

Note that by construction, we have $G(\lambda,\theta) = a(\lambda) c_-(\theta) + b(\lambda)c_+(\theta)$ for some polynomials $a$ and $b$ which are nonnegative on $[\lambda_-,\lambda_+]$. The formulas for the derivatives follow easily from this remark, the form of $a$ and $b$ and Lemma \[lem:convexInterpol\].

Set, for any $\lambda $ in $ [\lambda_-,\lambda_+]$, $T_\lambda = a(\lambda) T_- + b(\lambda) T_+$. The resulting set $T_\lambda$ is convex and has a positive curvature by Lemma \[l:curv\], and for $\lambda$ fixed $G(\lambda,\cdot)$ is the inverse of the Gauss map of $T_\lambda$, which constitutes a parametrization by normals of the boundary.

Assume that $\lambda < \lambda'$, using the monotonicity assumption ${(\mathcal{M})}$, we have for any $\theta,\theta'$, $$\begin{aligned}
                                \left\langle n(\theta'), G(\lambda,\theta) \right\rangle& \leq \sup_{y \in T_{\lambda}} \left\langle n(\theta'),y \right\rangle \\
                                & = \left\langle n(\theta'), G(\lambda,\theta') \right\rangle\\
                                &< \left\langle n(\theta'), G(\lambda',\theta') \right\rangle
                \end{aligned}$$ so that $G(\lambda,\theta) \neq G(\lambda',\theta')$. Furthermore, we have by definition of $G(\lambda',\theta')$ $$\begin{aligned}
                                G(\lambda,\theta) \in \bigcap_{\theta' \in {\mathbb{R}}/2\pi{\mathbb{Z}}} \left\{ y,\; \left\langle y, n(\theta')\right\rangle \leq \left\langle n(\theta'), G(\lambda',\theta') \right\rangle \right\} = T_{\lambda'},   
                \end{aligned}$$ where the equality follows from the convexity of $T_{\lambda'}$. By convexity and compactness, this entails that $T_\lambda = \mathrm{conv}({\mathrm{bd}\,}(T_\lambda)) \subset {\mathrm{int}\,}T_{\lambda'}$.

Let us show that the map $G$ is bijective, first consider proving surjectivity. Let $f$ be defined on $T_+ \setminus \mathrm{int}(T_-)$ through $$\begin{aligned}
                                f \colon x \mapsto \inf_{}\left\{ \lambda : \lambda\geq \lambda_-, \; x \in T_\lambda = a(\lambda) T_- + b(\lambda) T_+\right\}.
            \label{eq:defIntropol}
                \end{aligned}$$ Since $a(\lambda_+)=0,  b(\lambda_+)=1$ this function is well defined and by compactness and continuity the infimum is achieved. It must hold that $x $ belongs to $ {\mathrm{bd}\,}(T_{f(x)})$, indeed, if $f(x) = \lambda_-$, then $x $ belongs to $ {\mathrm{bd}\,}(T_-)$ and otherwise, if $x $ is in $ \mathrm{int}(T_{\lambda'})$ for $\lambda' > \lambda_-$, then $f(x) < \lambda'$. We deduce that $x$ is of the form $G(f(x),\theta)$ for a certain value of $\theta$, so that $G$ is surjective.

As for injectivity, we have already seen a first case, the monotonicity assumption ${(\mathcal{M})}$ ensures that $\lambda \neq \lambda'$ implies $G(\lambda,\theta) \neq G(\lambda',\theta')$ for any $\theta,\theta'$. Furthermore, we have the second case, for any $\lambda $ in $ [\lambda_-,\lambda_+]$ and any $\theta$, $G(\lambda,\theta) = \arg\max_{y\in T_\lambda} \left\langle y,n(\theta) \right\rangle$ so that $\theta \neq \theta'$ implies $G(\lambda,\theta) \neq G(\lambda,\theta')$. So in all cases, $(\lambda,\theta) \neq (\lambda',\theta')$ implies that $G(\lambda,\theta) \neq G(\lambda',\theta')$ and $G$ is injective.

Let us now show that the map $G$ is a local diffeomorphism by estimating its Jacobian map.

Since $0 \in \mathrm{int}(T_-)$, we have for any $\lambda,\theta$, $$\begin{aligned}
                                0 &< \sup_{y \in T_-} \left\langle y, n(\theta) \right\rangle \\                                
                                &= \left\langle  G(\lambda_-,\theta),n(\theta) \right\rangle \\                             
                                & \leq a(\lambda) \left\langle  c_-(\theta), n(\theta)\right\rangle + b(\lambda) \left\langle c_+(\theta), n(\theta)  \right\rangle,
                \end{aligned}$$ and both scalar products are positive so that $a(\lambda) + b(\lambda)  > 0$. Hence, for any $\theta $ in $ {\mathbb{R}}/ 2\pi {\mathbb{Z}}$, $$\begin{aligned}
                \label{eq:derivativeGtheta}
                                \frac{\partial G}{\partial \theta}(\lambda,\theta) = (a(\lambda) \rho_-(\theta) + b(\lambda) \rho_+ (\theta)) \tau(\theta),
                \end{aligned}$$ with $a(\lambda) \rho_-(\theta) + b(\lambda) \rho_+ (\theta) > 0$. Furthermore by assumption $\lambda\to  \max_{y \in T_\lambda} \left\langle y, n(\theta) \right\rangle=\left\langle  G(\lambda,\theta),n(\theta) \right\rangle$ has strictly positive derivative in $\lambda$, whence $$\begin{aligned}
                                \left\langle \frac{\partial G}{\partial \lambda} (\lambda,\theta), n(\theta)\right\rangle > 0.
                \end{aligned}$$ We deduce that for any fixed $\theta$, in the basis $(n(\theta), \tau(\theta))$, the Jacobian of $G$, denoted $J_G$, is triangular with positive diagonal entries. More precisely fix $\lambda, \theta$ and set $x = G(\lambda, \theta)$ such that $\lambda = f(x)$, $\theta = \theta(x)$. In the basis $(n(\theta), \tau(\theta))$, we deduce from that the Jacobian of $G$ is of the form $$\begin{aligned}
                    J_G(\lambda, \theta) = 
                    \begin{pmatrix}
                        \alpha &0\\
                        \gamma&\beta
                    \end{pmatrix},
                \end{aligned}$$ where $$\begin{aligned}
                        \alpha &= \left\langle \frac{\partial G}{\partial \lambda} (\lambda,\theta), n(\theta)\right\rangle>0  \\
                        \beta &= \left\langle \frac{\partial G}{\partial \theta} (\lambda,\theta), \tau(\theta)\right\rangle = a(\lambda) \rho_-(\theta) + b(\lambda) \rho_+ (\theta)   >0\\
                        \gamma &=  \left\langle \frac{\partial G}{\partial \lambda} (\lambda,\theta), \tau(\theta)\right\rangle.
                \end{aligned}$$ It is thus invertible and we have a local diffeomorphism. We deduce that $$\begin{aligned}
                    J_G(\lambda, \theta)^{-1} = 
                    \begin{pmatrix}
                        \alpha^{-1} &0\\
                        \frac{-\gamma}{\alpha\beta}&\beta^{-1}
                    \end{pmatrix}.
                \end{aligned}$$ We have $J_G(\lambda,\theta)^{-1} = J_{G^{-1}}(x)$ so that the first line is $\nabla f(x)$ and second line is $\nabla \theta(x)$, which proves the claimed expressions for gradients.

We also have $d n(\theta) / d\theta = \tau(\theta)$ so that $$\begin{aligned}
                    J_{n \circ \theta}(x) = \tau(\theta(x)) \nabla \theta(x)^T.
                \end{aligned}$$ Differentiating the gradient expression, we obtain ($\nabla$ denotes gradient with respect to $x$): $$\begin{aligned}
                    &\nabla^2 f(x)\\
                    =\quad& \frac{\partial}{\partial  x} \left(\frac{1}{\left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle} \, n(\theta(x))\right)\\
                    =\quad& \frac{1}{\left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle} \, J_{n\circ \theta}(x) \\
                    &+ n(\theta(x)) \nabla \left(\frac{1}{\left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle }\right)^T  \\
                    =\quad& \frac{1}{\left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle} \, J_{n\circ \theta}(x) \\
                    &- n(\theta(x)) \nabla \left(\left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle \right)^T \frac{1}{\left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle^2} \\
                    =\quad & \frac{1}{\left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle} \left( \tau(\theta) \nabla \theta(x)^T - \nabla f(x)  \nabla \left(\left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle \right)^T\right).
                \end{aligned}$$ We have $$\begin{aligned}
                     =\quad&\nabla \left(\left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle \right)^T \\
                     =\quad & \frac{\partial G}{\partial \lambda} (f(x), \theta(x))^T J_{n\circ\theta}(x) + n(\theta(x))^T J_{\frac{\partial G}{\partial \lambda}}(f(x),\theta(x)) J_{G^{-1}} (x)  \\
                     =\quad & \left\langle \frac{\partial G}{\partial \lambda}(f(x), \theta(x)), \tau(\theta(x)) \right\rangle \nabla \theta(x)^T + \left\langle n (\theta(x)), \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta)\right\rangle \nabla f(x)^T,
                \end{aligned}$$ where, for the last identity, we have used the fact that $$\begin{aligned}
                    n(\theta(x))^T J_{\frac{\partial G}{\partial \lambda}}(f(x),\theta(x)) &= \left\langle n (\theta(x)), \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta)\right\rangle n(\theta(x))^T + \left\langle n (\theta(x)), \frac{\partial^2 G}{\partial \lambda \partial \theta}(\lambda,\theta)\right\rangle \tau(\theta(x))^T\\
                    &= \left\langle n (\theta(x)), \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta)\right\rangle n(\theta(x))^T\\
                    n(\theta(x))^T J_{G^{-1}}(x) &= n(\theta(x))^T  J_G(f(x),\theta(x))^{-1}\\
                    &= \frac{1}{\left\langle \frac{\partial G}{\partial \lambda} (f(x),\theta(x)), n(\theta(x))\right\rangle} n(\theta(x))^T = \nabla f(x)^T.
                \end{aligned}$$ We deduce that $$\begin{aligned}
                    &\nabla^2 f(x)\\
                    =\quad & \frac{1}{\left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle}\\
                    &\Bigg( \tau(\theta(x)) \nabla \theta(x)^T - \left\langle \frac{\partial G}{\partial \lambda}(f(x), \theta(x)), \tau(\theta(x)) \right\rangle \nabla f(x) \nabla \theta(x)^T  \\
                    &-  \nabla f(x) \nabla f(x)^T \left\langle n (\theta(x)), \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta)\right\rangle \Bigg).
                \end{aligned}$$ We have $$\begin{aligned}
                   & \tau(\theta(x)) \nabla \theta(x)^T - \left\langle \frac{\partial G}{\partial \lambda}(f(x), \theta(x)), \tau(\theta(x)) \right\rangle \nabla f(x) \nabla \theta(x)^T \\
                    =\quad& \left(\tau(\theta(x)) - \frac{\left\langle \frac{\partial G}{\partial \lambda}(f(x), \theta(x)), \tau(\theta(x)) \right\rangle}{\left\langle \frac{\partial G}{\partial \lambda}(f(x), \theta(x)), n(\theta(x)) \right\rangle} n(\theta(x)) \right) \nabla \theta(x)^T\\
                    =\quad& (a(\lambda) \rho_-(\theta) + b(\lambda) \rho_+ (\theta)) \nabla \theta(x) \nabla \theta (x)^T
                \end{aligned}$$ So that we actually get $$\begin{aligned}
                    &\nabla^2 f(x) \left\langle \frac{\partial G}{\partial \lambda} (f(x), \theta(x)), n(\theta(x)) \right\rangle\\
                    =\quad &
                    (a(\lambda) \rho_-(\theta) + b(\lambda) \rho_+ (\theta)) \nabla \theta(x) \nabla \theta (x)^T -   \left\langle n (\theta(x)), \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta)\right\rangle \nabla f(x) \nabla f(x)^T\\
                    =\quad &\left\langle\frac{\partial G}{\partial \theta}(\lambda,\theta),\tau(\theta(x))\right\rangle
                     \nabla \theta(x) \nabla \theta (x)^T -   \left\langle n (\theta(x)), \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta)\right\rangle \nabla f(x) \nabla f(x)^T.
                \end{aligned}$$ This concludes the proof.

Smooth convex interpolation of smooth positively curved convex sequences
------------------------------------------------------------------------

In this section, we consider an indexing set $I$ with either $I = {\mathbb{N}}$ or $I = {\mathbb{Z}}$, and an increasing sequence of compact convex sets $\left( T_i \right)_{i\in I}$ such that for any $i $ in $ I$, the couple $T_+ := T_{i+1}$, $T_-: = T_i$ satisfies Assumption \[ass:curvature\]. In particular, for each $i $ in $ I$, $T_i$ is compact, convex with $C^k$ boundary and positive curvature. We denote by $c_i$ the corresponding parametrization by the normal, $T_i \subset \mathrm{int}\, T_{i+1}$. With no loss of generality we assume $\displaystyle 0 \in \cap_{i \in I} T_i$.

This is our main theoretical result.

\[th:smoothinterp\] Let $I = {\mathbb{N}}$ or $I = {\mathbb{Z}}$ and $\left( T_i \right)_{i\in I}$ such that for any $i \in I$, $T_i \subset {\mathbb{R}}^2$ and the couple $T_+ := T_{i+1}$, $T_-: = T_i$ satisfies Assumption \[ass:curvature\]. Then there exists a $C^k$ convex function $$\begin{aligned}
                                f \colon {\mathcal{T}}:= \mathrm{int}\left(\bigcup_{i\in I} T_i   \right) \mapsto {\mathbb{R}}\end{aligned}$$ such that

\(i) $T_i$ is a sublevel set of $f$ for all $i $ in $ I$.

\(ii) We have $$\begin{aligned}
                    \operatorname*{argmin}f= \begin{cases}
                    \bigcap_I T_i &\quad \text{ if } I = {\mathbb{Z}}\\
                    \{0\}&\quad \text{ if } I = {\mathbb{N}}.
                \end{cases}
                \end{aligned}$$

\(iii) $\nabla^2 f $ is positive definite on ${\mathcal{T}}\setminus \operatorname*{argmin}f$, and if $I = {\mathbb{N}}$, it is positive definite throughout ${\mathcal{T}}$.

$\:$\
[*Preconditionning.*]{} We have that $0 \in \cap_{i \in I } \mathrm{int}(T_i)$. Hence for any $i $ in $ I$ and $0 \leq \alpha <1$, $\alpha T_i \in \mathrm{int}(T_i)$. Furthermore, for $\alpha>0$, small enough, $(1 + \alpha)T_{i} \subset (1-\alpha) \mathrm{int}(T_{i+1})$. Set $\alpha_0$ such that $(1 + \alpha_0)T_0 \subset (1-\alpha_0) \mathrm{int} (T_1)$. By forward (and backward if $I = {\mathbb{Z}}$) induction, for all $i  $ in $ I$, we obtain $\alpha_i>0$ such that $$\begin{aligned}
                                (1 + \alpha_i) T_i &\subset (1-\alpha_{i}) \mathrm{int} (T_{i+1}).
                \end{aligned}$$ Setting for all $i $ in $ I$, $\epsilon_{i+1} = \min\{\alpha_i,\alpha_{i+1}\}$ ($\epsilon_0 = \alpha_0$ if $I = {\mathbb{N}}$), we have for all $i $ in $ I$ $$\begin{aligned}
                    (1 + \epsilon_i) T_i \subset (1 + \alpha_i) T_i &\subset (1-\alpha_{i}) \mathrm{int} (T_{i+1}) \subset(1-\epsilon_{i+1}) \mathrm{int} (T_{i+1}).
                \end{aligned}$$ For all $i $ in $ I$, we introduce $$\begin{aligned}
                S_{3i} &= T_i,\\
                S_{3i + 1} &= (1 + \epsilon_i) T_i\\
                S_{3i + 2} &= (1 - \epsilon_{i+1}) T_{i+1}.\end{aligned}$$

We have a new sequence of strictly increasing compact convex sets $\left( S_i\right)_{i \in I }$.

[*Value assignation.*]{} For each $i $ in $ I $, we set $$\begin{aligned}
                K_i = \max_{\|x\| = 1} \frac{\sigma_{S_{i+1}}(x)  - \sigma_{S_i}(x)}{ \sigma_{S_{i}}(x)  - \sigma_{S_{i-1}}(x)} \in \left( 0, + \infty \right).\end{aligned}$$ Note that for all $i $ in $ I$, $K_{3i} = 1$. We choose $\lambda_{1} = 2$, $\lambda_0 = 1$ and for all $i $ in $ I $, $$\label{value}
                \lambda_{i+1} = \lambda_i + K_i(\lambda_i - \lambda_{i-1}).$$ By construction, we have for all $i $ in $ I $ and all $\theta \in {\mathbb{R}}/ 2 \pi {\mathbb{Z}}$, $$\begin{aligned}
                \frac{\sigma_{S_{i+1}}(n(\theta))   - \sigma_{S_{i}}(n(\theta))}{\lambda_{i+1} - \lambda_{i}} \leq  \frac{\sigma_{S_{i}}(n(\theta)) - \sigma_{S_{i-1}}(n(\theta))  }{\lambda_{i} - \lambda_{i-1}}.\end{aligned}$$ If $I = {\mathbb{Z}}$, this entails $$\begin{aligned}
                0 < \lambda_i - \lambda_{i - 1} \leq \frac{\lambda_{1}- \lambda_0}{\sigma_{S_{1}}(n(\theta)) - \sigma_{S_{0}}(n(\theta))} (\sigma_{S_{i}}(n(\theta)) - \sigma_{S_{i-1}}(n(\theta))),\end{aligned}$$ and the right-hand side is summable over negative indices $i\leq 0$, so that $\lambda_{i} \to \underline{\lambda} \in {\mathbb{R}}$ as $i \to -\infty$. In all cases $(\lambda_i)_{i\in I}$ is an increasing sequence bounded from below.

[*Local interpolation.*]{} We fix $i $ in $ I $ and consider the function $G_i$ described in Proposition \[prop:diffeoMorphism\] with $T_+ = S_{3i + 3} = T_{i+1}$, $T_- = S_{3i} = T_i$, $\lambda_+ = \lambda_{3i+3}$, $\lambda_{1} = \lambda_{3i+2}$, $\lambda_0 = \lambda_{3i+1}$, $\lambda_- = \lambda_{3i}$, $e_0 = \epsilon_i$, $e_1 = \epsilon_{i+1}$. By linearity, we have for any $(\lambda,\theta) \in [\lambda_-,\lambda_+] \times {\mathbb{R}}/ 2\pi {\mathbb{Z}}$, $$\begin{aligned}
                \left\langle G_i(\lambda,\theta),n(\theta)\right\rangle = \tilde \gamma_{\Theta}(\lambda)\end{aligned}$$ where $ \tilde \gamma_{\Theta}$ is as in Lemma \[lem:convexInterpol\] with input data $q_0 = \left\langle c_{3i}(\theta), n(\theta)\right\rangle = \sigma_{S_{3i}}(n(\theta))$, $q_1 = \left\langle c_{3i+3}(\theta), n(\theta)\right\rangle = \sigma_{S_{3i+3}}(n(\theta))$, and $\lambda_-,\lambda_0,\lambda_1,\lambda_+,e_1,e_0$ as already described. This corresponds to the Bernstein approximation of the piecewise affine interpolation between the points $$\begin{aligned}
                &\left( \lambda_{3i}, \sigma_{S_{3i}}(n(\theta))\right)\nonumber\\
                &\left( \lambda_{3i+1}, \sigma_{S_{3i+1}}(n(\theta))\right),\nonumber \\
                &\left( \lambda_{3i+2}, \sigma_{S_{3i+2}}(n(\theta))\right),\nonumber \\
                &(\lambda_{3i+3}, \sigma_{3i+3}(n(\theta))),    \label{eq:piecewiseAffineInterpolExplicit}\end{aligned}$$ By construction of $\left( K_i \right)_{i \in I }$, we have for all $\theta$, $$\begin{aligned}
                0 < \frac{\sigma_{S_{3i+3}}(n(\theta)) - \sigma_{S_{3i+2}}(n(\theta))  }{\lambda_{3i+3} - \lambda_{3i + 2}} \leq    \frac{\sigma_{S_{3i+2}}(n(\theta)) - \sigma_{S_{3i+1}}(n(\theta))  }{\lambda_{3i+2} - \lambda_{3i + 1}} \leq \frac{\sigma_{S_{3i+1}}(n(\theta)) - \sigma_{S_{3i}}(n(\theta))  }{\lambda_{3i+1} - \lambda_{3i}}.\end{aligned}$$ Whence the affine interpolant between points in is strictly increasing and concave, and by using the shape preserving properties of Bernstein polynomials, $\left\langle G_i(\lambda,\theta),n(\theta)\right\rangle$ has strictly positive derivative. As a consequence $G_i$ is a diffeomorphism and its derivatives are as in Proposition \[prop:diffeoMorphism\]. Furthermore $$\begin{aligned}
                \lambda \mapsto \left\langle G_i(\lambda,\theta),n(\theta)\right\rangle  \end{aligned}$$ is a $C^k$ concave function of $\lambda$.

[*Global interpolation.*]{} Recall that $\underline{\lambda} = \inf_{i \in I} \lambda_i>-\infty$ and set $\bar{\lambda} = \sup_{i \in I} \lambda_i\in (-\infty,+\infty]$. For any $\lambda \in (\underline{\lambda},\bar{\lambda})$, there exists a unique $i_\lambda \in I$ such that $\lambda \in [\lambda_{3i_\lambda}, \lambda_{3i_\lambda+3})$. Define $$\begin{aligned}
    G \colon (\underline{\lambda}, \bar{\lambda}) \times {\mathbb{R}}/ 2\pi {\mathbb{Z}}&\mapsto {\mathbb{R}}^2 \\
    (\lambda,\theta) &\mapsto G_{i_\lambda} (\lambda, \theta).\end{aligned}$$ Fix $i$ in $I$. The boundary of $T_{i+1}$ is given by $G_{i+1}(\lambda_{3i+3}, {\mathbb{R}}/ 2 \pi {\mathbb{Z}}) = G_{i}(\lambda_{3i+3}, {\mathbb{R}}/ 2 \pi {\mathbb{Z}})$ with actually $$\label{coinc}
G_{i+1}(\lambda_{3i+3}, \theta) = G_{i}(\lambda_{3i+3}, \theta ), \mbox{ for all }\theta \mbox{ in }{\mathbb{R}}/ 2 \pi {\mathbb{Z}}.$$

Since $K_{3i} = 1$, we have $$\begin{aligned}
                \lambda_{3i+1} - \lambda_{3i} = \lambda_{3i} - \lambda_{3i-1}.\end{aligned}$$ The expressions of the derivatives in Proposition \[prop:diffeoMorphism\] and ensure that the derivatives of $G_{i+1}$ and $G_{i}$ agree on ${\lambda_{3i+3}} \times {\mathbb{R}}/ 2 \pi {\mathbb{Z}}$ up to order $k$. Hence $G$ is a local diffeomorphism. Bijectivity of each $G_i$ ensure that $G$ is also bijective and thus $G$ is a diffeomorphism. Furthermore $$\begin{aligned}
                \lambda \mapsto \left\langle G(\lambda,\theta),n(\theta)\right\rangle  \end{aligned}$$ is $C^k$ piecewise concave and thus concave.

[*Extending $G$.*]{} If $I = {\mathbb{N}}$, we may assume without loss of generality that $S_0 = B$ the Euclidean ball and $S_1 = 5/3 S_0$, which corresponds to $\epsilon_0 = 2/3$, eventually after adding a set in the list and rescaling. Let $\phi$ denote the function described in Lemma \[lem:interpolationAroundZero\] and $G_{-1}$ be described as in Lemma \[lem:diffGauge2\]. This allows to extend $G$ for $\lambda \in [0,1]$, $G$ is then $C^k$ on $(0,\bar{\lambda}) \times{\mathbb{R}}/ 2\pi {\mathbb{Z}}$ by using Lemma \[lem:diffGauge2\] and Proposition \[prop:diffeoMorphism\]. This does not affect the differentiability, monotonicity and concavity properties of $G$.

[*Defining the interpolant $f$.*]{} We assume without loss of generality that $\underline{\lambda} = 0$. We set $f$ to be the first component of the inverse of $G$ so that it is defined on $G^{-1}\left( (0,\bar{\lambda}) \times {\mathbb{R}}/ 2\pi {\mathbb{Z}}\right)$. We extend $f$ as follows:

-   $f(0) = 0$ if $I = {\mathbb{N}}$,

-   $f = 0$ on $\cap_{i \in I } T_i$ if $I = {\mathbb{Z}}$.

Since $G$ is $C^k$ and non-singular on $ (0,\bar{\lambda}) \times {\mathbb{R}}/ 2\pi {\mathbb{Z}}$, the inverse mapping theorem ensures that $f$ is $C^k$ on $\mathrm{int}({\mathcal{T}}) \setminus \arg\min_{{\mathcal{T}}} f$.

[*Convexity of $f$.*]{} For any $\theta $ in $ {\mathbb{R}}/ 2\pi {\mathbb{Z}}$, $$\begin{aligned}
(0, \bar{\lambda}) &\mapsto {\mathbb{R}}_+\\
                \lambda &\mapsto \sup_{z\in[f \leq \lambda]} n(\theta)^Tz \end{aligned}$$ is equal to $\left\langle G(\lambda,\theta),n(\theta)\right\rangle$ which is concave. It can be extended at $\lambda = 0$ by continuity. This preserves concavity hence, using Theorem \[th:crouzeix\], we have proved that $f$ is convex and $C^k$ on ${\mathcal{T}}\setminus \operatorname*{argmin}_{{\mathcal{T}}}f$.

[*Smoothness around the argmin and Hessian positivity.*]{} If $I = {\mathbb{N}}$, then the interpolant defined in Lemma \[lem:diffGauge2\] ensures that $f$ is proportional to the norm squared around $0$. Hence it is $C^k$ around $0$ with positive definite Hessian. We may compose $f$ with the function $g \colon t \mapsto \sqrt{t^2 + 1} + t$ which is increasing and has positive second derivative. This ensures that the resulting Hessian is positive definite outside $\operatorname*{argmin}f$ and thus everywhere since $$\begin{aligned}
                \nabla^2 g \circ f = g' \nabla^2 f + g'' \nabla f \nabla f^T\end{aligned}$$ is positive definite thanks to the expressions for the Hessian of $f$ in Proposition \[prop:diffeoMorphism\]

If $I = {\mathbb{Z}}$, we let all the derivatives of $f$ vanish around the solution set. The smoothing Lemma \[lem:CkSmoothing\] applies and provides a function $\phi$ with positive derivative on $(0,+\infty)$, such that $\phi \circ f$ is convex, $C^k$ with prescribed sublevel sets. Furthermore, we remark that $$\begin{aligned}
                \nabla^2 \phi \circ f = \phi' \nabla^2 f + \phi'' \nabla f \nabla f^T.\end{aligned}$$ We may compose $\phi \circ f$ with the function $g \colon t \mapsto \sqrt{t^2 + 1} + t$ which is increasing with positive second derivative, the expressions for the Hessian of $f$ in Proposition \[prop:diffeoMorphism\] ensure that the resulting Hessian is positive definite out of $\operatorname*{argmin}f$.

\[rem:alignedLevelSets\][In view of Remark \[rem:alignedInterpolation\], if we have $T_{i+1} = \alpha T_i$ for some $0<\alpha < 1$ and $i $ in $ {\mathbb{Z}}$, then the interpolated level sets between $T_i$ and $T_{i+1}$ are all of the form $s T_i$ for $\alpha \leq s \leq 1$.]{}

\[rem:strictConvexity\][ Recall that strict convexity of a differentiable function amounts to the injectivity of its gradient. In Theorem \[th:smoothinterp\] if there is a unique minimizer, then the invertibility of the Hessian outside argmin $f$ ensures that our interpolant is strictly convex (note that this is automatically the case if $I = {\mathbb{N}}$).]{}

Considerations on Legendre functions {#s:legendre}
------------------------------------

The following proposition provides some interpolant with additional properties as global Lipschitz continuity and finiteness properties for the dual function. At this stage these results appear as merely technical but they happen to be decisive in the construction of counterexamples involving Legendre functions. The properties of Legendre functions can be found in [@rockafellar1970convex Chapter 6]. We simply recall here that, given a convex body $C$ of ${\mathbb{R}}^p$, a convex function $h:C\to {\mathbb{R}}$ is [*Legendre*]{} if it is differentiable on ${\mathrm{int}\,}C$ and if $\nabla h$ defines a bijection from ${\mathrm{int}\,}C$ to $\nabla h({\mathrm{int}\,}C)$ with in addition $$\lim_{\begin{array}{l}
x\in {\mathrm{int}\,}C\\
x\to z\end{array}} \|\nabla h(x)\|=+\infty,$$ for all $z$ in ${\mathrm{bd}\,}C$. We also assume that $\mbox{epi}\,f:=\{(x,\lambda):f(x)\leq \lambda\}$ is closed in ${\mathbb{R}}^{p+1}$. The [*Legendre conjugate*]{} or [*dual function*]{} of $h$ is defined through $$h^*(z)=\sup\left\{\langle z,x\rangle -h(x):x \in C\right\},$$ for $z$ in ${\mathbb{R}}^p$, and its domain is $D:=\left\{z\in{\mathbb{R}}^p:h(z)<+\infty\right\}.$ The function $h^*$ is differentiable on the interior of $D$, and the inverse of $\nabla h:{\mathrm{int}\,}C \to {\mathrm{int}\,}D$ is $\nabla h^*:{\mathrm{int}\,}D \to {\mathrm{int}\,}C$.

We start with a simple technical lemma on the compactness of the domain of a Legendre function.

Let $h \colon {\mathbb{R}}^2 \mapsto {\mathbb{R}}$ be a globally Lipschitz continuous Legendre function, and set $D = \mathrm{int}(\mathrm{dom}(h^*))$ where $h^* \colon {\mathbb{R}}^2 \mapsto {\mathbb{R}}$ is the conjugate of $h$. For each $\lambda \geq \min_{{\mathbb{R}}^2}h$ let $\sigma_\lambda$ be the support function associated to the set $\left\{ z \in {\mathbb{R}}^2, \, h(z) \leq \lambda \right\}$. The following are equivalent

-   $h^*(x)\leq 0$ for all $x \in D$.

-   For all $y \in {\mathbb{R}}^2$, $\sigma_{h(y)}(\nabla h(y)) \leq h(y)$.

In both cases $h^*$ has compact domain. \[lem:legendreBounded\]

Let us establish beforehand the following formula $$\begin{aligned}
\label{e:h}
                                h^*(z) &= \sigma_{h(y)}(\nabla h(y))- h(y),
                \end{aligned}$$ with $y = \nabla h^* (z)$ and $y\in{\mathbb{R}}^2$. Since $h^*$ is Legendre, we have for all $y$ in $D$, $$\begin{aligned}
                                h^*(z) &= \sup_{y \in {\mathbb{R}}^2} \left\langle z, y  \right\rangle - h(y) =  \left\langle z, \nabla h^*(z) \right\rangle - h(\nabla h^*(z)).
                \end{aligned}$$ We have, setting $y = \nabla h^*(z)$ $$\begin{aligned}
                                \left\langle z, \nabla h^*(z) \right\rangle  = \left\langle \nabla h(y),y \right\rangle  = \sigma_{h(y)}(\nabla h(y))
                \end{aligned}$$ because $\nabla h(y)$ is normal to the sublevel set of $h$ which contains $y$ in its boundary. Hence we have $h^*(z) = \sigma_{h(y)}(\nabla h(y))- h(y)$ with $y = \nabla h^* (z)$, that is holds. Since $\nabla h^* \colon D \mapsto {\mathbb{R}}^2$ is a bijection, the equivalence follows. In this case the domain of $h^*$ is closed because $h^*$ is bounded and lower semicontinuous. The domain is also bounded by the Lipschitz continuity of $h$, whence compact.

\[th:globallyLipshitz\] Let $\left( S_i \right)_{i\in {\mathbb{N}}}$ be such that for any $i $ in $ I$, $T_- = S_{i}$, $T_+ = S_{i+1}$ satisfy Assumption \[ass:curvature\] and there exists a sequence $\left( \epsilon_i \right)_{i \in {\mathbb{N}}}$ in $(0,1)$ such that for all $i\geq 1$, $(1 - \epsilon_i)^{-1} S_{3i-1} = (1 +\epsilon_i)^{-1} S_{3i+1}=S_{3i}$.\
Assume in addition that, $$\begin{aligned}
                                \inf_{\|x\| = 1}& \sigma_{S_i}(x) - \sigma_{S_{i-1}}(x) =1 +  O\left(\frac{1}{i^3}\right)&\mbox{ (non degeneracy)}\label{nd}, \\
                                \sup_{\|x\|=1}& \left|\frac{\sigma_{S_{i+1}}(x) - \sigma_{S_i}(x)}{\sigma_{S_{i}}(x) - \sigma_{S_{i-1}}(x)} - 1\right| = O\left( \frac{1}{i^3} \right)& \mbox{ (moderate growth)}\label{mg}.
                \end{aligned}$$ Then there exists a convex $C^k$ function $h \colon {\mathbb{R}}^2 \mapsto {\mathbb{R}}$, such that

-   For all $i $ in $ {\mathbb{N}}$, $S_{3i}$ is a sublevel set of $h$,

-   $h$ has positive definite Hessian,

-   $h$ is globally Lipschitz continuous,

-   $h^*$ has a compact domain $D$ and is $C^k$ and strictly convex on $\mathrm{int}(D)$.

The construction of $h$ follows the exact same principle as that of Theorem \[th:smoothinterp\]. This ensures that the first two points are valid. Note that equation implies that the sets sequence grows by at least a fixed amount in each direction as $i$ grows. Hence we have ${\mathcal{T}}= {\mathbb{R}}^2$.\
*Global Lipschitz continuity of $h$:* The values of $h$ are defined through $$\begin{aligned}
                                \lambda_{i+1}-\lambda_{i}&=K_i(\lambda_{i}-\lambda_{i-1}),\;\forall i\in {\mathbb{N}}^*\\
                                K_i &= \max_{\|x\| = 1} \frac{\sigma_{S_{i+1}}(x)  - \sigma_{S_i}(x)}{ \sigma_{S_{i}}(x)  - \sigma_{S_{i-1}}(x)} \in \left( 0, + \infty \right),
                \end{aligned}$$ so that $K_i = 1 + O(1/i^3)$ thanks to equation . Note that the moderate growth assumption entails $$\begin{aligned}
                        \sup_{i \in {\mathbb{N}}} \sigma_{S_{i+1}}(x) - \sigma_{S_{i}}(x) = O\left( \prod_{i \in {\mathbb{N}}^*} K_i\right) = O(1). 
                        \label{eq:suportDiffBounded}
                \end{aligned}$$ For $i\geq 1$, one has $$\label{for}
    \lambda_{i+1}-\lambda_{i}=\prod_{1 \leq j \leq i} K_j(\lambda_1-\lambda_0).$$

On the other hand using the bounds , and the identity , there exists a constant $\kappa>0$ such that for all $i \geq 1$, all $\theta $ in $ {\mathbb{R}}/ 2 \pi {\mathbb{Z}}$, $$\begin{aligned}
                \frac{\sigma_{S_{i+1}}(n(\theta)) - \sigma_{S_{i}}(n(\theta))  }{\lambda_{i+1} - \lambda_{i}} = \frac{\left(\sigma_{S_{i+1}}(n(\theta)) - \sigma_{S_{i}}(n(\theta))\right)  }{\prod_{j=1}^i K_j (\lambda_{1} - \lambda_{0})} \geq \kappa>0.
                \label{eq:tempDerivLambda}\end{aligned}$$

By the interpolation properties described in Lemma \[lem:convexInterpol\] the function $\left\langle G(\lambda,\theta),n(\theta)\right\rangle$ constructed in Theorem \[th:smoothinterp\] has derivative with respect to $\lambda$ greater than $\kappa$. Recalling the expression of the gradient as given in Proposition \[prop:diffeoMorphism\] (and the concavity of $G$ with respect to $\lambda$), this shows that $$\|\nabla h(x)\|\leq \frac{1}{\kappa}$$ for all $x $ in $ {\mathbb{R}}^2$, and by the mean value theorem, $h$ is globally Lipschitz continuous on ${\mathbb{R}}^2$.

*Properties of the dual function:* $h$ is Legendre, its conjugate $h^*$ is therefore Legendre. From the definiteness of $\nabla^2 h$ and the fact that $\nabla h \colon {\mathbb{R}}^2 \mapsto \mathrm{int}(D)$ is a bijection, we deduce that $h^*$ is $C^k$ by the inverse mapping theorem. So the only property which we need to establish is that $h^*$ has a compact domain, in other words, using Lemma \[lem:legendreBounded\], it is sufficient to show that $\sup_{x \in \mathrm{int} D} h^*(x) \leq 0$.

Using the notation of the proof of Theorem \[th:smoothinterp\], we will show that it is possible to verify that, for all $\lambda,\theta$ in the domain of $G$ $$\begin{aligned}
\label{eq:toBeCheckedForBoundedness}
\frac{\left\langle n(\theta), G(\lambda,\theta)\right\rangle}{\left\langle \frac{\partial G}{\partial \lambda}(\lambda,\theta), n(\theta) \right\rangle} \leq \lambda.\end{aligned}$$ Equation is indeed the coordinate form of the characterization given in Lemma \[lem:legendreBounded\]. Let us observe that $$\begin{aligned}
\label{e:int}\frac{\partial}{\partial \lambda} \left( \left\langle n(\theta), G(\lambda,\theta)\right\rangle - \lambda \left\langle \frac{\partial G}{\partial \lambda}(\lambda,\theta), n(\theta) \right\rangle \right) = -\lambda \left\langle \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta), n(\theta) \right\rangle,\end{aligned}$$ and since $G$ is concave, the right hand side is positive.

Assume that we have proved that, $$\label{e:lim}
 \lambda \mapsto \sup_\theta -\lambda \left\langle \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta), n(\theta) \right\rangle$$ has finite integral as $\lambda \to \infty$. Since the function $$\label{e:funTheta}
 \theta \mapsto \lambda \left\langle \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta), n(\theta) \right\rangle$$ is continuous on ${\mathbb{R}}/2\pi{\mathbb{Z}}$ for any $\lambda$, Lebesgue dominated convergence theorem would ensure that $$\begin{aligned}
    \theta \mapsto \int_{\lambda \geq \lambda_0} -\lambda \left\langle \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta), n(\theta) \right\rangle d \lambda\end{aligned}$$ is continuous in $\theta$, so that: $$\begin{aligned}
                 \sup_{\theta} \left[\lim_{\lambda \to \infty} \left\langle n(\theta), G(\lambda,\theta)\right\rangle - \lambda \left\langle \frac{\partial G}{\partial \lambda}(\lambda,\theta), n(\theta) \right\rangle\right]<+\infty.\end{aligned}$$ Shifting values if necessary, we could assume that this upper bound is equal to zero to obtain equation . The latter being the condition required in Lemma \[lem:legendreBounded\], we would have reached a conclusion.

Let us therefore establish that is integrable over ${\mathbb{R}}_+$. Recall that $G$ is constructed using the Bernstein interpolation given in Lemma \[lem:convexInterpol\] between successive values of $\lambda$. As a result, for a fixed $\theta$, the function $\left\langle n(\theta), G(\lambda,\theta)\right\rangle $ is the interpolation of the piecewise affine function interpolating $$\begin{aligned}
                &\left( \lambda_{3i}, \sigma_{S_{3i}}(n(\theta))\right)\nonumber\\
                &\left( \lambda_{3i+1}, \sigma_{S_{3i+1}}(n(\theta))\right),\nonumber \\
                &\left( \lambda_{3i+2}, \sigma_{S_{3i+2}}(n(\theta))\right),\nonumber \\
                                &(\lambda_{3i+3}, \sigma_{3i+3}(n(\theta))),
                                \label{eq:piecewiseAffineInterpolExplicit2}\end{aligned}$$ as in equation . This interpolation is concave and increasing.

Assumption ensures that $K_j = 1 + O(1/j^3)$. Then $$\begin{aligned}
                \prod_{j=1}^m K_j = \bar{K} + O(1 / j^{2})\end{aligned}$$ where $\bar{K}$ is the finite, positive limit of the product (we can for example perform integral series comparison after taking the logarithm).

The recursion on the values writes for all $i\geq1$ $$\begin{aligned}
                \lambda_{i+1} = \lambda_i + K_i(\lambda_i - \lambda_{i-1}),\end{aligned}$$ so that $$\begin{aligned}
                \lambda_{i+1} - \lambda_i= (\lambda_1 - \lambda_0) \prod_{j=1}^i K_i = (\lambda_1 - \lambda_0) \bar{K} + O(1/i^{2}).\end{aligned}$$

This means that the gap between consecutive values tends to be constant. Thus by in Lemma \[lem:convexInterpol\], see also Remark \[rem:alignedInterpolation\], the degree of the Bernstein interpolants is bounded. Using this bound together with inequality , providing bounds for the derivatives of Bernstein’s polynomial, ensure that, for all $\lambda $ in $ [\lambda_{3i}, \lambda_{3i+3})$: $$\begin{aligned}
                &\left|\left\langle \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta), n(\theta) \right\rangle\right| \\
                =\,& O\left(  \max _{j = 3i+2, 3i+1}\left|  \frac{\sigma_{S_{j+1}}(n(\theta))   - \sigma_{S_{j}}(n(\theta))}{\lambda_{j+1} - \lambda_{j}} - \frac{\sigma_{S_{j}}(n(\theta)) - \sigma_{S_{j-1}}(n(\theta))  }{\lambda_{j} - \lambda_{j-1}}  \right| \right).\end{aligned}$$ Now for any $j = 3i+2, 3i+1$, $$\begin{aligned}
    &\left|  \frac{\sigma_{S_{j+1}}(n(\theta))   - \sigma_{S_{j}}(n(\theta))}{\lambda_{j+1} - \lambda_{j}} -    \frac{\sigma_{S_{j}}(n(\theta)) - \sigma_{S_{j-1}}(n(\theta))  }{\lambda_{j} - \lambda_{j-1}}  \right| \\
    =\,& \frac{1}{\lambda_{j+1} - \lambda_{j}} \left|  (\sigma_{S_{j+1}}(n(\theta))   - \sigma_{S_{j}}(n(\theta))) -    \frac{\lambda_{j+1} - \lambda_{j}}{\lambda_{j} - \lambda_{j-1}} (\sigma_{S_{j}}(n(\theta)) - \sigma_{S_{j-1}(n(\theta))  }) \right|\\
    =\,& \left(1 / ((\lambda_1 - \lambda_0) \bar{K}) + O(1/i^2)\right)\\
    &\times \left|  (\sigma_{S_{j+1}}(n(\theta))   - \sigma_{S_{j}}(n(\theta))) -   (1 + O(1/i^2)) (\sigma_{S_{j}}(n(\theta)) - \sigma_{S_{j-1}(n(\theta))  }) \right|\\
     =\,& \left(1 / ((\lambda_1 - \lambda_0) \bar{K}) + O(1/i^2)\right)\\
    &\times \left|  (\sigma_{S_{j+1}}(n(\theta))   - \sigma_{S_{j}}(n(\theta))) -    (\sigma_{S_{j}}(n(\theta)) - \sigma_{S_{j-1}(n(\theta))  }) \right|\end{aligned}$$ where the last identity follows from the triangle inequality because using $\sigma_{S_{j}}(n(\theta)) - \sigma_{S_{j-1}}(n(\theta)) = O(1)$ in . Hence $$\begin{aligned}
                &\left|\left\langle \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta), n(\theta) \right\rangle\right| \\
                =\,& 
                \left(1 / ((\lambda_1 - \lambda_0) \bar{K}) + O(1/i^2)\right)\\
                &\times O\left( \max _{j = 3i+2, 3i+1}\left|  \sigma_{S_{j+1}}(n(\theta))   - \sigma_{S_{j}}(n(\theta)) -   (\sigma_{S_{j}}(n(\theta)) - \sigma_{S_{j-1}}(n(\theta)) ) \right| \right)\\
                =\,& 
                \left(1 / ((\lambda_1 - \lambda_0) \bar{K}) + O(1/i^2)\right)\\
                &\times O\left(  \max _{j = 3i+2, 3i+1}\left|\sigma_{S_{j}}(n(\theta)) - \sigma_{S_{j-1}}(n(\theta))  \right| \times \left|  \frac{\sigma_{S_{j+1}}(n(\theta))   - \sigma_{S_{j}}(n(\theta))}{\sigma_{S_{j}}(n(\theta)) - \sigma_{S_{j-1}}(n(\theta)) } -   1 \right| \right)\\
                =\,& O(1/i^3),\end{aligned}$$ where the last inequality follows from and . Now as $i \to \infty$, $\lambda_{3i} \sim \lambda_{3i + 3} \sim i c$ for some constant $c>0$ and $$\begin{aligned}
                \sup_{\lambda \in [\lambda_{3i}, \lambda_{3i+3}], \theta \in [0, 2 \pi]} - \lambda \left\langle \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta), n(\theta) \right\rangle = O(1 / i^2)\end{aligned}$$ and $$\begin{aligned}
                \sup_{ \theta \in [0, 2 \pi]} - \lambda \left\langle \frac{\partial^2 G}{\partial \lambda^2}(\lambda,\theta), n(\theta) \right\rangle\end{aligned}$$ has finite integral as $\lambda \to \infty$. This implies and it concludes the proofs.

Smooth convex interpolation for sequences of polygons
=====================================================

Given a sequence of points $A_1,\ldots,A_n$, we denote by $A_1\ldots A_n$ the polygon obtained by joining successive points ending the loop with the segment $[A_n,A_1]$. In the sequel we consider mainly convex polygons, so that the vertices $A_1,\ldots,A_n$ are also the extreme points.

The purpose of this section is first to show that polygons can be approximated by smooth convex sets with prescribed normals under weak assumptions. Figure \[fig:illustrPolySmooth\] illustrates the result we would like to establish: given a target polygon with prescribed normals at its vertices, we wish to construct a smooth convex set interpolating the vertices with the desired normals and whose distance to the polygon is small.

Then given a sequence of nested polygons, we provide a smooth convex function which interpolates the polygons in the sense described just above.

Given a closed nonempty convex subset $S$ of ${\mathbb{R}}^p$ and $x$ in $S$, we recall that the [*normal cone to $S$ at $x$*]{} is $$N_S(x)=\left\{z\in {\mathbb{R}}^p:\langle z,y-x\rangle \leq 0,\forall y\in S\right\}.$$ Such vectors will often simply called normals (to $S$) at $x$.

Smooth approximations of polygons
---------------------------------

\[lem:approxSegment\] For any $r_-,r_+ > 0$, $t_- > 0$, $t_+<0$ and $\epsilon > 0$, $m \in {\mathbb{N}}$, $m \geq 3$, there exists a strictly concave polynomial function $p \colon [0,1] \mapsto [0,\epsilon]$ such that $$\begin{aligned}
                                p(0) &= 0 &p(1) &= 0\\
                                p'(0) &= t_-& p'(1) &= t_+\\
                                p''(0) &= -r_- &p''(1) &= -r_+.\\
                                p^{(q)}(0) & = 0 \quad  q \in\{3, \ldots, m\}.
                \end{aligned}$$

![Arrows designate the prescribed normals. We construct a strictly convex set with smooth boundary entirely contained in the auxiliary (blue) polygon. This set interpolates the normals and the distance to the original (red) polygon can be chosen arbitrarily small. The degree of smoothness of the boundary can be chosen arbitrarily high.[]{data-label="fig:illustrPolySmooth"}](illustrPolySmooth){width=".7\textwidth"}

Let us begin with preliminary remarks. Consider for any $a,b $ in $ {\mathbb{R}}$, the function $$\begin{aligned}
                                f \colon t \mapsto a(t-b)^2.
                \end{aligned}$$ For any $t $ in $ {\mathbb{R}}$, $q $ in $ {\mathbb{N}}$, $q > 2$, and any $c > 0$, we have $$\begin{aligned}
                                f(t + c) - f(t) &= \Delta^1_c f(t) = a c ( 2(t-b) + c)\nonumber\\
                                f(t + 2 c) - 2 f(t+c) + f(t) &= \Delta^2_c f(t) = a c ( 2 (t + c - b) + c - 2(t-b) - c) = 2ac^2\nonumber\\
                                \Delta_c^q f(t) &= 0
                                \label{eq:approxSegment00}
                \end{aligned}$$

[*Choosing the degree $d$ and constructing the polynomial.* ]{} For any $d $ in $ {\mathbb{N}}$, $d \geq 2m + 1$, we set $$\begin{aligned}
                                a_-(d) &= \frac{-dr_-}{2(d - 1)}<0&b_-(d) &= \frac{1}{2d} \left( 1 + 2 t_-\frac{d-1}{r_-} \right)> 0\nonumber\\
                                a_+(d) &= \frac{-dr_+}{2(d - 1)}<0&b_+(d) &= 1 + \frac{1}{2d} \left( -1 + 2t_+\frac{d-1}{r_+} \right)< 1,
                                \label{eq:approxSegment0}
                \end{aligned}$$ and define the functions $$\begin{aligned}
                                f_d \colon s &\mapsto
                                \begin{cases}
                                                a_-(d) (( s - b_-(d))^2 - b_-(d)^2) & \text{ if } s\leq b_-(d)\\
                                                - a_-(d) b_-(d)^2& \text{ if } s\geq b_-(d)
                                \end{cases}\\
                                g_d \colon t &\mapsto
                                \begin{cases}
                                                a_+(d) (( s - b_+(d))^2 - (1-b_+(d))^2) & \text{ if } s\geq b_+(d)\\
                                                - a_+(d) (1-b_+(d))^2& \text{ if } s\leq b_+(d).
                                \end{cases}
                \end{aligned}$$ Furthermore, we set $$\begin{aligned}
                                f \colon t &\mapsto
                                \begin{cases}
                                                \frac{r_-}{2} \left( \left( \frac{t_-}{r_-}  \right)^2- \left(s - \frac{t_-}{r_-}  \right)^2 \right)  & \text{ if } s\leq \frac{t_-}{r_-}\\
                                                \frac{r_-}{2} \left( \frac{t_-}{r_-}  \right)^2 & \text{ if } s\geq \frac{t_-}{r_-}
                                \end{cases}\\
                                g \colon t &\mapsto
                                \begin{cases}
                                                \frac{r_+}{2} \left(\left( \frac{t_+}{r_+}  \right)^2 -  \left(s -1 - \frac{t_+}{r_+}  \right)^2 \right)  & \text{ if } s\geq 1 + \frac{t_+}{r_+}\\
                                                \frac{r_+}{2} \left( \frac{t_+}{r_+}  \right)^2 & \text{ if } s\leq 1 + \frac{t_+}{r_+}.
                                \end{cases}
                \end{aligned}$$ Note that $b_-(d) \to t_- / r_-$, $b_+(d) \to 1 + t_+ / r_+$, $a_-(d) \to - r_-/2$ and $a_+(d) \to - r_+/2$ as $d \to \infty$ so that $f_d \to f$ and $g_d \to g$ uniformly on $[0,1]$. For any $d$, $f_d$ is concave increasing and $g_d$ is concave decreasing and all of them are Lipschitz continuous on $[0,1]$ with constants that do not depend on $d$. Note also that $f(0) = 0 < g(0)$ and $g(1) = 0 < f(1)$. We choose $d \geq 2m + 1$ such that $$\begin{aligned}
                                f_d\left( \frac{m}{d} \right) &\leq \min\left( \epsilon, g_d\left( \frac{m}{d} \right) \right)\nonumber\\
                                g_d\left(1 -  \frac{m}{d} \right) &\leq \min\left( \epsilon, f_d\left(1 -  \frac{m}{d} \right) \right).
                                \label{eq:approxSegment1}
                \end{aligned}$$ Such a $d$ always exists because in both cases, the left hand side converges to $0$ and the right hand side converges to a strictly positive term as $d$ tends to $\infty$. For such a $d$, we set $h \colon s \mapsto \min\left\{ f_d(s), g_d(s),\epsilon \right\}$. By construction, $h$ is concave, agrees with $f_d$ on $\left[ 0,\frac{m}{d} \right] \subset [0,1/2]$ and with $g_d$ on $\left[ 1- \frac{m}{d},1 \right] \subset [1/2,1]$. Using equation with $c = 1/d$, we deduce that $$\begin{aligned}
                                d (d-1) \Delta^2 h(0) &= d (d-1) \Delta^2 f_d(0) = d(d-1)\frac{2 a_-(d)}{d^2} = -r_-\\
                                d \Delta h(0) &= d \Delta f_d(0) =  a_-(d)\left( \frac{1}{d}  -2b_-(d)\right) = t_-\\
                                \Delta^{q}h(0) &= 0 = \Delta^{q}f_d(0) \quad \forall m \geq q \geq 3\\
                                d (d-1) \Delta^2 h\left(1 - \frac{2}{d}\right) &= d (d-1) \Delta^2 g_d\left(1 - \frac{2}{d}\right)= d(d-1)\frac{2 a_+(d)}{d^2} = -r_+\\
                                d \Delta h\left( 1 - \frac{1}{d} \right) &= d \Delta g_d\left( 1 - \frac{1}{d} \right) =  a_+(d) \left( \frac{-1}{d}  + 2 (1 - b_+(d))\right) = t_+\\
                                \Delta^{q}h\left( 1 - \frac{m}{d} \right) &= \Delta^{q}g_d\left( 1 - \frac{m}{d} \right) = 0 \quad \forall m \geq q \geq 3.
                \end{aligned}$$ From the concavity of $h$ and the derivative formula , we deduce that the polynomial $B_{h,d}$ satisfies the desired properties.

![Illustration of the approximation result of Lemma \[lem:approxSegment\] with $\epsilon = 0.1$, $m = 3$, $t_- =0.7$, $t_+ = -2.2$, $r_- = 2$ and $r_+ = 0.2$. The resulting polynomial is of degree 66. Numerical estimations of the first and second order derivatives at 0 and 1 match the required values up to 3 precision digits. The polynomial is strongly concave, however, this is barely visible because the strong concavity constant is extremely small. []{data-label="fig:illustrApproxSeg"}](illustrApproxSeg){width="\textwidth"}

We deduce the following result

Let $a > 0$, $r > 0$, $\epsilon > 0$, and an integer $m \geq 3$. Consider two unit vectors: $v_-$ with strictly positive entries and $v_+$ with first entry strictly positive and second entry strictly negative. Then there exists a $C^m$ curve $\gamma \colon [0,M] \mapsto {\mathbb{R}}^2$, such that

1.  $\|\gamma'\| = 1$.

2.  $\gamma(0) = (-a, 0) := A$ and $\gamma(1) = (0,a) : = B$.

3.  $\gamma'(0) = v_-$ and $\gamma'(1) = v_+$.

4.  $\|\gamma''(0)\| = \|\gamma''(-1)\| = r$.

5.  $\mathrm{det}(\gamma',\gamma'') < 0$ along the curve.

6.  $\gamma^{(q)}(0) = \gamma^{(q)}(1) = 0$ for any $3 \leq q \leq m$.

7.  ${\mathrm{dist}}(\gamma([0,M]), [A,B]) \leq \epsilon$.

\[lem:existenceCurve\]

Consider the graph of a polynomial as given in Lemma \[lem:approxSegment\] with $t_- = v_-[2] / v_-[1]> 0$, $t_+ = v_+[2] / v_+[1]< 0$ and $r_- = \frac{r}{2a} (1 + t_-^2)^{\frac{3}{2}}$, $r_+ = \frac{r}{2a} (1 + t_+^2)^{\frac{3}{2}}$ and $\epsilon/2a$ as an approximation parameter. This graph is parametrized by $t$. It is possible to reparametrize it by arclength to obtain a $C^m$ curve $\gamma_0$ whose tangents at $0$ is $T_-$, at $1$ is $T_+$, and whose curvature at $0$ and $1$ is $- \frac{r}{2a}$. Furthermore, $\gamma_0$ has strictly negative curvature whence item 5. Consider the affine transform: $x \mapsto 2a(x - 1/2)$, $y \mapsto 2a y$. This results in a $C^m$ curve $\gamma$, parametrized by arclength which satisfies the desired assumptions.

Let $S=A_1...A_n$ be a convex polygon. For each $i$, let $V_i$ be in $ N_S(A_i)$ such that the angle between $V_i$ and each of the two neighboring faces is within $\left( \frac{\pi}{2},\pi \right)$. Then for any $\epsilon > 0$ and any $m \geq 2$, there exists a compact convex set $C \subset {\mathbb{R}}^2$ such that

-   the boundary of $C$ is $C^m$ with non vanishing curvature,

-   $S \subset C$,

-   for any $i = 1,\ldots, n$, $A_i $ in $ {\mathrm{bd}\,}(C)$ and the normal cone to $C$ at $A_i$ is given by $V_i$,

-   $\max_{y \in C} {\mathrm{dist}}(y,S) \leq \epsilon$.

\[lem:polyGon\]

We assume without loss of generality that $A_1,\ldots,A_n$ are ordered clockwise. For $i = 1, \ldots, n-1$, and each segment $[A_i,A_{i+1}]$, we may perform a rotation and a translation to obtain $A_i = -(a,0)$ and $A_{i+1} = (a,0)$. Working in this coordinate system, using the angle condition on $V_i$, we may choose $v^-_i$, $v^+_{i+1}$ satisfying the hypotheses of Lemma \[lem:existenceCurve\] respectively orthogonal to $V_i$ and $V_{i+1}$. Choosing $r = 1$, we obtain $\gamma_i \colon [0,M_i] \mapsto {\mathbb{R}}^2$ as given by Lemma \[lem:existenceCurve\]. Rotation and translations affect only the direction of the derivatives of curves, not their length. Hence, it is possible to concatenate curves $\left( \gamma_i \right)_{i=1}^{n-1}$ and to preserve the $C^m$ properties of the resulting curve. At end-points, tangents and second order derivatives coincide while higher derivatives vanish. Furthermore the curvature has constant sign and does not vanish. We obtain a closed $C^m$ curve which defines a convex set which satisfies all the requirements of the lemma.

[Given any polygon, choosing normal vectors as given by the direction of the bissector of each angles ensure that the above assumptions are satisfied. Hence all our approximation results hold given polygon without specifying the choice of outer normals.]{} \[rem:noNormal\]

Smooth convex interpolants of polygonal sequences
-------------------------------------------------

[For $n\geq3$, let $A_1\ldots A_n $ be a convex polygon $S$ and $V_i $ be in $ N_S(V_i)$ for $i=1,\ldots,n$. We say that $\left( A_i,V_i \right)_{i=1}^n $ is [*interpolable*]{} if for each $i = 1,\ldots,n$, the angle between $V_i$ and each of the two neighboring faces of the polygon is in $\left( \frac{\pi}{2},\pi \right)$. The collection $\left( A_i,V_i \right)_{i=1}^n $ is called [*a polygon-normal pair*]{}.]{} \[def:interpolable\]

Let $I = {\mathbb{Z}}$ or $I = {\mathbb{N}}$. Let $\left( PN_i \right)_{i \in I}$ be a sequence of interpolable polygon-normal pairs. Setting for $i $ in $ I$, $PN_i = \left\{ \left( A_{j,i} \right)_{j=1}^{n_i}, \left( V_{j,i} \right)_{j=1}^{n_i} \right\}$ where $n_j $ is in $ {\mathbb{N}}$ and denoting by $T_i$ the polygon $A_{1,i}\ldots A_{n_i,i}$, we say that the sequence $\left( PN_i \right)_{i \in I}$ is strictly increasing if for all $i $ in $ I$, $T_i \subset \mathrm{int}(T_{i+1})$.

Let $\left( PN_i \right)_{i \in I}$ be a strictly increasing sequence of interpolable polygon-normal pairs. A sequence $(\epsilon_i)_{i \in I}$ in $(0,1)$ is said to be [*admissible*]{} if $0 \in \mathrm{int}(T_i)$ for each $i $ in $ I$ and $$\gamma T_i\subset{\mathrm{int}\,}T_{i+1}$$ for all $\gamma \in [1-\epsilon_i,1+\epsilon_i].$ We have the following corollary of Theorem \[th:smoothinterp\].

Let $I = {\mathbb{Z}}$ or $I = {\mathbb{N}}$. Let $\left( PN_i \right)_{i \in I}$ be a strictly increasing sequence of interpolable polygon-normal pairs and $(\epsilon_i)_{i \in I}$ be admissible. Set $\displaystyle\mathcal{T}:=\mathrm{int}\left(\cup_{i\in I}T_i  \right).$

Then for any $k $ in $ {\mathbb{N}}$, $k \geq 2$ there exists a $C^k$ convex function $f \colon \mathcal{T} \mapsto {\mathbb{R}}$, and an increasing sequence $(\lambda_i)_{i \in I}$, with $\inf_{i \in I} \lambda_i > -\infty$, such that for each $i $ in $ I$

-   $T_i \subset \left\{ x,\,f(x) \leq \lambda_i \right\}$.

-   ${\mathrm{dist}}(T_i, \left\{ x,\,f(x) \leq \lambda_i \right\}) \leq \epsilon_i$.

-   For each $i $ in $ I$, $j$ in $\{1,\ldots,n_i\}$, we have $f(A_{i,j}) = \lambda_i$ and $\nabla f(x)$ is colinear to $V_{i,j}$.

-   $\nabla^2 f$ is positive definite outside $\operatorname*{argmin}f$. When there is a unique minimizer then $\nabla^2f$ is positive definite throughout ${\mathcal{T}}$ [(]{}this is the case when $I = {\mathbb{N}}$ or when $I={\mathbb{Z}}$ and $\cap_{i \in I} T_i$ is a singleton[)]{}.

\[cor:polygonDecrease\]

We add two remarks which will be useful for directional convergence issues and the construction of Legendre functions:

-   If two consecutive elements of the sequence of interpolable polygon-normal pairs are homothetic with center $0$ in the interior of both polytopes, then the restriction of the resulting convex function to this convex ring can be constructed such that all the sublevel sets within this ring are homothetic with the same center.

-   If further conditions are imposed on the elements of a strictly increasing interpolable polygon-normal pair, then the resulting function can be constructed to be Legendre and globally Lipschitz continuous (that is, its Legendre conjugate has bounded support). This is a consequence of Proposition \[th:globallyLipshitz\] and will be detailled in the next section.

More on Legendre functions and a pathological function with polyhedral domain
-----------------------------------------------------------------------------

Using intensively polygonal interpolation, we build below a finite continuous Legendre function $h$ on an $\ell^\infty$ square with oscillating “mirror lines": $t\to \nabla h^*(\nabla h(x_0)+tc)$.

We start with the following preparation proposition related to the Legendre interpolation of Proposition \[th:globallyLipshitz\].

Let $\left( PN_i \right)_{i \in {\mathbb{N}}^*}$ be a strictly increasing sequence of interpolable polygon-normal pairs. Setting for $i$ in ${\mathbb{N}}^*$, $PN_i = \left\{ \left( A_{j,i} \right)_{j=1}^{n_i}, \left( V_{j,i} \right)_{j=1}^{n_i} \right\}$ where $n_j$ is in ${\mathbb{N}}^*$ and denoting by $T_i$ the polygon $A_{1,i}\ldots A_{n_i,i}$, we assume that $$T_i=3i P,\,\forall i \in {\mathbb{N}}^*,$$ where $P$ is a fixed polygon which contains the unit Euclidean disk.\
Then for any $l$ in ${\mathbb{N}}$, $l \geq 2$, there exists a strictly increasing sequence of sets $\left( S_i \right)_{i \in {\mathbb{N}}, \,i\geq 2}$, such that for $j \geq 1$,

-   $S_{3j}$ interpolates the normals of $PN_j$ in the sense of Lemma \[lem:polyGon\] with $\mathrm{dist}(S_{3j}, T_j) \leq 1 / (4(3j+2)^l)$

-   $S_{3j - 1} = \frac{3j - 1}{3_j} S_{3j}$

-   $S_{3j + 1} = \frac{3j + 1}{3_j} S_{3j}$

This sequence has the following properties

-   there exists $c>0$, such that for all $j $ in ${\mathbb{N}}$, $j \geq 3$ and for all unit vector $x$, $$\begin{aligned}
                        c\geq\sigma_{S_{j+1}}(x) - \sigma_{S_j}(x) \geq 1 - \frac{1}{(j+1)^l}.
                        \label{eq:interpolGauge1}
                    \end{aligned}$$

-   for all unit vector $x$, $$\begin{aligned}
                                    \left| \frac{\sigma_{S_{j+1}}(x)  - \sigma_{S_j}(x)}{ \sigma_{S_{j}}(x)  - \sigma_{S_{j-1}}(x)} -1\right| \leq \frac{1}{j^l},\:\forall j\geq 3.
                                    \label{eq:interpolGauge2}
                    \end{aligned}$$

\[lem:interpolGauge\]

Set for all $j$ in ${\mathbb{N}}^*$, $\delta_j = \frac{1}{4(3j+2)^l}$ and let $S_{3j}$ be the $\delta_j$ interpolant of $T_j=3jP$ as given by Lemma \[lem:polyGon\] so that ${\mathrm{dist}}(S_{3j}, 3jP) \leq \delta_j$. Since $P$ contains the unit ball, $$\label{inc} 3j P \subset S_{3j} \subset (3j + \delta_j) P.$$ Now set $$\begin{aligned}
                                S_{3j - 1} &= \frac{3j - 1}{3_j} S_{3j}\\
                                S_{3j + 1} &= \frac{3j + 1}{3_j} S_{3j}.
                \end{aligned}$$ For any $j$ in ${\mathbb{N}}^*$, it is clear that $S_{3j-1} \subset \mathrm{int}(S_{3j})$ and $S_{3j} \subset \mathrm{int}(S_{3j+1})$. Furthermore, by , we have $$\begin{aligned}
                                S_{3j+1} \subset \frac{3j+1}{3j} ( 3j + \delta_j) P \subset ((3j+1) + 2\delta_j) P \subset \mathrm{int}((3j+2) P) \subset \mathrm{int} (S_{3j+2})
                \end{aligned}$$ so that we indeed have a strictly increasing sequence of sets. We obtain from the construction, for any $j$ in ${\mathbb{N}}^*$, and any unit vector $x$, $$\begin{aligned}
                                \sigma_{S_{3j}}(x) - \sigma_{S_{3j-1}}(x) &=\sigma_{S_{3j+1}}(x) - \sigma_{S_{3j}}(x)\notag \\
                                &= \frac{1}{3j} \sigma_{S_{3j}}(x) \in \left[ \sigma_P(x), \left( 1 + \frac{\delta_j}{3j} \right) \sigma_P(x) \right] \subset \left[ 1, \left(1+\frac{\delta_j}     {3j} \right) \sigma_P(x) \right]\notag \\
                                \sigma_{S_{3j+2}}(x) - \sigma_{S_{3j+1}}(x) &\leq \sigma_P(x) (3j +2) \left( 1 + \frac{\delta_{j+1}}{3j+3} \right) - \sigma_P(x)(3j+1)\notag\\
                                \sigma_{S_{3j+2}}(x) - \sigma_{S_{3j+1}}(x) &\geq \sigma_P(x) (3j +2) - \sigma_P(x)(3j+1)\left( 1 + \frac{\delta_j}{3j} \right) \notag\\
                                & = \sigma_P(x) \left( 1 - \delta_j \frac{3j+1}{3j} \right)\notag\\ 
                                &\geq 1 -  \delta_j \frac{4}{3}  \geq 1 - \frac{1}{(3j+2)^l}\label{intermed}
                \end{aligned}$$ This proves . We deduce that for all $j$ in ${\mathbb{N}}^*,$ $$\begin{aligned}
                                \frac{\sigma_{S_{3j+1}}(x)  - \sigma_{S_{3j}}(x)}{ \sigma_{S_{3j}}(x)  - \sigma_{S_{3j-1}}(x)}& = 1, \mbox{ for all nonzero vector $x$,} \\
                                \max_{\|x\| = 1} \frac{\sigma_{S_{3j+2}}(x)  - \sigma_{S_{3j+1}}(x)}{ \sigma_{S_{3j+1}}(x)  - \sigma_{S_{3j}}(x)} &\leq (3j +2) \left( 1 + \frac{\delta_{j+1}}{3j+3} \right) - (3j+1)\\
                                &= 1 + \frac{3j+2}{3j+3} \delta_{j+1} \leq 1 + \frac{1}{(3j +1)^l},\\
                                \min_{\|x\| = 1} \frac{\sigma_{S_{3j+2}}(x)  - \sigma_{S_{3j+1}}(x)}{ \sigma_{S_{3j+1}}(x)  - \sigma_{S_{3j}}(x)} &\underset{\eqref{intermed}}{\geq} \frac{1 - \delta_j \frac{3j+1}{3j} }{\left( 1 + \frac{\delta_j}{3j} \right) }\\
                                &\geq \left( 1 - \delta_j \frac{3j+1}{3j}  \right) \left( 1 - \frac{\delta_j}{3j} \right) \\ 
                                & \geq 1 - \delta_j\left( \frac{3j+1}{3j} + \frac{1}{3j} \right) \geq 1 - \delta_j \frac{5}{3} \geq 1 - \frac{1}{(3j+1)^l}.
                \end{aligned}$$ Furthermore, using the fact that $t \mapsto \frac{1+t}{1-t}$ is increasing on $(-\infty,1)$ and the fact that $\delta_{j+1} \leq \delta_j$, $$\begin{aligned}
                                \max_{\|x\| = 1} \frac{\sigma_{S_{3j+3}}(x)  - \sigma_{S_{3j+2}}(x)}{ \sigma_{S_{3j+2}}(x)  - \sigma_{S_{3j+1}}(x)} &\underset{\eqref{intermed}}{\leq} \frac{1 + \frac{\delta_{j+1}}{3j + 3}}{1 - (3j+1) \frac{\delta_j}{3j}} \notag\\
                                &\leq\frac{1 + (3j+1)\frac{\delta_{j}}{3j}}{1 - (3j+1) \frac{\delta_j}{3j}} \notag\\
                                &\leq\frac{1 + \delta_j\frac{4}{3}}{1 - \delta_j \frac{4}{3}}\label{borntobewild}.
                \end{aligned}$$ Setting $s(t)= (1+t)/(1-t)$, we have, for all $t \leq 1/2$ $$\begin{aligned}
s'(t) &= \frac{2}{(1-t)^2}, \, s(0)=1\\
                    s''(t) &= \frac{4}{(1-t)^3} \leq 24,\, s'(0)=2.
                \end{aligned}$$ Thus $s(t) \leq 1 + 2 t + 12 t^2$ on $(-\infty,1/2]$. Since $\frac{4}{3} \delta_j \leq \frac{4}{75} \leq \frac{1}{2}$, we deduce from the previous remark and above: $$\begin{aligned}
                                    \max_{\|x\| = 1} \frac{\sigma_{S_{3j+3}}(x)  - \sigma_{S_{3j+2}}(x)}{ \sigma_{S_{3j+2}}(x)  - \sigma_{S_{3j+1}}(x)} &\leq 1 + \frac{8}{3} \delta_j + \frac{64}{3} \delta_j^2 = 1 + \delta_j\left( \frac{8}{3} + \frac{64}{3} \delta_j \right)\\
                                &\leq 1 + \delta_j\left(3 + \frac{64}{3 \times 25}\right) \leq 1 + 4 \delta_j = 1 + \frac{1}{(3j+2)^l}.
                \end{aligned}$$ Finally using again, $$\begin{aligned}
                                \min_{\|x\| = 1} \frac{\sigma_{S_{3j+3}}(x)  - \sigma_{S_{3j+2}}(x)}{ \sigma_{S_{3j+2}}(x)  - \sigma_{S_{3j+1}}(x)} &\geq \frac{1 }{(3j +2) \left( 1 + \frac{\delta_{j+1}}{3j+3} \right) - (3j+1) } \\
                                &= \frac{1 }{1 + \delta_{j+1}\frac{3j+2}{3j+3}  } \\
                                &\geq 1 - \delta_{j+1}\frac{3j+2}{3j+3} \geq 1 - \delta_{j+1} \geq 1 - \frac{1}{(3j+2)^l}. \end{aligned}$$ This proves the desired result.

Combining Lemma \[lem:interpolGauge\] and Proposition \[th:globallyLipshitz\], we obtain the following result.

Let $\left( PN_i \right)_{i \in {\mathbb{N}}^*}$ be a strictly increasing sequence of interpolable polygon-normal pairs. Set for $i$ in ${\mathbb{N}}^*$, $PN_i = \left\{ \left( A_{j,i} \right)_{j=1}^{n_i}, \left( V_{j,i} \right)_{j=1}^{n_i} \right\}$ where $n_i$ is in ${\mathbb{N}}^*$, denote by $T_i$ the polygon $A_{1,i}\ldots A_{n_i,i}$, and assume that for all $i $ in $ {\mathbb{N}}^*$, $T_i=3i P$ where $P$ is a fixed polygon which contains the unit disk in its interior.\
Then for any $k $ in $ {\mathbb{N}}$, $k \geq 2$ and all $l \geq 3$, there exists a $C^k$ globally Lipschitz continuous Legendre function, $h \colon {\mathbb{R}}^2 \mapsto {\mathbb{R}}$, and an increasing sequence $(\lambda_i)_{i \in  {\mathbb{N}}}$, such that for each $i $ in $ {\mathbb{N}}$:

-   $T_i \subset \left\{ x,\,h(x) \leq \lambda_i \right\}$,

-   ${\mathrm{dist}}(T_i, \left\{ x,\,h(x) \leq \lambda_i \right\}) \leq \frac{1}{4(3i+2)^l}$,

-   For any $x$ with $h(x) = \lambda_i$ and $\nabla h(x)$ is colinear to $V_i$ for each vertex $x$ of $T_i$,

-   $h$ has positive definite Hessian and is globally Lipschitz continuous,

-   $h^*$ has compact domain and is $C^k$ on the interior of its domain.

\[th:Legendre\]

The function $h^*$ constructed in Theorem \[th:Legendre\] has compact polygonal domain and is continuous on this domain. \[cor:legendreContinuity\]

Since $P$ is a polygon and contains the unit Euclidean disk, the gauge function of $3P$ is polyhedral with full domain ${\mathbb{R}}^2$, call it $\omega$. Denote by $P^{\circ}$ the polar of $P$. This is a polytope and since $\omega$ is the gauge of $P$, we actually have $\omega = \sigma_{P^{\circ}}$, the support function of the polar of $P$ [@schneider1993convex Theorem 1.7.6]. Hence the the convex conjugate of $\omega$ is the indicator of the polytope $P^\circ$ [@rockafellar1970convex Theorem 13.2].

It can be easily seen from the proof of Proposition \[th:globallyLipshitz\] that $\lambda_i = \alpha i + r_i$ with $r(i) = O(1)$ as $i \to \infty$. Without loss of generality, we may suppose that $\alpha = 1$ (this is a simple rescaling) so that there is a positive constant $c$ such that $|\lambda_i - i| \leq c$ for all $i$.

Let $h$ be given as in Theorem \[th:Legendre\], fix $i \geq 1$ and $x \in {\mathbb{R}}^2$ such that $\lambda_{i-1} \leq h(x) \leq \lambda_i$. We have in ${\mathbb{R}}^2$ $$\begin{aligned}
       \{y:\, \omega(y) \leq i-1 \} \subset \{y:\, h(y) \leq  \lambda_{i-1} \} \subset  \{y:\, h(y) \leq \lambda_{i} \} \subset \{y:\, \omega(y) \leq i+1 \} 
    \end{aligned}$$ and hence $$\begin{aligned}
        i - 1 \leq \omega(x) \leq i+1
    \end{aligned}$$ and we deduce that $$\begin{aligned}
        \omega(x) - 2 - c \leq i - 1 -c \leq  \lambda_{i-1} \leq h(x) \leq \lambda_i \leq i+c \leq \omega(x) + c + 1.
    \end{aligned}$$ Since $i$ was arbitrary, this shows that there exists a constant $C>0$ such that $|h(x) - \omega(x)|\leq C$ for all $x \in {\mathbb{R}}^2$. Recall that $z \mapsto \sup_{y \in {\mathbb{R}}^2} \left\langle y,z\right\rangle - \omega(y)$ is the indicator function of $P^\circ$, hence, $$\begin{aligned}
        z \in P^{\circ} & \Rightarrow \sup_{y \in {\mathbb{R}}^2} \left\langle y,z\right\rangle - \omega(y) = 0 \Rightarrow \sup_{y \in {\mathbb{R}}^2} \left\langle y,z\right\rangle - h(y) \leq C < +\infty \\
        z \not\in P^{\circ} & \Rightarrow \sup_{y \in {\mathbb{R}}^2} \left\langle y,z\right\rangle - \omega(y) = +\infty \Rightarrow \sup_{y \in {\mathbb{R}}^2} \left\langle y,z\right\rangle - h(y) = +\infty
    \end{aligned}$$ which shows that the domain of $h^*$ is actually $P^{\circ}$ which is a polytope. Now, $h^*$ is convex and lower semicontinuous on $P^{\circ}$, invoking the results of [@gale68convex], it is also upper semicontinuous on $B^*$ and finally it is continuous on $B^*$.

For any $\theta \in \left( \frac{-\pi}{4},\frac{\pi}{4} \right)$ there exists a Legendre function $h\colon {\mathbb{R}}^2 \mapsto {\mathbb{R}}$ whose domain is a closed square, continuous on this domain and $C^k$ on its interior, such that for all $i \in {\mathbb{N}}^*$, $\nabla h^*(i, 0)$ is proportional to $(\cos(\theta), (-1)^i \sin(\theta))$. \[cor:legendreOscillating\]

For $x=(u,v)$, set $\|x\|_1=|u|+|v|$, and let $P = \left\{ x \in {\mathbb{R}}^2,\, \|x\|_1 \leq 2 \right\}$. Let us construct a strictly increasing sequence of interpolable polygon-normal pairs $\left( PN_i \right)_{i \in {\mathbb{N}}^*}$ as follows, we fix $\theta \in \left( \frac{-\pi}{4},\frac{\pi}{4} \right)$ and set for all $i\in {\mathbb{N}}^*$ :

-   $T_i = 3i P$, the polygon associated to the $i$-th term $PN_i$ of the sequence,

-   except at the rightmost corner, consider the normals given by the canonical basis vectors and their opposite,

-   at the rightmost corner, $(6i,0)$, one chooses the normal given by the vector $$(\cos(\theta), (-1)^i \sin(\theta)).$$

We now invoke Theorem \[th:Legendre\] to obtain a Lipschitz continuous Legendre function, denoted $h^*$, with full domain having all the $T_i$ as sublevel sets and satisfying the hypotheses of the corollary. Rescaling by a factor $6$ and setting $h = h^{**}$ gives the result.

Counterexamples in continuous optimization
==========================================

We are now in position to apply our interpolation results to build counterexamples to classical problems in convex optimization. We worked on situations ranging from structural questions to qualitative behavior of algorithms and ODEs. Through 9 counterexamples we tried to cover a large spectrum but there are many more possibilities that are left for future research. Some example are constructed from decreasing sequences of convex sets, they can be interpolated using Theorem \[th:smoothinterp\] with $I = {\mathbb{Z}}$, indexing the sequence with negative indices and adding artificially additional sets for positive indices. Nonetheless we sometimes depart from the notations of the first sections and index these sequences by ${\mathbb{N}}$ even though they are decreasing for simplification purposes.

Kurdyka-Łojasiewicz inequality may not hold
-------------------------------------------

The following result is proved in [@bolte2010characterization], it was crucial to construct a $C^2$ convex function which does not satisfy Kurdyka-Łojasiewicz (KL) inequality.

[[@bolte2010characterization Lemma 35]]{} There exists a decreasing sequence of compact convex sets $\left( T_i \right)_{i\in {\mathbb{N}}}$ such that for any $i $ in $ {\mathbb{N}}$, $T_- = T_{i+1}$ and $T_+ = T_i$ satisfy Assumption \[ass:curvature\] and $$\begin{aligned}
                                \sum_{i=0}^{+\infty} {\mathrm{dist}}(T_i, T_{i+1}) = +\infty
                \end{aligned}$$ \[lem:tamsLoja\]

As a corollary, we improve the counterexample in [@bolte2010characterization] and provide a $C^k$ convex counterexamples for any $k\geq 2 $ in $ {\mathbb{N}}$.

There exists a $C^k$ convex function $f \colon {\mathbb{R}}^2 \mapsto {\mathbb{R}}$ which does not satisfy KL inequality. More precisely, for any $r > \inf f$ and $\varphi \colon [\inf f, r] \mapsto {\mathbb{R}}$ continuous and differentiable on $(\inf f, r)$ with $\varphi' > 0$ and $\varphi(\inf f) = 0$, we have $$\begin{aligned}
                    \inf \{\| \nabla (\varphi \circ f)(x) \|: \: x\in {\mathbb{R}}^2, \,\inf f < f(x)  < r \} = 0.
                \end{aligned}$$

Using [@schneider1993convex Theorem 1.8.13], each $T_i$ can be approximated up to arbitrary precision by a polygon. Hence we may assume that all $T_i$ are polygonal while preserving the property of Lemma \[lem:tamsLoja\] as well as Assumption \[ass:curvature\]. Furthermore, using Lemma \[lem:polyGon\] and Remark \[rem:noNormal\] each $T_i$ can in turn be approximated with arbitrary precision by a convex set with $C^k$ boundary and positive curvature. Hence we may also assume that all $T_i$ satisfy both the result of Lemma \[lem:tamsLoja\] and have $C^k$ boundary with nonvanishing curvature. Reversing the order of the sets and adding additional sets artificially, we are in the conditions of application of Theorem \[th:smoothinterp\] with $I = {\mathbb{Z}}$ and the resulting $f$ follows from the same argument as in [@bolte2010characterization Theorem 36].

Block coordinate descent may not converge
-----------------------------------------

![Illustration of the alternating minimization (resp. exact line search) example: on the left, the sublevel sets in gray and the corresponding alternating minimization (resp. exact line search) sequence in dashed lines. On the right the interpolating polygons together with their normal vectors as in Lemma \[lem:polyGon\].[]{data-label="fig:illustrAltMin"}](illustrAltMin1 "fig:"){width=".45\textwidth"}![Illustration of the alternating minimization (resp. exact line search) example: on the left, the sublevel sets in gray and the corresponding alternating minimization (resp. exact line search) sequence in dashed lines. On the right the interpolating polygons together with their normal vectors as in Lemma \[lem:polyGon\].[]{data-label="fig:illustrAltMin"}](illustrAltMin2 "fig:"){width=".4\textwidth"}

The following polygonal construction is illustrated in Figure \[fig:illustrAltMin\]. For any $n\geq2$ in ${\mathbb{N}}$, we set $$\begin{aligned}
                A_n &= \left(\frac{1}{4} + \frac{1}{n},   \frac{1}{4} + \frac{1}{n}\right)\\
                B_n &= \left( \frac{1}{4} + \frac{1}{2(
            n-1)} + \frac{1}{2n}, \frac{1}{4} + \frac{1}{2n} + \frac{1}{2(n+1)} \right)\\
                C_n &= \left(\frac{1}{4} + \frac{1}{n},  - \frac{1}{4} - \frac{1}{n}\right)\\
                D_n &= \left(-\frac{1}{4} - \frac{1}{n},  - \frac{1}{4} - \frac{1}{n}\right)\\
                E_n &= \left(-\frac{1}{4} - \frac{1}{n},  + \frac{1}{4} + \frac{1}{n}\right).\end{aligned}$$ This defines a convex polygon. We may choose the normals at $A_n, C_n, D_n, E_n$ to be bisectors of the corresponding corners and the normal at $B_n$ to be horizontal (see Figure \[fig:illustrAltMin\]). Rotating by an angle of $-\frac{n\pi}{2}$ and repeating the process indefinitely, we obtain the sequence of polygons depicted in Figure \[fig:illustrAltMin\]. It can be checked that the polygons form a strictly decreasing sequence of sets, as for $n>1$, the polygon $A_nB_nC_nD_nE_n$ is contained in the interior of the square $A_{n-1}C_{n-1}D_{n-1}E_{n-1}$. This fulfills the requirement of Corollary \[cor:polygonDecrease\].

There exists a $C^k$ convex function $f \colon {\mathbb{R}}^2 \mapsto {\mathbb{R}}$ and an initialization $x_0=(u_0,v_0)$ such that the recursion, for $i \geq 1$ $$\begin{aligned}
                                u_{i+1} &\in \operatorname*{argmin}_{u} f(u,v_i) \\
                                v_{i+1} &\in \operatorname*{argmin}_{v} f(u_{i+1},v)
                \end{aligned}$$ produces a non converging sequence $\displaystyle (x_i)_{i\in{\mathbb{N}}}=\left((u_i,v_i)\right)_{i\in{\mathbb{N}}}$. \[cor:altMin\]

We apply Corollary \[cor:polygonDecrease\] to the proposed decreasing sequence and by choosing $(u_0,v_0) = B_2$ for example. This requires to shift indices (start with $i = 2$) and use Theorem \[th:smoothinterp\] with $I = {\mathbb{Z}}$. Note that the optimality condition for partial minimization and the fact that level sets have nonvanishing curvature ensure that the partial minima are unique.

In the nonsmooth convex case cyclic minimization is known to fail to provide the infimum value, see e.g., [@auslender p. 94]. Smoothness is sufficient for establishing value convergence (see e.g. [@beck2013convergence; @wright2015coordinate] and references therein), whether it is enough or not for obtaining convergence was an open question. Our counterexample closes this question and shows that cyclic minimization does not yield converging sequences even for $C^k$ convex functions. This result also closes the question for the more general nonconvex case for which we are not aware of a nontrivial counterexample for convergence of alternating minimization. Let us mention however Powell’s example [@powell1973search] which shows that cyclic minimization with three blocks does not converge for smooth functions.

It would also be interesting to understand how our result may impact dual methods and counterexamples in that field, as for instance the recent three blocks counterexample in [@chen].

Gradient descent with exact line search may not converge
--------------------------------------------------------

Gradient descent with exact line search is governed by the recursion: $$x^+\in \operatorname*{argmin}\left\{f(y):y=x-t\nabla f(x),\,t\in{\mathbb{R}}\right\},$$ where $x$ is a point in the plane.

Observe that the step coincides with partial minimization when the gradient $\nabla f(x)$ is colinear to one of the axis of the canonical basis. From the previous section, we thus deduce the following.

There exists a $C^k$ convex function $f \colon {\mathbb{R}}^2 \mapsto {\mathbb{R}}$ and an initialization $z_0$ in the plane such that the recursion, for $i \geq 1$ $$\begin{aligned}
    x_{i+1} \in \operatorname*{argmin}\left\{f(y): y=x_i-t\nabla f(x_i),\, t\in{\mathbb{R}}\right\}
        \end{aligned}$$ produces a well defined non converging sequence $\left(x_i \right)_{i \in {\mathbb{N}}}$. \[cor:exactLineSearch\]

Convergence failure for gradient descent with exact line search is new up to our knowledge. Let us mention that despite non convergence, the constructed sequence satisfy sublinear convergence rates in function values [@beck2013convergence].

Tikhonov regularization path may have infinite length
-----------------------------------------------------

![Illustration of the Tikhonov regularization example, on the left in gray, polygons used to build the sublevel sets of the constructed $f$ and the corresponding solutions to for some values of $r$ (solutions are joined by dotted lines). On the right the normal to be chosen to apply Lemma \[lem:polyGon\] (for $n = 1$, see main text for details). The point $P$ represents $\bar x$, it sits on the $x$-axis and is constantly contained in the normal cone at $B_n$ for any $n \geq 1$.[]{data-label="fig:regTorralba"}](regTorralba1 "fig:"){width=".38\textwidth"}![Illustration of the Tikhonov regularization example, on the left in gray, polygons used to build the sublevel sets of the constructed $f$ and the corresponding solutions to for some values of $r$ (solutions are joined by dotted lines). On the right the normal to be chosen to apply Lemma \[lem:polyGon\] (for $n = 1$, see main text for details). The point $P$ represents $\bar x$, it sits on the $x$-axis and is constantly contained in the normal cone at $B_n$ for any $n \geq 1$.[]{data-label="fig:regTorralba"}](regTorralba2 "fig:"){width=".55\textwidth"}

Following [@torralba96], we consider for any $r > 0$ $$\begin{aligned}
                x(r) = \operatorname*{argmin}\left\{ f(x) + r \|x - \bar{x}\|_2^2:x\in {\mathbb{R}}^2\right\}
                \label{eq:tikhonovRegularization}\end{aligned}$$ where $f$ is $C^k$ convex and where $\bar x$ is any anchor point. We would like to show that the curve $r \mapsto x(r)$ may have infinite length. Torralba provided a counterexample in his PhD Thesis for [*continuous*]{} convex functions, see [@torralba96]. This work extends his result to smooth $C^k$ convex functions in ${\mathbb{R}}^2$.

For any $n $ in $ {\mathbb{N}}^*$, we set $$\begin{aligned}
                A_n &= \left(\frac{2}{n}, \frac{2}{n}\right)\\
                B_n &= \left(\frac{2}{n} + \frac{1}{n^2}, -\frac{1}{n}  \right)\\
                C_n &= \left(\frac{2}{n},  - \frac{2}{n}\right)\\
                D_n &= \left(- \frac{2}{n}, - \frac{2}{n}\right)\\
                E_n &= \left(- \frac{2}{n}, \frac{2}{n}\right).\end{aligned}$$ This is depicted in Figure \[fig:regTorralba\]. For all $n \geq 1$, denote by $M_n$ the point on the $x$ axis above $B_n$ and $N_n$, the intersection of the normal cone at $B_n$ and the $x$ axis. We have $$\begin{aligned}
    \frac{M_nN_n}{M_nB_n} = n \times M_nN_n= \frac{A'_nB_n}{A_nA'_n} = \frac{3/n}{1/n^2} = 3 n,\end{aligned}$$ so that for all $n \geq 1$, $M_nN_n = 3$ and $N_n = (3 + 2/n + 1 / n^2, 0)$. Choosing $P = (7,0)$, since for $n \geq 1$, $3 + 2/n + 1 / n^2 \leq 6 < 7$ , this shows that $P$ constantly belongs to the interior of the normal cone at $B_n$ for all $n \geq 1$. The sequence of level sets is constructed as in Figure \[fig:regTorralba\] by considering alternating symmetries with respect to the $x$-axis of the sequence of polygons above. It can be checked that the polygons form a strictly decreasing sequence of sets, as for $n>1$, the polygon $A_nB_nC_nD_nE_n$ is contained in the interior of the square $A_{n-1}C_{n-1}D_{n-1}E_{n-1}$. We choose the normal at $A_n, C_n, D_n, E_n$ to belong to the bisector at the corner and the normal at $B_n$ to be proportional to the vector $B_n P$. Applying Corollary \[cor:polygonDecrease\], we construct $f$ and choose $\bar{x} = P$ in to obtain the following:

There exists a $C^k$ strictly convex function $f \colon {\mathbb{R}}^2 \mapsto {\mathbb{R}}$ and $\bar{x} \in {\mathbb{R}}^2$ such that the curve $x((0,1))$ given by has infinite length. \[cor:homotopy\]

We apply Corollary \[cor:polygonDecrease\] with $I = {\mathbb{Z}}$ and revert the indices set to match the sequence that we have described. For any $n \geq 1$ there exists a value of $\lambda_n$ such that $f(B_n)= \lambda_n$ and $\nabla f(B_n)$ is colinear to the vector $B_n P$. Set $$\begin{aligned}
            r = \frac{\|\nabla f(B_n)\|}{2 B_nP}
        \end{aligned}$$ we have $\nabla f(B_n) + 2r(B_n - P) = 0$ which is the optimality condition in with $\bar{x} = P$. Hence we have shown that there exists a value of $r$ such that $B_n$ is the solution to . Since $n$ was arbitrary this is true for all $n$ and the curve $r \mapsto x(r)$ has to go through a sequence of points whose second coordinate is of the form $(-1)^n/n$ for all $n \geq 1$. Since this sequence is not absolutely summable, the curve has infinite length.

This result is in contrast with the definable case for which we have finite length by the monotonicity lemma, since the whole trajectory is definable and bounded.

Secants of gradient curves at infinity may not converge
-------------------------------------------------------

#### Thom’s gradient conjecture and Kurdyka-Mostowski-Parusinski’s theorem

A theorem of Łojasiewicz [@lojasiewicz1984trajectoires] asserts that bounded solutions to the gradient system $$\begin{aligned}
    \dot{x}(t) = - \nabla f(x(t))\end{aligned}$$ converge when $f$ is a real analytic potential. Thom conjectured in [@thom1989problemes] that this convergence should occur in a stronger form: trajectories converging to a given $\bar{x}$ should admit a tangent at infinity, that is $$\begin{aligned}
\label{secants}
    \frac{x(t) - \bar{x}}{\|x(t) - \bar{x}\|}\end{aligned}$$ should have a limit as $t \to \infty$. Lines passing through $\bar x$ and having as a slope are called [*secants of $x$ at $\bar x$*]{}. This conjecture was proved to be true in [@kurdyka2000proof]. In the convex world, it is well known that solutions to the gradient system converge for general potentials (this is a Féjer monotonicity argument due to Bruck); see also the original approaches by Manselli and Pucci [@Manselli] and Daniilidis et al. [@dan]. It is then natural to wonder whether this convergence satisfies higher order rigidity properties as in the analytic case. The answer turns out to be negative in general yielding a quite mysterious phase portrait.

#### Absence of tangential convergence for convex potentials

![$A = (-5,0)$, $B=\left( - \frac{5}{3} - \frac{25}{16}, \frac{10}{3} + \frac{5}{4} \right)$, $C=\left( \frac{-5}{2},5 \right)$, $D=\left( 0,5 \right)$, $E = \left( 5,0 \right)$, $F = \left( 0, -5 \right)$. All normals are chosen to be bisectors except $w$ which is parallel to the line $(DE)$. The vector $v$ is orthogonal to the segment $[BC]$. The point $C'$ is obtained by considering the intersection between the line $(Bw)$ (starting from $B$ with direction $w$), and the segment $[OC]$. The points $A'$, $B'$, $D'$, $E'$, $F'$ are obtained by performing a scaling of $A,B,D,E,F$ of a factor $\frac{OC'}{OC}$. The polygon $A''B''C''D''E''F''$ is $ABCDEF$ scaled by a factor $\frac{OC' + OC}{2 OC}$.[]{data-label="fig:thom"}](thom){width=".6\textwidth"}

The construction given in this paragraph is more complex than the previous ones, we start with a technical lemma which will be the basic building block for our counterexample.

Let $S$ be a convex set with $C^k$ boundary interpolating $ABCDEF$ in Figure \[fig:thom\] and let $g$ be the gauge function associated to $S$. The function $g$ is differentiable outside the origin. Consider any initialization $x_0 $ in $ [BC]$ with corresponding trajectory solution to the equation $$\begin{aligned}
                                \dot{x}(t) &= - \nabla g(x(t)), \, t\geq 0,\\
                                x(0) &= x_0.
                \end{aligned}$$ Set $\bar{t} = \sup_{x(t) \in OBC} t$, we have $\bar{t} < + \infty$ and $x(\bar{t}) $ in $ [CC']$. \[lem:techThom\]

The fact that $g$ is differentiable comes from the fact that its subgradient is uniquely determined by the normal cone to $S$ which has dimension one because of the smoothness of the boundary of $S$. Since $S$ is interpolating the polygon, we have $g(B) = g(C) = 1$. Furthermore, we have for all $t$, $\frac{d}{dt} g(x(t)) = -\|\nabla g(x(t))\|^2= -1$, thence $\bar{t} \leq 1 - g(C')$. By homogeneity, for any $x \neq 0$ and $s > 0$, $\nabla g (sx) = \nabla g(x)$. For any $x $ in $ [BC]$, by convexity $$\begin{aligned}
                                0 \leq \left\langle C-x,\nabla g(C) -  \nabla g(x) \right\rangle = \left\langle C-B,\nabla g(C) -  \nabla g(x)  \right\rangle \frac{\|C-x\|}{\|C-B\|},
                \end{aligned}$$ and therefore $$\begin{aligned}
                                -\left\langle C - B, \nabla g(x) \right\rangle \geq  -\left\langle C - B, \nabla g(C) \right\rangle  
                                \label{eq:thomGradient1}
                \end{aligned}$$ By homogeneity of $g$, is true for any $x$ in the triangle $OCB$ (different from $0$) and thus in the triangle $C'CB$ . Denote by $y$ the solution to the equation $$\begin{aligned}
                                \dot{y} &= - \nabla g(C)\\
                                y(0) &= B,
                \end{aligned}$$ which integrates to $y(t) = B - t w$ for all $t$. Equation ensures that for any $0\leq t \leq \bar{t}$ $$\begin{aligned}
                                \frac{d}{dt} (\left\langle C - B, x(t) \right\rangle) &\geq \frac{d}{dt} (\left\langle C - B, y(t) \right\rangle)
                \end{aligned}$$ Hence, we have for any $0\leq t \leq \bar{t}$, integrating on $[0,t]$ $$\begin{aligned}
                                \left\langle C - B, x(t) \right\rangle &\geq \left\langle C - B, y(t)  \right\rangle + \left\langle C - B, x_0 - B  \right\rangle \nonumber\\
                            &\geq \left\langle C - B, y(t)  \right\rangle.
                            \label{eq:thomGradient2}
                \end{aligned}$$ Furthermore, for all $x $ in $ [BC]$, we have $$\begin{aligned}
                    1 = \|\nabla g(x)\|^2 = \frac{1}{\|C - B\|^2}\left\langle C - B, \nabla g(x) \right\rangle^2 + \left\langle v, \nabla g(x) \right\rangle^2,
                \end{aligned}$$ because $v$ is orthogonal to $C-B$. The first term is maximal for $x = C$ and thus the second term is minimal for $x = C$, we have thus for all $x $ in $ [BC]$ $$\begin{aligned}
                    0 < \left\langle \nabla g(C), v \right\rangle = \left\langle - \nabla g(C), -v \right\rangle \leq \left\langle \nabla g(x), v \right\rangle =  \left\langle -\nabla g(x), -v \right\rangle   \leq 1.
                    \label{eq:thomGradient3}
                \end{aligned}$$ Equation holds for all $x$ in $OCB$ different from $O$ by homogeneity. We deduce that for all $0 \leq t \leq \bar{t}$, we have $$\begin{aligned}
                                \frac{d}{dt} (\left\langle -v, x(t) \right\rangle) &\geq \frac{d}{dt} (\left\langle -v, y(t) \right\rangle)
                \end{aligned}$$ and by integration $$\begin{aligned}
                            \left\langle -v, x(t) \right\rangle &\geq \left\langle -v, y(t)  \right\rangle + \left\langle -v, x_0 - B  \right\rangle \nonumber\\
                            &= \left\langle -v, y(t)  \right\rangle.
                            \label{eq:thomGradient4}
                \end{aligned}$$ Hence, in the coordinate system $(C - B, -v)$, which is orthogonal, for all $t $ in $ [0,\bar{t}]$, $x(t)$ has larger coordinates than $y(t)$.

The trajectory $y(t)$, of equation $t \mapsto B - t w$ is the line going from $B$ to $C'$. From equations and , we may write for all $t$ in $ [0,\bar{t}]$, $x(t) = y(t) + \alpha (t) (C-B) + \beta(t) (-v)$ where $\alpha$ and $\beta$ are positive functions. Since $y(t)$ belongs to the line $(BC')$, this shows that $x(t)$ has to be above this line for all $t \geq 0$, $t\leq \bar{t}$ and actually, $x(\bar{t}) $ in $ BCC'$. Hence at time $\bar{t}$, we have $x(\bar{t}) $ in $ [CC']$. This holds true because $x(\bar{t})$ is on the boundary of $OCB$ and on the boundary of $BCC'$. Hence either $x(\bar{t}) $ in $ [CC']$, either $x(\bar{t}) $ in $ [BC]$. Equation ensures that if $x(\bar{t}) $ in $ [BC]$ then $x(\bar{t}) = C$ which concludes the proof.

There exists a $C^k$ strictly convex function on ${\mathbb{R}}^2$ with a unique minimizer $\bar{x}$, such that any nonconstant solution to the gradient flow equation $$\begin{aligned}
                                \dot{x}(t) = -\nabla f(x(t))
                \end{aligned}$$ is such that $$\begin{aligned}
                                \frac{x(t) - \bar{x}}{\|x(t) - \bar{x}\|}
                \end{aligned}$$ does not have a limit as $t \to \infty$.\
The function $f$ has a positive definite Hessian everywhere except at $0$. \[cor:Thom\]

We assume without loss of generality that $\bar{x} = O$ is the origin. Writting $x(t) = (r(t),\theta(t))$ in polar coordinate, we will construct a function $f$ such that each solution to the ODE produces nonconverging trajectories $\theta(t)$.

We start with an interpolating set $S_0 = ABCDE$ as in Lemma \[lem:techThom\] and let $S_1 = A'B'C'D'E'$ be its scaled version as described in Figure \[fig:thom\].

Let $\alpha$ be the value of the angle $\widehat{BOC}$ and $m = \left\lceil \frac{2 \pi}{\alpha} \right\rceil + 1$. We have $$\begin{aligned}
                    \frac{2\pi}{m} < \alpha.
                \end{aligned}$$ To obtain $S_{2}$, we rotate $S_0$ by an angle $2\pi / m$, we denote $S'_0$ the resulting set. We rescale $S'_0$ by a factor $\beta $ in $ (0,1)$ so that $\beta S'_0$ lies in the interior of $S_1$. Call the resulting set $S_2$ and $S_3$ is obtained from $S_2$ exactly the same way as $S_1$ is obtained from $S_0$. We repeat the same process indefinitely to obtain a strictly decreasing sequence of $C^k$ sets. Note that for any $k $ in $ {\mathbb{N}}$, $S_{2km}$ and $S_{2km + 1}$ are homothetic to $S_0$.

We now invoke Corollary \[cor:polygonDecrease\] (with $I = {\mathbb{Z}}$ and revert the indices) to obtain a $C^k$ function $f$ with those prescribed level sets. Using Remark \[rem:alignedLevelSets\] it turns out that the level sets of $f$ between $S_0$ and $S_1$ are simple scalings of $S_0$. Hence the gradient curves of $f$ and those of the gauge function of $S$ are the same between $S_0$ and $S_1$, up to time reparametrization.

Using Lemma \[lem:techThom\] any trajectory crossing $[BC]$ in Figure \[fig:thom\], must also be crossing $[CC']$ and leave the triangle $BOC$ in finite time. The same statement holds after scaling the level sets and since for all $k $ in $ {\mathbb{N}}$, $S_{2km}$ and $S_{2km+1}$ are homothetic to $S_0$, this shows that no solution stays indefinitely in the triangle $BOC$.

Lemma \[lem:techThom\] still holds after rotations and by our construction, for any triangle $T$ obtained by rotating $BOC$ by a multiple of $2\pi/m$, no trajectory stays indefinitely within $T$. Since $2\pi / m < \alpha$, the union of these triangles $U$ contains $O$ in its interior.

Note first that any gradient curve converges to $\bar x$. Let us argue by contradiction and assume that there exists a continuous gradient curve $t \mapsto z(t)$ distinct from the stationary solution $\bar{x}$, such that $$\begin{aligned}
                    \frac{z(t) - \bar{x}}{\|z(t) - \bar{x}\|}
                \end{aligned}$$ converges. This exactly means that the angle $\theta(t)$ of the curve has a limit in $[0,2\pi)$ as $t$ goes to infinity. There is a rotation of $BOC$ by a multiple of $2\pi/m$ whose interior intersects the half line given by the direction $\theta$, call this triangle $T$. The directional convergence entails that there exists $t_0 \geq 0$ such that $z(t)$ belongs to $T$ for all $t \geq t_0$. Hence $z$ can not be a gradient curve. To complete the proof, we may add disks of increasing size to the list of sets to obtain a full domain function and invoke Theorem \[th:smoothinterp\] with $I = {\mathbb{Z}}$.

Newton’s flow may not converge
------------------------------

Given a twice differentiable convex function $f$, we define the open set $\Omega:=\{x\in{\mathbb{R}}^2:\nabla^2f \mbox{ is invertible}\}$ and we consider maximal solutions to the differential equation $$\begin{aligned}
     \dot{x}(t) = - \nabla^2 f(x(t))^{-1} \nabla f(x(t)),
     \label{eq:newton}\end{aligned}$$ on $\Omega$. This is the continuous counterpart of Newton’s method, it has been studied in [@aubin1984differential] and [@alvarez1998dynamical]. Let $x_0$ be in $\Omega$, there exists a unique maximal nontrivial interval $I$ containing $0$ and a unique solution $x$ to on $I$ with $x(0) = x_0$. Equation may be rewritten as $$\begin{aligned}
    \frac{d}{dt} \nabla f(x(t)) = - \nabla f(x(t))\end{aligned}$$ and thus for all $t $ in $ I$, we have $$\begin{aligned}
    \nabla f(x(t)) = e^{-t} \nabla f(x_0).
    \label{eq:newtonIntegrated}\end{aligned}$$ If we could ensure that $I = {\mathbb{R}}$ and $f$ has oscillating gradients close to its minimum, then entails that the direction of the gradient is constant along the solution, which requires oscillations in space to compensate for gradient oscillations.

For any $k\geq2$, there exists a $C^k$ convex coercive function $f \colon {\mathbb{R}}^2 \mapsto {\mathbb{R}}$ and an initial condition $x_0 $ in $ {\mathbb{R}}^2$ such that the solution to is bounded, defined on ${\mathbb{R}}$ and has at least two distinct accumulation points.\[cor:newton\]

The counterexample is sketched in Figure \[fig:illustrNewton\], the construction is the same as for Corollary \[cor:altMin\] but instead of doing quarter rotations, we use symmetry with respect to the first axis. We can then call for Corollary \[cor:polygonDecrease\] to construct the function $f$ and equation ensures that the solution interval is unbounded.

![Illustration of the continuous time Newton’s dynamics. On the left, the “skeletons" of the sublevel sets in gray and a sketch of the corresponding curve. On the right, the normals to be chosen in order to apply Lemma \[lem:polyGon\].[]{data-label="fig:illustrNewton"}](illustrNewton "fig:"){width=".45\textwidth"}![Illustration of the continuous time Newton’s dynamics. On the left, the “skeletons" of the sublevel sets in gray and a sketch of the corresponding curve. On the right, the normals to be chosen in order to apply Lemma \[lem:polyGon\].[]{data-label="fig:illustrNewton"}](illustrAltMin2 "fig:"){width=".4\textwidth"}

Bregman descent (mirror descent) may not converge
-------------------------------------------------

The mirror descent algorithm was introduced in [@newmirovsky1983problem] as an efficient method to solve constrained convex problems. In [@beck2003mirror], this method is shown to be equivalent to a projected subgradient method, using non-Euclidean projections. It plays an important role for some categories of constrained optimization problem; see e.g., [@bauschke2016descent] for recent developments and [@Dragomir19] for a surprising example.

Let us recall beforehand some definitions. Given a Legendre function $h$ with domain ${\mathrm{dom}\,}h$, define the [*Bregman distance*]{}[^5] associated to $h$ as $D_h(u,v)=h(u)-h(v)-\langle \nabla h(v),u-v\rangle$ where $u$ is in $ {\mathrm{dom}\,}h$ and $v$ is in the interior of ${\mathrm{dom}\,}h$.

Given a smooth convex function $f$ that we wish to minimize on $\overline{ {\mathrm{dom}\,}h}$, we consider the Bregman method $$x_{i+1}=\operatorname*{argmin}\left\{\langle \nabla f(x_i),u-x_i\rangle+\lambda D_h(u,x_i):u\in {\mathbb{R}}^2\right\},$$ where $x_0$ is in ${\mathrm{int}\,}{\mathrm{dom}\,}h$ and $\lambda>0$ is a step size. When the above iteration is well defined, e.g. when ${\mathrm{dom}\,}h$ is bounded, it writes: $$x_{i+1}=\nabla h^*\left(\nabla h(x_i)-\lambda \nabla f(x_i)\right).$$ In [@bauschke2016descent] the authors identified a generalized smoothness condition which confers good minimizing properties to the above method: $$\begin{aligned}
\label{hsmooth}Lh-f \mbox{ convex},\\ \lambda\in(0,L).\label{step}\end{aligned}$$

The corollary below shows that such an algorithm may not converge, even though we assume the cost to satisfy , the step to satisfy , and the Legendre function to have a compact domain.

There exists a Legendre function $h\colon D \mapsto {\mathbb{R}}$, defined on a closed square $D$, continuous on $D$, a vector $c$ in $ {\mathbb{R}}^2$, and $x_0$ in ${\mathbb{R}}^2$ such that the Bregman descent recursion $$\begin{aligned}
                    x_{i+1} = \nabla h^*\left( \nabla h(x_i) - c  \right),
                \end{aligned}$$ produces a bounded sequence $(x_i)_{i \in {\mathbb{N}}}$ which has at least two distinct accumulation points. \[cor:mirrorDescent\]

We fix $\theta \in \left( \frac{-\pi}{4},\frac{\pi}{4} \right)$, $\theta \neq 0$, and consider $h$ constructed in Corollary \[cor:legendreOscillating\] and choose $c = (-1,0)$. In this case the Bregman descent recursion writes for all $i$ in ${\mathbb{N}}$, $$\begin{aligned}
                                \nabla h (x_{i+1}) - \nabla h(x_i) = -c  
                \end{aligned}$$ so that we actually have $\nabla h (x_{i}) - \nabla h(x_{0}) = \nabla h (x_{i}) = - i c$ and thus $$x_i = \nabla h^*(-ic) = \nabla h^*(i,0).$$ By Corollary  \[cor:legendreOscillating\], we have for all $i \in {\mathbb{N}}$ that $\nabla h^*(i,0)$ proportional to $(\cos(\theta), (-1)^i \sin(\theta))$. Since the norm of the gradient of $h^*$ cannot vanish at infinitiy (no flat direction) and is bounded, this proves that the sequence $(x_i)_{i \in {\mathbb{N}}}$ has at least two accumulation points which is the desired result.

Central paths of Legendre barriers may not converge
---------------------------------------------------

Consider the problem $$\begin{aligned}
                \min_{x \in D} \left\langle c, x\right\rangle
                \label{eq:LP}\end{aligned}$$ where $D$ is a subset of ${\mathbb{R}}^2$ and $c$. Given a Legendre function $h$ on $D$, we introduce the $h$ central path through $$\begin{aligned}
                x(r) = \operatorname*{argmin}\left\{ \left\langle c, x\right\rangle + r h(x):x\in{\mathbb{R}}^2\right\}
                \label{eq:interiorLegendre}\end{aligned}$$ where $r>0$ is meant to tend to $0$. Central paths are one of the essential tools behind interior point methods, see e.g., [@NN; @Aus99] and references therein.

Note that the accumulation points of $x(r)$ as $r \to 0$, have to be in the the solution set of . It is even tempting to think that the convergence of the path to some specific minimizer could occur, as it is the case for many barriers, see e.g. [@Aus99]. We have however:

There exists a Legendre function $h\colon D \mapsto {\mathbb{R}}$, defined on a closed square $D$, continuous on $D$, a vector $c$ in $ {\mathbb{R}}^2$, such that the $h$ central path $r \mapsto x(r)$ has two distinct accumulation points. \[cor:interiorPoint\]

The optimality condition which characterizes $x(r)$ for any $r > 0$ writes, $$\begin{aligned}
                                x(r)= \nabla h^*\left( \frac{c}{r} \right),
                \end{aligned}$$ and the construction is the same as in Corollary \[cor:mirrorDescent\].

Hessian Riemannian gradient dynamics may not converge
-----------------------------------------------------

The construction of this paragraph is similar to the two previous paragraphs. Consider a $C^k$ ($k\geq 2$) Legendre function $h \colon D \mapsto {\mathbb{R}}$ and the continuous time dynamics $$\begin{aligned}
    \dot{x}(t) = - \nabla_H f(x(t)), \, t\geq 0,
    \label{eq:hessianRiemannian}\end{aligned}$$ where $H = \nabla^2 h$ is the Hessian of $h$ and $\nabla_H f = H^{-1}  \nabla f$ is the gradient of some differentiable function $f$ in the Riemannian metric induced by $H$ on ${\mathrm{int}\,}D$. Such dynamics were considered in [@BT03; @alvarez2004hessian].

We have the following result:

There exists a Legendre function $h\colon D \mapsto {\mathbb{R}}$, defined on a closed square $D$, continuous on $D$, a vector $c$ in $ {\mathbb{R}}^2$, and $x_0$ in ${\mathbb{R}}^2$ such that the solution to with $f=\langle c,\cdot\rangle$ has two distinct accumulation points. \[cor:hessianRiemannian\]

Equation may be rewritten $$\begin{aligned}
                                \frac{d}{dt} \nabla h(x(t))= - \nabla f(x(t)),
                \end{aligned}$$ so choosing $c = (-1,0)$, we have for all $t \in {\mathbb{R}}$, $\nabla h(x(t)) = \nabla h(x(0)) + (t,0) = (t,0)$ and the construction is the same as in Corollary \[cor:mirrorDescent\].

Appendix
========

There exists a $C^\infty$ strictly increasing concave function $\phi\colon [0,1] \mapsto [0,1]$ such that $$\begin{aligned}
                                \phi(t) &= \sqrt{2t/3} \quad \forall t \leq 1/6\\
                                \phi(1) &= 1 \\
                                \phi'(1) &= 2/3\\
                                \phi^{(m)}(1) &= 0, \quad \forall m \geq 2
                \end{aligned}$$ \[lem:interpolationAroundZero\]

Consider a $C^\infty$ function $g_0 \colon {\mathbb{R}}\mapsto [0,1]$ such that $g_0 = 1$ on $(-\infty,-1)$, $g_0 = 0$ on $(1, +\infty)$ (for example convoluting the step function with a smooth bump function). Set $g(t) = \frac{1}{2}\left( g_0(t) + 1 - g_0(-t) \right)$ we have that $g$ is $C^\infty$, $g = 1$ on $(-\infty,-1)$, $g = 0$ on $(1, +\infty)$ and $g(t) + g(-t) = 1$ for all $t$. We have $$\begin{aligned}
                                \int_{-1}^1 g(s) ds = 1\\
                                \int_{-1}^1 \left( \int_{-1}^t g(s)ds \right) dt = 1
                \end{aligned}$$ Set $\phi_0 \colon [-3,3] \mapsto {\mathbb{R}}$, such that $$\begin{aligned}
                                \phi_0(t) = \int_{-3}^t \left( \int_{-3}^r g(s) ds\right)dr.
                \end{aligned}$$ For all $r $ in $ [-3,3]$, we have $$\begin{aligned}
                                \int_{-3}^r g(s) ds = 
                                \begin{cases}
                                                r+3 & \text{ if } r \leq -1\\
                                                2 + \int_{-1}^r g(s)ds & \text{ if } -1 \leq r \leq 1\\
                                                3 & \text{ if } r \geq 1
                                \end{cases}
                \end{aligned}$$ and thus $$\begin{aligned}
                                \phi_0(t) = 
                                \begin{cases}
                                                \frac{t^2}{2} - 9/2 + 3(t+3) & \text{ if } t \leq -1\\
                                                2 + 2(t + 1) + \int_{-1}^t \left( \int_{-1}^r g(s)ds \right)dr& \text{ if } -1 \leq t \leq 1\\
                                                6 + 3(t-1) & \text{ if } 1 \geq t
                                \end{cases}
                \end{aligned}$$ and in particular $\phi_0(3) = 12$ and $\phi_0'(3) = 3$. Set $\phi_1(s) = \phi_0(6 s -3)/12$. $$\begin{aligned}
                                \phi_1(0) &= 0\\
                                \phi_1(t) &= \left(\frac{(6t-3)^2}{2} - 9/2 + 2(3t )\right)/12 =  3 t^2 / 2 =  \text{ if } t \leq 1/3\\
                                \phi_1(1) &= 1\\
                                \phi_1'(1) &= 3/2.
                \end{aligned}$$ $\phi_1$ is stricly increasing, let $\phi \colon [0,1] \mapsto [0,1]$ denote the inverse of $\phi_1$, we have $$\begin{aligned}
                                \phi(1) &= 1\\
                                \phi'(1) & = 2 / 3\\
                                \phi(t) &=  \sqrt{2t/3} \text{ if } t\leq 1/6.
                \end{aligned}$$

Consider any strictly increasing $C^k$ function $\phi \colon (0,2) \mapsto {\mathbb{R}}$ such that $\phi(1) = 1$ and $\phi^{(m)}(1) = 0$, $m = 2,\ldots k$. Then the function $$\begin{aligned}
                                G \colon (0,2) \times {\mathbb{R}}/ 2 \pi {\mathbb{Z}}& \mapsto {\mathbb{R}}^2 \\
                                (s,\theta) &\mapsto \phi(s) n(\theta)
                \end{aligned}$$ is diffeomorphism which satisfies for any $m=1 \ldots,k$ and $l =2,\ldots, k$, $$\begin{aligned}
                                &\frac{\partial^m G}{\partial \theta^m}(1,\theta) = n^{(m)}(\theta)\\
                                &\frac{\partial^{m+1} G}{\partial \lambda\partial \theta^m} (1,\theta) = \phi'(1)n^{(m)}(\theta) \\
                                &\frac{\partial^{l+m} G}{\partial \lambda^l \partial \theta^m }(\lambda_-,\theta) = 0.
                \end{aligned}$$ \[lem:diffGauge2\]

[**: Combinatorial Arbogast-Faà di Bruno Formula (from [@ma2009higher]).**]{} Let $g \colon {\mathbb{R}}\mapsto {\mathbb{R}}$ and $f \colon {\mathbb{R}}^p \mapsto [0, +\infty)$ be $C^k$ functions. Then we have for any $m \leq k$ and any indices $i_1,\ldots,i_m \in \left\{ 1,\ldots, p \right\}$. $$\begin{aligned}
                                \frac{\partial^m}{\prod_{l=1}^{m}\partial x_{i_l}} g \circ f(x) = \sum_{\pi \in {\mathcal{P}}} g^{(|\pi|)}(f(x)) \prod_{B \in \pi} \frac{\partial^{|B|} f}{\prod_{l \in B}\partial x_{i_l}}(x),
                \end{aligned}$$ where ${\mathcal{P}}$ denotes all partitions of $\left\{ 1,\ldots, m \right\}$, the product is over subsets of $\left\{ 1,\ldots,m \right\}$ given by the partition $\pi$ and $|\cdot|$ denotes the number of elements of a set. We rewrite this as follows $$\begin{aligned}
                                \frac{\partial^m}{\prod_{l=1}^{m}\partial x_{i_l}} g \circ f(x) = \sum_{k = 1}^m\sum_{\pi \in {\mathcal{P}}_k} g^{(k)}(f(x)) \prod_{B \in \pi} \frac{\partial^{|B|} f}{\prod_{l=1}^{m}\partial x_{i_l}}(x),
                \end{aligned}$$ where ${\mathcal{P}}_k$ denotes all partitions of size $k$ of $\left\{ 1,\ldots, m \right\}$. \[lem:faaDiBruno\]

[From [@bolte2010characterization Lemma 45]]{} Let $h $ in $ C^0\left( (0,r_0],{\mathbb{R}}_+^* \right)$ be an increasing function. Then there exists a function $\psi $ in $ C^\infty({\mathbb{R}},{\mathbb{R}}_+)$ such that $\psi = 0$ on, ${\mathbb{R}}_-$ and $0 < \psi(s) \leq h(s)$ for any $s$ in $(0,r_0]$ and $\psi$ is increasing on ${\mathbb{R}}$ \[lem:CinfiniteLowerBound\]

Let $D \subset {\mathbb{R}}^p$ be a nonempty compact convex set and $f \colon D \mapsto {\mathbb{R}}$ convex, continuous on $D$ and $C^k$ on $D \setminus \operatorname*{argmin}_{D} f$. Assume further that $\operatorname*{argmin}_D f \subset \mathrm{int}(D)$, $k \geq 1$, with $\min_D f = 0$. Then there exists $\phi \colon {\mathbb{R}}\mapsto {\mathbb{R}}_+$, $C^k$, convex and increasing with positive derivative on $(0,+\infty)$, such that $\phi\circ f$ is convex and $C^k$ on $D$. \[lem:CkSmoothing\]

By a simple translation, we may assume that $\min_D f = 0$ and $\max_D f = 1$. Any convex function is locally Lipschitz continuous on the interior of its domain so that $f$ is globally Lipschitz continuous on $D$ and its gradient is bounded. Hence, $f^2$ is $C^1$ and convex on $D$. We now proceed by recursion. For any $m =1,\ldots, k$, we let $Q_m$ denote the $m$-order tensor of partial derivatives of order $m$. Fix $m$ in $\{1,\ldots,k\}$. Assume that $f$ is $C^m$ throughout $D$ while it is $C^{m+1}$ on $D \setminus \arg\min_D f$. Note that all the derivatives up to order $m$ are bounded. We wish to prove that $f$ is globally $C^{m+1}$.

Consider the increasing function $$\begin{aligned}
                                h \colon (0,1] &\mapsto {\mathbb{R}}_+^*\\
                                s &\mapsto \frac{s}{1 + \sup_{s \leq f(x) \leq 1}\|Q_{m+1}(x)\|_{\infty}}
                \end{aligned}$$ and set $\psi$ as in Lemma \[lem:CinfiniteLowerBound\]. Recall that $\psi$ is $C^\infty$ and all its derivative vanish at $0$ and $\psi \leq h$ on $(0,1]$. Let $\phi$ denote the anti-derivative of $\psi$ such that $\phi(0) = 0$. $\phi$ is $C^\infty$ and convex increasing on ${\mathbb{R}}$ and, since its derivatives at $0$ vanish as well, one has, for any $q $ in $ {\mathbb{N}}$, $\phi^{(q)}(z) = o(z)$. Consider the function $\phi \circ f$. It is $C^m$ on $D$ and it has bounded derivatives up to order $m$. Furthermore, it is $C^{m+1}$ on $D \setminus \operatorname*{argmin}_D f$. Let $\bar{y} $ in $ \operatorname*{argmin}_D f$. If $\bar{y} $ in $ \mathrm{int}(\operatorname*{argmin}_D f)$, then $f$ and $\phi \,\circ f$ have derivatives of all order vanishing at $\bar{y}$. Assuming that $\bar{y} $ in $ \operatorname*{argmin}_D f\setminus \mathrm{int}(\operatorname*{argmin}_D f)$. By the induction assumption and Lemma \[lem:faaDiBruno\], we have for any indices $i_1,\ldots,i_m \in \left\{ 1,\ldots, p \right\}$ and any $h $ in $ {\mathbb{R}}^p$: $$\begin{aligned}
                                &\frac{\partial^m}{\prod_{l=1}^{m}\partial x_{i_l}} (\phi \circ f)(\bar{y} + z) - \frac{\partial^m}{\prod_{l=1}^{m}\partial x_{i_l}} (\phi \circ f)(\bar{y}) \\
                                =\;& \frac{\partial^m}{\prod_{l=1}^{m}\partial x_{i_l}}( \phi \circ f)(\bar{y} + z) \\
                                =\;&\sum_{q = 1}^{m}\sum_{\pi \in {\mathcal{P}}_q} \phi^{(q)}(f(\bar{y} + z)) \prod_{B \in \pi} \frac{\partial^{|B|} f}{\prod_{l=1}^{m}\partial x_{i_l}}(\bar{y} + z).
                \end{aligned}$$ All the derivatives of $f$ are of order less or equal to $m$ and thus remain bounded as $z \to 0$. Further more $f$ is Lipschitz continuous on $D$ so that $f(\bar{y} + z) = O(\|z\|)$ near $0$, and, for any $q $ in $ {\mathbb{N}}$, $\phi^{(q)}(f(\bar{y} + z)) = o(\|z\|)$. Hence $\phi \circ f$ has derivative of order $m+1$ at $\bar{y}$ and it is $0$.

Since $\operatorname*{argmin}_D f \subset \mathrm{int}(D)$, we may consider any sequence of point $(y_{j})_{j \in {\mathbb{N}}}$ in $D \setminus \operatorname*{argmin}_D f$ converging to $\bar{y}$. By Lemma \[lem:faaDiBruno\], we have for any indices $i_1,\ldots,i_{m+1} \in \left\{ 1,\ldots, p \right\}$, and any $j $ in $ {\mathbb{N}}$, $$\begin{aligned}
                                \frac{\partial^{(m+1)}}{\prod_{l=1}^{m+1}\partial x_{i_l}} (\phi \circ f)(y_j) &= \phi'(f(y_j)) \frac{\partial^{(m+1)} f}{\prod_{l=1}^{m}\partial x_{i_l}}(y_j) + \sum_{q = 2}^{m+1}\sum_{\pi \in \Pi_q} \phi^{(q)}(f(y_j)) \prod_{B \in \pi} \frac{\partial^{|B|} f}{\prod_{l=1}^{m}\partial x_{i_l}}(x)\\
                                &\leq h(f(y_j))\frac{\partial^{(m+1)} f}{\prod_{l=1}^{m}\partial x_{i_l}}(y_j) + \sum_{q = 2}^{m+1}\sum_{\pi \in \Pi_q} \phi^{(q)}(f(y_j)) \prod_{B \in \pi} \frac{\partial^{|B|} f}{\prod_{l=1}^{m}\partial x_{i_l}}(x)\\
                                &= f(y_j) \frac{\frac{\partial^{(m+1)} f}{\prod_{l=1}^{m}\partial x_{i_l}}(y_j)}{1 + \sup_{f(y_j) \leq f(x) \leq 1}\|Q_{m+1}(x)\|_{\infty}}  + O(f(y_j))\\
                                &= O(f(y_j)),
                \end{aligned}$$ where the inequality follows from the construction of $\phi$. The third step follows using the definition of $h$ and the fact that, for any $q \geq 2$,

1.  Each partition of $\left\{ 1,\ldots,m+1 \right\}$ of size $q$ contains subsets of size at most $m$. Thus in the product, the terms $\partial^{|B|} f$ correspond to bounded derivatives of $f$ by the induction hypothesis.

2.  $\phi^{(q)}(a) = o(a)$ as $a \to 0$.

The last step stems from the fact that the ratio has asbolute value less than $1$. This shows that the derivatives of order $m+1$ of $\phi \circ f$ are decreasing to $0$ as $j \to \infty$ and $\phi \circ f$ is actually $C^{m+1}$ and convex on $D$. The result follows by induction up to $m = k$ and by the fact that a composition of increasing convex functions is increasing and convex.

Let $p \colon {\mathbb{R}}_+ \mapsto {\mathbb{R}}_+$ be concave increasing and $C^1$ with $p' \geq c$ for some $c > 0$. Assume that there exists $ A > 0$ such that for all $x $ in $ {\mathbb{R}}_+$ $$\begin{aligned}
                    p(x) - x p'(x) \leq A.
                \end{aligned}$$ Then setting $a= A/c$, we have for all $x \geq a$, $$\begin{aligned}
                    p(x-a) - x p'(x-a) \leq 0
                \end{aligned}$$ \[lem:techConcaveAsymptotic\]

For all $x \geq a$, we have $$\begin{aligned}
        f(x-a) - (x-a)f'(x-a) \leq A,
    \end{aligned}$$ hence $$\begin{aligned}
        f(x-a) - xf'(x-a) \leq A - af'(x-a) \leq A - ac = 0.
    \end{aligned}$$

[**Acknowledgements.**]{} The authors acknowledge the support of AI Interdisciplinary Institute ANITI funding, through the French “Investing for the Future – PIA3” program under the Grant agreement n°ANR-19-PI3A-0004, Air Force Office of Scientific Research, Air Force Material Command, USAF, under grant numbers FA9550-19-1-7026, FA9550-18-1-0226, and ANR MasDol. J. Bolte acknowledges the support of ANR Chess, grant ANR-17-EURE-0010 and ANR OMS.

[apalike]{}

Alvarez, F., Bolte, J. and Brahic, O. (2004). Hessian Riemannian gradient flows in convex programming. SIAM journal on control and optimization, 43(2), 477–501.

Alvarez D. F. and Pérez C., J. M. (1998). A dynamical system associated with Newton’s method for parametric approximations of convex minimization problems. Applied Mathematics and Optimization, 38, 193–217.

Auslender, A., Optimisation Méthodes Numériques, Masson, Paris, France, 1976.

Auslender, A. (1999). Penalty and barrier methods: a unified framework. SIAM Journal on Optimization, 10(1), 211–230.

Aubin, J.-P. and Cellina, A. (1984). Differential inclusions: set-valued maps and viability theory. Springer.

Bauschke, H. H., Bolte, J. and Teboulle, M. (2016). A descent lemma beyond Lipschitz gradient continuity: first-order methods revisited and applications. Mathematics of Operations Research, 42(2), 330–348.

Bauschke, H. H., and Combettes, P. L. (2011). Convex analysis and monotone operator theory in Hilbert spaces (Vol. 408). New York: Springer.

Beck, A. (2017). First-order methods in optimization (Vol. 25). SIAM.

Beck, A. and Teboulle, M. (2003). Mirror descent and nonlinear projected subgradient methods for convex optimization. Operations Research Letters, 31(3), 167–175.

Beck, A. and Tetruashvili, L. (2013). On the convergence of block coordinate descent type methods. SIAM journal on Optimization, 23(4), 2037-2060.

Bertsekas, D. P., Scientific, A. (2015). Convex optimization algorithms. Belmont: Athena Scientific.

Bolte, J., Daniilidis, A., Ley, O., and Mazet, L. (2010). Characterizations of Łojasiewicz inequalities: subgradient flows, talweg, convexity. Transactions of the American Mathematical Society, 362(6).

Bolte, J., Nguyen, T. P., Peypouquet, J., and Suter, B. W. (2017). From error bounds to the complexity of first-order descent methods for convex functions. Mathematical Programming, 165(2), 471–507.

Bolte, J., and Teboulle, M. (2003). Barrier operators and associated gradient-like dynamical systems for constrained minimization problems. SIAM journal on control and optimization, 42(4), 1266-1292.

Borwein, J. M., Li, G., and Yao, L. (2014). Analysis of the convergence rate for the cyclic projection algorithm applied to basic semialgebraic convex sets. SIAM Journal on Optimization, 24(1), 498–527.

Boyd, S., and Vandenberghe, L. (2004). Convex optimization. Cambridge university press.

Chen, C., He, B., Ye, Y., and Yuan, X. (2016). The direct extension of ADMM for multi-block convex minimization problems is not necessarily convergent. Mathematical Programming, 155(1-2), 57–79.

Crouzeix, J.-P. (1980). Conditions for convexity of quasiconvex functions. Mathematics of Operations Research, 5(1), 120-125.

Daniilidis, A., Ley, O., and Sabourau, S. (2010). Asymptotic behaviour of self-contracted planar curves and gradient orbits of convex functions. Journal de matéhmatiques pures et appliquées, 94(2), 183–199.

Dragomir, R. A., Taylor, A., d’Aspremont, A., and Bolte, J. (2019). Optimal complexity and certification of Bregman first-order methods. arXiv preprint arXiv:1911.08510.

Fenchel, W., Convex Cones, Sets and Functions, Mimeographed lecture note, Princeton University, 1951.

de Finetti, B. (1949). Sulle stratificazioni convesse. Annali di Matematica Pura ed Applicata, 30(1), 173–183.

Gale, D., Klee, V. and Rockafellar, R. T. (1968). Convex functions on convex polytopes. Proceedings of the American Mathematical Society, 19(4), 867–873.

Golub, G. H., Hansen, P. C., and O’Leary, D. P. (1999). Tikhonov regularization and total least squares. SIAM Journal on Matrix Analysis and Applications, 21(1), 185–194.

Kannai, Y. (1977). Concavifiability and constructions of concave utility functions. Journal of mathematical Economics, 4(1), 1–56.

Kurdyka, K., Mostowski, T. and Parusinski, A. (2000). Proof of the gradient conjecture of R. Thom. Annals of Mathematics, 152(3), 763–792.

Łojasiewicz, S. (1984), Sur les trajectoires du gradient d’une fonction analytique. Seminari di Geometria, Bologna (1982/83), Universita’ degli Studi di Bologna, Bologna (1984), 115–117.

Lorentz, G. G. (1954). Bernstein polynomials. American Mathematical Soc.

Ma, T. W. (2009). Higher chain formula proved by combinatorics. The electronic journal of combinatorics, 16(1).

Manselli, P., and Pucci, C. (1991). Maximum length of steepest descent curves for quasi-convex functions. Geometriae Dedicata, 38(2), 211–227.

A. S. Nemirovsky and D. B. Yudin (1983). Problem Complexity and Method Efficiency in Optimization. Wiley-Interscience, New York.

Nesterov, Y. (2003). Lectures on convex optimization (Vol. 137). Springer.

Nesterov, Y., and Nemirovskii, A. (1994). Interior-point polynomial algorithms in convex programming (Vol. 13). Siam.

Powell, M. J. (1973). On search directions for minimization algorithms. Mathematical programming, 4(1), 193–201.

Schneider, R. (1993). Convex bodies: the Brunn-Minkowski theory (No. 151). Cambridge university press.

Torralba, D., Convergence épigraphique et changements d’échelle en analyse variationnelle et optimisation, Ph.D. Thesis (Université Montpellier 2, 1996).

Rockafellar, R. T. (1970). Convex analysis (Vol. 28). Princeton university press.

Thom, R. (1989). Problèmes rencontrés dans mon parcours mathématique : un bilan. Publications mathématiques de l’IHES, 70, 199–214.

Wright, S. J. (2015). Coordinate descent algorithms. Mathematical Programming, 151(1), 3–34.

[^1]: Toulouse School of Economics, Université Toulouse 1 Capitole, France.

[^2]: IRIT, ANITI, Université de Toulouse, CNRS. DEEL, IRT Saint Exupery, Toulouse, France.

[^3]: In the sense of sets inclusion the sequence being indexed on ${\mathbb{N}}$ or ${\mathbb{Z}}$

[^4]: By structural, we include homotopic deformations by mere summation

[^5]: It is actually not a proper distance.
