1 / 22

Chapter 9. Interior Point Methods

Chapter 9. Interior Point Methods. Three major variants Affine scaling algorithm - easy concept, good performance Potential reduction algorithm - poly time Path following algorithm - poly time, good performance, theoretically elegant . 9.4 The primal path following algorithm.

tarak
Download Presentation

Chapter 9. Interior Point Methods

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Chapter 9. Interior Point Methods • Three major variants • Affine scaling algorithm - easy concept, good performance • Potential reduction algorithm - poly time • Path following algorithm - poly time, good performance, theoretically elegant

  2. 9.4 The primal path following algorithm • min max Nonnegativity makes the problem difficult, hence use barrier function in the objective and consider unconstrained problem ( in the affine space • Barrier function: if for some Solve min s.t. (9.15) is strictly convex, hence has unique min point if min exists.

  3. ex) min , s. t. ,  min at 0 1

  4. min s.t. Let is optimal solution given . , when varies, is called the central path ( hence the name path following ) It can be shown that optimal solution to LP. When , is called the analytic center. • For dual problem, the barrier problem is max , s.t. (9.16) ( equivalent to min , minimizing convex function)

  5. Figure 9.4: The central path and the analytic center analytic center

  6. Results from nonlinear programming (NLP) min s.t., , : , all twice continuously differentiable ( gradient is given as a column vector) • Thm (Karush 1939, Kuhn-Tucker 1951, first order necessary optimality condition) If is a local minimum for (NLP) and some conditions (called constraint qualification) hold at , then there exist such that (1) (2) (3)

  7. Remark: • (2) is CS conditions and it implies that for non-active constraint . • (1) says is a nonnegative linear combination of for active constraints and (compare to strong duality theorem in p. 173 and its Figure ) • CS conditions for LP are KKT conditions • KKT conditions are necessary conditions for optimality, but it is also sufficient in some situations. One case is when objective function is convex and constraints are linear, which includes our barrier problem.

  8. Deriving KKT for barrier problem: • min , s.t. ( is row vector of , expressed as a column vector and , is the vector having 1 in all components.) Using (Lagrangian) multiplier for , we get (ignoring the sign of ) Note that and . () If we define , KKT becomes (), where .

  9. For dual barrier problem, min s.t. () ( is column vector of and is unit vector.) Using (Lagrangian) multiplier for , we get Now , hence we have the conditions , which is the same conditions we obtained from the primal barrier function.

  10. The conditions are given in the text as (9.17) where . Note that when , they are primal, dual feasibility and complementary slackness conditions. • Lemma 9.5: If and satisfy conditions (9.17), then they are optimal solutions to problems (9.15) and (9.16)

  11. Pf) Let and satisfy (9.17), and let be an arbitrary vector that satisfies and . Then attains min at . equality holds iff Hence for all feasible . In particular, is the unique optimal solution and . Similarly for and for dual barrier problem. 

  12. Primal path following algorithm • Starting from some and primal and dual feasible find solution of the barrier problem iteratively while . • To solve the barrier problem, we use quadratic approximation (2nd order Taylor expansion) of the barrier function and use the minimum of the approximate function as the next iterates. Taylor expansion is Also need to satisfy 

  13. Using KKT, solution to this problem is The duality gap is Hence stop the algorithm if Need a scheme to have an initial feasible solution (see text)

  14. The primal path following algorithm 1. (Initialization) Start with some primal and dual feasible and set . 2. (Optimality test) If , stop; else go to Step 3. 3. Let 4. (Computation of directions) Solve the linear system , , for and . 5. (Update of solutions) Let , , . 6. Let and go to Step 2.

  15. 9.5 The primal-dual path following algorithm • Find Newton directions both in the primal and dual space. Instead of finding min of quadratic approximation of barrier function, it finds the solution for KKT system. (9.26) System of nonlinear equations because of the last ones. • Let . Want such that We use first order Taylor approximation around Here is the Jacobian matrix whose -th element is given by

  16. Try to find that satisfies . We then set . The direction is called a Newton direction. Here F(z) is given by  This is equivalent to (9.28) (9.29) (9.30)

  17. Solution to the previous system is where Also limit the step length to ensure

  18. The primal-dual path following algorithm 1. (Initialization) Start with some feasible and set . 2. (Optimality test) If ,stop; else go to Step 3. 3. (Computation of Newton directions) Let , Solve the linear system (9.28) – (9.30) for and 4. (Find step lengths) Let

  19. (continued) 5. (Solution update) Update the solution vectors according to 6. Let and go to Step 2.

  20. Infeasible primal-dual path following methods • A variation of primal-dual path following. Starts from which is not necessarily feasible for either the primal or the dual, i.e. and/or . Iteration same as the primal-dual path following except feasibility not maintained in each iteration. Excellent performance.

  21. Self-dual method • Alternative method to find initial feasible solution w/o using big-M. Given an initial possibly infeasible point with and , consider the problem minimize subject to (9.33) where . This LP is self-dual. Note that is a feasible interior solution to (9.33)

  22. Since both the primal and dual are feasible, they have optimal solutions and the optimal value is 0. • Primal-dual path following method finds an optimal solution that satisfies ( satisfies strict complementarity ) • Can find optimal solution or determine unboundedness depending on the value of . (see Thm 9.8) • Running time : worst case : observed :

More Related