Calculating a function of a matrix with a real spectrum

Let T be a square matrix with a real spectrum, and let f be an analytic function. The problem of the approximate calculation of f(T) is discussed. Applying the Schur triangular decomposition and the reordering, one can assume that T is triangular and its diagonal entries t_ii are arranged in increasing order. To avoid calculations using the differences t_ii-t_jj with close (including equal) t_ii and t_jj, it is proposed to represent T in a block form and calculate the two main block diagonals using interpolating polynomials. The rest of the f(T) entries can be calculated using the Parlett recurrence algorithm. It is also proposed to perform scalar operations (such as the building of interpolating polynomials) with an enlarged number of decimal digits.

• 1 publication
• 2 publications
• 2 publications
research
11/28/2022

On the Rellich eigendecomposition of para-Hermitian matrices and the sign characteristics of *-palindromic matrix polynomials

We study the eigendecompositions of para-Hermitian matrices H(z), that i...
research
09/25/2019

A variant of Schur's product theorem and its applications

We show the following version of the Schur's product theorem. If M=(M_j,...
research
09/02/2021

Hankel Determinants of Certain Sequences Of Bernoulli Polynomials: A Direct Proof of an Inverse Matrix Entry from Statistics

We calculate the Hankel determinants of sequences of Bernoulli polynomia...
research
12/27/2022

Efficient method for calculating the eigenvalue of the Zakharov-Shabat system

In this paper, a direct method is proposed to calculate the eigenvalue o...
research
01/06/2022

Simple algorithm for GCD of polynomials

Based on the Bezout approach we propose a simple algorithm to determine ...
research
08/04/2021

An estimate of approximation of an analytic function of a matrix by a rational function

Let A be a square complex matrix; z_1, ..., z_N∈ℂ be arbitrary (possibly...
research
04/23/2022

New features in MathPartner 2021

We introduce new features in the MathPartner service that have recently ...