Parallel Solution of Large Lyapunov Equations*

In this paper two algorithms for the solution of large-order $(100 \leq n \leq 1000)$ Lyapunov equations $AX + XA^T + Q = 0$ are presented. First, a parallel version of the Hammarling algorithm for the solution of Lyapunov equations where the coefficient matrix A is large and dense is presented. Then a novel iterative parallel algorithm, called full-rank perturbed iteration (FRPI), is presented for the solution of Lyapunov equations where the matrix A is large and banded.

[1]  Dianne P. O'Leary,et al.  Data-flow algorithms for parallel matrix computation , 1985, CACM.

[2]  Christian H. Bischof,et al.  The WY representation for products of householder matrices , 1985, PPSC.

[3]  E. Wachspress.  Iterative solution of the Lyapunov matrix equation , 1988 .

[4]  A. S. Hodel,et al.  Homotopy methods for the solution of general modified algebraic Riccati equations , 1993 .

[5]  A. Hodel,et al.  Numerical Methods for the Solution of Large and Very Large, Sparse Lyapunov Equations , 1989 .

[6]  C. Loan,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix , 1978 .

[7]  A. Laub,et al.  Computation of system balancing transformations and other applications of simultaneous diagonalization algorithms , 1987 .

[8]  D. Bernstein,et al.  The optimal projection equations for fixed-order dynamic compensation , 1984 .

[9]  Michael T. Heath,et al.  Sparse Cholesky factorization on a local-memory multiprocessor , 1988 .

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

[11]  G. A. Geist,et al.  Stabilized Gaussian reduction of an arbitrary matrix to tridiagonal form , 1989 .

[12]  G. Golub,et al.  A Hessenberg-Schur method for the problem AX + XB= C , 1979 .

[13]  K. Glover.  All optimal Hankel-norm approximations of linear multivariable systems and their L, ∞ -error bounds† , 1984 .

[14]  Michael T. Heath,et al.  Parallel solution of triangular systems on distributed-memory multiprocessors , 1988 .

[15]  Emmanuel G. Collins,et al.  A homotopy algorithm for reduced order compensator design using the optimal projection equations , 1989, Proceedings of the 28th IEEE Conference on Decision and Control,.

[16]  D. Heller.  A Survey of Parallel Algorithms in Numerical Linear Algebra. , 1978 .

[17]  G. Hewer,et al.  The sensitivity of the stable Lyapunov equation , 1987, 26th IEEE Conference on Decision and Control.

[18]  Y. Saad,et al.  A hybrid Chebyshev Krylov subspace algorithm for solving nonsymmetric systems of linear equations , 1986 .

[19]  Lawrence S. Kroll.  Mathematica--A System for Doing Mathematics by Computer. , 1989 .

[20]  Mark J. Balas,et al.  Trends in large space structure control theory: Fondest hopes, wildest dreams , 1982 .

[21]  Edmond A. Jonckheere,et al.  A new set of invariants for linear systems--Application to reduced order compensator design , 1983 .