Fast verified computation for the solution of the T-congruence Sylvester equation (Q1756718)

From MaRDI portal
scientific article
Language Label Description Also known as
English
Fast verified computation for the solution of the T-congruence Sylvester equation
scientific article

    Statements

    Fast verified computation for the solution of the T-congruence Sylvester equation (English)
    0 references
    0 references
    21 December 2018
    0 references
    An algorithm is proposed to compute an interval matrix that contains the symmetric matrix $X$ that solves $AX+X^TB=C$ (all matrices are in $\mathbb{C}^{n\times n}$). It was shown in [\textit{M. Oozawa} et al., J. Comput. Appl. Math. 329, 51--56 (2018; Zbl 1373.15023)] how, under certain conditions, the problem can be transformed into a Lyapunov equation. Based on this result, conditions on the numerically computed generalized eigendecomposition of the pencil $(B,A^T)$ or $(A,B^T)$ are formulated to guarantee the existence of a unique solution. The algorithm finding the required interval matrix has complexity $O(n^3)$ assuming that the employed enclosure routines for products and inverses of matrices have complexity at most $O(n^3)$. The MATLAB code using INTLAB (interval arithmetic toolbox) is available online.
    0 references
    T-congruence Sylvester equation
    0 references
    discrete Lyapunov equation
    0 references
    verified numerical computation
    0 references
    0 references
    0 references

    Identifiers