Solving conic optimization problems with MOSEK
|
|
- Michael Horton
- 5 years ago
- Views:
Transcription
1 Solving conic optimization problems with MOSEK Erling D. Andersen MOSEK ApS, Fruebjergvej 3, Box 16, 2100 Copenhagen, Denmark. WWW: August 8,
2 2 / 42
3 Outline Outline MOSEK What is MOSEK. Conic problems MOSEK can solve. Interfacing with MOSEK. Algorithms employed in MOSEK. Implementation.. Summary and conclusions. 3 / 42
4 MOSEK Outline MOSEK MOSEK is a software package for large scale optimization. Version 1 released April Version 7 released May Linear and conic quadratic (+ mixed-integer). Conic quadratic optimization (+ mixed-integer). Semi-definite optimization starting from version 7. Convex(functional) optimization. C, JAVA,.NET and Python APIs. AMPL, AIMMS, GAMS, MATLAB, and R interfaces. Free for academic use. See 4 / 42
5 5 / 42
6 The primal problem The primal problem min nj=1 c j x j + n j=1 Cj, X j st li c n j=1 a ij x j + n j=1 Āij, X j u c i, i = 1,...,m, lj x x j u x j, j = 1,...,n, x K, X j 0, j = 1,..., n. Explanation: A,B := tr(a T B). x j is a scalar variable. Xj is a square matrix variable. K represents conic quadratic constraints. Xj 0 represents X j = X T j and X j is PSD. C j andāj are assumed to be symmetric. 6 / 42
7 The primal problem An example ofk: { x 1 x 27 x 6 andx 3 x 8 x 9 }. Models the quadratic cones. (Both normal and rotated variant are allowed.) A scalar variable can only be member of one sub cone. 7 / 42
8 The primal problem Employs a primal formulation. Semi-definite addition is an extension on the linear and conic quadratic case. Makes it easy to extend the software. LMIs (=dual formulation) C j i Ā ij y i can be formulated. Requires explicit slacks. Normally C j andāij are sparse (is exploited). Dense but low rank in C j andāij can be exploited. 8 / 42
9 9 / 42
10 What is available What is available The MATLAB toolbox interface Fusion: A new interface Optimizer API. Matrix orientated. Flexible and efficient. Efficiency is more important than ease of use. Available in C, Java,.NET and Python. R and MATLAB toolbox. Matrix orientated. Fusion. Deals with variables and constraints. Good compromise between efficiency and ease of use. Available for Java, MATLAB,.NET and Python. 10 / 42
11 What is available The MATLAB toolbox interface Fusion: A new interface CVX. Yalmip. 11 / 42
12 The MATLAB toolbox interface What is available The MATLAB toolbox interface Fusion: A new interface An example: min x 1 + st x 1 + x 2 +x , X x 1 x 2 2 +x 3 2, X 1 0., X 1, X 1 = 1, = 0.5, 12 / 42
13 Source code What is available The MATLAB toolbox interface Fusion: A new interface prob.c = [1, 0, 0]; prob.bardim = [3]; prob.barc.subj = [1, 1, 1, 1, 1]; prob.barc.subk = [1, 2, 2, 3, 3]; prob.barc.subl = [1, 1, 2, 2, 3]; prob.barc.val = [2.0, 1.0, 2.0, 1.0, 2.0]; subj: Variable index. subk,subl: Local index i.e. ( C j ) kl. Only the lower triangular part is specified. 13 / 42
14 What is available The MATLAB toolbox interface Fusion: A new interface prob.a = sparse([1, 2, 2],... [1, 2, 3],... [1, 1, 1], 2, 3); prob.bara.subi = [1, 1, 1, 2, 2, 2, 2, 2, 2]; prob.bara.subj = [1, 1, 1, 1, 1, 1, 1, 1, 1]; prob.bara.subk = [1, 2, 3, 1, 2, 3, 2, 3, 3]; prob.bara.subl = [1, 2, 3, 1, 1, 1, 2, 2, 3]; prob.bara.val = [1, 1, 1, 1, 1, 1, 1, 1, 1]; prob.blc = [1.0, 0.5]; prob.buc = [1.0, 0.5]; 14 / 42
15 What is available The MATLAB toolbox interface Fusion: A new interface [r, res] = mosekopt( symbcon ); prob.cones.type = [res.symbcon.msk_ct_quad]; prob.cones.sub = [1, 2, 3]; prob.cones.subptr = [1]; [r,res] = mosekopt( minimize,prob); X = zeros(3); X([1,2,3,5,6,9]) = res.sol.itr.barx; X = X + tril(x,-1) ; x = res.sol.itr.xx; 15 / 42
16 What is available The MATLAB toolbox interface Fusion: A new interface Pros: Easy to explain. Efficient. Cons: Problem must be serialized. Small problem changes requires error prone changes to the code. Possible solutions: CVX (MATLAB only currently). YALMIP (MATLAB only). 16 / 42
17 Fusion: A new interface What is available The MATLAB toolbox interface Fusion: A new interface The basic idea of Fusion: min st c T x A i x+b i K i i. A linear objective. Several linear constraints involving the variables can be formulated. K i is a convex set. K i cannot have an arbitrary structure e.g. must be symmetric cones. Several variables (scalar and symmetric matrix). Not just one big giant matrix. (Contrast with SeDuMi). 17 / 42
18 Example: The nearest correlation matrix What is available The MATLAB toolbox interface Fusion: A new interface Given a symmetric matrixathen the nearest correlation matrix is given by X = arg min A X F. X 0,diag(X)=e Conic formulation exploiting symmetry of A X: min t st z 2 t, svec(a X) = z, diag(x) = e, X 0. Simple formulation on paper, but cumbersome to formulate in standard form as in MOSEK, SeDuMi, SDPA, / 42
19 Python Fusion def svec(e): N = e.get_shape().dim(0) S = Matrix.sparse(N * (N+1) / 2, N * N, range(n * (N+1) / 2), [ (i+j*n) for j in xrange(n) for i in xrange(j,n) ], [ (1.0 if i == j else 2**(0.5)) for j in xrange(n) for i in xrange(j,n) ]) return Expr.mul(S,Expr.reshape( e, N * N )) def nearestcorr(a): M = Model("NearestCorrelation") N = len(a) # Setting up the variables X = M.variable("X",Domain.inPSDCone(N)) # t > z _2 tz = M.variable("tz", Domain.inQCone(N*(N+1)/2+1)) t = tz.index(0) z = tz.slice(1,n*(n+1)/2+1) # svec (A-X) = z M.constraint( Expr.sub(svec(Expr.sub(DenseMatrix(A),X)), z), Domain.equalsTo(0.0) ) # diag(x) = e for i in range(n): M.constraint( X.index(i,i), Domain.equalsTo(1.0) ) # Objective: Minimize t M.objective(ObjectiveSense.Minimize, t) M.solve() 19 / 42
20 Fusion summary What is available The MATLAB toolbox interface Fusion: A new interface Easy to add and delete variables and constraints. Not a modeling language but an easy to use API. Designed for efficiency yet still convenient to use. Available for Java, MATLAB,.NET, and Python. An alternative to Yalmip and CVX. 20 / 42
21 21 / 42
22 Simplified notation Simplified notation The homogeneous model The central path Notation Schur complement Termination Practical details Primal problem: where min c T x st Ax = b, x 0, x 0 represents a mixture of linear, conic quadratic and semi-definite variables. Dual problem: max b T y st A T y +s = c, s / 42
23 The homogeneous model Simplified notation The homogeneous model The central path Notation Schur complement Termination Practical details Simplified homogeneous model: Notes: s 0 κ + 0 A T c A 0 b c T b T 0 (x;τ),(s;κ) 0. (x,y,s,κ,τ) = 0 is a feasible solution. All solutions satisfies: x T s+κτ = 0. x y τ = / 42
24 Simplified notation The homogeneous model The central path Notation Schur complement Termination Practical details Ifτ > 0, κ = 0 then(x,y,s)/τ is a primal-dual optimal solution, Ax = bτ, A T y +s = cτ, x T s = 0. Ifκ > 0, τ = 0 then the primal and/or dual is infeasible, Ax = 0, A T y +s = 0, c T x b T y < 0, Primal is infeasible ifb T y > 0, dual is infeasible ifc T x < 0. Comments: Find a solution such thatτ > 0 or κ > 0 if it exists. It exist under suitable regularity conditions. 24 / 42
25 The central path Simplified notation The homogeneous model The central path Notation Schur complement Termination Practical details Let (x (0) ;τ (0) ),(s (0) ;κ (0) ) 0 then the central path is defined as: s 0 A T c x A T y (0) +s (0) cτ (0) 0 + A 0 b y = γ Ax (0) +bτ (0) κ c T b T 0 τ c T x (0) b T y (0) +κ (0) x s = γµ (0) e, τκ = γµ (0), whereγ [0,1],µ (0) = (x(0) ) T s (0) +κ (0) τ (0) n+1 > 0 and e = svec(i n1 ). svec(i np ). 25 / 42
26 Notation Simplified notation The homogeneous model The central path Notation Schur complement Termination Practical details For symmetric matricesu S n we define svec(u) = (U 11, 2U 21,..., 2U n1,u 22, 2U 32,..., 2U n2,...,u n with inverse operation smat(u) = u 1 u 2 / 2 u n / 2 u 2 / 2 u n+1 u 2n 1 / 2... u n / 2 u n 1 / 2 u n(n+1), and a symmetric product u v := (1/2)(smat(u)smat(v) + smat(v)smat(u)). 26 / 42
27 The search direction Simplified notation The homogeneous model The central path Notation Schur complement Termination Practical details Ad x +bd τ = (γ 1)( Ax (0) +bτ (0) ) A T d y +d s cd τ = (γ 1)(A T y (0) +s (0) cτ (0) ) c T d x b T d y +d κ = (γ 1)(c T x (0) b T y (0) +κ (0) ) d x +Π (0) d s = γµ (0) (s (0) ) 1 x (0) d τ +τ (0) /κ (0) κ = γµ (0) /τ (0) κ (0) Residuals reduction: A(x+αd x )+b(τ +αd τ ) = (1 α(1 γ))( Ax (0) +bτ (0) ). Complementarity: (x (0) +αd x ) T (s (0) +αd s )+(κ (0) +α κ)(τ (0) +α τ) = (1 α(1 γ))((x (0) ) T s (0) +κ (0) τ (0) ). }{{} <1 Polynomial complexity for suitable choice ofπ, γ and α. 27 / 42
28 The scaling matrix Π is Nesterov-Todd scaling in MOSEK. Simplified notation The homogeneous model The central path Notation Schur complement Termination Practical details 28 / 42
29 Solving for the search direction Simplified notation The homogeneous model The central path Notation Schur complement Termination Practical details After block-elimination, we get a2 2 system: AΠA T d y +(AΠc b)d τ =, (AΠc+b)d y +(c T Π (0) c+κ (0) /τ (0) )d τ =. We factoraπa T = LL T and eliminate d y = L T L 1 ( (AΠc b)d τ + ). Then ( ) (AΠc+b) T L T L 1 (AΠc b)+(c T Πc+κ/τ) d τ = An extra (insignificant) solve compared to an infeasible method. 29 / 42
30 Forming the Schur complement Simplified notation The homogeneous model The central path Notation Schur complement Termination Practical details Schur complement computed as a sum of outer products: AΠA T = = p k=1 a T a T 1p.. a T m1... a T mp P T k ( P k A :,k Π k A T :,kp T k Π 1... ) P k, Π p a a m1.. a 1p... a mp wherep k is a permutation of A :,k. Complexity depends on choice ofp k. We compute onlytril(aπa T ), so by nnz(p k A :,k ) 1 nnz(p k A :,k ) 2... nnz(p k A :,k ) p, get a good reduction in complexity. More general than SDPA. 30 / 42
31 Forming the Schur complement Simplified notation The homogeneous model The central path Notation Schur complement Termination Practical details GivenΠ k = R k R T k then a T ikπ k a jk = R T ka ik R k,r T ka jk R k Using the idea in SeDuMi then = R k R T ka jk R k R T k = ˆMM +M T ˆMT, i.e., a symmetric (low-rank) product. A ik,r k R T ka jk R k R T k Sparse case: Few elements of ˆMM +M T ˆMT needed. Dense case: Many elements of ˆMM +M T ˆMT needed.. May exploit low-rank structure in dense A ik later. 31 / 42
32 Termination Simplified notation The homogeneous model The central path Notation Schur complement Termination Practical details Optimal case: A x(k) +b ε τ (k) p (1+ b ), A T x(k) + s(k) c ε τ (k) τ (k) d (1+ c ), ( ) ( min (x (k) ) T s (k) +τ (k) κ (k), ct x (k) bt y (k) ε (τ (k) ) 2 τ (k) τ (k) g max 1, Primal infeasible case: ε i b T y (k) > b max(1, c ) A T y (k) +s (k) ), τ (k) ct x (k) Dual infeasible case: ε i c T x (k) > c max(1, b ) Ax (k) All εs are user specified. 32 / 42
33 Practical details Simplified notation The homogeneous model The central path Notation Schur complement Termination Practical details Employs Mehrotra s predictor-corrector technique. Employs a presolve which can reduce problem size and complexity significantly. Primarily for the linear part. Scale the problem to improve conditioning of a problem. Exploit sparse linear algebra when applicable. Exploit tuned dense BLAS when applicable. Fully parallelized for linear and conic quadratic cases. 33 / 42
34 34 / 42
35 Benchmark environment Benchmark environment Optimized problems Result Profile Windows server CPU: Intel XEON Cores. 3.5 GHz. MOSEK: Only 2 threads allowed. Test problems: Taken from: Exclude the smallest and largest with respect to running time. 35 / 42
36 Optimized problems Benchmark environment Optimized problems Result Profile Name # con. # cone # var. # mat. var. G40 mb G40mc G48mc buck butcher cancer cphil foot hand inc inc mater neu3g nonc rabmo reimer rendl ros rose sensor sensor shmup shmup swissroll taha1a taha1b trto vibra yalsdp / 42
37 Result Benchmark environment Optimized problems Result Profile Name P. obj. # sig. fig. # iter time(s) G40 mb e G40mc e G48mc e buck e butcher e cancer e cphil e foot e hand e inc e inc e mater e neu3g e nonc e rabmo e reimer e rendl e ros e rose e sensor e sensor e shmup e shmup e swissroll e taha1a e taha1b e trto e vibra e yalsdp e / 42
38 Profile Benchmark environment Optimized problems Result Profile In % Name Time Prslv Fac. set. Shur. Chol Upd. G40 mb G40mc G48mc buck butcher cancer cphil foot hand inc inc mater neu3g nonc rabmo reimer rendl ros rose sensor sensor shmup shmup swissroll taha1a taha1b trto vibra yalsdp / 42
39 39 / 42
40 Current status Future work References MOSEK has been extended to handel semi-definite problems. A new modeling tool Fusion has been made available. MOSEK can solve large semi-definite problems (usual restrictions apply). 40 / 42
41 Future work Future work References Add warmstart as discussed by Y. Ye. Improve the parallelization (Exploit Intel Phi 1 tera flop card). Implement stability enhancements for the scaling matrix suggested by J. Sturm [1]. Use an iterative method for solution of the Newton equation system to improve numerical stability (not speed). Improve the presolve for the conic part. Exploit dense but low rank structure in C andā. Fusion: C++ variant. (Requires garbage collection like feature). 41 / 42
42 References [1] J. F. Sturm. Avoiding numerical cancellation in the interior point method for solving semidefinite programs. Math. Programming, 95(1): , Future work References 42 / 42
Semidefinite Optimization using MOSEK
Semidefinite Optimization using MOSEK Joachim Dahl ISMP Berlin, August 23, 2012 http://www.mose.com Introduction Introduction Semidefinite optimization Algorithms Results and examples Summary MOSEK is
More informationConvex Optimization : Conic Versus Functional Form
Convex Optimization : Conic Versus Functional Form Erling D. Andersen MOSEK ApS, Fruebjergvej 3, Box 16, DK 2100 Copenhagen, Blog: http://erlingdandersen.blogspot.com Linkedin: http://dk.linkedin.com/in/edandersen
More informationOn implementing a primal-dual interior-point method for conic quadratic optimization
On implementing a primal-dual interior-point method for conic quadratic optimization E. D. Andersen, C. Roos, and T. Terlaky December 18, 2000 Abstract Conic quadratic optimization is the problem of minimizing
More informationResearch overview. Seminar September 4, Lehigh University Department of Industrial & Systems Engineering. Research overview.
Research overview Lehigh University Department of Industrial & Systems Engineering COR@L Seminar September 4, 2008 1 Duality without regularity condition Duality in non-exact arithmetic 2 interior point
More information15. Conic optimization
L. Vandenberghe EE236C (Spring 216) 15. Conic optimization conic linear program examples modeling duality 15-1 Generalized (conic) inequalities Conic inequality: a constraint x K where K is a convex cone
More informationThe Ongoing Development of CSDP
The Ongoing Development of CSDP Brian Borchers Department of Mathematics New Mexico Tech Socorro, NM 87801 borchers@nmt.edu Joseph Young Department of Mathematics New Mexico Tech (Now at Rice University)
More informationIII. Applications in convex optimization
III. Applications in convex optimization nonsymmetric interior-point methods partial separability and decomposition partial separability first order methods interior-point methods Conic linear optimization
More informationWarmstarting the Homogeneous and Self-Dual Interior Point Method for Linear and Conic Quadratic Problems
Mathematical Programming Computation manuscript No. (will be inserted by the editor) Warmstarting the Homogeneous and Self-Dual Interior Point Method for Linear and Conic Quadratic Problems Anders Skajaa
More informationSemidefinite Programming, Combinatorial Optimization and Real Algebraic Geometry
Semidefinite Programming, Combinatorial Optimization and Real Algebraic Geometry assoc. prof., Ph.D. 1 1 UNM - Faculty of information studies Edinburgh, 16. September 2014 Outline Introduction Definition
More informationIMPLEMENTATION OF INTERIOR POINT METHODS
IMPLEMENTATION OF INTERIOR POINT METHODS IMPLEMENTATION OF INTERIOR POINT METHODS FOR SECOND ORDER CONIC OPTIMIZATION By Bixiang Wang, Ph.D. A Thesis Submitted to the School of Graduate Studies in Partial
More information18. Primal-dual interior-point methods
L. Vandenberghe EE236C (Spring 213-14) 18. Primal-dual interior-point methods primal-dual central path equations infeasible primal-dual method primal-dual method for self-dual embedding 18-1 Symmetric
More informationInfeasible Primal-Dual (Path-Following) Interior-Point Methods for Semidefinite Programming*
Infeasible Primal-Dual (Path-Following) Interior-Point Methods for Semidefinite Programming* Notes for CAAM 564, Spring 2012 Dept of CAAM, Rice University Outline (1) Introduction (2) Formulation & a complexity
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 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 informationIntroduction to Semidefinite Programs
Introduction to Semidefinite Programs Masakazu Kojima, Tokyo Institute of Technology Semidefinite Programming and Its Application January, 2006 Institute for Mathematical Sciences National University of
More informationSPARSE SECOND ORDER CONE PROGRAMMING FORMULATIONS FOR CONVEX OPTIMIZATION PROBLEMS
Journal of the Operations Research Society of Japan 2008, Vol. 51, No. 3, 241-264 SPARSE SECOND ORDER CONE PROGRAMMING FORMULATIONS FOR CONVEX OPTIMIZATION PROBLEMS Kazuhiro Kobayashi Sunyoung Kim Masakazu
More informationLecture 5. Theorems of Alternatives and Self-Dual Embedding
IE 8534 1 Lecture 5. Theorems of Alternatives and Self-Dual Embedding IE 8534 2 A system of linear equations may not have a solution. It is well known that either Ax = c has a solution, or A T y = 0, c
More informationFast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma
Fast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 8 September 2003 European Union RTN Summer School on Multi-Agent
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 informationInfeasible Primal-Dual (Path-Following) Interior-Point Methods for Semidefinite Programming*
Infeasible Primal-Dual (Path-Following) Interior-Point Methods for Semidefinite Programming* Yin Zhang Dept of CAAM, Rice University Outline (1) Introduction (2) Formulation & a complexity theorem (3)
More informationOn the solution of large-scale SDP problems by the modified barrier method using iterative solvers
Mathematical Programming manuscript No. (will be inserted by the editor) Michal Kočvara Michael Stingl On the solution of large-scale SDP problems by the modified barrier method using iterative solvers
More informationA direct formulation for sparse PCA using semidefinite programming
A direct formulation for sparse PCA using semidefinite programming A. d Aspremont, L. El Ghaoui, M. Jordan, G. Lanckriet ORFE, Princeton University & EECS, U.C. Berkeley A. d Aspremont, INFORMS, Denver,
More informationSecond-order cone programming
Outline Second-order cone programming, PhD Lehigh University Department of Industrial and Systems Engineering February 10, 2009 Outline 1 Basic properties Spectral decomposition The cone of squares The
More informationPENNON A Generalized Augmented Lagrangian Method for Convex NLP and SDP p.1/39
PENNON A Generalized Augmented Lagrangian Method for Convex NLP and SDP Michal Kočvara Institute of Information Theory and Automation Academy of Sciences of the Czech Republic and Czech Technical University
More informationOn the solution of large-scale SDP problems by the modified barrier method using iterative solvers
Mathematical Programming manuscript No. (will be inserted by the editor) Michal Kočvara Michael Stingl On the solution of large-scale SDP problems by the modified barrier method using iterative solvers
More informationDistributed Control of Connected Vehicles and Fast ADMM for Sparse SDPs
Distributed Control of Connected Vehicles and Fast ADMM for Sparse SDPs Yang Zheng Department of Engineering Science, University of Oxford Homepage: http://users.ox.ac.uk/~ball4503/index.html Talk at Cranfield
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 informationOn the solution of large-scale SDP problems by the modified barrier method using iterative solvers
Mathematical Programming manuscript No. (will be inserted by the editor) Michal Kočvara Michael Stingl On the solution of large-scale SDP problems by the modified barrier method using iterative solvers
More informationLecture 18: Optimization Programming
Fall, 2016 Outline Unconstrained Optimization 1 Unconstrained Optimization 2 Equality-constrained Optimization Inequality-constrained Optimization Mixture-constrained Optimization 3 Quadratic Programming
More informationLecture 17: Primal-dual interior-point methods part II
10-725/36-725: Convex Optimization Spring 2015 Lecture 17: Primal-dual interior-point methods part II Lecturer: Javier Pena Scribes: Pinchao Zhang, Wei Ma Note: LaTeX template courtesy of UC Berkeley EECS
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 informationSolving conic optimization problems via self-dual embedding and facial reduction: a unified approach
Solving conic optimization problems via self-dual embedding and facial reduction: a unified approach Frank Permenter Henrik A. Friberg Erling D. Andersen August 18, 216 Abstract We establish connections
More informationProblem structure in semidefinite programs arising in control and signal processing
Problem structure in semidefinite programs arising in control and signal processing Lieven Vandenberghe Electrical Engineering Department, UCLA Joint work with: Mehrdad Nouralishahi, Tae Roh Semidefinite
More informationA NEW SECOND-ORDER CONE PROGRAMMING RELAXATION FOR MAX-CUT PROBLEMS
Journal of the Operations Research Society of Japan 2003, Vol. 46, No. 2, 164-177 2003 The Operations Research Society of Japan A NEW SECOND-ORDER CONE PROGRAMMING RELAXATION FOR MAX-CUT PROBLEMS Masakazu
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 informationResearch Reports on Mathematical and Computing Sciences
ISSN 1342-284 Research Reports on Mathematical and Computing Sciences Exploiting Sparsity in Linear and Nonlinear Matrix Inequalities via Positive Semidefinite Matrix Completion Sunyoung Kim, Masakazu
More information4. Algebra and Duality
4-1 Algebra and Duality P. Parrilo and S. Lall, CDC 2003 2003.12.07.01 4. Algebra and Duality Example: non-convex polynomial optimization Weak duality and duality gap The dual is not intrinsic The cone
More informationThe Graph Realization Problem
The Graph Realization Problem via Semi-Definite Programming A. Y. Alfakih alfakih@uwindsor.ca Mathematics and Statistics University of Windsor The Graph Realization Problem p.1/21 The Graph Realization
More informationHow to generate weakly infeasible semidefinite programs via Lasserre s relaxations for polynomial optimization
CS-11-01 How to generate weakly infeasible semidefinite programs via Lasserre s relaxations for polynomial optimization Hayato Waki Department of Computer Science, The University of Electro-Communications
More informationFast ADMM for Sum of Squares Programs Using Partial Orthogonality
Fast ADMM for Sum of Squares Programs Using Partial Orthogonality Antonis Papachristodoulou Department of Engineering Science University of Oxford www.eng.ox.ac.uk/control/sysos antonis@eng.ox.ac.uk with
More informationApproximate Farkas Lemmas in Convex Optimization
Approximate Farkas Lemmas in Convex Optimization Imre McMaster University Advanced Optimization Lab AdvOL Graduate Student Seminar October 25, 2004 1 Exact Farkas Lemma Motivation 2 3 Future plans The
More informationEasy and not so easy multifacility location problems... (In 20 minutes.)
Easy and not so easy multifacility location problems... (In 20 minutes.) MINLP 2014 Pittsburgh, June 2014 Justo Puerto Institute of Mathematics (IMUS) Universidad de Sevilla Outline 1 Introduction (In
More informationA direct formulation for sparse PCA using semidefinite programming
A direct formulation for sparse PCA using semidefinite programming A. d Aspremont, L. El Ghaoui, M. Jordan, G. Lanckriet ORFE, Princeton University & EECS, U.C. Berkeley Available online at www.princeton.edu/~aspremon
More informationA semidefinite relaxation scheme for quadratically constrained quadratic problems with an additional linear constraint
Iranian Journal of Operations Research Vol. 2, No. 2, 20, pp. 29-34 A semidefinite relaxation scheme for quadratically constrained quadratic problems with an additional linear constraint M. Salahi Semidefinite
More informationSolving large Semidefinite Programs - Part 1 and 2
Solving large Semidefinite Programs - Part 1 and 2 Franz Rendl http://www.math.uni-klu.ac.at Alpen-Adria-Universität Klagenfurt Austria F. Rendl, Singapore workshop 2006 p.1/34 Overview Limits of Interior
More information6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC
6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC 2003 2003.09.02.10 6. The Positivstellensatz Basic semialgebraic sets Semialgebraic sets Tarski-Seidenberg and quantifier elimination Feasibility
More informationConic Linear Optimization and its Dual. yyye
Conic Linear Optimization and Appl. MS&E314 Lecture Note #04 1 Conic Linear Optimization and its Dual Yinyu Ye Department of Management Science and Engineering Stanford University Stanford, CA 94305, U.S.A.
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 informationSecond-Order Cone Program (SOCP) Detection and Transformation Algorithms for Optimization Software
and Second-Order Cone Program () and Algorithms for Optimization Software Jared Erickson JaredErickson2012@u.northwestern.edu Robert 4er@northwestern.edu Northwestern University INFORMS Annual Meeting,
More informationWhat s New in Active-Set Methods for Nonlinear Optimization?
What s New in Active-Set Methods for Nonlinear Optimization? Philip E. Gill Advances in Numerical Computation, Manchester University, July 5, 2011 A Workshop in Honor of Sven Hammarling UCSD Center for
More informationSection Notes 9. IP: Cutting Planes. Applied Math 121. Week of April 12, 2010
Section Notes 9 IP: Cutting Planes Applied Math 121 Week of April 12, 2010 Goals for the week understand what a strong formulations is. be familiar with the cutting planes algorithm and the types of cuts
More informationContinuous Optimisation, Chpt 9: Semidefinite Problems
Continuous Optimisation, Chpt 9: Semidefinite Problems Peter J.C. Dickinson DMMP, University of Twente p.j.c.dickinson@utwente.nl http://dickinson.website/teaching/2016co.html version: 21/11/16 Monday
More informationTamás Terlaky George N. and Soteria Kledaras 87 Endowed Chair Professor. Chair, Department of Industrial and Systems Engineering Lehigh University
5th SJOM Bejing, 2011 Cone Linear Optimization (CLO) From LO, SOCO and SDO Towards Mixed-Integer CLO Tamás Terlaky George N. and Soteria Kledaras 87 Endowed Chair Professor. Chair, Department of Industrial
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 informationImplementation and Evaluation of SDPA 6.0 (SemiDefinite Programming Algorithm 6.0)
Research Reports on Mathematical and Computing Sciences Series B : Operations Research Department of Mathematical and Computing Sciences Tokyo Institute of Technology 2-12-1 Oh-Okayama, Meguro-ku, Tokyo
More informationLinear Algebra and its Applications
Linear Algebra and its Applications 433 (2010) 1101 1109 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: www.elsevier.com/locate/laa Minimal condition number
More informationA SUFFICIENTLY EXACT INEXACT NEWTON STEP BASED ON REUSING MATRIX INFORMATION
A SUFFICIENTLY EXACT INEXACT NEWTON STEP BASED ON REUSING MATRIX INFORMATION Anders FORSGREN Technical Report TRITA-MAT-2009-OS7 Department of Mathematics Royal Institute of Technology November 2009 Abstract
More informationTamás Terlaky George N. and Soteria Kledaras 87 Endowed Chair Professor. Chair, Department of Industrial and Systems Engineering Lehigh University
BME - 2011 Cone Linear Optimization (CLO) From LO, SOCO and SDO Towards Mixed-Integer CLO Tamás Terlaky George N. and Soteria Kledaras 87 Endowed Chair Professor. Chair, Department of Industrial and Systems
More informationOn Two Measures of Problem Instance Complexity and their Correlation with the Performance of SeDuMi on Second-Order Cone Problems
2016 Springer International Publishing AG. Part of Springer Nature. http://dx.doi.org/10.1007/s10589-005-3911-0 On Two Measures of Problem Instance Complexity and their Correlation with the Performance
More informationIntroduction to Semidefinite Programming I: Basic properties a
Introduction to Semidefinite Programming I: Basic properties and variations on the Goemans-Williamson approximation algorithm for max-cut MFO seminar on Semidefinite Programming May 30, 2010 Semidefinite
More informationSMCP Documentation. Release Martin S. Andersen and Lieven Vandenberghe
SMCP Documentation Release 0.4.5 Martin S. Andersen and Lieven Vandenberghe May 18, 2018 Contents 1 Current release 3 2 Future releases 5 3 Availability 7 4 Authors 9 5 Feedback and bug reports 11 5.1
More informationA Homogeneous Interior-Point Algorithm for Nonsymmetric Convex Conic Optimization
Mathematical Programming manuscript No. (will be inserted by the editor) A Homogeneous Interior-Point Algorithm for Nonsymmetric Convex Conic Optimization Anders Skajaa Yinyu Ye Received: date / Accepted:
More informationConic Linear Programming. Yinyu Ye
Conic Linear Programming Yinyu Ye December 2004, revised January 2015 i ii Preface This monograph is developed for MS&E 314, Conic Linear Programming, which I am teaching at Stanford. Information, lecture
More informationThe Q Method for Second-Order Cone Programming
The Q Method for Second-Order Cone Programming Yu Xia Farid Alizadeh July 5, 005 Key words. Second-order cone programming, infeasible interior point method, the Q method Abstract We develop the Q method
More informationCourse Outline. FRTN10 Multivariable Control, Lecture 13. General idea for Lectures Lecture 13 Outline. Example 1 (Doyle Stein, 1979)
Course Outline FRTN Multivariable Control, Lecture Automatic Control LTH, 6 L-L Specifications, models and loop-shaping by hand L6-L8 Limitations on achievable performance L9-L Controller optimization:
More informationAdvances in Convex Optimization: Theory, Algorithms, and Applications
Advances in Convex Optimization: Theory, Algorithms, and Applications Stephen Boyd Electrical Engineering Department Stanford University (joint work with Lieven Vandenberghe, UCLA) ISIT 02 ISIT 02 Lausanne
More informationA CONIC DANTZIG-WOLFE DECOMPOSITION APPROACH FOR LARGE SCALE SEMIDEFINITE PROGRAMMING
A CONIC DANTZIG-WOLFE DECOMPOSITION APPROACH FOR LARGE SCALE SEMIDEFINITE PROGRAMMING Kartik Krishnan Advanced Optimization Laboratory McMaster University Joint work with Gema Plaza Martinez and Tamás
More informationSparse Matrix Theory and Semidefinite Optimization
Sparse Matrix Theory and Semidefinite Optimization Lieven Vandenberghe Department of Electrical Engineering University of California, Los Angeles Joint work with Martin S. Andersen and Yifan Sun Third
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 19: Stabilization via LMIs Optimization Optimization can be posed in functional form: min x F objective function : inequality constraints
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 informationEE 227A: Convex Optimization and Applications October 14, 2008
EE 227A: Convex Optimization and Applications October 14, 2008 Lecture 13: SDP Duality Lecturer: Laurent El Ghaoui Reading assignment: Chapter 5 of BV. 13.1 Direct approach 13.1.1 Primal problem Consider
More informationSpherical Euclidean Distance Embedding of a Graph
Spherical Euclidean Distance Embedding of a Graph Hou-Duo Qi University of Southampton Presented at Isaac Newton Institute Polynomial Optimization August 9, 2013 Spherical Embedding Problem The Problem:
More informationConstraint Reduction for Linear Programs with Many Constraints
Constraint Reduction for Linear Programs with Many Constraints André L. Tits Institute for Systems Research and Department of Electrical and Computer Engineering University of Maryland, College Park Pierre-Antoine
More information12. Interior-point methods
12. Interior-point methods Convex Optimization Boyd & Vandenberghe inequality constrained minimization logarithmic barrier function and central path barrier method feasibility and phase I methods complexity
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 informationOptimization Tools in an Uncertain Environment
Optimization Tools in an Uncertain Environment Michael C. Ferris University of Wisconsin, Madison Uncertainty Workshop, Chicago: July 21, 2008 Michael Ferris (University of Wisconsin) Stochastic optimization
More informationEE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 4
EE/ACM 150 - Applications of Convex Optimization in Signal Processing and Communications Lecture 4 Andre Tkacenko Signal Processing Research Group Jet Propulsion Laboratory April 12, 2012 Andre Tkacenko
More information1 A HOMOGENIZED CUTTING PLANE
OPTIMIZATION METHODS AND APPLICATIONS Published by KLUWER 1999 1 A HOMOGENIZED CUTTING PLANE METHOD TO SOLVE THE CONVEX FEASIBILITY PROBLEM E.D.Andersen,J.E.Mitchell,C.Roos,andT.Terlaky. Faculty of ITS,
More informationm i=1 c ix i i=1 F ix i F 0, X O.
What is SDP? for a beginner of SDP Copyright C 2005 SDPA Project 1 Introduction This note is a short course for SemiDefinite Programming s SDP beginners. SDP has various applications in, for example, control
More informationPOLYNOMIAL OPTIMIZATION WITH SUMS-OF-SQUARES INTERPOLANTS
POLYNOMIAL OPTIMIZATION WITH SUMS-OF-SQUARES INTERPOLANTS Sercan Yıldız syildiz@samsi.info in collaboration with Dávid Papp (NCSU) OPT Transition Workshop May 02, 2017 OUTLINE Polynomial optimization and
More informationA Parametric Simplex Algorithm for Linear Vector Optimization Problems
A Parametric Simplex Algorithm for Linear Vector Optimization Problems Birgit Rudloff Firdevs Ulus Robert Vanderbei July 9, 2015 Abstract In this paper, a parametric simplex algorithm for solving linear
More informationCroatian Operational Research Review (CRORR), Vol. 3, 2012
126 127 128 129 130 131 132 133 REFERENCES [BM03] S. Burer and R.D.C. Monteiro (2003), A nonlinear programming algorithm for solving semidefinite programs via low-rank factorization, Mathematical Programming
More informationSolving the normal equations system arising from interior po. linear programming by iterative methods
Solving the normal equations system arising from interior point methods for linear programming by iterative methods Aurelio Oliveira - aurelio@ime.unicamp.br IMECC - UNICAMP April - 2015 Partners Interior
More informationRobust and Optimal Control, Spring 2015
Robust and Optimal Control, Spring 2015 Instructor: Prof. Masayuki Fujita (S5-303B) G. Sum of Squares (SOS) G.1 SOS Program: SOS/PSD and SDP G.2 Duality, valid ineqalities and Cone G.3 Feasibility/Optimization
More informationParallel implementation of primal-dual interior-point methods for semidefinite programs
Parallel implementation of primal-dual interior-point methods for semidefinite programs Masakazu Kojima, Kazuhide Nakata Katsuki Fujisawa and Makoto Yamashita 3rd Annual McMaster Optimization Conference:
More information9. Numerical linear algebra background
Convex Optimization Boyd & Vandenberghe 9. Numerical linear algebra background matrix structure and algorithm complexity solving linear equations with factored matrices LU, Cholesky, LDL T factorization
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 informationPart 4: Active-set methods for linearly constrained optimization. Nick Gould (RAL)
Part 4: Active-set methods for linearly constrained optimization Nick Gould RAL fx subject to Ax b Part C course on continuoue optimization LINEARLY CONSTRAINED MINIMIZATION fx subject to Ax { } b where
More informationThe moment-lp and moment-sos approaches
The moment-lp and moment-sos approaches LAAS-CNRS and Institute of Mathematics, Toulouse, France CIRM, November 2013 Semidefinite Programming Why polynomial optimization? LP- and SDP- CERTIFICATES of POSITIVITY
More informationPrimal-Dual Geometry of Level Sets and their Explanatory Value of the Practical Performance of Interior-Point Methods for Conic Optimization
Primal-Dual Geometry of Level Sets and their Explanatory Value of the Practical Performance of Interior-Point Methods for Conic Optimization Robert M. Freund M.I.T. June, 2010 from papers in SIOPT, Mathematics
More informationSVM May 2007 DOE-PI Dianne P. O Leary c 2007
SVM May 2007 DOE-PI Dianne P. O Leary c 2007 1 Speeding the Training of Support Vector Machines and Solution of Quadratic Programs Dianne P. O Leary Computer Science Dept. and Institute for Advanced Computer
More informationAgenda. Interior Point Methods. 1 Barrier functions. 2 Analytic center. 3 Central path. 4 Barrier method. 5 Primal-dual path following algorithms
Agenda Interior Point Methods 1 Barrier functions 2 Analytic center 3 Central path 4 Barrier method 5 Primal-dual path following algorithms 6 Nesterov Todd scaling 7 Complexity analysis Interior point
More informationl p -Norm Constrained Quadratic Programming: Conic Approximation Methods
OUTLINE l p -Norm Constrained Quadratic Programming: Conic Approximation Methods Wenxun Xing Department of Mathematical Sciences Tsinghua University, Beijing Email: wxing@math.tsinghua.edu.cn OUTLINE OUTLINE
More informationContinuous Optimisation, Chpt 9: Semidefinite Optimisation
Continuous Optimisation, Chpt 9: Semidefinite Optimisation Peter J.C. Dickinson DMMP, University of Twente p.j.c.dickinson@utwente.nl http://dickinson.website/teaching/2017co.html version: 28/11/17 Monday
More informationThe Simplest Semidefinite Programs are Trivial
The Simplest Semidefinite Programs are Trivial Robert J. Vanderbei Bing Yang Program in Statistics & Operations Research Princeton University Princeton, NJ 08544 January 10, 1994 Technical Report SOR-93-12
More informationModule 04 Optimization Problems KKT Conditions & Solvers
Module 04 Optimization Problems KKT Conditions & Solvers Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html September
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 informationSolving linear equations with Gaussian Elimination (I)
Term Projects Solving linear equations with Gaussian Elimination The QR Algorithm for Symmetric Eigenvalue Problem The QR Algorithm for The SVD Quasi-Newton Methods Solving linear equations with Gaussian
More informationA Quick Tour of Linear Algebra and Optimization for Machine Learning
A Quick Tour of Linear Algebra and Optimization for Machine Learning Masoud Farivar January 8, 2015 1 / 28 Outline of Part I: Review of Basic Linear Algebra Matrices and Vectors Matrix Multiplication Operators
More informationThere are several approaches to solve UBQP, we will now briefly discuss some of them:
3 Related Work There are several approaches to solve UBQP, we will now briefly discuss some of them: Since some of them are actually algorithms for the Max Cut problem (MC), it is necessary to show their
More information