ScalingStacks

3.5. The piecewise affine case [02M8]

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.5. The piecewise affine case

The Legendre-Fenchel duality for piecewise affine concave functions can be described in combinatorial terms. Moreover, some technical issues of the general theory disappear when dealing with piecewise affine concave functions on convex polyhedra and uniform limits of such functions.

Definition 3.58.

Let C⊂NℝC\subset N_{\mathbb{R}} be a convex polyhedron. A function f:C→ℝf\colon C\to\mathbb{R} is piecewise affine if there a finite cover of CC by closed subsets such that the restriction of ff to each of these subsets is an affine function. A concave function f:Nℝ→ℝ¯f\colon N_{\mathbb{R}}\to{\underline{\mathbb{R}}} is said to be piecewise affine if dom⁡(f){\operatorname{dom}}(f) is a convex polyhedron and the restriction f|dom⁡(f)f|_{{\operatorname{dom}}(f)} piecewise affine.

Lemma 3.59.

Let ff be a piecewise affine function defined on a convex polyhedron C⊂NℝC\subset N_{\mathbb{R}}. Then there exists a polyhedral complex Π\Pi in CC such that the restriction of ff to each polyhedron of Π\Pi is an affine function.

Proof.

This is an easy consequence of the max-min representation of piecewise affine functions in [Ovc02]. ∎

Definition 3.60.

Let CC be a convex polyhedron, Π\Pi a polyhedral complex in CC and f:C→ℝf\colon C\to\mathbb{R} a piecewise affine function. We say that Π\Pi and ff are compatible if ff is affine on each polyhedron of Π\Pi. Alternatively, we say that ff is a piecewise affine function on Π\Pi. If the function ff is concave, it is said to be strictly concave on Π\Pi if Π=Π⁡(f)\Pi=\Pi(f). The polyhedral complex Π\Pi is said to be regular if there exists a concave piecewise affine function ff such that Π=Π⁡(f)\Pi=\Pi(f).

As was the case for convex polyhedra, piecewise affine concave functions can be described in two dual ways, which we refer as the H-representation and the V-representation. For the H-representation, we consider a convex polyhedron

Λ=⋂1≤j≤k{u∈Nℝ∣⟨aj,u⟩+αj≥0}\Lambda=\bigcap_{1\leq j\leq k}\{u\in N_{\mathbb{R}}\mid\langle a_{j},u\rangle+\alpha_{j}\geq 0\}

as in (3.4) and a set of affine equations {(aj,αj)}k+1≤j≤l⊂Mℝ×ℝ\{(a_{j},\alpha_{j})\}_{k+1\leq j\leq l}\subset M_{\mathbb{R}}\times\mathbb{R}. We then define a concave function on NℝN_{\mathbb{R}} as

(3.61) f⁡(u)=mink+1≤j≤l⁡(⟨aj,u⟩+αj) for ​u∈Λf(u)=\min_{k+1\leq j\leq l}(\langle a_{j},u\rangle+\alpha_{j})\quad\text{ for }u\in\Lambda

and f⁡(u)=−∞f(u)=-\infty for u∉Λu\notin\Lambda. With this representation, the recession function of ff is given by

rec⁡(f)​(u)=mink+1≤j≤l⁡⟨aj,u⟩, for ​u∈rec⁡(Λ)\operatorname{rec}(f)(u)=\min_{k+1\leq j\leq l}\langle a_{j},u\rangle,\quad\text{ for }u\in\operatorname{rec}(\Lambda)

and rec⁡(f)​(u)=−∞\operatorname{rec}(f)(u)=-\infty for u∉rec⁡(Λ)u\notin\operatorname{rec}(\Lambda). In particular,

(3.62) dom(rec(f))=rec(dom(f)),stab(rec(f))=stab(f).{\operatorname{dom}}(\operatorname{rec}(f))=\operatorname{rec}({\operatorname{dom}}(f)),\quad\operatorname{stab}(\operatorname{rec}(f))=\operatorname{stab}(f).

For the V-representation, we consider a polyhedron

Λ′=cone⁡(b1,…,bk)+conv⁡(bk+1,…,bl)\Lambda^{\prime}=\operatorname{cone}(b_{1},\dots,b_{k})+\operatorname{conv}(b_{k+1},\dots,b_{l})

as in (3.5), a set of slopes {βj}1≤j≤k⊂ℝ\{\beta_{j}\}_{1\leq j\leq k}\subset\mathbb{R} and a set of values {βj}k+1≤j≤l⊂ℝ\{\beta_{j}\}_{k+1\leq j\leq l}\subset\mathbb{R}. We then define a concave function on MℝM_{\mathbb{R}} as

(3.63) g(u)=sup{∑j=1lλjβj|λj≥0,∑j=k+1lλj=1,∑j=1lλjbj=u}.g(u)=\sup\bigg\{\sum_{j=1}^{l}\lambda_{j}\beta_{j}\bigg|\ \lambda_{j}\geq 0,\ \sum_{j=k+1}^{l}\lambda_{j}=1,\ \sum_{j=1}^{l}\lambda_{j}b_{j}=u\bigg\}.

With this second representation, we obtain the recession function as

rec(g)(u)=sup{∑j=1kλjβj|λj≥0,∑j=1kλjbj=u}.\operatorname{rec}(g)(u)=\sup\bigg\{\sum_{j=1}^{k}\lambda_{j}\beta_{j}\bigg|\ \lambda_{j}\geq 0,\ \sum_{j=1}^{k}\lambda_{j}b_{j}=u\bigg\}.

As we have already mentioned, the Legendre-Fenchel duality of piecewise affine concave functions can be described in combinatorial terms.

Proposition 3.64.

Let Λ\Lambda be a polyhedron in NℝN_{\mathbb{R}} and ff a piecewise affine concave function with dom⁡(f)=Λ{\operatorname{dom}}(f)=\Lambda given as

Λ\displaystyle\Lambda =⋂1≤j≤k{u∈Nℝ∣⟨aj,u⟩+αj≥0},\displaystyle=\bigcap_{1\leq j\leq k}\{u\in N_{\mathbb{R}}\mid\langle a_{j},u\rangle+\alpha_{j}\geq 0\},
f⁡(u)\displaystyle f(u) =mink+1≤j≤l⁡(⟨aj,u⟩+αj)for ​u∈Λ\displaystyle=\min_{k+1\leq j\leq l}(\langle a_{j},u\rangle+\alpha_{j})\quad\text{for }u\in\Lambda

with aj∈Mℝa_{j}\in M_{\mathbb{R}} and αj∈ℝ\alpha_{j}\in\mathbb{R}. Then

stab⁡(f)\displaystyle\operatorname{stab}(f) =cone⁡(a1,…,ak)+conv⁡(ak+1,…,al),\displaystyle=\operatorname{cone}(a_{1},\dots,a_{k})+\operatorname{conv}(a_{k+1},\dots,a_{l}),
f∨​(x)\displaystyle f^{\vee}(x) =sup{∑j=1l−λjαj|λj≥0,∑j=k+1lλj=1,∑j=1lλjaj=x} for x∈stab(f).\displaystyle=\sup\bigg\{\sum_{j=1}^{l}-\lambda_{j}\alpha_{j}\bigg|\ \lambda_{j}\geq 0,\sum_{j=k+1}^{l}\lambda_{j}=1,\ \sum_{j=1}^{l}\lambda_{j}a_{j}=x\bigg\}\ \text{ for }x\in\operatorname{stab}(f).
Proof.

This is proved in [Roc70, pp. 172-174]. ∎

Example 3.65.

Let Λ\Lambda be a convex polyhedron in NℝN_{\mathbb{R}}. Then both the indicator function ιΛ\iota_{\Lambda} and the support function ΨΛ\Psi_{\Lambda} are concave and piecewise affine. We have ΨΛ∨=ιΛ\Psi_{\Lambda}^{\vee}=\iota_{\Lambda}. In particular, if we fix an isomorphism Nℝ≃ℝnN_{\mathbb{R}}\simeq\mathbb{R}^{n}, the function

ΨΔn:Nℝ⟶ℝ,(u1,…,un)⟼min⁡{0,u1,…,un}\Psi_{\Delta^{n}}\colon N_{\mathbb{R}}\longrightarrow\mathbb{R},\quad(u_{1},\dots,u_{n})\longmapsto\min\{0,u_{1},\dots,u_{n}\}

is the support function of the standard simplex Δn=conv⁡(𝟎,e1∨,…,en∨)⊂Mℝ\Delta^{n}=\operatorname{conv}(\boldsymbol{0},e^{\vee}_{1},\dots,e^{\vee}_{n})\subset M_{\mathbb{R}}, where {e1,…,en}\{e_{1},\dots,e_{n}\} is the standard basis of ℝn\mathbb{R}^{n} and {e1∨,…,en∨}\{e_{1}^{\vee},\dots,e_{n}^{\vee}\} is the dual basis. Hence, stab⁡(ΨΔn)=Δn\operatorname{stab}(\Psi_{\Delta^{n}})=\Delta^{n} and ΨΔn∨=ιΔn\Psi_{\Delta^{n}}^{\vee}=\iota_{\Delta^{n}}.

Let Λ\Lambda be a polyhedron in NℝN_{\mathbb{R}} and ff a piecewise affine concave function with dom⁡(f)=Λ{\operatorname{dom}}(f)=\Lambda. Then dom⁡(∂f)=Λ{\operatorname{dom}}(\partial f)=\Lambda and Π⁡(f)\Pi(f) and Π⁡(f∨)\Pi(f^{\vee}) are convex decompositions of Λ\Lambda and of Λ′:=stab⁡(f)\Lambda^{\prime}:=\operatorname{stab}(f) respectively. By Theorem 3.33, the Legendre-Fenchel correspondence

ℒ​f:Π⁡(f)⟶Π⁡(f∨){\mathcal{L}}f\colon\Pi(f)\longrightarrow\Pi(f^{\vee})

is a duality in the sense of Definition 3.32. However in the polyhedral case, these decompositions are dual in a stronger sense. We need to introduce some more definitions before we can properly state this duality.

Definition 3.66.

Let Λ\Lambda be a polyhedron and KK a face of Λ\Lambda. The angle of Λ\Lambda at KK is defined as

∠(K,Λ)={t(u−v)∣u∈Λ,v∈K,t≥0}.\angle(K,\Lambda)=\{t(u-v)\mid u\in\Lambda,v\in K,t\geq 0\}.

It is a polyhedral cone.

Definition 3.67.

The dual of a convex cone σ⊂Nℝ\sigma\subset N_{\mathbb{R}} is defined as

σ∨={x∈Mℝ∣⟨x,u⟩≥0​ for all ​u∈σ}.\sigma^{\vee}=\{x\in M_{\mathbb{R}}\mid\langle x,u\rangle\geq 0\text{ for all }u\in\sigma\}.

This is a convex closed cone.

If σ\sigma is a convex closed cone, then σ∨⁣∨=σ\sigma^{\vee\vee}=\sigma. For a piecewise affine concave function ff on NℝN_{\mathbb{R}}, by Proposition 3.64 we have

rec⁡(dom⁡(f))∨=rec⁡(stab⁡(f)).\operatorname{rec}({\operatorname{dom}}(f))^{\vee}=\operatorname{rec}(\operatorname{stab}(f)).
Definition 3.68.

Let C,C′C,C^{\prime} be convex polyhedra in NℝN_{\mathbb{R}} and MℝM_{\mathbb{R}}, respectively, and Π,Π′\Pi,\Pi^{\prime} polyhedral complexes in CC and C′C^{\prime}, respectively. We say that Π\Pi and Π′\Pi^{\prime} are dual polyhedral complexes if there is a bijective map Π→Π′,Λ↦Λ∗\Pi\to\Pi^{\prime},\Lambda\mapsto\Lambda^{*} such that

  1. (1)

    for all Λ,K∈Π\Lambda,K\in\Pi, the inclusion K⊂ΛK\subset\Lambda hols if and only if K∗⊃Λ∗K^{*}\supset\Lambda^{*};

  2. (2)

    for all Λ,K∈Π\Lambda,K\in\Pi, if K⊂ΛK\subset\Lambda, then ∠⁡(Λ∗,K∗)=∠​(K,Λ)∨\angle(\Lambda^{*},K^{*})=\angle(K,\Lambda)^{\vee}.

For Λ∈Π\Lambda\in\Pi, the angle ∠⁡(Λ,Λ)\angle(\Lambda,\Lambda) is the linear subspace generated by differences of points in Λ\Lambda. Condition (2) above implies that ∠⁡(Λ,Λ)\angle(\Lambda,\Lambda) and ∠⁡(Λ∗,Λ∗)\angle(\Lambda^{*},\Lambda^{*}) are orthogonal. In particular, dim(Λ)+dim(Λ∗)=n\dim(\Lambda)+\dim(\Lambda^{*})=n.

Proposition 3.69.

Let ff be a piecewise affine concave function with Λ=dom⁡(f)\Lambda={\operatorname{dom}}(f) and Λ′=stab⁡(f)\Lambda^{\prime}=\operatorname{stab}(f). Then Π⁡(f)\Pi(f) and Π⁡(f∨)\Pi(f^{\vee}) are polyhedral complexes in Λ\Lambda and Λ′\Lambda^{\prime} respectively. Moreover, they are dual of each other. In particular, the vertices of Π⁡(f)\Pi(f) are in bijection with the polyhedra of Π⁡(f∨)\Pi(f^{\vee}) of maximal dimension.

Proof.

This is proved in [PR04, Proposition 1]. ∎

Example 3.70.

Consider the standard simplex Δn\Delta^{n} of Example 3.65. Its indicator function induces the standard polyhedral complex in Δn\Delta^{n} consisting of the collection of its faces. The dual of ιΔn\iota_{\Delta^{n}}, the support function ΨΔn\Psi_{\Delta^{n}}, induces a fan ΣΔn:=Π⁡(ΨΔn)\Sigma_{\Delta^{n}}:=\Pi(\Psi_{\Delta^{n}}) of NℝN_{\mathbb{R}}. The duality between these polyhedral complexes can be made explicit as

Π⁡(ιΔn)⟶ΣΔn,F⟼∠​(F,Δn)∨.\Pi(\iota_{\Delta^{n}})\longrightarrow\Sigma_{\Delta^{n}},\quad F\longmapsto\angle(F,\Delta^{n})^{\vee}.
Example 3.71.

The previous example can be generalized to an arbitrary polytope Δ⊂Mℝ\Delta\subset M_{\mathbb{R}} . The indicator function ιΔ\iota_{\Delta} induces the standard decomposition of Δ\Delta into its faces and dually, the support function ΨΔ\Psi_{\Delta} induces a polyhedral complex ΣΔ:=Π⁡(ΨΔ)\Sigma_{\Delta}:=\Pi(\Psi_{\Delta}) made of cones. If Δ\Delta is of maximal dimension, then ΣΔ\Sigma_{\Delta} is a fan.

The faces of Δ\Delta are in one-to-one correspondence with the cones of ΣΔ\Sigma_{\Delta} through the Legendre-Fenchel correspondence. For a face FF of Δ\Delta, its corresponding cone is

σF:=F∗={u∈Nℝ∣⟨u,x−y⟩≥0 for all x∈Δ,y∈F}.\sigma_{F}:=F^{*}=\{u\in N_{\mathbb{R}}\mid\langle u,x-y\rangle\geq 0\mbox{ for all }x\in\Delta,y\in F\}.

Reciprocally, to each cone σ\sigma corresponds a face of Δ\Delta of complementary dimension

Fσ:=σ∗={x∈Δ∣⟨x,u⟩=ΨΔ(u) for all u∈σ}.F_{\sigma}:=\sigma^{*}=\{x\in\Delta\mid\langle x,u\rangle=\Psi_{\Delta}(u)\mbox{ for all }u\in\sigma\}.

On a cone σ∈Σ\sigma\in\Sigma, the function ΨΔ\Psi_{\Delta} is defined by any vector mσm_{\sigma} in the affine space aff⁡(Fσ)\operatorname{aff}(F_{\sigma}). The cone σ\sigma is normal to FσF_{\sigma}.

For piecewise affine concave functions, the operations of taking the recession function and the associated polyhedral convex commute with each other.

Proposition 3.72.

Let ff be a piecewise affine concave function on NℝN_{\mathbb{R}}. Then

Π⁡(rec⁡(f))=rec⁡(Π⁡(f)).\Pi(\operatorname{rec}(f))=\operatorname{rec}(\Pi(f)).
Proof.

Let Pf​(u,x)=f⁡(u)+f∨​(x)−⟨u,x⟩P_{f}(u,x)=f(u)+f^{\vee}(x)-\left<u,x\right> be the function introduced in (3.20). For each x∈stab⁡(f)x\in\operatorname{stab}(f) write Pf,x​(u)=P⁡(u,x)P_{f,x}(u)=P(u,x). Let CxC_{x} be as in Definition 3.23. By Lemma 3.24,

Cx={u∈dom⁡(f)∣Pf,x​(u)=0}.C_{x}=\{u\in{\operatorname{dom}}(f)\mid P_{f,x}(u)=0\}.

Write P′​(v)=rec⁡(f)​(v)−⟨u,x⟩P^{\prime}(v)=\operatorname{rec}(f)(v)-\left<u,x\right>. Then P′=rec⁡(Pf,x)P^{\prime}=\operatorname{rec}(P_{f,x}).

We claim that, for each x∈stab⁡(f)x\in\operatorname{stab}(f),

rec⁡(Cx)={v∈dom⁡(rec⁡(f))∣P′​(v)=0}.\operatorname{rec}(C_{x})=\{v\in{\operatorname{dom}}(\operatorname{rec}(f))\mid P^{\prime}(v)=0\}.

Let v∈rec⁡(Cx)v\in\operatorname{rec}(C_{x}). Clearly v∈dom⁡(rec⁡(f))v\in{\operatorname{dom}}(\operatorname{rec}(f)) and, since x∈stab⁡(f)x\in\operatorname{stab}(f), the set CxC_{x} is non-empty. Let u0∈Cxu_{0}\in C_{x}. Then, for each λ>0\lambda>0, u0+λ​v∈Cxu_{0}+\lambda v\in C_{x}. Therefore,

P′​(v)=limλ→∞Pf,x​(u0+λ​v)−Pf,x​(u0)λ=0.P^{\prime}(v)=\lim_{\lambda\to\infty}\frac{P_{f,x}(u_{0}+\lambda v)-P_{f,x}(u_{0})}{\lambda}=0.

Conversely, let v∈dom⁡(rec⁡(f))v\in{\operatorname{dom}}(\operatorname{rec}(f)) satisfying P′​(v)=0P^{\prime}(v)=0 and u∈Cxu\in C_{x}. On the one hand, by the properties of the function PfP_{f}, we have Pf,x​(u+v)≤0P_{f,x}(u+v)\leq 0. On the other hand, since P′=rec⁡(Pf,x)P^{\prime}=\operatorname{rec}(P_{f,x}),

Pf,x​(u+v)−Pf,x​(u)≥P′​(v)=0.P_{f,x}(u+v)-P_{f,x}(u)\geq P^{\prime}(v)=0.

Thus Pf,x​(u+v)≥Pf,x​(u)=0P_{f,x}(u+v)\geq P_{f,x}(u)=0 and finally Pf,x​(u+v)=0P_{f,x}(u+v)=0. This implies that, if u∈Cxu\in C_{x} then u+v∈Cxu+v\in C_{x}, showing v∈rec⁡(Cx)v\in\operatorname{rec}(C_{x}). Hence the claim is proved.

By definition Π⁡(f)={Cx}x∈stab⁡(f)\Pi(f)=\{C_{x}\}_{x\in\operatorname{stab}(f)}. Hence rec⁡(Π⁡(f))={rec⁡(Cx)}x∈stab⁡(f)\operatorname{rec}(\Pi(f))=\{\operatorname{rec}(C_{x})\}_{x\in\operatorname{stab}(f)}. For each x∈stab⁡(rec⁡(f))x\in\operatorname{stab}(\operatorname{rec}(f)), write

Cx′={v∈dom⁡(rec⁡(f))∣P′​(v)=0}.C^{\prime}_{x}=\{v\in{\operatorname{dom}}(\operatorname{rec}(f))\mid P^{\prime}(v)=0\}.

Then Π⁡(rec⁡(f))={Cx′}x∈stab⁡(rec⁡(f))\Pi(\operatorname{rec}(f))=\{C^{\prime}_{x}\}_{x\in\operatorname{stab}(\operatorname{rec}(f))}. The result follows from the previous claim and the fact that stab⁡(f)=stab⁡(rec⁡(f))\operatorname{stab}(f)=\operatorname{stab}(\operatorname{rec}(f)) by (3.62). ∎

Now we want to study the compatibility of Legendre-Fenchel duality and integral and rational structures. Let N≃ℤnN\simeq\mathbb{Z}^{n} be a lattice of rank nn such that Nℝ=N⊗ℝN_{\mathbb{R}}=N\otimes\mathbb{R}. Set M=N∨=Hom⁡(N,ℤ)M=N^{\vee}=\operatorname{Hom}(N,\mathbb{Z}) for its dual lattice, so Mℝ=M⊗ℝM_{\mathbb{R}}=M\otimes\mathbb{R}. We also set Nℚ=N⊗ℚN_{\mathbb{Q}}=N\otimes\mathbb{Q} and Mℚ=M⊗ℚM_{\mathbb{Q}}=M\otimes\mathbb{Q}.

Definition 3.73.

A piecewise affine concave function ff on NℝN_{\mathbb{R}} is an H-lattice (respectively, a V-lattice) concave function if it has an H-representation (respectively, a V-representation) with integral coefficients. We say that ff is a rational piecewise affine concave function if it has an H-representation (or equivalently, a V-representation) with rational coefficients.

Observe that the domain of a V-lattice concave function is a lattice polyhedron, whereas the domain of an H-lattice concave function is a rational polyhedron.

Remark 3.74.

The notion of H-lattice concave functions defined on the whole NℝN_{\mathbb{R}} coincides with the notion of tropical Laurent polynomials over the integers, that is, the elements of the group semi-algebra ℤtrop​[N]\mathbb{Z}_{\text{trop}}[N], where the arithmetic operations of the base semi-ring ℤtrop=(ℤ,⊕,⊙)\mathbb{Z}_{\text{trop}}=(\mathbb{Z},\oplus,\odot) are defined as x⊕y=min⁡(x,y)x\oplus y=\min(x,y) and x⊙y=x+yx\odot y=x+y.

Proposition 3.75.

Let ff be a piecewise affine concave function on NℝN_{\mathbb{R}}.

  1. (1)

    ff is an H-lattice concave function (respectively, a rational piecewise affine concave function) if and only if f∨f^{\vee} is a V-lattice concave function (respectively, a rational piecewise affine concave function).

  2. (2)

    rec⁡(f)\operatorname{rec}(f) is an H-lattice concave function if and only if stab⁡(f)\operatorname{stab}(f) is a lattice polyhedron.

Proof.

This follows easily from Proposition 3.64. ∎

Example 3.76.

If Δ\Delta is a lattice polytope, its indicator function is a V-lattice function, its support function ΨΔ\Psi_{\Delta} is an H-lattice function and, when Δ\Delta has maximal dimension, the fan ΣΔ\Sigma_{\Delta} is a rational fan. In particular, if the isomorphism N≃ℝnN\simeq\mathbb{R}^{n} of Example 3.65 is given by the choice of an integral basis e1,…,ene_{1},\dots,e_{n} of NN, then Δn\Delta^{n} is a lattice polytope, the function ΨΔn\Psi_{\Delta^{n}} is an H-lattice concave function and ΣΔn\Sigma_{\Delta^{n}} is a rational fan. If we write e0=−∑i=1neie_{0}=-\sum_{i=1}^{n}e_{i}, this is the fan generated by the vectors e0,e1,…,ene_{0},e_{1},\dots,e_{n} in the sense that each cone of ΣΔn\Sigma_{\Delta^{n}} is the cone generated by a strict subset of the above set of vectors. Figure 3 illustrates the case n=2n=2.

Δ 2 σ 1 σ 0 σ 2 ↦ ( u 1 , u 2 ) 0 ↦ ( u 1 , u 2 ) - u 2 ↦ ( u 1 , u 2 ) - u 1
Figure 3. The standard simplex Δ2\Delta^{2}, its associated fan and support function

Let Λ\Lambda and Λ′\Lambda^{\prime} be polyhedra in NℝN_{\mathbb{R}} and in MℝM_{\mathbb{R}}, respectively. We set 𝒫⁡(Λ,Λ′)\mathscr{P}(\Lambda,\Lambda^{\prime}) for the space of piecewise affine concave functions with effective domain Λ\Lambda and stability set Λ′\Lambda^{\prime}. We also set 𝒫¯​(Λ,Λ′){\overline{\mathscr{P}}}(\Lambda,\Lambda^{\prime}) for the closure of this space with respect to uniform convergence. We set

𝒫⁡(Λ)=⋃Λ′𝒫⁡(Λ,Λ′),𝒫¯​(Λ)=⋃Λ′𝒫¯​(Λ,Λ′)\mathscr{P}(\Lambda)=\bigcup_{\Lambda^{\prime}}\mathscr{P}(\Lambda,\Lambda^{\prime}),\quad{\overline{\mathscr{P}}}(\Lambda)=\bigcup_{\Lambda^{\prime}}{\overline{\mathscr{P}}}(\Lambda,\Lambda^{\prime})

for the space of piecewise affine concave functions with effective domain Λ\Lambda and for its closure with respect to uniform convergence, respectively. We also set

𝒫=⋃Λ,Λ′𝒫⁡(Λ,Λ′),𝒫¯=⋃Λ,Λ′𝒫¯​(Λ,Λ′).\mathscr{P}=\bigcup_{\Lambda,\Lambda^{\prime}}\mathscr{P}(\Lambda,\Lambda^{\prime}),\quad{\overline{\mathscr{P}}}=\bigcup_{\Lambda,\Lambda^{\prime}}{\overline{\mathscr{P}}}(\Lambda,\Lambda^{\prime}).

When we need to specify the vector space NℝN_{\mathbb{R}} we will denote it as a subindex as in 𝒫Nℝ\mathscr{P}_{N_{\mathbb{R}}} or 𝒫¯Nℝ{\overline{\mathscr{P}}}_{N_{\mathbb{R}}}.

The following propositions contain the basic properties of the Legendre-Fenchel duality acting on 𝒫¯{\overline{\mathscr{P}}}. The elements in 𝒫¯{\overline{\mathscr{P}}} are continuous functions on polyhedra. In particular, they are closed concave functions. Observe that when working with uniform limits of piecewise affine concave functions, the technical issues in §3.2 disappear.

Proposition 3.77.

The concave piecewise affine functions and their uniform limits satisfy the following properties.

  1. (1)

    Let f∈𝒫¯Nℝf\in{\overline{\mathscr{P}}}_{N_{\mathbb{R}}}. Then f∨⁣∨=ff^{\vee\vee}=f.

  2. (2)

    If f∈𝒫⁡(Λ,Λ′)f\in\mathscr{P}(\Lambda,\Lambda^{\prime}) (respectively f∈𝒫¯​(Λ,Λ′)f\in{\overline{\mathscr{P}}}(\Lambda,\Lambda^{\prime})) then f∨∈𝒫⁡(Λ′,Λ)f^{\vee}\in\mathscr{P}(\Lambda^{\prime},\Lambda) (respectively f∨∈𝒫¯​(Λ′,Λ)f^{\vee}\in{\overline{\mathscr{P}}}(\Lambda^{\prime},\Lambda)).

  3. (3)

    If f∈𝒫¯​(Λ)f\in{\overline{\mathscr{P}}}(\Lambda) then dom⁡(rec⁡(f))=rec⁡(Λ){\operatorname{dom}}(\operatorname{rec}(f))=\operatorname{rec}(\Lambda).

  4. (4)

    Let fi∈𝒫⁡(Λi,Λi′)f_{i}\in\mathscr{P}(\Lambda_{i},\Lambda_{i}^{\prime}) (respectively fi∈𝒫¯​(Λi,Λi′)f_{i}\in{\overline{\mathscr{P}}}(\Lambda_{i},\Lambda_{i}^{\prime})), i=1,2i=1,2, with Λ1∩Λ2≠∅\Lambda_{1}\cap\Lambda_{2}\not=\emptyset. Then f1+f2∈𝒫⁡(Λ1∩Λ2,Λ1′+Λ2′)f_{1}+f_{2}\in\mathscr{P}(\Lambda_{1}\cap\Lambda_{2},\Lambda_{1}^{\prime}+\Lambda_{2}^{\prime}) (respectively f1+f2∈𝒫¯​(Λ1∩Λ2,Λ1′+Λ2′)f_{1}+f_{2}\in{\overline{\mathscr{P}}}(\Lambda_{1}\cap\Lambda_{2},\Lambda_{1}^{\prime}+\Lambda_{2}^{\prime})) and (f1+f2)∨=f1∨⊞f2∨(f_{1}+f_{2})^{\vee}=f_{1}^{\vee}\boxplus f_{2}^{\vee}.

  5. (5)

    Let fi∈𝒫⁡(Λi,Λi′)f_{i}\in\mathscr{P}(\Lambda_{i},\Lambda_{i}^{\prime}) (respectively fi∈𝒫¯​(Λi,Λi′)f_{i}\in{\overline{\mathscr{P}}}(\Lambda_{i},\Lambda_{i}^{\prime})), i=1,2i=1,2, with Λ1′∩Λ2′≠∅\Lambda_{1}^{\prime}\cap\Lambda_{2}^{\prime}\not=\emptyset. Then f1⊞f2∈𝒫⁡(Λ1+Λ2,Λ1′∩Λ2′)f_{1}\boxplus f_{2}\in\mathscr{P}(\Lambda_{1}+\Lambda_{2},\Lambda_{1}^{\prime}\cap\Lambda_{2}^{\prime}) (respectively f1+f2∈𝒫¯​(Λ1+Λ2,Λ1′∩Λ2′)f_{1}+f_{2}\in{\overline{\mathscr{P}}}(\Lambda_{1}+\Lambda_{2},\Lambda_{1}^{\prime}\cap\Lambda_{2}^{\prime})) and (f1⊞f2)∨=f1∨+f2∨(f_{1}\boxplus f_{2})^{\vee}=f_{1}^{\vee}+f_{2}^{\vee}.

  6. (6)

    Let (fi)i≥1⊂𝒫¯(f_{i})_{i\geq 1}\subset{\overline{\mathscr{P}}} be a sequence converging uniformly to a function ff. Then f∈𝒫¯f\in{\overline{\mathscr{P}}}.

Proof.

All the statements follow, either directly from the definition, or propositions 3.64 and 3.18. ∎

Proposition 3.78.

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 f∈𝒫Nℝf\in\mathscr{P}_{N_{\mathbb{R}}} (respectively f∈𝒫¯Nℝf\in{\overline{\mathscr{P}}}_{N_{\mathbb{R}}}) with dom⁡(f)∩im⁡(A)≠∅{\operatorname{dom}}(f)\cap\operatorname{im}(A)\not=\emptyset and g∈𝒫Qℝg\in\mathscr{P}_{Q_{\mathbb{R}}} (respectively g∈𝒫¯Qℝg\in{\overline{\mathscr{P}}}_{Q_{\mathbb{R}}}) such that stab⁡(g)∩im⁡(H∨)≠∅\operatorname{stab}(g)\cap\operatorname{im}(H^{\vee})\not=\emptyset. Then A∗​f∈𝒫QℝA^{\ast}f\in{\mathscr{P}}_{Q_{\mathbb{R}}} (respectively A∗​f∈𝒫¯QℝA^{*}f\in{\overline{\mathscr{P}}}_{Q_{\mathbb{R}}}) and A∗​g∈𝒫NℝA_{*}g\in\mathscr{P}_{N_{\mathbb{R}}} (respectively A∗​g∈𝒫¯NℝA_{*}g\in{\overline{\mathscr{P}}}_{N_{\mathbb{R}}}). Moreover,

  1. (1)

    stab⁡(A∗​f)=H∨​(stab⁡(f))\operatorname{stab}(A^{\ast}f)=H^{\vee}(\operatorname{stab}(f)), (A∗​f)∨=(H∨)∗​(f∨−u0)(A^{\ast}f)^{\vee}=(H^{\vee})_{\ast}(f^{\vee}-u_{0}) and, for all y∈stab⁡(A∗​f)y\in\operatorname{stab}(A^{\ast}f),

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

    stab⁡(A∗​g)=(H∨)−1​(stab⁡(g))\operatorname{stab}(A_{\ast}g)=(H^{\vee})^{-1}(\operatorname{stab}(g)), (A∗​g)∨=(H∨)∗​(g∨)+u0(A_{\ast}g)^{\vee}=(H^{\vee})^{\ast}(g^{\vee})+u_{0} and, for all u∈dom⁡(A∗​g)u\in{\operatorname{dom}}(A_{*}g),

    A∗​g​(u)=maxv∈A−1​(u)⁡g⁡(v).A_{*}g(u)=\max_{v\in A^{-1}(u)}g(v).
Proof.

These statements follow either from Proposition 3.46 or from [Roc70, Corollary 19.3.1]. ∎

We will be concerned mainly with functions in 𝒫¯{\overline{\mathscr{P}}} whose effective domain is either a polytope or the whole space NℝN_{\mathbb{R}}. These are the kind of functions that arise when considering proper toric varieties. The functions in 𝒫⁡(Nℝ)\mathscr{P}(N_{\mathbb{R}}) can be realized as the inverse image of the support function of the standard simplex, while the functions of 𝒫⁡(Δ)\mathscr{P}(\Delta) can be realized as direct images of the indicator function of the standard simplex.

Lemma 3.79.

Let f∈𝒫⁡(Nℝ)f\in\mathscr{P}(N_{\mathbb{R}}) and let f⁡(u)=min0≤i≤r⁡(ai​(u)+αi)f(u)=\min_{0\leq i\leq r}(a_{i}(u)+\alpha_{i}) be an H-representation of ff. Write 𝛂=(αi−α0)i=1,…,r\boldsymbol{\alpha}=(\alpha_{i}-\alpha_{0})_{i=1,\dots,r}, and consider the linear map H:Nℝ→ℝrH\colon N_{\mathbb{R}}\to\mathbb{R}^{r} given by H⁡(u)=(ai​(u)−a0​(u))i=1,…,rH(u)=(a_{i}(u)-a_{0}(u))_{i=1,\dots,r} and the affine map A=H+𝛂.A=H+\boldsymbol{\alpha}. Then

  1. (1)

    f=A∗​ΨΔr+a0+α0;f=A^{\ast}\Psi_{\Delta^{r}}+a_{0}+\alpha_{0};

  2. (2)

    f∨=τa0​(H∨)∗​(ιΔr−𝜶)−α0.f^{\vee}=\tau_{a_{0}}(H^{\vee})_{\ast}(\iota_{\Delta^{r}}-\boldsymbol{\alpha})-\alpha_{0}.

This second function can be alternatively described as the function which parameterizes the upper envelope of the extended polytope

conv⁡((a1,−α1),…,(al,−αl))⊂Mℝ×ℝ.\operatorname{conv}((a_{1},-\alpha_{1}),\dots,(a_{l},-\alpha_{l}))\subset M_{\mathbb{R}}\times\mathbb{R}.
Proof.

Statement (1) follows from the explicit description of ΨΔr\Psi_{\Delta^{r}} in Example 3.70. Statement (2) follows from Proposition 3.78. The last statement is a consequence of Proposition 3.64. ∎

The next proposition characterizes the elements of 𝒫¯​(Nℝ){\overline{\mathscr{P}}}(N_{\mathbb{R}}) and 𝒫¯​(Δ){\overline{\mathscr{P}}}(\Delta) for a polytope Δ\Delta.

Proposition 3.80.

Let Δ\Delta be a convex polytope of MℝM_{\mathbb{R}}.

  1. (1)

    The space 𝒫¯​(Δ,Nℝ){\overline{\mathscr{P}}}(\Delta,N_{\mathbb{R}}) agrees with the space of all continuous concave functions on Δ\Delta.

  2. (2)

    A concave function ff belongs to 𝒫¯​(Nℝ,Δ){\overline{\mathscr{P}}}(N_{\mathbb{R}},\Delta) if and only if dom⁡(f)=Nℝ{\operatorname{dom}}(f)=N_{\mathbb{R}} and |f−ΨΔ||f-\Psi_{\Delta}| is bounded.

Proof.

We start by proving (1). By the properties of uniform convergence, it is clear that any element of 𝒫¯​(Δ,Nℝ){\overline{\mathscr{P}}}(\Delta,N_{\mathbb{R}}) is concave and continuous. Conversely, a continuous function ff on Δ\Delta is uniformly continuous because Δ\Delta is compact. Therefore, given ε>0\varepsilon>0 there is a δ>0\delta>0 such that |f⁡(u)−f⁡(v)|<ε|f(u)-f(v)|<\varepsilon for all u,v∈Δu,v\in\Delta such that ‖u−v‖<δ\|u-v\|<\delta. By compactness, we can find a triangulation Δ=⋃iΔi\Delta=\bigcup_{i}\Delta_{i} with diam⁡(Δi)<δ\operatorname{diam}(\Delta_{i})<\delta. Let {bj}j\{b_{j}\}_{j} be the vertices of this triangulation and consider the function g∈𝒫⁡(Δ,Nℝ)g\in\mathscr{P}(\Delta,N_{\mathbb{R}}) defined as

g(u)=sup{∑j=1lλjf(bj)|λj≥0,∑jλj=1∑jλjaj=x}.g(u)=\sup\bigg\{\sum_{j=1}^{l}\lambda_{j}f(b_{j})\bigg|\ \lambda_{j}\geq 0,\sum_{j}\lambda_{j}=1\sum_{j}\lambda_{j}a_{j}=x\bigg\}.

For u∈Δu\in\Delta, let bj0,…,bjnb_{j_{0}},\dots,b_{j_{n}} denote the vertices of an element of the triangulation containing uu. We write u=λj0​uj0+⋯+λjn​ujnu=\lambda_{j_{0}}u_{j_{0}}+\dots+\lambda_{j_{n}}u_{j_{n}} for some λji≥0\lambda_{j_{i}}\geq 0 and λj0+⋯+λjn=1\lambda_{j_{0}}+\dots+\lambda_{j_{n}}=1. By concavity, we have

f⁡(u)≥g⁡(u)≥∑k=0nλjk​f​(ujk)≥f⁡(u)−ε,f(u)\geq g(u)\geq\sum_{k=0}^{n}\lambda_{j_{k}}f(u_{j_{k}})\geq f(u)-\varepsilon,

which shows that any continuous function on Δ\Delta can be arbitrarily approximated by elements of 𝒫⁡(Δ,Nℝ)\mathscr{P}(\Delta,N_{\mathbb{R}}).

We now prove (2). Let f∈𝒫¯​(Nℝ,Δ)f\in{\overline{\mathscr{P}}}(N_{\mathbb{R}},\Delta). By definition, for each ε>0\varepsilon>0 we can find a function g∈𝒫⁡(Nℝ,Δ)g\in\mathscr{P}(N_{\mathbb{R}},\Delta) with sup|f−g|≤ε\sup|f-g|\leq\varepsilon. In particular, |f−g||f-g| is bounded. Furthermore, rec⁡(g)=ΨΔ\operatorname{rec}(g)=\Psi_{\Delta} and |g−rec⁡(g)||g-\operatorname{rec}(g)| is bounded because g∈𝒫⁡(Nℝ)g\in\mathscr{P}(N_{\mathbb{R}}). Hence dom⁡(f)=dom⁡(g)=Nℝ{\operatorname{dom}}(f)={\operatorname{dom}}(g)=N_{\mathbb{R}} and |f−ΨΔ||f-\Psi_{\Delta}| is bounded.

Conversely, let ff be a concave function such that dom⁡(f)=Nℝ{\operatorname{dom}}(f)=N_{\mathbb{R}} and |f−ΨΔ||f-\Psi_{\Delta}| is bounded. Then stab⁡(f)=stab⁡(ΨΔ)=Δ\operatorname{stab}(f)=\operatorname{stab}(\Psi_{\Delta})=\Delta and f∨f^{\vee} is a continuous concave function on Δ\Delta. Hence we can apply (1) to f∨f^{\vee} to obtain functions gi∈𝒫⁡(Δ,Nℝ)g_{i}\in\mathscr{P}(\Delta,N_{\mathbb{R}}) approaching f∨f^{\vee} uniformly. We conclude that the functions gi∨∈𝒫⁡(Nℝ,Δ)g_{i}^{\vee}\in\mathscr{P}(N_{\mathbb{R}},\Delta) approach ff uniformly and so f∈𝒫¯​(Nℝ,Δ)f\in{\overline{\mathscr{P}}}(N_{\mathbb{R}},\Delta). ∎

Proposition 3.81.

Let Δ\Delta be a lattice polytope of MℝM_{\mathbb{R}}. Then the subset of rational piecewise affine concave functions in 𝒫¯​(Δ,Nℝ){\overline{\mathscr{P}}}(\Delta,N_{\mathbb{R}}) (respectively, in 𝒫¯​(Nℝ,Δ){\overline{\mathscr{P}}}(N_{\mathbb{R}},\Delta)) is dense with respect to uniform convergence.

Proof.

This follows from Proposition 3.80 and the density of rational numbers. ∎

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