Abstract: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.

Computing the Lambert W function in arbitrary-precision complex interval arithmetic

New approximations to the principal real-valued branch of the Lambert W-function

Arb: Efficient Arbitrary-Precision Midpoint-Radius Interval Arithmetic

Complex branches of a generalised Lambert $W$ function arising from $p,q$--binomial coefficients

Numerical integration in arbitrary-precision ball arithmetic

Numerical Evaluation of Elliptic Functions, Elliptic Integrals and Modular Forms

Polyarc bounded complex interval arithmetic

Series representations of all branches of the Lambert W function

A Robust Approximation to a Lambert-Type Function

Winding Numbers, Unwinding Numbers, and the Lambert W Function

Rigorous Roundoff Error Analysis of Probabilistic Floating-Point Computations

Efficient implementation of elementary functions in the medium-precision range

Modified Affine Arithmetic is More Accurate Than Centered Interval Arithmetic or Affine Arithmetic

Faster arbitrary-precision dot product and matrix multiplication

A Probabilistic Approach to Floating-Point Arithmetic

Improved Laguerre Spectral Methods with Less Round-off Errors and Better Stability

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

Acknowledgment of priority: Usage of the Lambert W function in statistics

Computation of polynomial and rational approximations in complex domains by the -method

On Wirsing's problem in small exact degree

Effective Quadratic Error Bounds for Floating-Point Algorithms Computing the Hypotenuse Function