An increasing rank Riemannian method for generalized Lyapunov equations
Fuente:
arXiv
Saved in:
| Main Authors: | , |
|---|---|
| Format: | Preprint |
| Published: |
2023
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866916113739153408 |
|---|---|
| author | Huang, Zhenwei Huang, Wen |
| author_facet | Huang, Zhenwei Huang, Wen |
| contents | In this paper, we consider finding a low-rank approximation to the solution of a large-scale generalized Lyapunov matrix equation in the form of $A X M + M X A = C$, where $A$ and $M$ are symmetric positive definite matrices. An algorithm called an Increasing Rank Riemannian Method for Generalized Lyapunov Equation (IRRLyap) is proposed by merging the increasing rank technique and Riemannian optimization techniques on the quotient manifold $\mathbb{R}_*^{n \times p} / \mathcal{O}_p$. To efficiently solve the optimization problem on $\mathbb{R}_*^{n \times p} / \mathcal{O}_p$, a line-search-based Riemannian inexact Newton's method is developed with its global convergence and local superlinear convergence rate guaranteed. Moreover, we derive a preconditioner which takes $M \neq I$ into consideration. Numerical experiments show that the proposed Riemannian inexact Newton's method and preconditioner has superior performance and IRRLyap is preferable compared to the tested state-of-the-art methods when the lowest rank solution is desired. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2308_00213 |
| institution | arXiv |
| publishDate | 2023 |
| record_format | arxiv |
| spellingShingle | An increasing rank Riemannian method for generalized Lyapunov equations Huang, Zhenwei Huang, Wen Optimization and Control In this paper, we consider finding a low-rank approximation to the solution of a large-scale generalized Lyapunov matrix equation in the form of $A X M + M X A = C$, where $A$ and $M$ are symmetric positive definite matrices. An algorithm called an Increasing Rank Riemannian Method for Generalized Lyapunov Equation (IRRLyap) is proposed by merging the increasing rank technique and Riemannian optimization techniques on the quotient manifold $\mathbb{R}_*^{n \times p} / \mathcal{O}_p$. To efficiently solve the optimization problem on $\mathbb{R}_*^{n \times p} / \mathcal{O}_p$, a line-search-based Riemannian inexact Newton's method is developed with its global convergence and local superlinear convergence rate guaranteed. Moreover, we derive a preconditioner which takes $M \neq I$ into consideration. Numerical experiments show that the proposed Riemannian inexact Newton's method and preconditioner has superior performance and IRRLyap is preferable compared to the tested state-of-the-art methods when the lowest rank solution is desired. |
| title | An increasing rank Riemannian method for generalized Lyapunov equations |
| topic | Optimization and Control |
| url | https://arxiv.org/abs/2308.00213 |