Details

Autor(en) / Beteiligte

• VANDEREYCKEN, Bart
• VANDEWALLE, Stefan

Titel

A RIEMANNIAN OPTIMIZATION APPROACH FOR COMPUTING LOW-RANK SOLUTIONS OF LYAPUNOV EQUATIONS

Ist Teil von

• SIAM journal on matrix analysis and applications, 2010, Vol.31 (5), p.2553-2579

Ort / Verlag

Philadelphia, PA: Society for Industrial and Applied Mathematics

Erscheinungsjahr

2010

Link zum Volltext

Quelle

EBSCOhost Business Source Ultimate

Beschreibungen/Notizen

• We propose a new framework based on optimization on manifolds to approximate the solution of a Lyapunov matrix equation by a low-rank matrix. The method minimizes the error on the Riemannian manifold of symmetric positive semidefinite matrices of fixed rank. We detail how objects from differential geometry, like the Riemannian gradient and Hessian, can be efficiently computed for this manifold. As a minimization algorithm we use the Riemannian trust-region method of [P.-A. Absil, C. Baker, and K. Gallivan, Found. Comput. Math., 7 (2007), pp. 303-330] based on a second-order model of the objective function on the manifold. Together with an efficient preconditioner, this method can find low-rank solutions with very little memory. We illustrate our results with numerical examples.