Computational Krylov-based methods for large-scale differential Sylvester matrix problems
From MaRDI portal
Publication:6288738
DOI10.1002/NLA.2187arXiv1707.02078WikidataQ129940420 ScholiaQ129940420MaRDI QIDQ6288738FDOQ6288738
Authors: M. Hached, Khalide Jbilou
Publication date: 7 July 2017
Abstract: In the present paper, we propose Krylov-based methods for solving large-scale differential Sylvester matrix equations having a low rank constant term. We present two new approaches for solving such differential matrix equations. The first approach is based on the integral expression of the exact solution and a Krylov method for the computation of the exponential of a matrix times a block of vectors. In the second approach, we first project the initial problem onto a block (or extended block) Krylov subspace and get a low-dimensional differential Sylvester matrix equation. The latter problem is then solved by some integration numerical methods such as BDF or Rosenbrock method and the obtained solution is used to build the low rank approximate solution of the original problem. We give some new theoretical results such as a simple expression of the residual norm and upper bounds for the norm of the error. Some numerical experiments are given in order to compare the two approaches.
Matrix equations and identities (15A24) Numerical methods for initial value problems involving ordinary differential equations (65L05) Numerical methods for matrix equations (65F45)
This page was built for publication: Computational Krylov-based methods for large-scale differential Sylvester matrix problems
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6288738)