Consistency and efficient solution of the Sylvester equation for {star operator}-congruence

Abstract

In this paper, the matrix equation AX + X{star operator}B = C is considered, where the matrices A and B have sizes m × n and n × m, respectively, the size of the unknown X is n × m, and the operator (·){star operator} denotes either the transpose or the conjugate transpose of a matrix. In the first part of the paper, necessary and sufficient conditions for the existence and uniqueness of solutions are reviewed. These conditions were obtained previously by Wimmer [H.K. Wimmer. Roth's theorems for matrix equations with symmetry constraints. Linear Algebra Appl., 199:357- 362, 1994.], by Byers and Kressner [R. Byers and D. Kressner. Structured condition numbers for invariant subspaces. SIAM J. Matrix Anal. Appl., 28:326-347, 2006.], and by Kressner, Schr̈oder and Watkins [D. Kressner, C. Schröder, and D.S. Watkins. Implicit QR algorithms for palindromic and even eigenvalue problems. Numer. Algorithms, 51:209-238, 2009.]. This review generalizes to fields of characteristic different from two the existence condition that Wimmer originally proved for the complex field. In the second part, an algorithm is developed, in the real or complex square case m = n, to solve the equation in O(n3) flops when the solution is unique. This algorithm is based on the generalized Schur decomposition of the matrix pencil A - λB{star operator}. The equation AX + X{star operator}B = C is connected with palindromic eigenvalue problems and, as a consequence, the square complex case has attracted recently the attention of several authors.