ScalingStacks

3.3. Operations on concave functions and duality [02LH]

Original official author HTML, exact retained edition. Historical TeX conversion verdicts remain unchanged. Cited-edition alignment and mathematical self-containment are not assessed.

Complete original source context · Original author HTML

3.3. Operations on concave functions and duality

In this section we consider the basic operations on concave functions and their interplay with the Legendre-Fenchel duality.

Let f1f_{1} and f2f_{2} be two concave functions such that their stability sets are not disjoint. Their sup-convolution is the function

f1⊞f2:Mℝ⟶ℝ¯,v⟼supu1+u2=v(f1​(u1)+f2​(u2)).f_{1}\boxplus f_{2}\colon M_{\mathbb{R}}\longrightarrow{\underline{\mathbb{R}}},\quad v\longmapsto\sup_{u_{1}+u_{2}=v}(f_{1}(u_{1})+f_{2}(u_{2})).

This is a concave function whose effective domain is the Minkowski sum dom⁡(f1)+dom⁡(f2){\operatorname{dom}}(f_{1})+{\operatorname{dom}}(f_{2}). This operation is associative and commutative whenever the terms are defined.

The operations of pointwise addition and sup-convolution are dual to each other. When working with general concave functions, there are some technical issues in this duality that will disappear when considering uniform limits of piecewise affine concave functions.

Proposition 3.38.

Let f1,…,flf_{1},\dots,f_{l} be concave functions.

  1. (1)

    If stab⁡(f1)∩⋯∩stab⁡(fl)≠∅\operatorname{stab}(f_{1})\cap\dots\cap\operatorname{stab}(f_{l})\not=\emptyset, then

    (f1⊞⋯⊞fl)∨=f1∨+⋯+fl∨.(f_{1}\boxplus\dots\boxplus f_{l})^{\vee}=f_{1}^{\vee}+\dots+f_{l}^{\vee}.
  2. (2)

    If dom⁡(f1)∩⋯∩dom⁡(fl)≠∅{\operatorname{dom}}(f_{1})\cap\dots\cap{\operatorname{dom}}(f_{l})\not=\emptyset, then

    (cl⁡(f1)+⋯+cl⁡(fl))∨=cl⁡(f1∨⊞⋯⊞fl∨).({\operatorname{cl}}(f_{1})+\dots+{\operatorname{cl}}(f_{l}))^{\vee}={\operatorname{cl}}(f_{1}^{\vee}\boxplus\dots\boxplus f_{l}^{\vee}).
  3. (3)

    If ri⁡(dom⁡(f1))∩⋯∩ri⁡(dom⁡(fl))≠∅\operatorname{ri}({\operatorname{dom}}(f_{1}))\cap\dots\cap\operatorname{ri}({\operatorname{dom}}(f_{l}))\not=\emptyset, then

    (f1+⋯+fl)∨=f1∨⊞⋯⊞fl∨.(f_{1}+\dots+f_{l})^{\vee}=f_{1}^{\vee}\boxplus\dots\boxplus f_{l}^{\vee}.
Proof.

This is proved in [Roc70, Theorem 16.4]. ∎

Remark 3.39.

When some of the fif_{i}, say f1,…,fkf_{1},\dots,f_{k}, are piecewise affine, the statement (3) of the previous proposition holds under the weaker hypothesis [Roc70, Theorem 20.1]

dom⁡(f1)∩⋯∩dom⁡(fk)∩ri⁡(dom⁡(fk+1))∩⋯∩ri⁡(dom⁡(fl))≠∅.{\operatorname{dom}}(f_{1})\cap\dots\cap{\operatorname{dom}}(f_{k})\cap\operatorname{ri}({\operatorname{dom}}(f_{k+1}))\cap\dots\cap\operatorname{ri}({\operatorname{dom}}(f_{l}))\not=\emptyset.

Let ff be a concave function. For λ>0\lambda>0, the left and right scalar multiplication of ff by λ\lambda are the functions defined, for u∈Nℝu\in N_{\mathbb{R}}, by (λ​f)​(u)=λ​f​(u)(\lambda f)(u)=\lambda f(u) and (f​λ)​(u)=λ​f​(u/λ)(f\lambda)(u)=\lambda f(u/\lambda) respectively. For a point u0∈Nℝu_{0}\in N_{\mathbb{R}}, the translate of ff by u0u_{0} is the concave function defined as (τu0​f)​(u)=f⁡(u−u0)(\tau_{u_{0}}f)(u)=f(u-u_{0}) for u∈Nℝu\in N_{\mathbb{R}}.

Proposition 3.40.

Let ff be a concave function on NℝN_{\mathbb{R}}, λ>0\lambda>0, u0∈Nℝu_{0}\in N_{\mathbb{R}} and x0∈Mℝx_{0}\in M_{\mathbb{R}}. Then

  1. (1)

    dom⁡(λ​f)=dom⁡(f){\operatorname{dom}}(\lambda f)={\operatorname{dom}}(f), stab⁡(λ​f)=λ​stab⁡(f)\operatorname{stab}(\lambda f)=\lambda\operatorname{stab}(f) and (λ​f)∨=f∨​λ(\lambda f)^{\vee}=f^{\vee}\lambda;

  2. (2)

    dom⁡(f​λ)=λ​dom⁡(f){\operatorname{dom}}(f\lambda)=\lambda{\operatorname{dom}}(f), stab⁡(f​λ)=stab⁡(f)\operatorname{stab}(f\lambda)=\operatorname{stab}(f) and (f​λ)∨=λ​f∨(f\lambda)^{\vee}=\lambda f^{\vee};

  3. (3)

    dom⁡(τu0​f)=dom⁡(f)+u0{\operatorname{dom}}(\tau_{u_{0}}f)={\operatorname{dom}}(f)+u_{0}, stab⁡(τu0​f)=stab⁡(f)\operatorname{stab}(\tau_{u_{0}}f)=\operatorname{stab}(f) and (τu0​f)∨=f∨+u0(\tau_{u_{0}}f)^{\vee}=f^{\vee}+u_{0};

  4. (4)

    dom⁡(f+x0)=dom⁡(f){\operatorname{dom}}(f+x_{0})={\operatorname{dom}}(f), stab⁡(f+x0)=stab⁡(f)+x0\operatorname{stab}(f+x_{0})=\operatorname{stab}(f)+x_{0} and (f+x0)∨=τx0​f∨(f+x_{0})^{\vee}=\tau_{x_{0}}f^{\vee}.

Proof.

This follows easily from the definitions. ∎

We next consider direct and inverse images of concave functions by affine maps. Let QℝQ_{\mathbb{R}} be a another finite dimensional real vector space and set Pℝ=Qℝ∨P_{\mathbb{R}}=Q_{\mathbb{R}}^{\vee} for its dual space. For a linear map H:Qℝ→NℝH\colon Q_{\mathbb{R}}\to N_{\mathbb{R}} we denote by H∨:Mℝ→PℝH^{\vee}\colon M_{\mathbb{R}}\to P_{\mathbb{R}} the dual map. We need the following lemma in order to properly define direct images.

Lemma 3.41.

Let H:Qℝ→NℝH\colon Q_{\mathbb{R}}\to N_{\mathbb{R}} be a linear map and gg a concave function on QℝQ_{\mathbb{R}}. If stab⁡(g)∩im⁡(H∨)≠∅\operatorname{stab}(g)\cap\operatorname{im}(H^{\vee})\not=\emptyset then, for all u∈Nℝu\in N_{\mathbb{R}},

supv∈H−1​(u)g⁡(v)<∞.\sup_{v\in H^{-1}(u)}g(v)<\infty.
Proof.

Let x∈Mℝx\in M_{\mathbb{R}} such that H∨​(x)∈stab⁡(g)H^{\vee}(x)\in\operatorname{stab}(g). By the definition of the stability set, supv∈Qℝ(g⁡(v)−⟨H∨​(x),v⟩)<∞\sup_{v\in Q_{\mathbb{R}}}(g(v)-\langle H^{\vee}(x),v\rangle)<\infty. Thus, for any u∈Nℝu\in N_{\mathbb{R}},

supv∈Qℝ(g⁡(v)−⟨H∨​(x),v⟩)\displaystyle\sup_{v\in Q_{\mathbb{R}}}(g(v)-\langle H^{\vee}(x),v\rangle) =supv∈Qℝ(g⁡(v)−⟨x,H⁡(v)⟩)\displaystyle=\sup_{v\in Q_{\mathbb{R}}}(g(v)-\langle x,H(v)\rangle)
≥supv∈H−1​(u)(g⁡(v)−⟨x,H⁡(v)⟩)=supv∈H−1​(u)g⁡(v)−⟨x,u⟩\displaystyle\geq\sup_{v\in H^{-1}(u)}(g(v)-\langle x,H(v)\rangle)=\sup_{v\in H^{-1}(u)}g(v)-\langle x,u\rangle

and so supv∈H−1​(u)g⁡(v)\sup_{v\in H^{-1}(u)}g(v) is bounded above, as stated. ∎

Definition 3.42.

Let A:Qℝ→NℝA\colon Q_{\mathbb{R}}\to N_{\mathbb{R}} be an affine map defined as A=H+u0A=H+u_{0} for a linear map HH and a point u0∈Nℝu_{0}\in N_{\mathbb{R}}. Let ff be a concave function on NℝN_{\mathbb{R}} such that dom⁡(f)∩im⁡(A)≠∅{\operatorname{dom}}(f)\cap\operatorname{im}(A)\not=\emptyset and gg a concave function on QℝQ_{\mathbb{R}} such that stab⁡(g)∩im⁡(H∨)≠∅\operatorname{stab}(g)\cap\operatorname{im}(H^{\vee})\not=\emptyset. Then the inverse image of ff by AA is defined as

A∗​f:Qℝ⟶ℝ,v⟼f∘A⁡(v),A^{\ast}f\colon Q_{\mathbb{R}}\longrightarrow\mathbb{R},\quad v\longmapsto f\circ A(v),

and the direct image of gg by AA is defined as

A∗​g:Nℝ⟶ℝ,u⟼supv∈A−1​(u)g⁡(v).A_{\ast}g\colon N_{\mathbb{R}}\longrightarrow\mathbb{R},\quad u\longmapsto\sup_{v\in A^{-1}(u)}g(v).

It is easy to see that the inverse image A∗​fA^{\ast}f is concave with effective domain dom⁡(A∗​f)=A−1​(dom⁡(f)){\operatorname{dom}}(A^{\ast}f)=A^{-1}({\operatorname{dom}}(f)). Similarly, the direct image A∗​gA_{\ast}g is concave with effective domain dom⁡(A∗​g)=A⁡(dom⁡(g)){\operatorname{dom}}(A_{\ast}g)=A({\operatorname{dom}}(g)), thanks to Lemma 3.41.

The inverse image of a closed function is also closed. In contrast, the direct image of a closed function is not necessarily closed: consider for instance the indicator function ιC\iota_{C} of the set C={(x,y)∈ℝ2∣xy≥1,x>0}C=\{(x,y)\in\mathbb{R}^{2}\mid xy\geq 1,x>0\}, which is a closed concave function. Let A:ℝ2→ℝA\colon\mathbb{R}^{2}\to\mathbb{R} be the first projection. Then A∗​ιCA_{\ast}\iota_{C} is the indicator function of the subset ℝ>0\mathbb{R}_{>0}, which is not a closed concave function.

We now turn to the behaviour of the sup-differential with respect to the basic operations. A first important property is the additivity.

Proposition 3.43.

For each i=1,…,li=1,\dots,l, let fif_{i} be a concave function and λi>0\lambda_{i}>0 a real number. Then

  1. (1)

    ∂(∑iλi​fi)⊃∑iλi​∂(fi)\partial\left(\sum_{i}\lambda_{i}f_{i}\right)\supset\sum_{i}\lambda_{i}\partial(f_{i});

  2. (2)

    if ri⁡(dom⁡(f1))∩⋯∩ri⁡(dom⁡(fl))≠∅\operatorname{ri}({\operatorname{dom}}(f_{1}))\cap\dots\cap\operatorname{ri}({\operatorname{dom}}(f_{l}))\neq\emptyset, then

    (3.44) ∂(∑iλi​fi)=∑iλi​∂(fi).\partial\bigg(\sum_{i}\lambda_{i}f_{i}\bigg)=\sum_{i}\lambda_{i}\partial(f_{i}).
Proof.

This is [Roc70, Theorem 23.8]. ∎

As in Remark 3.39, if f1,…,fkf_{1},\dots,f_{k} are piecewise affine, then (3.44) holds under the weaker hypothesis

dom⁡(f1)∩⋯∩dom⁡(fk)∩ri⁡(dom⁡(fk+1))∩⋯∩ri⁡(dom⁡(fl))≠∅.{\operatorname{dom}}(f_{1})\cap\dots\cap{\operatorname{dom}}(f_{k})\cap\operatorname{ri}({\operatorname{dom}}(f_{k+1}))\cap\dots\cap\operatorname{ri}({\operatorname{dom}}(f_{l}))\neq\emptyset.

The following result gives the behaviour of the sup-differential with respect to linear maps

Proposition 3.45.

Let H:Qℝ→NℝH\colon Q_{\mathbb{R}}\to N_{\mathbb{R}} be a linear map, u0∈Nℝu_{0}\in N_{\mathbb{R}} and A=H+u0A=H+u_{0} the associated affine map. Let ff be a concave function on NℝN_{\mathbb{R}}, then

  1. (1)

    ∂(A∗​f)​(v)⊃H∨​∂f⁡(A​v)\partial(A^{*}f)(v)\supset H^{\vee}\partial f(Av) for all v∈Qℝv\in Q_{\mathbb{R}};

  2. (2)

    if either ri⁡(dom⁡(f))∩im⁡(A)≠∅\operatorname{ri}({\operatorname{dom}}(f))\cap\operatorname{im}(A)\neq\emptyset or ff is piecewise affine and dom⁡(f)∩im⁡(A)≠∅{\operatorname{dom}}(f)\cap\operatorname{im}(A)\neq\emptyset, then for all v∈Qℝv\in Q_{\mathbb{R}} we have

    ∂(A∗​f)​(v)=H∨​∂f⁡(A​v).\partial(A^{*}f)(v)=H^{\vee}\partial f(Av).
Proof.

The linear case u0=0u_{0}=0 is [Roc70, Theorem 23.9]. The general case follows from the linear case and the commutativity of the sup-differential and the translation. ∎

We summarize the behaviour of direct and inverse images of affine maps with respect to the Legendre-Fenchel duality.

Proposition 3.46.

Let A:Qℝ→NℝA\colon Q_{\mathbb{R}}\to N_{\mathbb{R}} be an affine map defined as A=H+u0A=H+u_{0} for a linear map HH and a point u0∈Nℝu_{0}\in N_{\mathbb{R}}. Let ff be a concave function on NℝN_{\mathbb{R}} such that dom⁡(f)∩im⁡(A)≠∅{\operatorname{dom}}(f)\cap\operatorname{im}(A)\not=\emptyset and gg a concave function on QℝQ_{\mathbb{R}} such that stab⁡(g)∩im⁡(H∨)≠∅\operatorname{stab}(g)\cap\operatorname{im}(H^{\vee})\not=\emptyset. Then

  1. (1)

    stab⁡(A∗​g)=(H∨)−1​(stab⁡(g))\operatorname{stab}(A_{\ast}g)=(H^{\vee})^{-1}(\operatorname{stab}(g)) and

    (A∗​g)∨=(H∨)∗​(g∨)+u0;(A_{\ast}g)^{\vee}=(H^{\vee})^{\ast}(g^{\vee})+u_{0};
  2. (2)

    H∨​(stab⁡(f))⊂stab⁡(A∗​f)⊂H∨​(stab⁡(f))¯H^{\vee}(\operatorname{stab}(f))\subset\operatorname{stab}(A^{\ast}f)\subset{\overline{H^{\vee}(\operatorname{stab}(f))}} and

    (A∗​cl⁡(f))∨=cl⁡((H∨)∗​(f∨−u0));(A^{\ast}{\operatorname{cl}}(f))^{\vee}={\operatorname{cl}}((H^{\vee})_{\ast}(f^{\vee}-u_{0}));
  3. (3)

    if ri⁡(dom⁡(f))∩im⁡(A)≠∅\operatorname{ri}({\operatorname{dom}}(f))\cap\operatorname{im}(A)\not=\emptyset then stab⁡(A∗​f)=H∨​(stab⁡(f))\operatorname{stab}(A^{\ast}f)=H^{\vee}(\operatorname{stab}(f)) and, for all yy in this set,

    (A∗​f)∨​(y)=(H∨)∗​(f∨−u0)​(y)=maxx∈(H∨)−1​(y)⁡(f∨​(x)−⟨x,u0⟩).(A^{\ast}f)^{\vee}(y)=(H^{\vee})_{\ast}(f^{\vee}-u_{0})(y)=\max_{x\in(H^{\vee})^{-1}(y)}(f^{\vee}(x)-\langle x,u_{0}\rangle).

    Moreover, for y∈ri⁡(stab⁡(A∗​f))y\in\operatorname{ri}(\operatorname{stab}(A^{\ast}f)), a point x∈(H∨)−1​(y)x\in(H^{\vee})^{-1}(y) realizes this maximum if and only if x∈∂f⁡(A​v)x\in\partial f(Av) for a v∈Qℝv\in Q_{\mathbb{R}} such that y∈∂(A∗​f)​(v)y\in\partial(A^{*}f)(v).

Observe that the last assertion in the above proposition can be also expressed as

(3.47) (A∗​f)∨​(∂(A∗​f)​(v))=f∨​(∂f⁡(A​v))−⟨∂f⁡(A​v),u0⟩.(A^{\ast}f)^{\vee}(\partial(A^{*}f)(v))=f^{\vee}(\partial f(Av))-\langle\partial f(Av),u_{0}\rangle.
Proof.

By Proposition 3.40(3,4),

A∗​(f)=(H+u0)∗​(f)=H∗​(τ−u0​f),A∗​g=(H+u0)∗​g=τu0​(H∗​g).A^{\ast}(f)=(H+{u_{0}})^{\ast}(f)=H^{\ast}(\tau_{-u_{0}}f),\quad A_{\ast}g=(H+{u_{0}})_{\ast}g=\tau_{u_{0}}(H_{\ast}g).

Then, except for the last assertion, the result follows by combining this with the case when AA is a linear map, treated in [Roc70, Theorem 16.3].

To prove the last assertion of the proposition, we first note that the concave function

(f∨−u0)|(H∨)−1​(y)(f^{\vee}-u_{0})|_{(H^{\vee})^{-1}(y)}

attains its maximum at a point xx if and only if its sup-differential at xx contains 00. We fix a point x0x_{0} in (H∨)−1​(y)(H^{\vee})^{-1}(y) and we consider the affine inclusion

ι:Ker⁡(H∨)↪Mℝ,z↦z+x0.\iota\colon\operatorname{Ker}(H^{\vee})\hookrightarrow M_{\mathbb{R}},\quad z\mapsto z+x_{0}.

We denote by ι∨:Nℝ→Nℝ/im⁡(H)\iota^{\vee}\colon N_{\mathbb{R}}\to N_{\mathbb{R}}/\operatorname{im}(H) the dual of the linear part of ι\iota. Set F=ι∗​(f∨−u0)F=\iota^{*}(f^{\vee}-u_{0}), then for z∈Ker⁡(H∨)z\in\operatorname{Ker}(H^{\vee}), by Proposition 3.45, we have

∂F⁡(z)=ι∨​(∂f∨​(z+x0)−u0)\partial F(z)=\iota^{\vee}(\partial f^{\vee}(z+x_{0})-u_{0})

and so 0∈∂F⁡(z)0\in\partial F(z) if and only if ∂f∨​(z+x0)∩im⁡(A)≠∅\partial f^{\vee}(z+x_{0})\cap\operatorname{im}(A)\not=\emptyset. Hence x=z+x0x=z+x_{0} realizes the maximum if and only if x∈∂f⁡(A​v)x\in\partial f(Av) for some v∈Qℝv\in Q_{\mathbb{R}} such that y∈∂(A∗​f)​(v)y\in\partial(A^{*}f)(v), as stated. ∎

In particular, the operations of direct and inverse image of linear maps are dual to each other. In the notation of Proposition 3.46 and assuming for simplicity ri⁡(dom⁡(f))∩im⁡(H)≠∅\operatorname{ri}({\operatorname{dom}}(f))\cap\operatorname{im}(H)\not=\emptyset, we have

(H∗​g)∨=(H∨)∗​(g∨),(H∗​f)∨=(H∨)∗​(f∨),(H_{\ast}g)^{\vee}=(H^{\vee})^{\ast}(g^{\vee}),\quad(H^{\ast}f)^{\vee}=(H^{\vee})_{\ast}(f^{\vee}),

while the stability sets relate by stab⁡(H∗​g)=(H∨)−1​(stab⁡(g))\operatorname{stab}(H_{\ast}g)=(H^{\vee})^{-1}(\operatorname{stab}(g)) and stab⁡(H∗​f)=H∨​(stab⁡(f))\operatorname{stab}(H^{\ast}f)=H^{\vee}(\operatorname{stab}(f)).

The last concept we recall in this section is the notion of recession of a concave function.

Definition 3.48.

The recession function of a concave function f:Nℝ→ℝ¯f\colon N_{\mathbb{R}}\to{\underline{\mathbb{R}}}, denoted rec⁡(f)\operatorname{rec}(f), is the function

rec⁡(f):Nℝ⟶ℝ¯,u⟼infv∈dom⁡(f)(f⁡(u+v)−f⁡(v)).\operatorname{rec}(f)\colon N_{\mathbb{R}}\longrightarrow{\underline{\mathbb{R}}},\quad u\longmapsto\inf_{v\in{\operatorname{dom}}(f)}(f(u+v)-f(v)).

This is a concave conical function. If ff is closed, its recession function can be defined as the limit

(3.49) rec⁡(f)​(u)=limλ→∞λ−1​f​(v0+λ​u)\operatorname{rec}(f)(u)=\lim_{\lambda\to\infty}\lambda^{-1}f(v_{0}+\lambda u)

for any v0∈dom⁡(f)v_{0}\in{\operatorname{dom}}(f) [Roc70, Theorem 8.5].

It is clear from the definition that dom⁡(rec⁡(f))⊂rec⁡(dom⁡(f)){\operatorname{dom}}(\operatorname{rec}(f))\subset\operatorname{rec}({\operatorname{dom}}(f)). The equality does not hold in general, as can be seen by considering the concave function ℝ→ℝ\mathbb{R}\to\mathbb{R}, u↦−exp⁡(u)u\mapsto-\exp(u).

If ff is closed then the function rec⁡(f)\operatorname{rec}(f) is closed [Roc70, Theorem 8.5]. Hence it is natural to regard recession functions as support functions.

Proposition 3.50.

Let ff be a concave function. Then rec⁡(f∨)\operatorname{rec}(f^{\vee}) is the support function of dom⁡(f){\operatorname{dom}}(f). If ff is closed, then rec⁡(f)\operatorname{rec}(f) is the support function of stab⁡(f)\operatorname{stab}(f).

Proof.

This is [Roc70, Theorem 13.3]. ∎

Original mathematics by the credited authors. Source-backed reader collection; mathematical self-containment is not assessed.