Fast Interpolation-based Globality Certificates for Computing Kreiss Constants and the Distance to Uncontrollability

From MaRDI portal
Publication:6326423

DOI10.1137/20M1358955arXiv1910.01069MaRDI QIDQ6326423FDOQ6326423

Tim Mitchell

Publication date: 2 October 2019

Abstract: We propose a new approach to computing global minimizers of singular value functions in two real variables. Specifically, we present new algorithms to compute the Kreiss constant of a matrix and the distance to uncontrollability of a linear control system, both to arbitrary accuracy. Previous state-of-the-art methods for these two quantities rely on 2D level-set tests that are based on solving large eigenvalue problems. Consequently, these methods are costly, i.e., mathcalO(n6) work using dense eigensolvers, and often multiple tests are needed before convergence. Divide-and-conquer techniques have been proposed that reduce the work complexity to mathcalO(n4) on average and mathcalO(n5) in the worst case, but these variants are nevertheless still very expensive and can be numerically unreliable. In contrast, our new interpolation-based globality certificates perform level-set tests by building interpolant approximations to certain one-variable continuous functions that are both relatively cheap and numerically robust to evaluate. Our new approach has a mathcalO(kn3) work complexity and uses mathcalO(n2) memory, where k is the number of function evaluations necessary to build the interpolants. Not only is this interpolation process mostly "embarrassingly parallel," but also low-fidelity approximations typically suffice for all but the final interpolant, which must be built to high accuracy. Even without taking advantage of the aforementioned parallelism, k is sufficiently small that our new approach is generally orders of magnitude faster than the previous state-of-the-art.












This page was built for publication: Fast Interpolation-based Globality Certificates for Computing Kreiss Constants and the Distance to Uncontrollability

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6326423)