A primal affine-scaling algorithm for constrained convex programs

The affine-scaling algorithm was initially developed for linear programming problems. Its extension to problems with a nonlinear objective performs at each iteration a scaling followed by a line search along the steepest descent direction. In this paper we prove that any accumulation point generated by this algorithm when applied to a convex function is an … Read more

The Trust Region Subproblem and Semidefinite Programming

The trust region subproblem (the minimization of a quadratic objective subject to one quadratic constraint and denoted TRS) has many applications in diverse areas, e.g. function minimization, sequential quadratic programming, regularization, ridge regression, and discrete optimization. In particular, it determines the step in trust region algorithms for function minimization. Trust region algorithms are popular for … Read more

Computation of Minimum Volume Covering Ellipsoids

We present a practical algorithm for computing the minimum volume n-dimensional ellipsoid that must contain m given points a_1, …, a_m \in R^n. This convex constrained problem arises in a variety of applied computational settings, particularly in data mining and robust statistics. Its structure makes it particularly amenable to solution by interior-point methods, and it … Read more

Limiting behavior of the central path in semidefinite optimization

It was recently shown that, unlike in linear optimization, the central path in semidefinite optimization (SDO) does not converge to the analytic center of the optimal set in general. In this paper we analyze the limiting behavior of the central path to explain this unexpected phenomenon. This is done by deriving a new necessary and … Read more

Characterization of the limit point of the central path in semidefinite programming

In linear programming, the central path is known to converge to the analytic center of the set of optimal solutions. Recently, it has been shown that this is not necessarily true for linear semidefinite programming in the absence of strict complementarity. The present paper deals with the formulation of a convex problem whose solution defines … Read more

SOSTOOLS: Sum of Squares Optimization Toolbox for MATLAB — User’s Guide

SOSTOOLS is a free MATLAB toolbox for formulating and solving sum of squares (SOS) optimization programs. It uses a simple notation and a flexible and intuitive high-level user interface to specify the SOS programs. Currently these are solved using SeDuMi, a well-known semidefinite programming solver, while SOSTOOLS handles internally all the necessary reformulations and data … Read more

Relating Homogeneous Cones and Positive Definite Cones via hBcalgebras

$T$-algebras are non-associative algebras defined by Vinberg in the early 1960’s for the purpose of studying homogeneous cones. Vinberg defined a cone $K(\mathcal A)$ for each $T$-algebra $\mathcal A$ and proved that every homogeneous cone is isomorphic to one such $K(\mathcal A)$. We relate each $T$-algebra $\mathcal A$ with a space of linear operators in … Read more

PENNON – A Code for Convex Nonlinear and Semidefinite Programming

We introduce a computer program PENNON for the solution of problems of convex Nonlinear and Semidefinite Programming (NLP-SDP). The algorithm used in PENNON is a generalized version of the Augmented Lagrangian method, originally introduced by Ben-Tal and Zibulevsky for convex NLP problems. We present generalization of this algorithm to convex NLP-SDP problems, as implemented in … Read more

A primal-dual symmetric relaxation for homogeneous conic systems

We address the feasibility of the pair of alternative conic systems of constraints Ax = 0, x \in C, and -A^T y \in C^*, defined by an m by n matrix A and a cone C in the n-dimensional Euclidean space. We reformulate this pair of conic systems as a primal-dual pair of conic programs. … Read more

A new class of potential affine algorithms for linear convex programming

We obtain a new class of primal affine algorithms for the linearly constrained convex programming. It is constructed from a family of metrics generated the r power, r>=1, of the diagonal iterate vector matrix. We obtain the so called weak convergence. That class contains, as particular cases, the multiplicative Eggermont algorithm for the minimization of … Read more