Low rank solution of Lyapunov equations

Jing-Rebecca Li 1 Jacob White
1 POEMS - Propagation des Ondes : Étude Mathématique et Simulation
Inria Saclay - Ile de France, ENSTA ParisTech UMA - Unité de Mathématiques Appliquées, CNRS - Centre National de la Recherche Scientifique : UMR7231
Abstract : This paper presents the Cholesky factor--alternating direction implicit (CF--ADI) algorithm, which generates a low-rank approximation to the solution X of the Lyapunov equation AX+XAT = -BBT. The coefficient matrix A is assumed to be large, and the rank of the right-hand side -BBT is assumed to be much smaller than the size of A. The CF--ADI algorithm requires only matrix-vector products and matrix-vector solves by shifts of A. Hence, it enables one to take advantage of any sparsity or structure in A. This paper also discusses the approximation of the dominant invariant subspace of the solution X. We characterize a group of spanning sets for the range of X. A connection is made between the approximation of the dominant invariant subspace of X and the generation of various low-order Krylov and rational Krylov subspaces. It is shown by numerical examples that the rational Krylov subspace generated by the CF--ADI algorithm, where the shifts are obtained as the solution of a rational minimax problem, often gives the best approximation to the dominant invariant subspace of X. Copyright © 2004 Society for Industrial and Applied Mathematics
Type de document :
Article dans une revue
SIAM Review, Society for Industrial and Applied Mathematics, 2004, 46 (4), pp.693-713. <10.1137/S0036144504443389>
Liste complète des métadonnées

https://hal-ensta.archives-ouvertes.fr/hal-00781126
Contributeur : Aurélien Arnoux <>
Soumis le : vendredi 25 avril 2014 - 14:30:12
Dernière modification le : jeudi 23 mars 2017 - 01:02:27

Identifiants

Collections

Citation

Jing-Rebecca Li, Jacob White. Low rank solution of Lyapunov equations. SIAM Review, Society for Industrial and Applied Mathematics, 2004, 46 (4), pp.693-713. <10.1137/S0036144504443389>. <hal-00781126v2>

Partager

Métriques

Consultations de la notice

200