FFakultätakultät fürfür MMathematikathematik uundnd IInformatiknformatik
Preprint 2016-06
Susanne Franke, Patrick Mehlitz, Maria Pilecka Optimality conditions for the simple convex bilevel programming problem in Banach spaces
ISSN 1433-9307
Susanne Franke, Patrick Mehlitz, Maria Pilecka
Optimality conditions for the simple convex bilevel programming problem in Banach spaces
TU Bergakademie Freiberg
Fakultät für Mathematik und Informatik
Prüferstraße 9
09599 FREIBERG
http://tu-freiberg.de/fakult1
ISSN 1433 – 9307
Herausgeber: Dekan der Fakultät für Mathematik und Informatik
Herstellung: Medienzentrum der TU Bergakademie Freiberg TU Bergakademie Freiberg Preprint
Optimality conditions for the simple convex bilevel programming problem in Banach spaces
Susanne Franke, Patrick Mehlitz, Maria Pilecka
Abstract The simple convex bilevel programming problem is a convex minimization problem whose feasible set is the solution set of another convex optimization problem. Such problems appear frequently when searching for the projection of a certain point onto the solution set of another program. Due to the nature of the problem, Slater’s constraint qualification generally fails to hold at any feasible point. Hence, one has to formulate weaker constraint qualifications or stationarity notions in order to state optimality conditions. In this paper, we use two different single-level reformulations of the problem, the optimal value and the Karush-Kuhn-Tucker approach, to derive optimality conditions for the original program. Since all these considerations are carried out in Banach spaces, the results are not limited to standard optimization problems in Rn. On the road, we introduce and discuss a certain concept of M-stationarity for mathematical programs with complementarity constraints in Banach spaces. Keywords Bilevel Programming · Convex Programming · Constraint Qualifications · Mathematical Program with Complementarity Constraints · Programming in Banach Spaces Mathematics Subject Classification (2000) 46N10 · 49K27 · 90C25 · 90C33
Dedicated to Professor Stephan Dempe on the occasion of his 60th birthday.
1 Introduction
Bilevel programming problems were first formulated by Stackelberg in his investigation of market econ- omy in 1934 [49], whereas the first mathematical model can be found in [10]. Since then, bilevel problems have been studied thoroughly from both a theoretical point of view and with regard to applications as well as solution strategies, see the monographs [6,12,46]. In order to handle a bilevel programming prob- lem mathematically, it is usually transformed into an optimization problem with only one level. Two such possibilities which we will use in our paper are the optimal value reformulation, introduced in [39], and the Karush-Kuhn-Tucker (KKT) approach, see for example [34]. Instead of the standard bilevel programming problem which has one variable for the upper and one for the lower level, we investigate the special (simple) case where both levels share the same variable. Hence, consider the simple convex bilevel programming problem (SCBPP for short) which is given as stated below: F (x) → min (SCBPP) x ∈ Ψ := Argmin{f(y) | y ∈ Θ}.
Here, F, f : X → R are convex, continuous, and sufficiently smooth functionals of a Banach space X , while Θ ⊆ X is a nonempty, closed, and convex set. Note that these assumptions guarantee that the set Ψ is convex, i.e. (SCBPP) is a convex optimization problem. However, the derivation of optimality
S. Franke, P. Mehlitz (corresponding author), and M. Pilecka Technische Universit¨atBergakademie Freiberg, Faculty of Mathematics and Computer Science, 09596 Freiberg, Germany E-mail: {susanne.franke,mehlitz,pilecka}@math.tu-freiberg.de, the work of the first author has been supported by the Deutsche Forschungsgemeinschaft, grant DE-650/7-1 2 Susanne Franke, Patrick Mehlitz, Maria Pilecka
conditions is a challenging problem since most of the standard constraint qualifications generally fail to hold at any feasible point of (SCBPP). Problems of this type arise frequently when a best point (in a certain sense) among the optimal solutions of a convex optimization problem has to be found. Note that searching for an efficient solution of a convex biobjective optimization problems can also be reformulated as (SCBPP). However, solving (SCBPP) provides only special efficient solutions of such a vector optimization problem. Regarding the finite-dimensional situation, (SCBPP) was first introduced and investigated in [48] where an application of this problem is given as well. Further discussions and its relationship to standard bilevel programming problems can be found in [13]. In order to stay as close as possible to the terminology of bilevel programming, we call
f(x) → min x ∈ Θ
the lower level problem of (SCBPP). If the lower level feasible set Θ is described by conic constraints, it is not difficult to show, see Theorem 5.5, that (SCBPP) is equivalent to a certain mathematical program with complementarity constraints (MPCC for short) in Banach spaces, recently introduced and studied in [36, 51, 53]. These problems in finite-dimensional spaces have been the subject of theoretical research for many years. Since the usual KKT conditions for such problems are very restrictive, several different stationarity conditions have been defined, see [18, 41, 45, 60]. Furthermore, since some of the well-known constraint qualifications are violated at every feasible point of such problems, other regularity conditions for MPCCs have been introduced, cf. [17, 19]. The numerical approach to MPCCs is developed as well, see [29] and references therein. For the systematic and comprehensive study on this topic, we refer the interested reader to [34]. Moreover, a study of MPCCs in very general settings can also be found in [38]. Due to the framework of the considered problem in Banach spaces, it is also connected to works on optimal control. Hierarchical programming in function spaces is the subject of discussion in [8,36,55,56]. The special case where the solution of a variational inequality in Sobolev spaces has to be controlled was investigated in [22, 27, 28, 31, 37, 50]. It turns out that the strong stationarity conditions of the surrogate problem are, in general, too strong to hold at the local minimizers of (SCBPP), see Section5. Hence, weaker stationarity notions have to be considered. Therefore, in this paper, we want to consider a generalized concept of Mordukhovich- stationarity (M-stationarity for short) which is known to be weaker than strong stationarity for common finite-dimensional MPCCs. We organized this article as follows: In Section2, we introduce the notation and preliminary results we use throughout the paper. Section3 is dedicated to the study of constraint qualifications for optimization problems in Banach spaces. Afterwards, MPCCs in Banach spaces are discussed in Section4. Here, we recall the notions of weak and strong stationarity already known from [36,51]. Furthermore, we introduce a concept of M-stationarity which is shown to be reasonable in the case where the cone which generates the complementarity constraint is polyhedral. Finally, we consider (SCBPP) in Section5. The lower level optimal objective value and the lower KKT conditions of (SCBPP) are used to state two single-level surrogate problems which are exploited to derive optimality conditions for (SCBPP). The aforementioned KKT approach leads to necessary optimality conditions of M-stationarity-type.
2 Notation and preliminary results
2.1 Notation
n n,+ Let us start this section by introducing the notation we exploit in this paper. We use N, R, R, R , R0 , n×m R , and Sn to denote the natural numbers (without zero), the real numbers, the extended real line, the set of all real vectors with n components, the cone of all vectors from Rn possessing nonnegative components, the set of all real matrices with n rows and m columns, and the set of all symmetric matrices n×n + 1,+ n×m from R , respectively. Furthermore, we stipulate R0 := R0 . For an arbitrary matrix Q ∈ R and |I|×m an index set I ⊆ {1, . . . , n}, QI ∈ R denotes the matrix composed of the rows of Q whose indices n×m n×n come from I. Furthermore, O ∈ R and I ∈ R represent the zero matrix and the identity matrix Optimality conditions for the simple convex bilevel programming problem in Banach spaces 3
+ − of appropriate dimension, respectively. By Sn and Sn we denote the cone of all positive and negative semidefinite matrices from Sn, respectively. The forthcoming definitions and results in this section are based on [5, 22, 38, 40]. Let X be a real Banach space with norm k·kX and zero vector 0, and let A ⊆ X be nonempty. We denote by lin(A), cone(A), conv(A), int(A), and cl(A) the smallest subspace of X containing A, the smallest convex cone containing A, the convex hull of A, the interior of A, and the closure of A, respectively. The indicator function δA : X → R of the set A is defined by ( 0 if x ∈ A, ∀x ∈ X : δA(x) := +∞ if x∈ / A.
The (topological) dual space of X is denoted by X ?, and h·, ·i: X × X ? → R is the corresponding dual pairing. Now fix some x ∈ A. Then we define the radial cone, the tangent (or Bouligand) cone, and the weak tangent cone to A at x as stated below:
RA(x) := {d ∈ X | ∃α > 0 ∀t ∈ (0, α): x + td ∈ A}, TA(x) := {d ∈ X | ∃{dk} ⊆ X ∃{tk} ⊆ R: dk → d, tk & 0, x + tkdk ∈ A ∀k ∈ N}, w TA (x) := {d ∈ X | ∃{dk} ⊆ X ∃{tk} ⊆ R: dk * d, tk & 0, x + tkdk ∈ A ∀k ∈ N}. Here, → and * denote the norm and the weak convergence in X , respectively. Furthermore, we will use ? ? → to denote the weak-?-convergence in X . Observe that the cone TA(x) is always a closed set which satisfies w RA(x) ⊆ TA(x) ⊆ TA (x). ε ε Let BX be the closed unit ball of X , while for any ε > 0 and x ∈ X , we use UX (x) and BX (x) to denote the open and closed ε-ball around x, respectively. Let B ⊆ X ? be nonempty. Then cl?(B) denotes the weak-?-closure of B, i.e. the closure with respect to (w.r.t.) the weak-?-topology. We define the polar cones and annihilators of the sets A and B by
◦ ∗ ? ∗ ∗ ∗ A := {x ∈ X | ∀x ∈ A: hx, x i ≤ 0},B◦ := {x ∈ X | ∀x ∈ B : hx, x i ≤ 0}, ⊥ ∗ ? ∗ ∗ ∗ A := {x ∈ X | ∀x ∈ A: hx, x i = 0},B⊥ := {x ∈ X | ∀x ∈ B : hx, x i = 0},
∼ ?? ◦ respectively. In the case where X is reflexive, i.e. if X = X is valid, we can exploit B◦ = B and ⊥ B⊥ = B . For the purpose of simplicity, we omit curly brackets when considering singletons, i.e. we set lin(x) := lin({x}) and x⊥ := {x}⊥ for any x ∈ X , and similar definitions shall hold for all x∗ ∈ X ?. Now X and Y may be arbitrary Banach spaces again. Then the product space X × Y is a Banach space, too, when, e.g., equipped with the sum norm induced by k·kX and k·kY . In the case where Y = X 2 2 holds, we use X := X × X and the components of any x ∈ X are addressed by x1, x2 ∈ X . We exploit a similar notation for all n ∈ N satisfying n ≥ 3 as well as arbitrary subsets of A ⊆ X , i.e., An denotes the Cartesian product of order n of A. For a closed, convex set C ⊆ X and a fixed point x ∈ C, we obtain
RC (x) = cone(C − {x}) and TC (x) = cl(RC (x)),
∗ ◦ ∗ respectively. Moreover, if x ∈ TC (x) is chosen, then the critical cone to C w.r.t. x and x is defined by
∗ ∗ ⊥ KC (x, x ) := TC (x) ∩ (x ) .
In the case that C is additionally a cone, we easily obtain
◦ ◦ ⊥ RC (x) = C + lin(x), TC (x) = cl(C + lin(x)), TC (x) = C ∩ x .
∗ ∗ ◦ The (not necessarily conic) set C is called polyhedric w.r.t. (x, x ), where x ∈ C and x ∈ TC (x) hold, if ∗ ⊥ ∗ cl(RC (x) ∩ (x ) ) = KC (x, x ) is satisfied. Moreover, C is called polyhedric if it is polyhedric w.r.t. all points (x, x∗) satisfying x ∈ C ∗ ◦ and x ∈ TC (x) . The concept of polyhedricity was introduced in [22,37] and was recently studied in [52]. 4 Susanne Franke, Patrick Mehlitz, Maria Pilecka
∗ ∗ ? Recall that C is called polyhedral if there exist functionals x1, . . . , xn ∈ X and scalars β1, . . . , βn ∈ R such that C possesses the representation
∗ C = {x ∈ X | hx, xi i ≤ βi ∀i ∈ {1, . . . , n}}. Since the radial cone to a polyhedral set at an arbitrary point is already closed, polyhedral sets are polyhedric everywhere, see Lemma 2.3 as well. For a Banach space X and a closed (but not necessarily convex) set D ⊆ X satisfying x ∈ D, we define the set of all Fr´echet σ-normals to D at x as stated below: ∗ σ ∗ ? hy − x, x i NbD(x) := x ∈ X lim sup ≤ σ . y→x, y∈D ky − xkX
0 The closed, convex cone NbD(¯x) := NbD(¯x) is called the Fr´echet normal cone to D atx ¯. If X is re- w ◦ flexive, NbD(x) = TD (x) follows from [38, Theorem 1.10]. Next, we introduce the basic (or limiting, Mordukhovich) normal cone to D at x by
( ∗ ? ) ∃{σk} ⊆ R ∃{xk} ⊆ D ∃{xk} ⊆ X : N (x) := x∗ ∈ X ? . D ∗ ? ∗ ∗ σ σk & 0, xk → x,¯ xk → x , xk ∈ NbD(xk) ∀k ∈ N From [38, Theorem 2.35], we can fix σ to be zero in the above definition as long as X is an Asplund space, i.e. a Banach space whose separable subspaces possess separable duals. Note that any reflexive Banach space possesses the Asplund property. Furthermore, we can replace the weak-?-convergence in X ? by the weak convergence in the setting of reflexive Banach spaces. In the case of the set D being convex, we have ◦ ND(x) = NbD(x) = TD(x) .
We call D sequentially normally compact (SNC for short) at x if for any sequences {σk} ⊆ R, {xk} ⊆ D, ∗ ? ∗ σk ∗ ? and {xk} ⊆ X such that xk ∈ NbD (xk) is valid for all k ∈ N and σk & 0, xk → x, as well as xk → 0 ∗ hold, xk → 0 is satisfied. Again, we can fix σ to zero in the Asplund setting. Obviously, any subset of a finite-dimensional Banach space is SNC at any of its points. On the other hand, a singleton in X is SNC if and only if X is finite-dimensional, see [38, Theorem 1.21]. Let U, V, and W be arbitrary Banach spaces. By L[U, V] we denote the space of bounded linear operators mapping from U to V. For any operator F ∈ L[U, V], the operator F? ∈ L[V?, U ?] denotes the adjoint of F. The identical mapping of U, denoted by IU , is always a continuous operator as long as U is equipped with the same norm in its definition and image space. For another operator G ∈ L[U, W], the linear operator (F, G) ∈ L[U, V × W] is defined as stated below: ∀u ∈ U :(F, G)[u] := (F[u], G[u]). Let Γ : U → 2V be an arbitrary set-valued mapping. Then gph Γ := {(u, v) ∈ U×V | v ∈ Γ (u)} denotes its graph. Recall that Γ is called Lipschitz-like at some point (¯u, v¯) ∈ gph Γ if there are constants ε > 0, ρ > 0, and L > 0 such that
0 ε ρ 0 0 ∀u, u ∈ UU (¯u): Γ (u) ∩ UV (¯v) ⊆ Γ (u ) + Lku − u kU BV is satisfied. If the above condition holds whenever u0 :=u ¯ is fixed, then Γ is called calm at (¯u, v¯). Obviously, calmness is implied by the Lipschitz-like property. Note that the latter condition is also known as pseudo Lipschitz or Aubin property, see [12, 25]. ∗ ? U ? Assume that U and V are reflexive Banach spaces. The normal coderivative DN Γ (¯u, v¯): V → 2 is defined by ∗ ? ∗ ∗ ∗ ? ∗ ∗ ∀v ∈ V : DN Γ (¯u, v¯)(v ) = {u ∈ U | (u , −v ) ∈ NgphΓ (¯u, v¯)} . We call Γ partially sequentially normally compact (PSNC for short) at the point (¯u, v¯) ∈ gph Γ if for any ∗ ∗ ? ? ∗ ∗ sequences {(uk, vk)} ⊆ gph Γ and {(uk, vk)} ⊆ U × V which satisfy (uk, vk) → (¯u, v¯), uk * 0, vk → 0, ∗ ∗ ∗ and (uk, vk) ∈ Nbgph Γ (uk, vk) for any k ∈ N, we obtain uk → 0. Note that if Γ is PSNC at (¯u, v¯), satisfies ∗ DN Γ (¯u, v¯)(0) = {0}, and possesses a closed graph in a neighborhood of (¯u, v¯), then Γ is Lipschitz-like at the latter point, see [38, Theorem 4.10] where a weaker condition w.r.t. the so-called mixed coderivative has been formulated. Optimality conditions for the simple convex bilevel programming problem in Banach spaces 5
We drop the reflexivity assumption on U and V. Let φ: U → R be a convex functional andu ¯ ∈ U be a point where |φ(¯u)| < ∞ is satisfied. Then
∂φ(¯u) := {u∗ ∈ U ? | φ(u) ≥ φ(¯u) + hu − u,¯ u∗i ∀u ∈ U}
denotes the subdifferential of φ atu ¯. For a nonempty, closed, and convex set S ⊆ U, the functional δS is ◦ convex, and it is easy to see that for any pointu ¯ ∈ S, we have ∂δS(¯u) = TS(¯u) . For some convex cone K ⊆ V, we say that a mapping ψ : U → V is K-convex if the following condition holds:
∀u, u0 ∈ U ∀σ ∈ (0, 1): − σψ(u) − (1 − σ)ψ(u0) + ψ(σu + (1 − σ)u0) ∈ K.
Note that for a K-convex mapping, the set {u ∈ U | ψ(u) ∈ K} is convex. Suppose that θ : U → V is a twice Fr´echet differentiable mapping andu ¯ ∈ U is fixed. Then the linear operators θ0(¯u) ∈ L[U, V] and θ(2)(¯u) ∈ L[U, L[U, V]] denote its first and second order Fr´echet derivative atu ¯, respectively. For any ξ ∈ V?, we use hθ(2)(¯u), ξi to represent the linear operator in L[U, U ?] defined below: (2) (2) ? (2) ∀du ∈ U : hθ (¯u), ξi[du] := θ (¯u)[du] [ξ] = ξ ◦ θ (¯u)[du] . Here, ◦ denotes the composition of mappings. For some bounded domain Ω ⊆ Rd, C(Ω) denotes the Banach space of all scalar continuous functions on cl(Ω) equipped with the maximum norm. We identify Ω with the measure space (Ω,Σ, | · |) where | · | denotes the Lebesgue measure and Σ represents the σ-algebra of all Lebesgue-measurable subsets of Ω. Fix some constant p ∈ [1, ∞]. Then Lp(Ω) denotes the Banach space of (equivalence classes of) measurable functions on Ω equipped with the norm 1 Z p p p ∀u ∈ L (Ω): kukLp(Ω) := |u(ω)| dω Ω in the case p ∈ [1, ∞) and
∞ ∀u ∈ L (Ω): kukL∞(Ω) := inf sup |u(ω)|, N∈Σ, |N|=0 ω∈Ω\N
0 otherwise. For p ∈ [1, ∞), the dual space of Lp(Ω) is Lp (Ω) where p0 ∈ (1, ∞] denotes the so-called conjugate coefficient of p characterized via 1/p + 1/p0 = 1. The corresponding dual pairing is given by
0 Z ∀u ∈ Lp(Ω) ∀v ∈ Lp (Ω): hu, vi := u(ω)v(ω)dω. Ω
p If p ∈ (1, ∞) is valid, L (Ω) is reflexive. For any set S ∈ Σ, χS : Ω → R denotes the characteristic function of S given by ( 1 if ω ∈ S ∀ω ∈ Ω : χS(ω) := 0 if ω∈ / S.
p 1/p Clearly, we have χS ∈ L (Ω) for all p ∈ [1, ∞]. Moreover, kχSkLp(Ω) = |S| is valid for any p ∈ [1, ∞). For some p ∈ (1, ∞), W 1,p(Ω) denotes the Sobolev space of all first order weakly differentiable functions p p from L (Ω) whose weak derivatives ∂ω1 , . . . , ∂ωd belong to L (Ω) as well. We equip this space with the norm defined below:
1 d ! p 1,p p X p 1,p ∀u ∈ W (Ω): kukW (Ω) := kukLp(Ω) + k∂ωi ukLp(Ω) . i=1 A detailed introduction to these and more general function spaces can be found in [3].
2.2 Preliminary results
In this section, we present all supplementary results which are used in this paper. 6 Susanne Franke, Patrick Mehlitz, Maria Pilecka
2.2.1 Some facts on bounded linear operators
The first part of the following result is called Generalized Farkas Lemma and is stated in a more abstract setting, e.g., in [20, Theorem 1, Lemma 3].
Lemma 2.1 Let X , Y be Banach spaces, and suppose that K ⊆ Y is a nonempty, closed, convex cone. Assume that φ: X → R is a positively homogeneous and convex functional while A ∈ L[X , Y] is a bounded linear operator.
1. For any ξ ∈ X ?, the following assertions are equivalent: (a) ξ ∈ cl? ∂φ(0) + A?[K◦], (b) ∀x ∈ X : A[x] ∈ K =⇒ hx, ξi ≤ φ(x).
Particularly, setting φ(·) := δTS (¯x)(·) (indicator function of the tangent cone to S at x¯) for some closed, convex set S ⊆ X and x¯ ∈ S, we have
? ◦ ? ◦ ◦ cl TS(¯x) + A [K ] = {x ∈ TS(¯x) | A[x] ∈ K} .
2. Suppose that A[X ] − K = Y is satisfied. Then we have
A?[K◦] = {x ∈ X | A[x] ∈ K}◦.
Note that if X is reflexive, then we can replace the weak-?-closure in Lemma 2.1 by the common closure since the weak and weak-?-topology in X ? are equivalent, and convex sets are closed if and only if they are weakly closed by Mazur’s theorem. Indeed, the first statement of Lemma 2.1 is a generalization of the famous Farkas Lemma.
n m m×n n m,+ Remark 2.1 Choosing X = R , Y = R , A ∈ R , b ∈ R , K := R0 , as well as ξ = 0 and defining φ: Rn → R by φ(x) := b>x for any x ∈ Rn, Lemma 2.1 yields that the assertions
1. ∃y ∈ Rm : A>y = b, y ≥ 0, 2. ∀x ∈ Rn : Ax ≥ 0 =⇒ b>x ≥ 0 are equivalent. Hence, precisely one of the following systems possesses a solution:
( ( A>y = b Ax ≥ 0 y ≥ 0 b>x < 0.
Moreover, we can characterize the adjoint of a dense range operator simply by applying Lemma 2.1.
Remark 2.2 Let X and Y be reflexive Banach spaces, while A ∈ L[X , Y] is an arbitrary operator. Then A? is injective if and only if A possesses a dense range.
Lemma 2.2 Let A ∈ L[X , Y] be an injective, bounded, linear operator with a closed range between Banach spaces X and Y. Then there is a constant c > 0 satisfying ckxkX ≤ kA[x]kY for all x ∈ X .
Proof Clearly, due to its closedness, A[X ] is a Banach space as well. Hence, we can define a bijection A˜ ∈ L[X , A[X ]] by means of A˜[x] := A[x] for any x ∈ X . Then A˜−1 is a bounded linear operator as well, −1 i.e., there is a constantc ˜ > 0 which satisfies kA˜ [y]kX ≤ c˜kykY for all y ∈ A[X ]. For arbitrary x ∈ X , we obtain
−1 kxkX = kA˜ [A[x]]kX ≤ c˜kA[x]kY .
1 Hence, the statement of the lemma follows by choosing c := c˜. ut Optimality conditions for the simple convex bilevel programming problem in Banach spaces 7
2.2.2 On polyhedral cones
∗ ∗ ? Let Z be a reflexive Banach space and {z1 , . . . , zk} ⊆ Z be a set of linear independent functionals. We consider the closed, convex, polyhedral cone
∗ K := {z ∈ Z | hz, zi i ≤ 0 ∀i ∈ {1, . . . , k}}. (1)
As mentioned earlier, this set is polyhedric everywhere. Although this result is well-known, we provide a short proof for the sake of completeness.
Lemma 2.3 Let the polyhedral cone K be given as stated in (1). Then K is polyhedric.
Proof The assertion clearly follows if we can show RK (¯u) = TK (¯u) for arbitraryu ¯ ∈ K. Therefore, ∗ chooseu ¯ ∈ K arbitrarily, define I(¯u) := {i ∈ {1, . . . , k} | hu,¯ zi i = 0}, and observe
∗ RK (¯u) = K + lin(¯u) = {z + κu¯ ∈ Z | κ ∈ R, hz, zi i ≤ 0 ∀i ∈ {1, . . . , k}} ∗ = {d ∈ Z | ∃κ ∈ R: hd − κu,¯ zi i ≤ 0 ∀i ∈ {1, . . . , k}} ∗ ∗ ∗ = {d ∈ Z | ∃κ ∈ R: hd, zi i ≤ 0 ∀i ∈ I(¯u), hd, zi i ≤ κhu,¯ zi i ∀i ∈ {1, . . . , k}\ I(¯u)} ∗ = {d ∈ Z | hd, zi i ≤ 0 ∀i ∈ I(¯u)}.
Obviously, RK (¯u) is closed and, hence, it coincides with TK (¯u). This completes the proof. ut
Later we will define a generalized notion of M-stationarity for MPCCs where the cone which induces the complementarity is given by (1). Before we formulate a lemma which summarizes the results obtained in [24], we introduce some necessary notation. From the definition of the polar cone,
n o ◦ Pk ∗ ? K = i=1αizi ∈ Z αi ≥ 0 ∀i ∈ {1, . . . , k}
follows easily. Let us analyze the complementarity set
F := {(u, v) ∈ K × K◦ | hu, vi = 0}.
Choosing (¯u, v¯) ∈ F, we define index sets I(¯u, v¯) and J(¯u, v¯) by
I(¯u, v¯) := {i ∈ {1, . . . , k} | hu,¯ z∗i = 0}, i (2) J(¯u, v¯) := {i ∈ I(¯u, v¯) | α¯i > 0}.
Pk ∗ Here, we exploited the unique representationv ¯ = i=1 α¯izi for someα ¯1,..., α¯k ≥ 0. For any index sets J := J(¯u, v¯) ⊆ P ⊆ Q ⊆ I(¯u, v¯) =: I, we define
∗ ∗ CQ,P := cone({zi | i ∈ Q \ P }) + lin({zi | i ∈ P }), ∗ ∗ ⊥ DQ,P := {z ∈ Z | hz, zi i ≤ 0 ∀i ∈ Q \ P } ∩ {zi | i ∈ P } .
Lemma 2.4 Using the above notations, the following relations hold:
[ NbF (¯u, v¯) = CI,J × DI,J , NF (¯u, v¯) = CQ,P × DQ,P . J⊆P ⊆Q⊆I 8 Susanne Franke, Patrick Mehlitz, Maria Pilecka
2.2.3 On the SNC property of the nonnegative cone in different function spaces
Let Ω ⊆ Rd be a bounded domain. In this section, we want to comment on the SNC property of the cone of nonnegative functions in different function spaces defined on Ω. Therefore, we provide the following three lemmata which describe the situation w.r.t. continuous functions, functions from Lebesgue spaces, and functions from certain Sobolev spaces.
Lemma 2.5 The closed, convex cone
+ C(Ω)0 := {u ∈ C(Ω) | u(ω) ≥ 0 for all ω ∈ Ω}
is SNC everywhere.
+ Proof For some function u ∈ C(Ω), we obtain max{u; 0}, max{−u; 0} ∈ C(Ω)0 . This shows the relation + + lin C(Ω)0 = C(Ω). Noting that int C(Ω)0 6= ∅ is valid, the lemma’s assertion follows from [38, Theorem 1.21]. ut
Lemma 2.6 For any p ∈ [1, ∞), the closed, convex cone
p + p L (Ω)0 := {u ∈ L (Ω) | u(ω) ≥ 0 f.a.a. ω ∈ Ω}
∞ + is nowhere SNC. On the other hand, the closed, convex cone L (Ω)0 is SNC everywhere.
∞ + Proof For p = ∞, we can adapt the proof of Lemma 2.5 since L (Ω)0 possesses a nonempty interior and its span already equals L∞(Ω). Thus, let us fix p ∈ (1, ∞), its conjugate coefficient p0 ∈ (1, ∞) (i.e. the number satisfying 1/p+1/p0 = 1), p + Ω and some arbitraryu ¯ ∈ L (Ω)0 . Clearly, since Ω is open, we find a sequence {Ωk} ⊆ 2 of measurable
sets with positive measure such that |Ωk| & 0 holds true. Let us define uk := (1 − χΩk )¯u. We easily see p uk → u¯ in L (Ω) from the dominated convergence theorem, see [47, Theorem 5.2.2]. Let us define
1 p0 ∀ω ∈ Ω : ηk(ω) := −|Ωk| χΩk (ω).
This leads to