Source author record

Biel Roig-Solvas

Biel Roig-Solvas appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.OC

Catalog footprint

What is connected

3works

1topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A Globally Convergent LP and SOCP-based algorithm for Semidefinite Programming

Semidefinite programs (SDP) are one of the most versatile frameworks in numerical optimization, serving as generalizations of many conic programs and as relaxations of NP-hard combinatorial problems. Their main drawback is their computational and memory complexity, which sets a practical limit to the size of problems solvable by off-the-shelf SDP solvers. To circumvent this fact, many algorithms have been proposed to exploit the structure of particular problems and increase the scalability of SDPs for those problem instances. Progress has been less steep, however, for general-case SDPs. In this paper, motivated by earlier results by Ahmadi and Hall, we show that a general SDP can be solved to $ε$-optimality, in polynomial time, by performing a sequence of less computationally demanding Linear or Second Order Cone programs. In addition, we provide a bound on the number of iterations required to achieve $ε$-optimality. These results are illustrated using random SDPs and well-known problems from the SDPLib dataset.

preprint2022arXiv

Euclidean Distance Bounds for LMI Analytic Centers using a Novel Bound on the Lambert Function

Linear matrix inequalities (LMIs) are ubiquitous in modern control theory, as well as in a variety of other fields in science and engineering. Their analytic centers, i.e. the maximum determinant elements of the feasible set spanned by these LMIs, are the solution of many well-known problems in statistics, communications, geometry, and control and can be approximated to arbitrary precision by semidefinite programs (SDPs). The quality of these approximations is measured with respect to the difference in log-determinant of both the exact and the approximate solution to these SDPs, a quantity that follows directly from the duality theory of semidefinite programming. However, in many applications the relevant parameters are functions of the entries of the LMI argument $X$. In these cases it is of interest to develop metrics that quantify the quality of approximate solutions based on the error on these parameters, something that the log-determinant error fails to capture due to the non-linear interaction of all the matrix entries. In this work we develop upper bounds on the Frobenius norm error between suboptimal solutions $X_f$ and the exact optimizer $X_*$ of maximum determinant problems, a metric that provides a direct translation to the entry-wise error of $X$ and thus to the relevant parameters of the application. We show that these bounds can be expressed through the use of the Lambert function $W(x)$, i.e. the solution of the equation $W(x) e^{W(x)}= x$, and derive novel bounds for one of its branches to generate efficient closed-form bounds on the Euclidean distance to the LMI analytic center. Finally, we test the quality of these bounds numerically in the context of interior point methods termination criteria.

preprint2019arXiv

Chordal Decomposition in Rank Minimized Semidefinite Programs with Applications to Subspace Clustering

Semidefinite programs (SDPs) often arise in relaxations of some NP-hard problems, and if the solution of the SDP obeys certain rank constraints, the relaxation will be tight. Decomposition methods based on chordal sparsity have already been applied to speed up the solution of sparse SDPs, but methods for dealing with rank constraints are underdeveloped. This paper leverages a minimum rank completion result to decompose the rank constraint on a single large matrix into multiple rank constraints on a set of smaller matrices. The re-weighted heuristic is used as a proxy for rank, and the specific form of the heuristic preserves the sparsity pattern between iterations. Implementations of rank-minimized SDPs through interior-point and first-order algorithms are discussed. The problem of subspace clustering is used to demonstrate the computational improvement of the proposed method.