The NEOS Server for Optimization: Version 4 and Beyond

We describe developments associated with Version 4 of the NEOS Server and note that these developments have led to an exponential growth in the number of job submissions. We also provide an overview of some of the research and educational uses for the NEOS Server and discuss future research challenges. Citation Preprint ANL/MCS-P947-0202, Argonne National … 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

Semidefinite Programming in the Space of Partial Positive Semidefinite Matrices

We build upon the work of Fukuda et al.\ \cite{FuKoMuNa01} and Nakata et al.\ \cite{NaFuFuKoMu01}, in which the theory of partial positive semidefinite matrices has been applied to the semidefinite programming (SDP) problem as a technique for exploiting sparsity in the data. In contrast to their work, which improves an existing algorithm that is based … 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

Mesh Topology Design in Overlay Virtual Private Networks

We study the mesh topology design problem in overlay virtual private networks. Given a set of customer nodes and associated traffic matrix, tunnels that are connecting node pairs through a public service provider network subject to degree constraints are determined so as to minimize total multihopped traffic. Valid inequalities strengthening the LP relaxation and a … Read more

Exact Solutions of Some Nonconvex Quadratic Optimization Problems via SDP and SOCP Relaxations

We show that SDP (semidefinite programming) and SOCP (second order cone programming) relaxations provide exact optimal solutions for a class of nonconvex quadratic optimization problems. It is a generalization of the results by S.~Zhang for a subclass of quadratic maximization problems that have nonnegative off-diagonal coefficient matrices of objective quadratic functions and diagonal coefficient matrices … Read more

Computing All Nonsingular Solutions of Cyclic-n Polynomial Using Polyhedral Homotopy Continuation Methods

All isolated solutions of the cyclic-n polynomial equations are not known for larger dimensions than 11. We exploit two types of symmetric structures in the cyclic-n polynomial to compute all isolated nonsingular solutions of the equations efficiently by the polyhedral homotopy continuation method and to verify the correctness of the generated approximate solutions. Numerical results … Read more

Solving Stability Problems on a Superclass of Interval Graphs

We introduce a graph invariant, called thinness, and show that a maximum weighted stable set on a graph $G(V, E)$ with thinness $k$ may be found in $O(\frac{|V|}{k})^k$-time, if a certain representation is given. We show that a graph has thinness 1 if and only if it is an interval graph, while a graph with … Read more