English

Computing the Lyapunov operator \varphi-functions, with an application to matrix-valued exponential integrators

Numerical Analysis 2022-04-28 v1 Numerical Analysis

Abstract

In this paper, we develop efficient and accurate evaluation for the Lyapunov operator function φl(LA)[Q],\varphi_l(\mathcal{L}_A)[Q], where φl()\varphi_l(\cdot) is the function related to the exponential, LA\mathcal{L}_A is a Lyapunov operator and QQ is a symmetric and full-rank matrix. An important application of the algorithm is to the matrix-valued exponential integrators for matrix differential equations such as differential Lyapunov equations and differential Riccati equations. The method is exploited by using the modified scaling and squaring procedure combined with the truncated Taylor series. A quasi-backward error analysis is presented to determine the value of the scaling parameter and the degree of the Taylor approximation. Numerical experiments show that the algorithm performs well in both accuracy and efficiency.

Keywords

Cite

@article{arxiv.2204.12976,
  title  = {Computing the Lyapunov operator \varphi-functions, with an application to matrix-valued exponential integrators},
  author = {Dongping Li and Yue Zhang and Xiuying Zhang},
  journal= {arXiv preprint arXiv:2204.12976},
  year   = {2022}
}

Comments

15 pages

R2 v1 2026-06-24T11:00:25.259Z