Characterizing the universal rigidity of generic frameworks

A framework is a graph and a map from its vertices to E^d (for some d). A framework is universally rigid if any framework in any dimension with the same graph and edge lengths is a Euclidean image of it. We show that a generic universally rigid framework has a positive semi-definite stress matrix of maximal rank. Connelly showed that the existence of such a positive semi-definite stress matrix is sufficient for universal rigidity, so this provides a characterization of universal rigidity for generic frameworks. We also extend our argument to give a new result on the genericity of strict complementarity in semidefinite programming.Comment: 18 pages, v2: updates throughout; v3: published versio

Conic convex programming and self-dual embedding

How to initialize an algorithm to solve an optimization problem is of great theoretical and practical importance. In the simplex method for linear programming this issue is resolved by either the two-phase approach or using the so-called big M technique. In the interior point method, there is a more elegant way to deal with the initialization problem, viz. the self-dual embedding technique proposed by Ye, Todd and Mizuno (30). For linear programming this technique makes it possible to identify an optimal solution or conclude the problem to be infeasible/unbounded by solving its embedded self-dual problem. The embedded self-dual problem has a trivial initial solution and has the same structure as the original problem. Hence, it eliminates the need to consider the initialization problem at all. In this paper, we extend this approach to solve general conic convex programming, including semidefinite programming. Since a nonlinear conic convex programming problem may lack the so-called strict complementarity property, it causes difficulties in identifying solutions for the original problem, based on solutions for the embedded self-dual system. We provide numerous examples from semidefinite programming to illustrate various possibilities which have no analogue in the linear programming case

Computational analysis of real-time convex optimization for control systems

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Aeronautics and Astronautics, 2000.Includes bibliographical references (p. 177-189).Computational analysis is fundamental for certification of all real-time control software. Nevertheless, analysis of on-line optimization for control has received little attention to date. On-line software must pass rigorous standards in reliability, requiring that any embedded optimization algorithm possess predictable behavior and bounded run-time guarantees. This thesis examines the problem of certifying control systems which utilize real-time optimization. A general convex programming framework is used, to which primal-dual path-following algorithms are applied. The set of all optimization problem instances which may arise in an on-line procedure is characterized as a compact parametric set of convex programming problems. A method is given for checking the feasibility and well-posedness of this compact set of problems, providing certification that every problem instance has a solution and can be solved in finite time. The thesis then proposes several algorithm initialization methods, considering the fixed and time-varying constraint cases separately. Computational bounds are provided for both cases. In the event that the computational requirements cannot be met, several alternatives to on-line optimization are suggested. Of course, these alternatives must provide feasible solutions with minimal real-time computational overhead. Beyond this requirement, these methods approximate the optimal solution as well as possible. The methods explored include robust table look-up, functional approximation of the solution set, and ellipsoidal approximation of the constraint set. The final part of this thesis examines the coupled behavior of a receding horizon control scheme for constrained linear systems and real-time optimization. The driving requirement is to maintain closed-loop stability, feasibility and well-posedness of the optimal control problem, and bounded iterations for the optimization algorithm. A detailed analysis provides sufficient conditions for meeting these requirements. A realistic example of a small autonomous air vehicle is furnished, showing how a receding horizon control law using real-time optimization can be certified.by Lawrence Kent McGovern.Ph.D

A Riemannian low-rank method for optimization over semidefinite matrices with block-diagonal constraints

We propose a new algorithm to solve optimization problems of the form $\min f(X)$ for a smooth function $f$ under the constraints that $X$ is positive semidefinite and the diagonal blocks of $X$ are small identity matrices. Such problems often arise as the result of relaxing a rank constraint (lifting). In particular, many estimation tasks involving phases, rotations, orthonormal bases or permutations fit in this framework, and so do certain relaxations of combinatorial problems such as Max-Cut. The proposed algorithm exploits the facts that (1) such formulations admit low-rank solutions, and (2) their rank-restricted versions are smooth optimization problems on a Riemannian manifold. Combining insights from both the Riemannian and the convex geometries of the problem, we characterize when second-order critical points of the smooth problem reveal KKT points of the semidefinite problem. We compare against state of the art, mature software and find that, on certain interesting problem instances, what we call the staircase method is orders of magnitude faster, is more accurate and scales better. Code is available.Comment: 37 pages, 3 figure