Numerical solution of generalized Lyapunov equations

Two efficient methods for solving generalized Lyapunov equations and their implementations in FORTRAN 77 are presented. The first one is a generalization of the Bartels–Stewart method and the second is an extension of Hammarling's method to generalized Lyapunov equations. Our LAPACK based subroutines are implemented in a quite flexible way. They can handle the transposed equations and provide scaling to avoid overflow in the solution. Moreover, the Bartels–Stewart subroutine offers the optional estimation of the separation and the reciprocal condition number. A brief description of both algorithms is given. The performance of the software is demonstrated by numerical experiments.

[1]  Sven Hammarling.  Numerical solution of the discrete-time, convergent, non-negative definite Lyapunov equation , 1991 .

[2]  E. M. Wright.  Solution of the equation $ze^z = a$ , 1959 .

[3]  Alan J. Laub,et al.  Solution of the Sylvester matrix equation AXBT + CXDT = E , 1992, TOMS.

[4]  Jack Dongarra,et al.  LINPACK Users' Guide , 1987 .

[5]  Peter Lancaster,et al.  The theory of matrices , 1969 .

[6]  K. Chu.  The solution of the matrix equations AXB−CXD=E AND (YA−DZ,YC−BZ)=(E,F) , 1987 .

[7]  B. A. J. Cockle,et al.  Solution of an equation , 1843 .

[8]  S. Hammarling.  Numerical Solution of the Stable, Non-negative Definite Lyapunov Equation , 1982 .

[9]  J. Snyders,et al.  On Nonnegative Solutions of the Equation $AD + DA' = - C$ , 1970 .

[10]  Antony Jameson,et al.  Solution of the Equation $AX + XB = C$ by Inversion of an $M \times M$ or $N \times N$ Matrix , 1968 .

[11]  Ed Anderson,et al.  LAPACK Users' Guide , 1995 .

[12]  Richard H. Bartels,et al.  Algorithm 432 [C2]: Solution of the matrix equation AX + XB = C [F4] , 1972, Commun. ACM.

[13]  Alan J. Laub,et al.  Algorithm 705; a FORTRAN-77 software package for solving the Sylvester matrix equation AXBT + CXDT = E , 1992, TOMS.

[14]  B. Kågström,et al.  Generalized Schur methods with condition estimators for solving the generalized Sylvester equation , 1989 .

[15]  Nicholas J. Higham,et al.  Perturbation theory and backward error forAX−XB=C , 1993 .

[16]  W. Bowen,et al.  Philadelphia , 1892 .

[17]  F. R. Gantmakher.  The Theory of Matrices , 1984 .

[18]  Jack J. Dongarra,et al.  Matrix Eigensystem Routines — EISPACK Guide Extension , 1977, Lecture Notes in Computer Science.

[19]  Gene H. Golub,et al.  Matrix computations , 1983 .

[20]  R. Byers.  A LINPACK-style condition estimator for the equation AX-XB^{T} = C , 1984 .

[21]  Nicholas J. Higham,et al.  FORTRAN codes for estimating the one-norm of a real or complex matrix, with applications to condition estimation , 1988, TOMS.

[22]  Bo Kågström,et al.  LAPACK-style algorithms and software for solving the generalized Sylvester equation and estimating the separation between regular matrix pairs , 1994, TOMS.