An O(nL) Infeasible-Interior-Point Algorithm for Linear Programming arxiv: v2 [math.oc] 29 Jun 2015
|
|
- Shavonne Taylor
- 5 years ago
- Views:
Transcription
1 An O(nL) Infeasible-Interior-Point Algorithm for Linear Programming arxiv: v [math.oc] 9 Jun 015 Yuagang Yang and Makoto Yamashita September 8, 018 Abstract In this paper, we propose an arc-search infeasible-interior-point algorithm. We show that this algorithm is polynomial and the polynomial bound is O(nL) which is at least as good as the best existing bound for infeasible-interior-point algorithms for linear programming. Keywords: polynomial algorithm, infeasible-interior-point method, linear programming. US NRC, Office of Research, 1 Church Street, Rockville, yaguang.yang@verizon.net. Department of Mathematical and Computing Sciences, Tokyo Institute of Technology. Makoto.Yamashita@is.titech.ac.jp. 1
2 1 Introduction Since Klee and Minty [1] showed that a simplex method for linear programming is not a polynomial algorithm, polynomial complexity bound has become a popular metric to measure the efficiency of optimization algorithms. Searching for polynomial algorithms for linear programming was a major research area of optimization between 1980 s and 1990 s after Khachiyan [] announced the first polynomial algorithm for linear programming. Although Khachiyan s algorithm was shown to be much less efficient in practice than the simplex method [3], Karmarkar s interior-point method [4] demonstrated the possibility of existence of efficient polynomial algorithms. For feasible starting point, people quickly established polynomial bounds for various interior-point algorithms [5, 6, 7, 8, 9, 10]. The lowest bound of these algorithms is O( nl) which has not been improved for more than two decades. To have an efficient implementation for interior-point algorithms, Mehrotra [11] and Lustig et. al. [1] realized that higher-order method and infeasible starting point are two necessary improvements. However, algorithms with either one of these features had poorer complexity bounds than O( nl). Monteiro, Adler, and Resende [13] showed that a higher-order algorithm starting from a feasible point has the polynomial bound O(nL). For infeasible-interior-point method, Zhang [14], Mizuno [15], and Miao [16] established polynomiality for several different algorithms(none of them is a higher-order algorithm). The best complexity bound O(nL) for infeasible interior-point methods has not been changed since the eary of 1990 s. Recently, Yang[17, 18] showed that for a higher-order interior-point method starting from a feasible point, the polynomial bound can be improved to O( nl) by using an arc-search method. Very recently, Yang et. al. [19] used the same idea and proposed a polynomial arc-search infeasible-interior-point algorihtm with a complexity bound of O(n 5 4L). In this paper, we show that for higher-order infeasible-interior-point method using arc-search, the polynomial bound can be improved to O(nL), which is a bound at least as good as the best bound of existing infeasible-interior-point algorithms. The remainder of the paper is organized as follows. Section describes the problem. Section 3 provides an infeasible-predictor-corrector algorithm. Section 4 proves its polynomiality. Section 5 summarizes the conclusions. Problem Descriptions The standard form of linear programming in this paper is given as follows: min c T x, subject to Ax = b, x 0, (1) where A R m n, b R m, c R n are given, and x R n is the vector to be optimized. Associated with the linear programming is the dual programming that is also presented in the standard form: max b T y, subject to A T y+s = c, s 0, ()
3 where dual variable vector y R m, and dual slack vector s R n. We use S to denote the set of all the optimal solutions (x,y,s ) of (1) and (). It is well known that x R n is an optimal solution of (1) if and only if x, y, and s satisfy the following KKT conditions Ax = b, A T y+s = c, (x,s) 0, x i s i = 0, i = 1,...,n. (3a) (3b) (3c) (3d) To simplify the notation, we will denote Hadamard (element-wise) product of two vectors x and s by x s, the element-wise division of the two vectors by s 1 x, or x s 1, or x s if min s i > 0, the Euclidean norm of x by x, the infinite norm of x by x, the identity matrix of any dimension by I, the vector of all ones with appropriate dimension by e, the block column vectors, for example, [x T,s T ] T by (x,s). For x R n, we will denote a related diagonal matrix by X R n n whose diagonal elements are the components of the vector x. Finally, we define an initial vector point of a sequence by x 0, an initial scalar point of a sequence by µ 0, the vector point after the kth iteration by x k, the scalar point after the kth iteration by µ k. Let r k b = Axk b, r k c = AT y k +s k c. (4a) (4b) Given a strictly positive current point (x k,s k ) > 0, the infeasible-predictor-corrector algorithm is to find the solution of (1) approximately along a curve C(t) defined by the following system Ax(t) b = tr k b A T y(t)+s(t) c = tr k c x(t) s(t) = tx k s k (x(t),s(t)) > 0, (5a) (5b) (5c) (5d) where t (0,1]. As t 0, (x(t),y(t),s(t)) appraches the solution of (1). Since C(t) is not easy to obtain, we will use an ellipse E [0] in the n+m dimensional space to approximate the curve defined by (5), where E is given by E = {(x(α),y(α),s(α)) : (x(α),y(α),s(α)) = acos(α)+ bsin(α)+ c}, (6) a R n+m and b R n+m are the axes of the ellipse, and c R n+m is the center of the ellipse. Taking the derivatives of (5) gives A 0 0 ẋ r k 0 A T b I ẏ = r k S k 0 X k c, (7) ṡ x k s k 3
4 A A T I S k 0 X k ẍ ÿ s = 0 0 ẋ ṡ. (8) We require the ellipse to pass the same point (x k,y k,s k ) on C(t) and to have the same derivatives given by (7) and (8). The ellipse is given in [17, 18] as Theorem.1 Let (x(α),y(α),s(α)) be an arc defined by (6) passing through a point (x,y,s) E C(t), and its first and second derivatives at (x,y,s) be (ẋ,ẏ,ṡ) and (ẍ,ÿ, s) which are defined by (7) and (8). Then, the ellipse approximation of (5) is given by x(α) = x ẋsin(α)+ẍ(1 cos(α)). (9) y(α) = y ẏsin(α)+ÿ(1 cos(α)). (10) s(α) = s ṡsin(α)+ s(1 cos(α)). (11) 3 Infeasible predictor-corrector algorithm We denote the duality measure by and define the set of neighborhood by µ = xt s n, (1) N(θ) := {(x,s) (x,s) > 0, x s µe θµ}. (13) The proposed algorithm searches an optimizer along the ellipse while staying inside N(θ). Algorithm Data: A, b, c, θ (0, + ], ǫ > 0, initial point (x0,y 0,s 0 ) N(θ). for iteration k = 1,,... Step 1: If stop. Otherwise continue. r k b = Axk b ǫ, r k c = AT y k +s k c ǫ, µ k ǫ, (14a) (14b) (14c) (x k,s k ) > 0. (14d) Step : Solve the linear systems of equations (7) and (8) to get (ẋ,ẏ,ṡ) and (ẍ,ÿ, s). 4
5 Step 3: Find the smallest positive ᾱ (0,π/] such that for all α (0,ᾱ], (x(α),s(α)) > 0 and (x(α) s(α)) (1 sin(α))µ k e θ(1 sin(α))µ k. (15) Set (to simplify the notation, we use α in stead of ᾱ in the rest of the paper) (x(α),y(α),s(α)) = (x k,y k,s k ) (ẋ,ẏ,ṡ)sin(α)+(ẍ,ÿ, s)(1 cos(α)). (16) Step 4: Calculate ( x, y, s) by solving A 0 0 x 0 0 A T I y = 0 S(α) 0 X(α) s (1 sin(α))µ k e x(α) s(α) Update and end (for). (17) (x k+1,y k+1,s k+1 ) = (x(α),y(α),s(α))+( x, y, s) (18) Step 5: Set k +1 k. Go back to Step 1. s k+1 µ k+1 = xk+1t. (19) n In the rest of this section, we will show (1) r k b 0, rk c 0, and µ k 0; () there exist α (0,π/] such that (x(α),s(α)) > 0 and (15) holds; (3) (x k,s k ) N(θ). It is easy to show that r k b, rk c, and µ k decrease at the same rate in every iteration. Lemma 3.1 r k+1 b = r k b (1 sin(α)), rk+1 c = r k c (1 sin(α)), µ k+1 = µ k (1 sin(α)). (0) Proof: Using (4), (18), (16), (8), and (7), we have r k+1 b r k b = A(x k+1 x k ) = A(x(α)+ x x k ) = A(x k ẋsin(α) x k ) = Aẋsin(α) = r k b sin(α). This shows the first relation. The second relation follows a similar derivation. From (17), it holds that ( x) T s = ( x) T ( A T y) = (A x) T y = 0. Using (18), we have x k+1t s k+1 = (x(α)+ x) T (s(α)+ s) = x(α) T s(α)+x(α) T s+s(α) T x = x(α) T s(α)+(1 sin(α))µ k n x(α) T s(α) = (1 sin(α))µ k n. Dividing both sides by n proves the last relation. 5
6 Clearly, if sin(α) = 1 (α = π ), we will find the optimal solution (allowing some x i = 0 and/or s j = 0) in one step, which is rarely the case. Therefore, from now on, we assume α (0, π ). We will use the following lemma of [16]. Lemma 3. Let ( x, s) be given by (17). Then 1 x s 4 (X(α)S(α)) (x(α) s(α) µk+1 e). (1) Theorem 3.1 Let (x(α),y(α),s(α)) and (x k,s k ) N(θ). Then, for all k 0 (i) there is an α > 0, such that (x(α),s(α)) > 0 and (15) holds. (ii) if θ 1 +, then (xk+1,s k+1 ) N(θ) for all the iterations. Proof: Using 1 cos(α) 1 cos (α) = sin (α) and (x k,y k,s k ) N(θ), we have x(α) s(α) (1 sin(α))µ k e = (x k s k µ k e)(1 sin(α))+(ẍ s ẋ ṡ)(1 cos(α)) (ẍ ṡ+ẋ s)sin(α)(1 cos(α)) θµ k (1 sin(α))+( ẍ s + ẋ ṡ )sin 4 (α) +( ẍ ṡ + ẋ s )sin 3 (α). () Clearly, if ( ) ( ) + ẋ ṡ q(α) := ẍ s sin 4 + ẍ ṡ (α)+ ẋ s sin 3 (α)+θµ k sin(α) θµ k 0, (3) then, (15) holds. Indeed, since q(0) = θµ < 0, by continuity, there exist α > 0 such that (3) holds. This shows that (15) holds. From (15), we have x i (α)s i (α) (1 θ)(1 sin(α))µ k > 0, θ [0,0.5) and α [0,π/). This shows (x(α), s(α)) > 0. Therefore, we finish part (i). Furthermore, from Lemma 3.1, (15) is now equivalent to x(α) s(α) µ k+1 e θµ k+1. Using (18), (17), Lemmas 3.1 and 3., and part (i) of this theorem, we have x k+1 s k+1 µ k+1 e = (x(α)+ x) (s(α)+ s) µ k+1 e 1 = x s 4 (X(α)S(α)) (x(α) s(α) µk+1 e) x(α) s(α) µ k+1 e 4 min i x i (α)s i (α) (θ) µ k+1 4(1 θ)µ k+1 θ (1 θ) µ k+1. (4) 6
7 Itiseasytocheck thatforθ , θ 1 θ holds, therefore, forθ (1 θ) +, we have x k+1 s k+1 µ k+1 e θµ k+1. Wenow show that (x k+1,s k+1 ) > 0. Let x k+1 (t) = x(α)+t x ands k+1 (t) = s(α)+t s. Then, x k+1 (0) = x(α) and x k+1 (1) = x k+1. Since x k+1 (t) s k+1 (t) = (x(α)+t x) (s(α)+t s) = x(α) s(α)+t(x(α) s+s(α) x)+t x s, using (17), (15), (4), and the assumption that θ 1 +, we have x k+1 (t) s k+1 (t) µ k+1 e = (1 t)(x(α) s(α) µ k+1 e)+t x s θ (1 t)θµ k+1 +t 1 θ µ k+1 ((1 t)+t )θµ k+1 := f(t)θµ k+1. (5) The function f(t) is a monotonical decreasing function of t [0,1], and f(0) =. This proves x k+1 (t) s k+1 (t) µ k+1 e θµ k+1. Therefore, x k+1 i (t)s k+1 i (t) (1 θ)µ k+1 > 0 for all t [0,1], which means (x k+1 s k+1 ) > 0. This finishes the proof of part (ii). This theorem indicates that the proposed algorithm is well-defined. 4 Polynomiality The analysis follows similar ideas in many existing literatures, such as [16, ]. Let the initial point be selected to satisfy (x 0,s 0 ) N(θ), x ρx 0, s ρs 0, (x,y,s ) S, (6) where ρ 1. Let ω f and ω o be the quality of the initial point which are the distances from feasibility and optimility given by ω f = min x,y,s {max{ (X0 ) 1 (x x 0 ), (S 0 ) 1 (s s 0 ) } Ax = b, A T y+s = c}. (7) and {max{x T ω 0 s 0 s Tx 0 = min x,y,s x 0T s0, x 0T s 0,1} (x,y,s ) S}. (8) Let ωp r and ωr d be the ratios of the feasibility and the total complementarity defined by ωp r = Ax0 b, x 0T s 0 (9a) ωd r = AT y 0 +s 0 c. x 0T s 0 (9b) 7
8 In view of Lemma 3.1, we have that Ax k b = ω r px kt s k, A T y k +s k c = ω r d xkt s k. (30a) (30b) Invoking Lemma 3.3 of [] for λ p = λ d = ξ = 1 and (7), we have the following two Lemmas [16]. Lemma 4.1 Let (ẋ,ṡ) be defined by (7), and D k = (X k ) 1 (S k ) 1. Then max{ (D k ) 1 ẋ, (D k )ṡ } (x k s k ) 1 +ω f (1+ω o (x k ) T s k ). (31) min i (x k i sk i )1 Lemma 4. Let (x 0,s 0 ) be defined by (6). Then ω f ρ, ω o ρ. (3) This leads to the following lemma. Lemma 4.3 Let (ẋ,ṡ) be defined by (7). Then, there exists a positive constant C 0, independent of n, such that max{ (D k ) 1 ẋ, (D k )ṡ } C 0 n(xk ) T s k. (33) Proof: First, it is easy to see (x k s k ) 1 = i x k i sk i = (x k ) T s k. (34) Since (x k,s k ) N(θ), we have min i (x k is k i) (1 θ)µ k = (1 θ) (xk ) T s k. Therefore, (x k ) T s k min i (x k i sk i )1 n(x k ) T s k n (1 θ). (35) Substituting (34) and (35) into (31) and using Lemma 4. prove (33) with C 0 = 1 + ρ(1+ρ) (1 θ) 1+ ωf (1+ω o ) (1 θ). 8
9 From Lemma 4.3, we can establish several useful inequalities. The following simple facts will be used several times. Let u and v be two vectors, then u v = ( )( ) (u i v i ) u i vi = u v. (36) i i If u and v satisfy u T v = 0, then, max{ u, v } u + v = u+v, (37) i and (see [3, Lemma 5.3]) u v 3 u+v. (38) Lemma 4.4 Let (ẋ,ṡ) and (ẍ, s) be defined by (7) and (8), respectively. Then, there exist positive constants C 1, C, C 3, and C 4, independent of n, such that ẋ ṡ C 1 n µ k, (39) ẍ s C n 4 µ k, (40) max{ (D k ) 1 ẍ, (D k ) s } C 3 n µ k, (41) max{ ẍ ṡ, ẋ s } C 4 n 3 µ k (4) Proof: First, using (36) and Lemma 4.3, we have ẋ ṡ = (D k ) 1 ẋ (D k )ṡ (D k ) 1 ẋ (D k )ṡ C 0n(x k ) T s k := C 1 n µ k. (43) Second, using (38), (8), (39), and (34), we have ẍ s = (D k ) 1 ẍ (D k ) s 3 (D k ) 1 ẍ+(d k ) s 3 (XS) 1 (ẋ ṡ) n ( ) = 1 ẋi ṡ i n = 1 (ẋ i ṡ i ) xi si x i=1 i=1 i s i n 1 i=1 (ẋ iṡ i ) min i=1,...,n x i s i 1 ẋ ṡ 1 C1n 4 µ k (1 θ)µ k (1 θ)µ k C = 1 1 n4 µ k 1 θ := C n 4 µ k. (44) Third, using (37), (8), and (39), we have = max{ (D k ) 1 ẍ, (D k ) s } (D k ) 1 ẍ+(d k ) s (XS) 1 4C (ẋ ṡ) 1 n4 µ k := C3 1 θ n4 µ k. (45) 9
10 Taking square root on both sides proves (41). Finally, using (36), (41), and Lemma 4.3, we have Similarly, we can show This finishes the proof. ẍ ṡ = (D k ) 1 ẍ (D k )ṡ (D k ) 1 ẍ (D k )ṡ (C 3 n µ k )(C 0 n µ k ) := C 4 n 3 µ k. (46) ẋ s C 4 n 3 µ k. (47) Now we are ready to estimate a conservative bound for sin(α). Lemma 4.5 Let (x k,y k,s k ) be generated by Algorithm 3.1. Then, sin(α) obtained in Step 3 satisfies the following inequality. where C = max{1,c 1 3 4,(C 1 +C ) 1 4}. sin(α) θ Cn, (48) Proof: Let sin(α) = θ. In view of (3) and Lemma 4.4, we have Cn q(α) µ k ((C 1 +C )n 4 sin 4 (α)+c 4 n 3 sin 3 (α)+θsin(α) θ) := µ k p(α) ( (C1 +C )θ 4 µ k + C ) 4θ 3 + θ 16C 4 8C 3 Cn θ ( ) θ 4 µ k 16 + θ3 4 + θ θ 0. Since p(α) is a monotonic function of sin(α), for all sin(α) θ, the above inequalities Cn hold (the last inequality holds because of θ 1). Therefore, for all sin(α) θ, the Cn inequality (15) holds. This finishes the proof. Remark 4.1 It is worthwhile to point out that the constant C depends on C 0 which depends on ρ, but ρ is an unknown before we find the solution. Also, we can always find a better steplength sin(α) by solving the quartic q(α) = 0 and the calculation of the roots for a quartic polynomial is deterministic, negligible, and independent to n [4, 5]. Following the standard argument developed in [3], we have the main theorem. Theorem 4.1 Let (x k,y k,s k ) be generated by Algorithm 3.1 with an initial point given by (6). For any ǫ > 0, the algorithm will terminate with (x k,y k,s k ) satisfying (14) in at most O(nL) iterations, where L = max{ln((x 0 ) T s 0 /ǫ),ln( r 0 b /ǫ),ln( r 0 c /ǫ)}. Proof: In view of Lemma 3.1, r k b, rk c, and µ k decrease at the same rate (1 sin(α)) in every iteration. Using the Lemma 4.5 and [3, Theorem 3.] proves the claim. 10
11 5 Conclusions We proposed an infeasible-interior-point algorithm that searches the optimizer along an ellipse that approximates the central-path. We showed that the proposed algorithm is polynomial and that the polynomial bound is at least as good as the best existing bound for infeasible-interior-point algorithms for linear programming. References [1] V. Klee and G. Minty, How good is the simplex algorithm? In: Shisha, O., (eds.) Inequalities, Vol. III, pp , Academic Press, (197). [] L. Khachiyan, A polynomial algorithm in linear programming, Doklady Akademiia Nauk SSSR, 4, (1979). [3] R. Bland, D. Goldfarb, and M. Todd, The ellipsoid method: A survey, Operations Research, 9, (1981). [4] N. Karmarkar, A new polynomial-time algorithm for linear programming, Combinatorica, 4, (1984). [5] J. Renegar, A polynomial-time algorithm based on Newton s method, for linear programming, Mathematical Programming, 40, (1988). [6] M. Kojima, S. Mizuno, and A. Yoshise, A polynomial-time algorithm for a class of linear complementarity problem, Mathematical Programming, 44, 1-6 (1989). [7] M. Kojima, S. Mizuno, and A. Yoshise, A primal-dual interior point algorithm for linear programming, in Progress in Mathematical Programming: Interior-point and Related Methods, N. Megiddo, ed., Springer-Verlag, New York, [8] C. Gonzaga, Polynomial affine algorithms for linear programming, Mathematical Programming, 49, 6-1 (1990). [9] Y. Ye, An O(n 3 L) potential reduction algorithm for linear programming, Mathematical Programming, 50, (1991). [10] S. Mizuno, M. Todd, and Y. Ye, On adaptive step primal-dual interior-point algorithms for linear programming, Mathematics of Operations Research, 18, (1993). [11] S. Mehrotra, On the implementation of a primal-dual interior point method, SIAM Journal on Optimization,, (199). [1] I. Lustig, R. Marsten, and D. Shannon, On implementing Mehrotra s predictorcorrector interior-point method for linear programming, SIAM Journal on Optimization,, (199). 11
12 [13] R. Monteiro, I. Adler, and M. Resende, A polynomial-time primal-dual affine scaling algorithm for linear and convex quadratic programming and its power series extension, Mathematics of Operations Research, 15, (1990). [14] Y. Zhang, On the convergence of a class of infeasible-interior-point methods for the horizontal linear Complementarity problem, SIAM Journal on Optimization, 4, 08-7 (1994). [15] S. Mizuno, Polynomiality of the Kojima-Megiddo-Mizuno infeasible interior point algorithm for linear programming, Mathematical Programming, 67, (1994) [16] J. Miao, Two infeasible interior-point predictor-corrector algorithms for linear programming, SIAM Journal on Optimization, 6, (1996). [17] Y. Yang, A Polynomial Arc-Search Interior-Point Algorithm for Linear Programming, Journal of Optimization Theory and Applications, 158, (013). [18] Y. Yang, A polynomial arc-search interior-point algorithm for convex quadratic programming, European Journal of Operational Research, 15, 5-38 (011). [19] X. Yang, Y. Zhang, and H. Liu, A wide neighborhood infeasible-interior-point method with arc-search for linear programming, Journal of Applied Mathematics and Computing, 1-17 (015). [0] M. P. Do Carmo, Differential Geometry of Curves and Surfaces, Prentice-Hall, New Jersey, (1976). [1] S. Mizuno, A new polynomial time method for a linear complemenetarity problem, Mathematical Programming 56, (199). [] M. Kojima, Basic lemmas in polynomial-time infeasible-interior-point methods for linear programming, Annals of Operations Research, 6, 1-8 (1996). [3] S. Wright, Primal-Dual Interior-Point Methods, SIAM, Philadelphia, (1997). [4] S.L. Shmakov, A universal method of solving quartic equations, International Journal of Pure and Applied Mathematics, 71 (), (011). [5] Y. Yang and Z. Zhou, An analytic solution to Wahba s problem, Aerospace Science and Technology, 30, (013) 1
A WIDE NEIGHBORHOOD PRIMAL-DUAL INTERIOR-POINT ALGORITHM WITH ARC-SEARCH FOR LINEAR COMPLEMENTARITY PROBLEMS 1. INTRODUCTION
J Nonlinear Funct Anal 08 (08), Article ID 3 https://doiorg/0395/jnfa083 A WIDE NEIGHBORHOOD PRIMAL-DUAL INTERIOR-POINT ALGORITHM WITH ARC-SEARCH FOR LINEAR COMPLEMENTARITY PROBLEMS BEIBEI YUAN, MINGWANG
More informationEnlarging neighborhoods of interior-point algorithms for linear programming via least values of proximity measure functions
Enlarging neighborhoods of interior-point algorithms for linear programming via least values of proximity measure functions Y B Zhao Abstract It is well known that a wide-neighborhood interior-point algorithm
More informationA Generalized Homogeneous and Self-Dual Algorithm. for Linear Programming. February 1994 (revised December 1994)
A Generalized Homogeneous and Self-Dual Algorithm for Linear Programming Xiaojie Xu Yinyu Ye y February 994 (revised December 994) Abstract: A generalized homogeneous and self-dual (HSD) infeasible-interior-point
More informationLecture 10. Primal-Dual Interior Point Method for LP
IE 8534 1 Lecture 10. Primal-Dual Interior Point Method for LP IE 8534 2 Consider a linear program (P ) minimize c T x subject to Ax = b x 0 and its dual (D) maximize b T y subject to A T y + s = c s 0.
More informationAn Infeasible Interior-Point Algorithm with full-newton Step for Linear Optimization
An Infeasible Interior-Point Algorithm with full-newton Step for Linear Optimization H. Mansouri M. Zangiabadi Y. Bai C. Roos Department of Mathematical Science, Shahrekord University, P.O. Box 115, Shahrekord,
More informationOn Superlinear Convergence of Infeasible Interior-Point Algorithms for Linearly Constrained Convex Programs *
Computational Optimization and Applications, 8, 245 262 (1997) c 1997 Kluwer Academic Publishers. Manufactured in The Netherlands. On Superlinear Convergence of Infeasible Interior-Point Algorithms for
More informationInterior Point Methods for Mathematical Programming
Interior Point Methods for Mathematical Programming Clóvis C. Gonzaga Federal University of Santa Catarina, Florianópolis, Brazil EURO - 2013 Roma Our heroes Cauchy Newton Lagrange Early results Unconstrained
More informationInterior Point Methods in Mathematical Programming
Interior Point Methods in Mathematical Programming Clóvis C. Gonzaga Federal University of Santa Catarina, Brazil Journées en l honneur de Pierre Huard Paris, novembre 2008 01 00 11 00 000 000 000 000
More informationCCO Commun. Comb. Optim.
Communications in Combinatorics and Optimization Vol. 3 No., 08 pp.5-70 DOI: 0.049/CCO.08.580.038 CCO Commun. Comb. Optim. An infeasible interior-point method for the P -matrix linear complementarity problem
More informationA FULL-NEWTON STEP INFEASIBLE-INTERIOR-POINT ALGORITHM COMPLEMENTARITY PROBLEMS
Yugoslav Journal of Operations Research 25 (205), Number, 57 72 DOI: 0.2298/YJOR3055034A A FULL-NEWTON STEP INFEASIBLE-INTERIOR-POINT ALGORITHM FOR P (κ)-horizontal LINEAR COMPLEMENTARITY PROBLEMS Soodabeh
More informationOptimization: Then and Now
Optimization: Then and Now Optimization: Then and Now Optimization: Then and Now Why would a dynamicist be interested in linear programming? Linear Programming (LP) max c T x s.t. Ax b αi T x b i for i
More informationOn Mehrotra-Type Predictor-Corrector Algorithms
On Mehrotra-Type Predictor-Corrector Algorithms M. Salahi, J. Peng, T. Terlaky April 7, 005 Abstract In this paper we discuss the polynomiality of Mehrotra-type predictor-corrector algorithms. We consider
More informationA New Class of Polynomial Primal-Dual Methods for Linear and Semidefinite Optimization
A New Class of Polynomial Primal-Dual Methods for Linear and Semidefinite Optimization Jiming Peng Cornelis Roos Tamás Terlaky August 8, 000 Faculty of Information Technology and Systems, Delft University
More informationA path following interior-point algorithm for semidefinite optimization problem based on new kernel function. djeffal
Journal of Mathematical Modeling Vol. 4, No., 206, pp. 35-58 JMM A path following interior-point algorithm for semidefinite optimization problem based on new kernel function El Amir Djeffal a and Lakhdar
More informationInterior-Point Methods
Interior-Point Methods Stephen Wright University of Wisconsin-Madison Simons, Berkeley, August, 2017 Wright (UW-Madison) Interior-Point Methods August 2017 1 / 48 Outline Introduction: Problems and Fundamentals
More informationA Full Newton Step Infeasible Interior Point Algorithm for Linear Optimization
A Full Newton Step Infeasible Interior Point Algorithm for Linear Optimization Kees Roos e-mail: C.Roos@tudelft.nl URL: http://www.isa.ewi.tudelft.nl/ roos 37th Annual Iranian Mathematics Conference Tabriz,
More informationCS711008Z Algorithm Design and Analysis
CS711008Z Algorithm Design and Analysis Lecture 8 Linear programming: interior point method Dongbo Bu Institute of Computing Technology Chinese Academy of Sciences, Beijing, China 1 / 31 Outline Brief
More informationA full-newton step feasible interior-point algorithm for P (κ)-lcp based on a new search direction
Croatian Operational Research Review 77 CRORR 706), 77 90 A full-newton step feasible interior-point algorithm for P κ)-lcp based on a new search direction Behrouz Kheirfam, and Masoumeh Haghighi Department
More informationResearch Note. A New Infeasible Interior-Point Algorithm with Full Nesterov-Todd Step for Semi-Definite Optimization
Iranian Journal of Operations Research Vol. 4, No. 1, 2013, pp. 88-107 Research Note A New Infeasible Interior-Point Algorithm with Full Nesterov-Todd Step for Semi-Definite Optimization B. Kheirfam We
More informationAn Infeasible Interior Point Method for the Monotone Linear Complementarity Problem
Int. Journal of Math. Analysis, Vol. 1, 2007, no. 17, 841-849 An Infeasible Interior Point Method for the Monotone Linear Complementarity Problem Z. Kebbiche 1 and A. Keraghel Department of Mathematics,
More informationNumerical Optimization
Linear Programming - Interior Point Methods Computer Science and Automation Indian Institute of Science Bangalore 560 012, India. NPTEL Course on Example 1 Computational Complexity of Simplex Algorithm
More informationA full-newton step infeasible interior-point algorithm for linear complementarity problems based on a kernel function
Algorithmic Operations Research Vol7 03) 03 0 A full-newton step infeasible interior-point algorithm for linear complementarity problems based on a kernel function B Kheirfam a a Department of Mathematics,
More informationPrimal-Dual Interior-Point Methods. Javier Peña Convex Optimization /36-725
Primal-Dual Interior-Point Methods Javier Peña Convex Optimization 10-725/36-725 Last time: duality revisited Consider the problem min x subject to f(x) Ax = b h(x) 0 Lagrangian L(x, u, v) = f(x) + u T
More informationA full-newton step infeasible interior-point algorithm for linear programming based on a kernel function
A full-newton step infeasible interior-point algorithm for linear programming based on a kernel function Zhongyi Liu, Wenyu Sun Abstract This paper proposes an infeasible interior-point algorithm with
More informationA Strongly Polynomial Simplex Method for Totally Unimodular LP
A Strongly Polynomial Simplex Method for Totally Unimodular LP Shinji Mizuno July 19, 2014 Abstract Kitahara and Mizuno get new bounds for the number of distinct solutions generated by the simplex method
More informationChapter 14 Linear Programming: Interior-Point Methods
Chapter 14 Linear Programming: Interior-Point Methods In the 1980s it was discovered that many large linear programs could be solved efficiently by formulating them as nonlinear problems and solving them
More informationA tight iteration-complexity upper bound for the MTY predictor-corrector algorithm via redundant Klee-Minty cubes
A tight iteration-complexity upper bound for the MTY predictor-corrector algorithm via redundant Klee-Minty cubes Murat Mut Tamás Terlaky Department of Industrial and Systems Engineering Lehigh University
More informationChapter 6 Interior-Point Approach to Linear Programming
Chapter 6 Interior-Point Approach to Linear Programming Objectives: Introduce Basic Ideas of Interior-Point Methods. Motivate further research and applications. Slide#1 Linear Programming Problem Minimize
More informationc 2005 Society for Industrial and Applied Mathematics
SIAM J. OPTIM. Vol. 15, No. 4, pp. 1147 1154 c 2005 Society for Industrial and Applied Mathematics A NOTE ON THE LOCAL CONVERGENCE OF A PREDICTOR-CORRECTOR INTERIOR-POINT ALGORITHM FOR THE SEMIDEFINITE
More informationInterior Point Methods for Linear Programming: Motivation & Theory
School of Mathematics T H E U N I V E R S I T Y O H F E D I N B U R G Interior Point Methods for Linear Programming: Motivation & Theory Jacek Gondzio Email: J.Gondzio@ed.ac.uk URL: http://www.maths.ed.ac.uk/~gondzio
More informationA Redundant Klee-Minty Construction with All the Redundant Constraints Touching the Feasible Region
A Redundant Klee-Minty Construction with All the Redundant Constraints Touching the Feasible Region Eissa Nematollahi Tamás Terlaky January 5, 2008 Abstract By introducing some redundant Klee-Minty constructions,
More informationA Second-Order Path-Following Algorithm for Unconstrained Convex Optimization
A Second-Order Path-Following Algorithm for Unconstrained Convex Optimization Yinyu Ye Department is Management Science & Engineering and Institute of Computational & Mathematical Engineering Stanford
More informationInterior Point Methods. We ll discuss linear programming first, followed by three nonlinear problems. Algorithms for Linear Programming Problems
AMSC 607 / CMSC 764 Advanced Numerical Optimization Fall 2008 UNIT 3: Constrained Optimization PART 4: Introduction to Interior Point Methods Dianne P. O Leary c 2008 Interior Point Methods We ll discuss
More informationA PREDICTOR-CORRECTOR PATH-FOLLOWING ALGORITHM FOR SYMMETRIC OPTIMIZATION BASED ON DARVAY'S TECHNIQUE
Yugoslav Journal of Operations Research 24 (2014) Number 1, 35-51 DOI: 10.2298/YJOR120904016K A PREDICTOR-CORRECTOR PATH-FOLLOWING ALGORITHM FOR SYMMETRIC OPTIMIZATION BASED ON DARVAY'S TECHNIQUE BEHROUZ
More informationPolynomiality of Linear Programming
Chapter 10 Polynomiality of Linear Programming In the previous section, we presented the Simplex Method. This method turns out to be very efficient for solving linear programmes in practice. While it is
More informationSemidefinite Programming
Chapter 2 Semidefinite Programming 2.0.1 Semi-definite programming (SDP) Given C M n, A i M n, i = 1, 2,..., m, and b R m, the semi-definite programming problem is to find a matrix X M n for the optimization
More informationFollowing The Central Trajectory Using The Monomial Method Rather Than Newton's Method
Following The Central Trajectory Using The Monomial Method Rather Than Newton's Method Yi-Chih Hsieh and Dennis L. Bricer Department of Industrial Engineering The University of Iowa Iowa City, IA 52242
More informationinforms DOI /moor.xxxx.xxxx c 200x INFORMS
MATHEMATICS OF OPERATIONS RESEARCH Vol. xx, No. x, Xxxxxxx 2x, pp. xxx xxx ISSN 364-765X EISSN 526-547 x xxx xxx informs DOI.287/moor.xxxx.xxxx c 2x INFORMS A New Complexity Result on Solving the Markov
More informationA Simpler and Tighter Redundant Klee-Minty Construction
A Simpler and Tighter Redundant Klee-Minty Construction Eissa Nematollahi Tamás Terlaky October 19, 2006 Abstract By introducing redundant Klee-Minty examples, we have previously shown that the central
More informationfrom the primal-dual interior-point algorithm (Megiddo [16], Kojima, Mizuno, and Yoshise
1. Introduction The primal-dual infeasible-interior-point algorithm which we will discuss has stemmed from the primal-dual interior-point algorithm (Megiddo [16], Kojima, Mizuno, and Yoshise [7], Monteiro
More informationOn well definedness of the Central Path
On well definedness of the Central Path L.M.Graña Drummond B. F. Svaiter IMPA-Instituto de Matemática Pura e Aplicada Estrada Dona Castorina 110, Jardim Botânico, Rio de Janeiro-RJ CEP 22460-320 Brasil
More information1. Introduction A number of recent papers have attempted to analyze the probabilistic behavior of interior point algorithms for linear programming. Ye
Probabilistic Analysis of an Infeasible-Interior-Point Algorithm for Linear Programming Kurt M. Anstreicher 1, Jun Ji 2, Florian A. Potra 3, and Yinyu Ye 4 Final Revision June, 1998 Abstract We consider
More informationAN INTERIOR POINT METHOD, BASED ON RANK-ONE UPDATES, Jos F. Sturm 1 and Shuzhong Zhang 2. Erasmus University Rotterdam ABSTRACT
October 13, 1995. Revised November 1996. AN INTERIOR POINT METHOD, BASED ON RANK-ONE UPDATES, FOR LINEAR PROGRAMMING Jos F. Sturm 1 Shuzhong Zhang Report 9546/A, Econometric Institute Erasmus University
More informationA SECOND ORDER MEHROTRA-TYPE PREDICTOR-CORRECTOR ALGORITHM FOR SEMIDEFINITE OPTIMIZATION
J Syst Sci Complex (01) 5: 1108 111 A SECOND ORDER MEHROTRA-TYPE PREDICTOR-CORRECTOR ALGORITHM FOR SEMIDEFINITE OPTIMIZATION Mingwang ZHANG DOI: 10.1007/s1144-01-0317-9 Received: 3 December 010 / Revised:
More informationOn Generalized Primal-Dual Interior-Point Methods with Non-uniform Complementarity Perturbations for Quadratic Programming
On Generalized Primal-Dual Interior-Point Methods with Non-uniform Complementarity Perturbations for Quadratic Programming Altuğ Bitlislioğlu and Colin N. Jones Abstract This technical note discusses convergence
More informationA Second Full-Newton Step O(n) Infeasible Interior-Point Algorithm for Linear Optimization
A Second Full-Newton Step On Infeasible Interior-Point Algorithm for Linear Optimization H. Mansouri C. Roos August 1, 005 July 1, 005 Department of Electrical Engineering, Mathematics and Computer Science,
More informationA new primal-dual path-following method for convex quadratic programming
Volume 5, N., pp. 97 0, 006 Copyright 006 SBMAC ISSN 00-805 www.scielo.br/cam A new primal-dual path-following method for convex quadratic programming MOHAMED ACHACHE Département de Mathématiques, Faculté
More informationPenalty and Barrier Methods General classical constrained minimization problem minimize f(x) subject to g(x) 0 h(x) =0 Penalty methods are motivated by the desire to use unconstrained optimization techniques
More informationPrimal-dual relationship between Levenberg-Marquardt and central trajectories for linearly constrained convex optimization
Primal-dual relationship between Levenberg-Marquardt and central trajectories for linearly constrained convex optimization Roger Behling a, Clovis Gonzaga b and Gabriel Haeser c March 21, 2013 a Department
More informationAn E cient A ne-scaling Algorithm for Hyperbolic Programming
An E cient A ne-scaling Algorithm for Hyperbolic Programming Jim Renegar joint work with Mutiara Sondjaja 1 Euclidean space A homogeneous polynomial p : E!R is hyperbolic if there is a vector e 2E such
More informationPrimal-Dual Interior-Point Methods. Ryan Tibshirani Convex Optimization /36-725
Primal-Dual Interior-Point Methods Ryan Tibshirani Convex Optimization 10-725/36-725 Given the problem Last time: barrier method min x subject to f(x) h i (x) 0, i = 1,... m Ax = b where f, h i, i = 1,...
More informationPredictor-corrector methods for sufficient linear complementarity problems in a wide neighborhood of the central path
Copyright information to be inserted by the Publishers Predictor-corrector methods for sufficient linear complementarity problems in a wide neighborhood of the central path Florian A. Potra and Xing Liu
More informationImproved Full-Newton-Step Infeasible Interior- Point Method for Linear Complementarity Problems
Georgia Southern University Digital Commons@Georgia Southern Mathematical Sciences Faculty Publications Mathematical Sciences, Department of 4-2016 Improved Full-Newton-Step Infeasible Interior- Point
More informationConvergence Analysis of Inexact Infeasible Interior Point Method. for Linear Optimization
Convergence Analysis of Inexact Infeasible Interior Point Method for Linear Optimization Ghussoun Al-Jeiroudi Jacek Gondzio School of Mathematics The University of Edinburgh Mayfield Road, Edinburgh EH9
More information2 The SDP problem and preliminary discussion
Int. J. Contemp. Math. Sciences, Vol. 6, 2011, no. 25, 1231-1236 Polynomial Convergence of Predictor-Corrector Algorithms for SDP Based on the KSH Family of Directions Feixiang Chen College of Mathematics
More informationA Full-Newton Step O(n) Infeasible Interior-Point Algorithm for Linear Optimization
A Full-Newton Step On) Infeasible Interior-Point Algorithm for Linear Optimization C. Roos March 4, 005 February 19, 005 February 5, 005 Faculty of Electrical Engineering, Computer Science and Mathematics
More informationImproved Full-Newton Step O(nL) Infeasible Interior-Point Method for Linear Optimization
J Optim Theory Appl 2010) 145: 271 288 DOI 10.1007/s10957-009-9634-0 Improved Full-Newton Step OnL) Infeasible Interior-Point Method for Linear Optimization G. Gu H. Mansouri M. Zangiabadi Y.Q. Bai C.
More informationCurvature as a Complexity Bound in Interior-Point Methods
Lehigh University Lehigh Preserve Theses and Dissertations 2014 Curvature as a Complexity Bound in Interior-Point Methods Murat Mut Lehigh University Follow this and additional works at: http://preserve.lehigh.edu/etd
More informationLocal Self-concordance of Barrier Functions Based on Kernel-functions
Iranian Journal of Operations Research Vol. 3, No. 2, 2012, pp. 1-23 Local Self-concordance of Barrier Functions Based on Kernel-functions Y.Q. Bai 1, G. Lesaja 2, H. Mansouri 3, C. Roos *,4, M. Zangiabadi
More informationPrimal-Dual Interior-Point Methods. Ryan Tibshirani Convex Optimization
Primal-Dual Interior-Point Methods Ryan Tibshirani Convex Optimization 10-725 Given the problem Last time: barrier method min x subject to f(x) h i (x) 0, i = 1,... m Ax = b where f, h i, i = 1,... m are
More informationPrimal-Dual Interior-Point Methods by Stephen Wright List of errors and typos, last updated December 12, 1999.
Primal-Dual Interior-Point Methods by Stephen Wright List of errors and typos, last updated December 12, 1999. 1. page xviii, lines 1 and 3: (x, λ, x) should be (x, λ, s) (on both lines) 2. page 6, line
More informationNew Infeasible Interior Point Algorithm Based on Monomial Method
New Infeasible Interior Point Algorithm Based on Monomial Method Yi-Chih Hsieh and Dennis L. Bricer Department of Industrial Engineering The University of Iowa, Iowa City, IA 52242 USA (January, 1995)
More informationLecture 14 Barrier method
L. Vandenberghe EE236A (Fall 2013-14) Lecture 14 Barrier method centering problem Newton decrement local convergence of Newton method short-step barrier method global convergence of Newton method predictor-corrector
More informationAn interior-point gradient method for large-scale totally nonnegative least squares problems
An interior-point gradient method for large-scale totally nonnegative least squares problems Michael Merritt and Yin Zhang Technical Report TR04-08 Department of Computational and Applied Mathematics Rice
More informationInterior Point Methods for LP
11.1 Interior Point Methods for LP Katta G. Murty, IOE 510, LP, U. Of Michigan, Ann Arbor, Winter 1997. Simplex Method - A Boundary Method: Starting at an extreme point of the feasible set, the simplex
More informationIntroduction to integer programming II
Introduction to integer programming II Martin Branda Charles University in Prague Faculty of Mathematics and Physics Department of Probability and Mathematical Statistics Computational Aspects of Optimization
More informationProperties of a Simple Variant of Klee-Minty s LP and Their Proof
Properties of a Simple Variant of Klee-Minty s LP and Their Proof Tomonari Kitahara and Shinji Mizuno December 28, 2 Abstract Kitahara and Mizuno (2) presents a simple variant of Klee- Minty s LP, which
More informationA priori bounds on the condition numbers in interior-point methods
A priori bounds on the condition numbers in interior-point methods Florian Jarre, Mathematisches Institut, Heinrich-Heine Universität Düsseldorf, Germany. Abstract Interior-point methods are known to be
More information4TE3/6TE3. Algorithms for. Continuous Optimization
4TE3/6TE3 Algorithms for Continuous Optimization (Algorithms for Constrained Nonlinear Optimization Problems) Tamás TERLAKY Computing and Software McMaster University Hamilton, November 2005 terlaky@mcmaster.ca
More informationPrimal-dual IPM with Asymmetric Barrier
Primal-dual IPM with Asymmetric Barrier Yurii Nesterov, CORE/INMA (UCL) September 29, 2008 (IFOR, ETHZ) Yu. Nesterov Primal-dual IPM with Asymmetric Barrier 1/28 Outline 1 Symmetric and asymmetric barriers
More informationNew Interior Point Algorithms in Linear Programming
AMO - Advanced Modeling and Optimization, Volume 5, Number 1, 2003 New Interior Point Algorithms in Linear Programming Zsolt Darvay Abstract In this paper the abstract of the thesis New Interior Point
More informationNonsymmetric potential-reduction methods for general cones
CORE DISCUSSION PAPER 2006/34 Nonsymmetric potential-reduction methods for general cones Yu. Nesterov March 28, 2006 Abstract In this paper we propose two new nonsymmetric primal-dual potential-reduction
More informationPrimal-Dual Interior-Point Methods for Linear Programming based on Newton s Method
Primal-Dual Interior-Point Methods for Linear Programming based on Newton s Method Robert M. Freund March, 2004 2004 Massachusetts Institute of Technology. The Problem The logarithmic barrier approach
More informationA PRIMAL-DUAL INTERIOR POINT ALGORITHM FOR CONVEX QUADRATIC PROGRAMS. 1. Introduction Consider the quadratic program (PQ) in standard format:
STUDIA UNIV. BABEŞ BOLYAI, INFORMATICA, Volume LVII, Number 1, 01 A PRIMAL-DUAL INTERIOR POINT ALGORITHM FOR CONVEX QUADRATIC PROGRAMS MOHAMED ACHACHE AND MOUFIDA GOUTALI Abstract. In this paper, we propose
More informationDetecting Infeasibility in Infeasible-Interior-Point Methods for Optimization
Detecting Infeasibility in Infeasible-Interior-Point Methods for Optimization M. J. Todd January 16, 2003 Abstract We study interior-point methods for optimization problems in the case of infeasibility
More informationInterior-point algorithm for linear optimization based on a new trigonometric kernel function
Accepted Manuscript Interior-point algorithm for linear optimization based on a new trigonometric kernel function Xin Li, Mingwang Zhang PII: S0-0- DOI: http://dx.doi.org/./j.orl.0.0.0 Reference: OPERES
More informationLOWER BOUNDS FOR THE MAXIMUM NUMBER OF SOLUTIONS GENERATED BY THE SIMPLEX METHOD
Journal of the Operations Research Society of Japan Vol 54, No 4, December 2011, pp 191 200 c The Operations Research Society of Japan LOWER BOUNDS FOR THE MAXIMUM NUMBER OF SOLUTIONS GENERATED BY THE
More informationA Constraint-Reduced Variant of Mehrotra s Predictor-Corrector Algorithm
A Constraint-Reduced Variant of Mehrotra s Predictor-Corrector Algorithm Luke B. Winternitz, Stacey O. Nicholls, André L. Tits, Dianne P. O Leary September 24, 2007 Abstract Consider linear programs in
More informationCorrector-predictor methods for monotone linear complementarity problems in a wide neighborhood of the central path
Mathematical Programming manuscript No. will be inserted by the editor) Florian A. Potra Corrector-predictor methods for monotone linear complementarity problems in a wide neighborhood of the central path
More informationFull-Newton-Step Interior-Point Method for the Linear Complementarity Problems
Georgia Southern University Digital Commons@Georgia Southern Electronic Theses and Dissertations Graduate Studies, Jack N. Averitt College of Summer 2011 Full-Newton-Step Interior-Point Method for the
More informationAn interior-point trust-region polynomial algorithm for convex programming
An interior-point trust-region polynomial algorithm for convex programming Ye LU and Ya-xiang YUAN Abstract. An interior-point trust-region algorithm is proposed for minimization of a convex quadratic
More informationInterior-Point Methods for Linear Optimization
Interior-Point Methods for Linear Optimization Robert M. Freund and Jorge Vera March, 204 c 204 Robert M. Freund and Jorge Vera. All rights reserved. Linear Optimization with a Logarithmic Barrier Function
More informationCOMPARATIVE STUDY BETWEEN LEMKE S METHOD AND THE INTERIOR POINT METHOD FOR THE MONOTONE LINEAR COMPLEMENTARY PROBLEM
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume LIII, Number 3, September 2008 COMPARATIVE STUDY BETWEEN LEMKE S METHOD AND THE INTERIOR POINT METHOD FOR THE MONOTONE LINEAR COMPLEMENTARY PROBLEM ADNAN
More informationarxiv: v2 [math.oc] 11 Jan 2018
A one-phase interior point method for nonconvex optimization Oliver Hinder, Yinyu Ye January 12, 2018 arxiv:1801.03072v2 [math.oc] 11 Jan 2018 Abstract The work of Wächter and Biegler [40] suggests that
More informationA polynomial time interior point path following algorithm for LCP based on Chen Harker Kanzow smoothing techniques
Math. Program. 86: 9 03 (999) Springer-Verlag 999 Digital Object Identifier (DOI) 0.007/s007990056a Song Xu James V. Burke A polynomial time interior point path following algorithm for LCP based on Chen
More informationNew stopping criteria for detecting infeasibility in conic optimization
Optimization Letters manuscript No. (will be inserted by the editor) New stopping criteria for detecting infeasibility in conic optimization Imre Pólik Tamás Terlaky Received: March 21, 2008/ Accepted:
More informationLecture 1. 1 Conic programming. MA 796S: Convex Optimization and Interior Point Methods October 8, Consider the conic program. min.
MA 796S: Convex Optimization and Interior Point Methods October 8, 2007 Lecture 1 Lecturer: Kartik Sivaramakrishnan Scribe: Kartik Sivaramakrishnan 1 Conic programming Consider the conic program min s.t.
More informationImproved Full-Newton-Step Infeasible Interior- Point Method for Linear Complementarity Problems
Georgia Southern University Digital Commons@Georgia Southern Electronic Theses & Dissertations Graduate Studies, Jack N. Averitt College of Summer 2015 Improved Full-Newton-Step Infeasible Interior- Point
More informationLecture 15: October 15
10-725: Optimization Fall 2012 Lecturer: Barnabas Poczos Lecture 15: October 15 Scribes: Christian Kroer, Fanyi Xiao Note: LaTeX template courtesy of UC Berkeley EECS dept. Disclaimer: These notes have
More informationSolving Obstacle Problems by Using a New Interior Point Algorithm. Abstract
Solving Obstacle Problems by Using a New Interior Point Algorithm Yi-Chih Hsieh Department of Industrial Engineering National Yunlin Polytechnic Institute Huwei, Yunlin 6308 Taiwan and Dennis L. Bricer
More informationThe Q Method for Symmetric Cone Programmin
The Q Method for Symmetric Cone Programming The Q Method for Symmetric Cone Programmin Farid Alizadeh and Yu Xia alizadeh@rutcor.rutgers.edu, xiay@optlab.mcma Large Scale Nonlinear and Semidefinite Progra
More informationarxiv: v1 [math.oc] 21 Jan 2019
STATUS DETERMINATION BY INTERIOR-POINT METHODS FOR CONVEX OPTIMIZATION PROBLEMS IN DOMAIN-DRIVEN FORM MEHDI KARIMI AND LEVENT TUNÇEL arxiv:1901.07084v1 [math.oc] 21 Jan 2019 Abstract. We study the geometry
More information10 Numerical methods for constrained problems
10 Numerical methods for constrained problems min s.t. f(x) h(x) = 0 (l), g(x) 0 (m), x X The algorithms can be roughly divided the following way: ˆ primal methods: find descent direction keeping inside
More informationRachid Benouahboun 1 and Abdelatif Mansouri 1
RAIRO Operations Research RAIRO Oper. Res. 39 25 3 33 DOI:.5/ro:252 AN INTERIOR POINT ALGORITHM FOR CONVEX QUADRATIC PROGRAMMING WITH STRICT EQUILIBRIUM CONSTRAINTS Rachid Benouahboun and Abdelatif Mansouri
More informationLecture 3 Interior Point Methods and Nonlinear Optimization
Lecture 3 Interior Point Methods and Nonlinear Optimization Robert J. Vanderbei April 16, 2012 Machine Learning Summer School La Palma http://www.princeton.edu/ rvdb Example: Basis Pursuit Denoising L
More informationIMPLEMENTING THE NEW SELF-REGULAR PROXIMITY BASED IPMS
IMPLEMENTING THE NEW SELF-REGULAR PROXIMITY BASED IPMS IMPLEMENTING THE NEW SELF-REGULAR PROXIMITY BASED IPMS By Xiaohang Zhu A thesis submitted to the School of Graduate Studies in Partial Fulfillment
More informationMotivating examples Introduction to algorithms Simplex algorithm. On a particular example General algorithm. Duality An application to game theory
Instructor: Shengyu Zhang 1 LP Motivating examples Introduction to algorithms Simplex algorithm On a particular example General algorithm Duality An application to game theory 2 Example 1: profit maximization
More informationLinear Programming: Simplex
Linear Programming: Simplex Stephen J. Wright 1 2 Computer Sciences Department, University of Wisconsin-Madison. IMA, August 2016 Stephen Wright (UW-Madison) Linear Programming: Simplex IMA, August 2016
More informationFull Newton step polynomial time methods for LO based on locally self concordant barrier functions
Full Newton step polynomial time methods for LO based on locally self concordant barrier functions (work in progress) Kees Roos and Hossein Mansouri e-mail: [C.Roos,H.Mansouri]@ewi.tudelft.nl URL: http://www.isa.ewi.tudelft.nl/
More informationInterior Point Methods for Nonlinear Optimization
Interior Point Methods for Nonlinear Optimization Imre Pólik 1 and Tamás Terlaky 2 1 School of Computational Engineering and Science, McMaster University, Hamilton, Ontario, Canada, imre@polik.net 2 School
More information