A tangential low-rank ADI method for solving indefinite Lyapunov equations

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Smith, Rudi, Werner, Steffen W. R.
Natura: Preprint
Pubblicazione: 2025
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866908694290104320
author Smith, Rudi
Werner, Steffen W. R.
author_facet Smith, Rudi
Werner, Steffen W. R.
contents Continuous-time algebraic Lyapunov equations have become an essential tool in various applications. In the case of large-scale sparse coefficient matrices and indefinite constant terms, indefinite low-rank factorizations have successfully been used to allow methods like the alternating direction implicit (ADI) iteration to efficiently compute accurate approximations to the solution of the Lyapunov equation. However, classical block-type approaches quickly increase in computational costs when the rank of the constant term grows. In this paper, we propose a novel tangential reformulation of the ADI iteration that allows for the efficient construction of low-rank approximations to the solution of Lyapunov equations with indefinite right-hand sides even in the case of constant terms with higher ranks. We provide adaptive methods for the selection of the corresponding ADI parameters, namely shifts and tangential directions, which allow for the automatic application of the method to any relevant problem setting. The effectiveness of the developed algorithms is illustrated by several numerical examples.
format Preprint
id arxiv_https___arxiv_org_abs_2512_04983
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle A tangential low-rank ADI method for solving indefinite Lyapunov equations
Smith, Rudi
Werner, Steffen W. R.
Numerical Analysis
Optimization and Control
15A24, 65F45, 65F55, 65H10, 93A15
Continuous-time algebraic Lyapunov equations have become an essential tool in various applications. In the case of large-scale sparse coefficient matrices and indefinite constant terms, indefinite low-rank factorizations have successfully been used to allow methods like the alternating direction implicit (ADI) iteration to efficiently compute accurate approximations to the solution of the Lyapunov equation. However, classical block-type approaches quickly increase in computational costs when the rank of the constant term grows. In this paper, we propose a novel tangential reformulation of the ADI iteration that allows for the efficient construction of low-rank approximations to the solution of Lyapunov equations with indefinite right-hand sides even in the case of constant terms with higher ranks. We provide adaptive methods for the selection of the corresponding ADI parameters, namely shifts and tangential directions, which allow for the automatic application of the method to any relevant problem setting. The effectiveness of the developed algorithms is illustrated by several numerical examples.
title A tangential low-rank ADI method for solving indefinite Lyapunov equations
topic Numerical Analysis
Optimization and Control
15A24, 65F45, 65F55, 65H10, 93A15
url https://arxiv.org/abs/2512.04983