The MOSEK Python optimizer API manual Version 7.0 (Revision 141)

Optimizer API for Python - Documentation - Mosek Optimizer API for Python - Documentation - Mosek

25.11.2015 Views

128 CHAPTER 10. PROBLEM FORMULATION AND SOLUTIONS • The quadratic cone: ⎧ ⎨ Q n = ⎩ x ∈ n∑ Rn : x 1 ≥ √ j=2 x 2 j ⎫ ⎬ ⎭ . • The rotated quadratic cone: ⎧ ⎨ Q r n = ⎩ x ∈ Rn : 2x 1 x 2 ≥ ⎫ n∑ ⎬ x 2 j, x 1 ≥ 0, x 2 ≥ 0 ⎭ . j=3 Although these cones may seem to provide only limited expressive power they can be used to model a wide range of problems as demonstrated in [7]. 10.2.1 Duality for conic quadratic optimization The dual problem corresponding to the conic quadratic optimization problem (10.6) is given by maximize (l c ) T s c l − (u c ) T s c u + (l x ) T s x l − (u x ) T s x u + c f subject to A T y + s x l − s x u + s x n = c, − y + s c l − s c u = 0, s c l , s c u, s x l , s x u ≥ 0, s x n ∈ C ∗ , where the dual cone C ∗ is a Cartesian product of the cones (10.7) C ∗ = C ∗ 1× · · · ×C ∗ p, where each C ∗ t is the dual cone of C t . For the cone types MOSEK can handle, the relation between the primal and dual cone is given as follows: • The R n set: • The quadratic cone: C t = R nt ⇔ C ∗ t = {s ∈ R nt : s = 0} . ⎧ ⎨ C t = Q nt ⇔ Ct ∗ = Q nt = ⎩ s ∈ ∑n t Rnt : s 1 ≥ √ j=2 s 2 j ⎫ ⎬ ⎭ . • The rotated quadratic cone: ⎧ ⎨ C t = Q r n t ⇔ Ct ∗ = Q r n t = ⎩ s ∈ Rnt : 2s 1 s 2 ≥ ∑n t j=3 ⎫ ⎬ s 2 j, s 1 ≥ 0, s 2 ≥ 0 ⎭ .

10.2. CONIC QUADRATIC OPTIMIZATION 129 Please note that the dual problem of the dual problem is identical to the original primal problem. 10.2.2 Infeasibility for conic quadratic optimization In case MOSEK finds a problem to be infeasible it reports a certificate of the infeasibility. This works exactly as for linear problems (see Section 10.1.2). 10.2.2.1 Primal infeasible problems If the problem (10.6) is infeasible, MOSEK will report a certificate of primal infeasibility: The dual solution reported is the certificate of infeasibility, and the primal solution is undefined. A certificate of primal infeasibility is a feasible solution to the problem such that the objective value is strictly positive. maximize (l c ) T s c l − (u c ) T s c u + (l x ) T s x l − (u x ) T s x u subject to A T y + s x l − s x u + s x n = 0, − y + s c l − s c u = 0, s c l , s c u, s x l , s x u ≥ 0, s x n ∈ C ∗ , (10.8) 10.2.2.2 Dual infeasible problems If the problem (10.7) is infeasible, MOSEK will report a certificate of dual infeasibility: The primal solution reported is the certificate of infeasibility, and the dual solution is undefined. A certificate of dual infeasibility is a feasible solution to the problem minimize subject to ˆlc ˆlx c T x ≤ Ax ≤ û c , ≤ x ≤ û x , x ∈ C, (10.9) where and { 0 if l c ˆlc i = i > −∞, − ∞ otherwise, { and û c 0 if u c i := i < ∞, ∞ otherwise, { 0 if l x ˆlx j = j > −∞, − ∞ otherwise, such that the objective value is strictly negative. { 0 if u and û x x j := j < ∞, ∞ otherwise,

128 CHAPTER 10. PROBLEM FORMULATION AND SOLUTIONS<br />

• <strong>The</strong> quadratic cone:<br />

⎧<br />

⎨<br />

Q n =<br />

⎩ x ∈ n∑<br />

Rn : x 1 ≥ √<br />

j=2<br />

x 2 j<br />

⎫<br />

⎬<br />

⎭ .<br />

• <strong>The</strong> rotated quadratic cone:<br />

⎧<br />

⎨<br />

Q r n =<br />

⎩ x ∈ Rn : 2x 1 x 2 ≥<br />

⎫<br />

n∑<br />

⎬<br />

x 2 j, x 1 ≥ 0, x 2 ≥ 0<br />

⎭ .<br />

j=3<br />

Although these cones may seem to provide only limited expressive power they can be used to model a<br />

wide range of problems as demonstrated in [7].<br />

10.2.1 Duality for conic quadratic optimization<br />

<strong>The</strong> dual problem corresponding to the conic quadratic optimization problem (10.6) is given by<br />

maximize (l c ) T s c l − (u c ) T s c u + (l x ) T s x l − (u x ) T s x u + c f<br />

subject to A T y + s x l − s x u + s x n = c,<br />

− y + s c l − s c u = 0,<br />

s c l , s c u, s x l , s x u ≥ 0,<br />

s x n ∈ C ∗ ,<br />

where the dual cone C ∗ is a Cartesian product of the cones<br />

(10.7)<br />

C ∗ = C ∗ 1× · · · ×C ∗ p,<br />

where each C ∗ t is the dual cone of C t . For the cone types <strong>MOSEK</strong> can handle, the relation between the<br />

primal and dual cone is given as follows:<br />

• <strong>The</strong> R n set:<br />

• <strong>The</strong> quadratic cone:<br />

C t = R nt ⇔ C ∗ t = {s ∈ R nt : s = 0} .<br />

⎧<br />

⎨<br />

C t = Q nt ⇔ Ct ∗ = Q nt =<br />

⎩ s ∈ ∑n t<br />

Rnt : s 1 ≥ √<br />

j=2<br />

s 2 j<br />

⎫<br />

⎬<br />

⎭ .<br />

• <strong>The</strong> rotated quadratic cone:<br />

⎧<br />

⎨<br />

C t = Q r n t<br />

⇔ Ct ∗ = Q r n t<br />

=<br />

⎩ s ∈ Rnt : 2s 1 s 2 ≥<br />

∑n t<br />

j=3<br />

⎫<br />

⎬<br />

s 2 j, s 1 ≥ 0, s 2 ≥ 0<br />

⎭ .

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!