Direct methods for matrix Sylvester and Lyapunov equations

Abstract

We revisit the two standard dense methods for matrix Sylvester and Lyapunov equations: the Bartels-Stewart method forA1X+XA2+D=0and Hammarling's method forAX+XAT+BBT=0withAstable. We construct three schemes for solving the unitarily reduced quasitriangular systems. We also construct a new rank-1 updating scheme in Hammarling's method. This new scheme is able to accommodate aBwith more columns than rows as well as the usual case of aBwith more rows than columns, while Hammarling's original scheme needs to separate these two cases. We compared all of our schemes with the Matlab Sylvester and Lyapunov solverlyap.m; the results show that our schemes are much more efficient. We also compare our schemes with the Lyapunov solversllyapin the currently possibly the most efficient control library package SLICOT; numerical results show our scheme to be competitive.