65 research outputs found
Index Reduction for Differential-Algebraic Equations with Mixed Matrices
Differential-algebraic equations (DAEs) are widely used for modeling of
dynamical systems. The difficulty in solving numerically a DAE is measured by
its differentiation index. For highly accurate simulation of dynamical systems,
it is important to convert high-index DAEs into low-index DAEs. Most of
existing simulation software packages for dynamical systems are equipped with
an index-reduction algorithm given by Mattsson and S\"{o}derlind.
Unfortunately, this algorithm fails if there are numerical cancellations.
These numerical cancellations are often caused by accurate constants in
structural equations. Distinguishing those accurate constants from generic
parameters that represent physical quantities, Murota and Iri introduced the
notion of a mixed matrix as a mathematical tool for faithful model description
in structural approach to systems analysis. For DAEs described with the use of
mixed matrices, efficient algorithms to compute the index have been developed
by exploiting matroid theory.
This paper presents an index-reduction algorithm for linear DAEs whose
coefficient matrices are mixed matrices, i.e., linear DAEs containing physical
quantities as parameters. Our algorithm detects numerical cancellations between
accurate constants, and transforms a DAE into an equivalent DAE to which
Mattsson--S\"{o}derlind's index-reduction algorithm is applicable. Our
algorithm is based on the combinatorial relaxation approach, which is a
framework to solve a linear algebraic problem by iteratively relaxing it into
an efficiently solvable combinatorial optimization problem. The algorithm does
not rely on symbolic manipulations but on fast combinatorial algorithms on
graphs and matroids. Furthermore, we provide an improved algorithm under an
assumption based on dimensional analysis of dynamical systems.Comment: A preliminary version of this paper is to appear in Proceedings of
the Eighth SIAM Workshop on Combinatorial Scientific Computing, Bergen,
Norway, June 201
A Semidefinite Hierarchy for Containment of Spectrahedra
A spectrahedron is the positivity region of a linear matrix pencil and thus
the feasible set of a semidefinite program. We propose and study a hierarchy of
sufficient semidefinite conditions to certify the containment of a
spectrahedron in another one. This approach comes from applying a moment
relaxation to a suitable polynomial optimization formulation. The hierarchical
criterion is stronger than a solitary semidefinite criterion discussed earlier
by Helton, Klep, and McCullough as well as by the authors. Moreover, several
exactness results for the solitary criterion can be brought forward to the
hierarchical approach. The hierarchy also applies to the (equivalent) question
of checking whether a map between matrix (sub-)spaces is positive. In this
context, the solitary criterion checks whether the map is completely positive,
and thus our results provide a hierarchy between positivity and complete
positivity.Comment: 24 pages, 2 figures; minor corrections; to appear in SIAM J. Opti
The Multiplicative Compound of a Matrix Pencil with Applications to Difference-Algebraic Equations
The multiplicative and additive compounds of a matrix have important
applications in geometry, linear algebra, and dynamical systems described by
difference equations and by ordinary differential equations. Here, we introduce
a generalization of the multiplicative compound to matrix pencils. We analyze
the properties of this new compound and describe several applications to the
analysis of discrete-time dynamical systems described by difference-algebraic
equations
Combinatorial Optimization
Combinatorial Optimization is a very active field that benefits from bringing together ideas from different areas, e.g., graph theory and combinatorics, matroids and submodularity, connectivity and network flows, approximation algorithms and mathematical programming, discrete and computational geometry, discrete and continuous problems, algebraic and geometric methods, and applications. We continued the long tradition of triannual Oberwolfach workshops, bringing together the best researchers from the above areas, discovering new connections, and establishing new and deepening existing international collaborations
Real Algebraic Geometry With A View Toward Systems Control and Free Positivity
New interactions between real algebraic geometry, convex optimization and free non-commutative geometry have recently emerged, and have been the subject of numerous international meetings. The aim of the workshop was to bring together experts, as well as young researchers, to investigate current key questions at the interface of these fields, and to explore emerging interdisciplinary applications
Real Algebraic Geometry With a View Toward Moment Problems and Optimization
Continuing the tradition initiated in MFO workshop held in 2014, the aim of this workshop was to foster the interaction between real algebraic geometry, operator theory, optimization, and algorithms for systems control. A particular emphasis was given to moment problems through an interesting dialogue between researchers working on these problems in finite and infinite dimensional settings, from which emerged new challenges and interdisciplinary applications
- …