Abstract:The doubling algorithms are very efficient iterative methods for computing the unique minimal nonnegative solution to an M-matrix algebraic Riccati equation (MARE). They are globally and quadratically convergent, except for MARE in the critical case where convergence is linear with the linear rate 1 / 2. However, the initialization phase and the doubling iteration kernel of any doubling algorithm involve inverting nonsingular M-matrices. In particular, for MARE in the critical case, the M-matrices in the doubling iteration kernel, although nonsingular, move towards singular M-matrices at convergence. These inversions are causes of concerns on entrywise relative accuracy of the eventually computed minimal nonnegative solution. Fortunately, a nonsingular M-matrix can be inverted by the GTH-like algorithm of Alfa et al. (Math Comp 71:217-236, 2002) to almost full entrywise relative accuracy, provided a triplet representation of the matrix is known. Recently, Nguyen and Poloni (Numer Math 130(4):763-792, 2015) discovered a way to construct triplet representations in a cancellation-free manner for all involved M-matrices in the doubling iteration kernel, for a special class of MAREs arising from Markov-modulated fluid queues. In this paper, we extend Nguyen's and Poloni's work to all MAREs by also devising a way to construct the triplet representations cancellation-free. Our construction, however, is not a straightforward extension of theirs. It is made possible by an introduction of novel recursively computable auxiliary nonnegative vectors. As the second contribution, we propose an entrywise relative residual for an approximate solution. The residual has an appealing feature of being able to reveal the entrywise relative accuracies of all entries, large and small, of the approximation. This is in marked contrast to the usual legacy normalized residual which reflects relative accuracies of large entries well but not so much those of very tiny entries. Numerical examples are presented to demonstrate and confirm our claims.

Entrywise relative perturbation theory for nonsingular M-matrices and applications

Deviation Inequalities on the Spectral Norm of Random Matrices with Dependent Entries

Entrywise perturbation theory for diagonally dominant M-matrices with applications

Some Sufficient Conditions for Nonsingular H-Matrices

Entrywise relative perturbation bounds for exponentials of essentially non-negative matrices

Accurate solutions of M -matrix Sylvester equations

Accurate solutions of M -matrix algebraic Riccati equations

Computing the Smallest Eigenvalue of an M-Matrix.

A positivity preserving inexact Noda iteration for computing the smallest eigenpair of a large irreducible M -matrix

Relative perturbation bounds for the eigenvalues of diagonalizable and singular matrices – Application of perturbation theory for simple invariant subspaces

On a perturbation analysis of Higham squared maximum Gaussian elimination growth matrices

Highly Accurate Doubling Algorithms for M -Matrix Algebraic Riccati Equations

Perturbation of the Drazin Inverse for Matrices with Equal Eigenprojections at Zero

Accurate computation of the smallest eigenvalue of a diagonally dominant M-matrix.

On the perturbation and subproper splittings for the generalized inverse AT,S(2) of rectangular matrix A

Solutions and perturbation analysis of the matrix equation X - \sum_{i=1}^m A_i^* X^{-1} A_i = Q

Perturbation analysis of the matrix equation X - \sum_{i=1}^m A_i^* X^{p_i} A_i = Q

Entrywise perturbation theory for rate matrices of GI/M/1 type Markov chains

Preconditioning Technique for Symmetric M-matrices

Perturbation analysis of the Hermitian positive definite solution of the matrix equation ?*=

Solutions and improved perturbation analysis for the matrix equation X-A^{*}X^{-p}A=Q (p>0)