The MOSEK Python optimizer API manual Version 7.0 (Revision 141)
Optimizer API for Python - Documentation - Mosek Optimizer API for Python - Documentation - Mosek
130 CHAPTER 10. PROBLEM FORMULATION AND SOLUTIONS 10.3 Semidefinite optimization Semidefinite optimization is an extension of conic quadratic optimization (see Section 10.2) allowing positive semidefinite matrix variables to be used in addition to the usual scalar variables. A semidefinite optimization problem can be written as minimize subject to l c i ≤ n−1 ∑ ∑p−1 〈 〉 c j x j + Cj , X j + c f j=0 n−1 j=0 p−1 ∑ ∑ 〈 〉 a ij x j + Aij , X j j=0 j=0 ≤ u c i, i = 0, . . . , m − 1 lj x ≤ x j ≤ u x j , j = 0, . . . , n − 1 x ∈ C, X j ∈ S r + j , j = 0, . . . , p − 1 (10.10) where the problem has p symmetric positive semidefinite variables X j ∈ S r + j of dimension r j with symmetric coefficient matrices C j ∈ S rj and A i,j ∈ S rj . We use standard notation for the matrix inner product, i.e., for U, V ∈ R m×n we have 〈U, V 〉 := m−1 ∑ i=0 n−1 ∑ U ij V ij . With semidefinite optimization we can model a wide range of problems as demonstrated in [7]. j=0 10.3.1 Duality for semidefinite optimization The dual problem corresponding to the semidefinite optimization problem (10.10) is given by maximize subject to (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 c − A T y + s x u − s x l = s x n, m∑ C j − y i A ij = S j , j = 0, . . . , p − 1 i=0 s c l − s c u = y, s c l , s c u, s x l , s x u ≥ 0, s x n ∈ C ∗ , S j ∈ S r + j , j = 0, . . . , p − 1 (10.11) where A ∈ R m×n , A ij = a ij , which is similar to the dual problem for conic quadratic optimization (see Section 10.7), except for the addition of dual constraints m∑ (C j − y i A ij ) ∈ S r + j . Note that the dual of the dual problem is identical to the original primal problem. i=0
10.3. SEMIDEFINITE OPTIMIZATION 131 10.3.2 Infeasibility for semidefinite 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.3.2.1 Primal infeasible problems If the problem (10.10) is infeasible, MOSEK will report a certificate of primal infeasibility: The dual solution reported is a certificate of infeasibility, and the primal solution is undefined. A certificate of primal infeasibility is a feasible solution to the problem maximize subject to (l c ) T s c l − (u c ) T s c u + (l x ) T s x l − (u x ) T s x u A T y + s x l − s x u + s x n = 0, m−1 ∑ y i A ij + S j = 0, j = 0, . . . , p − 1 i=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 ∗ , S j ∈ S r + j , j = 0, . . . , p − 1 such that the objective value is strictly positive. (10.12) 10.3.2.2 Dual infeasible problems If the problem (10.11) 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 where minimize n−1 ∑ ∑p−1 〈 〉 c j x j + Cj , X j j=0 j=0 p−1 subject to ˆlc i ≤ ∑ ∑ 〈 〉 a ij x j + Aij , X j j=1 ˆlx j=0 ≤ û c i, i = 0, . . . , m − 1 ≤ x ≤ û x , x ∈ C, X j ∈ S + r j , j = 0, . . . , p − 1 (10.13) and { 0 if l c ˆlc i = i > −∞, − ∞ otherwise, { 0 if l x ˆlx j = j > −∞, − ∞ otherwise, { and û c 0 if u c i := i < ∞, ∞ otherwise, { 0 if u and û x x j := j < ∞, ∞ otherwise,
- Page 101 and 102: 6.1. SEPARABLE CONVEX (SCOPT) INTER
- Page 103 and 104: 6.1. SEPARABLE CONVEX (SCOPT) INTER
- Page 105 and 106: 6.1. SEPARABLE CONVEX (SCOPT) INTER
- Page 107 and 108: Chapter 7 Advanced API tutorial Thi
- Page 109 and 110: 7.1. THE PROGRESS CALL-BACK 87 71 p
- Page 111 and 112: 7.2. SOLVING LINEAR SYSTEMS INVOLVI
- Page 113 and 114: 7.2. SOLVING LINEAR SYSTEMS INVOLVI
- Page 115 and 116: 7.2. SOLVING LINEAR SYSTEMS INVOLVI
- Page 117 and 118: 7.2. SOLVING LINEAR SYSTEMS INVOLVI
- Page 119 and 120: Chapter 8 A case study 8.1 Portfoli
- Page 121 and 122: 8.1. PORTFOLIO OPTIMIZATION 99 e T
- Page 123 and 124: 8.1. PORTFOLIO OPTIMIZATION 101 is
- Page 125 and 126: 8.1. PORTFOLIO OPTIMIZATION 103 15
- Page 127 and 128: 8.1. PORTFOLIO OPTIMIZATION 105 63
- Page 129 and 130: 8.1. PORTFOLIO OPTIMIZATION 107 8.1
- Page 131 and 132: 8.1. PORTFOLIO OPTIMIZATION 109 110
- Page 133 and 134: 8.1. PORTFOLIO OPTIMIZATION 111 z j
- Page 135 and 136: 8.1. PORTFOLIO OPTIMIZATION 113 Var
- Page 137 and 138: 8.1. PORTFOLIO OPTIMIZATION 115 56
- Page 139 and 140: 8.1. PORTFOLIO OPTIMIZATION 117 172
- Page 141 and 142: Chapter 9 Usage guidelines The purp
- Page 143 and 144: 9.3. WRITING TASK DATA TO A FILE 12
- Page 145 and 146: Chapter 10 Problem formulation and
- Page 147 and 148: 10.1. LINEAR OPTIMIZATION 125 be a
- Page 149 and 150: 10.2. CONIC QUADRATIC OPTIMIZATION
- Page 151: 10.2. CONIC QUADRATIC OPTIMIZATION
- Page 155 and 156: 10.4. QUADRATIC AND QUADRATICALLY C
- Page 157 and 158: Chapter 11 The optimizers for conti
- Page 159 and 160: 11.1. HOW AN OPTIMIZER WORKS 137 11
- Page 161 and 162: 11.2. LINEAR OPTIMIZATION 139 11.2.
- Page 163 and 164: 11.2. LINEAR OPTIMIZATION 141 Whene
- Page 165 and 166: 11.2. LINEAR OPTIMIZATION 143 11.2.
- Page 167 and 168: 11.2. LINEAR OPTIMIZATION 145 • R
- Page 169 and 170: 11.5. NONLINEAR CONVEX OPTIMIZATION
- Page 171 and 172: 11.6. SOLVING PROBLEMS IN PARALLEL
- Page 173 and 174: 11.6. SOLVING PROBLEMS IN PARALLEL
- Page 175 and 176: 11.6. SOLVING PROBLEMS IN PARALLEL
- Page 177 and 178: Chapter 12 The optimizers for mixed
- Page 179 and 180: 12.3. THE MIXED-INTEGER CONIC OPTIM
- Page 181 and 182: 12.5. TERMINATION CRITERION 159 •
- Page 183 and 184: 12.7. UNDERSTANDING SOLUTION QUALIT
- Page 185 and 186: Chapter 13 The analyzers 13.1 The p
- Page 187 and 188: 13.1. THE PROBLEM ANALYZER 165 Cons
- Page 189 and 190: 13.2. ANALYZING INFEASIBLE PROBLEMS
- Page 191 and 192: 13.2. ANALYZING INFEASIBLE PROBLEMS
- Page 193 and 194: 13.2. ANALYZING INFEASIBLE PROBLEMS
- Page 195 and 196: 13.2. ANALYZING INFEASIBLE PROBLEMS
- Page 197 and 198: 13.2. ANALYZING INFEASIBLE PROBLEMS
- Page 199 and 200: Chapter 14 Primal feasibility repai
- Page 201 and 202: 14.2. AUTOMATIC REPAIR 179 One way
10.3. SEMIDEFINITE OPTIMIZATION 131<br />
10.3.2 Infeasibility for semidefinite optimization<br />
In case <strong>MOSEK</strong> finds a problem to be infeasible it reports a certificate of the infeasibility. This works<br />
exactly as for linear problems (see Section 10.1.2).<br />
10.3.2.1 Primal infeasible problems<br />
If the problem (10.10) is infeasible, <strong>MOSEK</strong> will report a certificate of primal infeasibility: <strong>The</strong> dual<br />
solution reported is a certificate of infeasibility, and the primal solution is undefined.<br />
A certificate of primal infeasibility is a feasible solution to the problem<br />
maximize<br />
subject to<br />
(l c ) T s c l − (u c ) T s c u + (l x ) T s x l − (u x ) T s x u<br />
A T y + s x l − s x u + s x n = 0,<br />
m−1<br />
∑<br />
y i A ij + S j = 0, j = 0, . . . , p − 1<br />
i=0<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 ∗ , S j ∈ S r + j<br />
, j = 0, . . . , p − 1<br />
such that the objective value is strictly positive.<br />
(10.12)<br />
10.3.2.2 Dual infeasible problems<br />
If the problem (10.11) is infeasible, <strong>MOSEK</strong> will report a certificate of dual infeasibility: <strong>The</strong> primal<br />
solution reported is the certificate of infeasibility, and the dual solution is undefined.<br />
A certificate of dual infeasibility is a feasible solution to the problem<br />
where<br />
minimize<br />
n−1<br />
∑ ∑p−1<br />
〈 〉<br />
c j x j + Cj , X j<br />
j=0<br />
j=0<br />
p−1<br />
subject to ˆlc i ≤ ∑ ∑ 〈 〉<br />
a ij x j + Aij , X j<br />
j=1<br />
ˆlx<br />
j=0<br />
≤ û c i, i = 0, . . . , m − 1<br />
≤ x ≤ û x ,<br />
x ∈ C, X j ∈ S + r j<br />
, j = 0, . . . , p − 1<br />
(10.13)<br />
and<br />
{ 0 if l<br />
c<br />
ˆlc i = i > −∞,<br />
− ∞ otherwise,<br />
{ 0 if l<br />
x<br />
ˆlx j = j > −∞,<br />
− ∞ otherwise,<br />
{<br />
and û c 0 if u<br />
c<br />
i :=<br />
i < ∞,<br />
∞ otherwise,<br />
{ 0 if u<br />
and û x x<br />
j :=<br />
j < ∞,<br />
∞ otherwise,