The MOSEK command line tool Version 7.0 (Revision 141)

The MOSEK command line tool. Version 7.0 ... - Documentation The MOSEK command line tool. Version 7.0 ... - Documentation

25.11.2015 Views

38 CHAPTER 5. THE OPTIMIZERS FOR CONTINUOUS PROBLEMS A x∗ τ ∗ = b, A T y∗ τ ∗ + s∗ = c, τ∗ − c T x∗ y∗ + bT τ ∗ τ ∗ = 0, x ∗ , s ∗ , τ ∗ , κ ∗ ≥ 0. This shows that x∗ τ is a primal optimal solution and ( y∗ ∗ as the optimal interior-point solution since is a primal-dual optimal solution. On other hand, if κ ∗ > 0 then (x, y, s) = τ , s∗ ∗ τ∗ ( ) x ∗ τ ∗ , y∗ τ ∗ , s∗ τ∗ ) is a dual optimal solution; this is reported This implies that at least one of Ax ∗ = 0, A T y ∗ + s ∗ = 0, − c T x ∗ + b T y ∗ = κ ∗ , x ∗ , s ∗ , τ ∗ , κ ∗ ≥ 0. or − c T x ∗ > 0 (5.3) b T y ∗ > 0 (5.4) is satisfied. If (5.3) is satisfied then x ∗ is a certificate of dual infeasibility, whereas if (5.4) is satisfied then y ∗ is a certificate of dual infeasibility. In summary, by computing an appropriate solution to the homogeneous model, all information required for a solution to the original problem is obtained. A solution to the homogeneous model can be computed using a primal-dual interior-point algorithm [5]. 5.2.2.1 Interior-point termination criterion For efficiency reasons it is not practical to solve the homogeneous model exactly. Hence, an exact optimal solution or an exact infeasibility certificate cannot be computed and a reasonable termination criterion has to be employed. In every iteration, k, of the interior-point algorithm a trial solution to homogeneous model is generated where (x k , y k , s k , τ k , κ k )

5.2. LINEAR OPTIMIZATION 39 Whenever the trial solution satisfies the criterion yk AT ∥ ( ∣ ∣ (x k ) T s k ∣∣∣ min (τ k ) 2 , c T x k τ k − bT y k ) ∣∣∣ τ k the interior-point optimizer is terminated and x k , s k , τ k , κ k > 0. ∥ ∥∥∥ ∥ Axk τ k − b ∞ ≤ ɛ p (1 + ‖b‖ ∞ ), ∥ τ k + sk ∥∥∥ τ k − c ∞ ≤ ɛ d (1 + ‖c‖ ∞ ), and ( ≤ ɛ g max 1, min(∣ ∣c T x k∣ ∣ , ∣ ∣b T y k∣ ) ∣ ) τ k , (5.5) (x k , y k , s k ) τ k is reported as the primal-dual optimal solution. The interpretation of (5.5) is that the optimizer is terminated if • xk τ k • is approximately primal feasible, ( ) y k , sk τ k τ k is approximately dual feasible, and • the duality gap is almost zero. On the other hand, if the trial solution satisfies −ɛ i c T x k > ‖c‖ ∞ ∥ ∥ Ax k ∞ max(1, ‖b‖ ∞ ) then the problem is declared dual infeasible and x k is reported as a certificate of dual infeasibility. The motivation for this stopping criterion is as follows: First assume that ∥ ∥ Ax k ∥ ∥ ∞ = 0 ; then x k is an exact certificate of dual infeasibility. Next assume that this is not the case, i.e. and define ∥ Ax k ∥ ∥ ∞ > 0, It is easy to verify that ¯x := ɛ i max(1, ‖b‖ ∞ ) ‖Ax k ‖ ∞ ‖c‖ ∞ x k . ‖A¯x‖ ∞ = ɛ i max(1, ‖b‖ ∞ ) ‖c‖ ∞ and − c T ¯x > 1, which shows ¯x is an approximate certificate of dual infeasibility where ɛ i controls the quality of the approximation. A smaller value means a better approximation.

5.2. LINEAR OPTIMIZATION 39<br />

Whenever the trial solution satisfies the criterion<br />

yk<br />

AT<br />

∥<br />

( ∣ ∣ (x k ) T s k ∣∣∣<br />

min<br />

(τ k ) 2 , c T x k<br />

τ k − bT y k ) ∣∣∣<br />

τ k<br />

the interior-point optimizer is terminated and<br />

x k , s k , τ k , κ k > 0.<br />

∥ ∥∥∥ ∥ Axk τ k − b ∞ ≤ ɛ p (1 + ‖b‖ ∞ ),<br />

∥<br />

τ k + sk ∥∥∥<br />

τ k − c ∞ ≤ ɛ d (1 + ‖c‖ ∞ ), and<br />

(<br />

≤<br />

ɛ g max<br />

1, min(∣ ∣c T x k∣ ∣ , ∣ ∣b T y k∣ )<br />

∣ )<br />

τ k ,<br />

(5.5)<br />

(x k , y k , s k )<br />

τ k<br />

is reported as the primal-dual optimal solution. <strong>The</strong> interpretation of (5.5) is that the optimizer is<br />

terminated if<br />

• xk<br />

τ k<br />

•<br />

is approximately primal feasible,<br />

( )<br />

y<br />

k<br />

, sk<br />

τ k τ k<br />

is approximately dual feasible, and<br />

• the duality gap is almost zero.<br />

On the other hand, if the trial solution satisfies<br />

−ɛ i c T x k ><br />

‖c‖ ∞<br />

∥ ∥ Ax<br />

k ∞<br />

max(1, ‖b‖ ∞ )<br />

then the problem is declared dual infeasible and x k is reported as a certificate of dual infeasibility.<br />

<strong>The</strong> motivation for this stopping criterion is as follows: First assume that ∥ ∥ Ax<br />

k ∥ ∥ ∞ = 0 ; then x k is<br />

an exact certificate of dual infeasibility. Next assume that this is not the case, i.e.<br />

and define<br />

∥ Ax<br />

k ∥ ∥ ∞ > 0,<br />

It is easy to verify that<br />

¯x := ɛ i<br />

max(1, ‖b‖ ∞ )<br />

‖Ax k ‖ ∞ ‖c‖ ∞<br />

x k .<br />

‖A¯x‖ ∞ = ɛ i<br />

max(1, ‖b‖ ∞ )<br />

‖c‖ ∞<br />

and − c T ¯x > 1,<br />

which shows ¯x is an approximate certificate of dual infeasibility where ɛ i controls the quality of the<br />

approximation. A smaller value means a better approximation.

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

Saved successfully!

Ooh no, something went wrong!