Substitution algorithms for rational matrix equations

Fasi, Massimiliano and Iannazzo, Bruno (2019) Substitution algorithms for rational matrix equations. [MIMS Preprint]

Text faia19.pdf Download (463kB)

Abstract

We study equations of the form $r(X) = A$, where $r$ is a rational function and $A$ and $X$ are square matrices of the same size. We develop two techniques for solving these equations by inverting, through a substitution strategy, two schemes for the evaluation of rational functions of matrices. For block triangular matrices, the new methods yield the same computational cost as the evaluation schemes from which they are obtained. A general equation is reduced to block upper triangular form by exploiting the Schur decomposition of the given matrix. For real data, using the real Schur decomposition, the algorithms compute the real solutions using only real arithmetic, and numerical experiments show that our implementations are superior, in terms of both accuracy and speed, to existing alternatives. Finally, we discuss how these techniques can be used as building blocks in algorithms for computing functions of matrices defined through matrix equation of the type $f(X) = A$, where $f$ is a primary matrix function.

Item Type:	MIMS Preprint
Subjects:	MSC 2010, the AMS's Mathematics Subject Classification > 15 Linear and multilinear algebra; matrix theory MSC 2010, the AMS's Mathematics Subject Classification > 65 Numerical analysis
Depositing User:	Mr Massimiliano Fasi
Date Deposited:	12 May 2019 20:50
Last Modified:	12 May 2019 20:50
URI:	https://eprints.maths.manchester.ac.uk/id/eprint/2708

Actions (login required)

View Item