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.

Iterative and doubling algorithms for Riccati-type matrix equations: a comparative introduction

Two Iterative Algorithms for Stochastic Algebraic Riccati Matrix Equations.

Performance Enhancement of Doubling Algorithms for a Class of Complex Nonsymmetric Algebraic Riccati Equations

Solving Large-Scale Continuous-Time Algebraic Riccati Equations by Doubling

A Novel Iterative Algorithm for Solving Coupled Riccati Equations

Solving Large-scale Discrete-time Algebraic Riccati Equations by Doubling

A New Relaxative Algorithm for Coupled Riccati Matrix Equations

Highly accurate decoupled doubling algorithm for large-scale M-matrix algebraic Riccati equations

Convergence Analysis of Structure-Preserving Doubling Algorithms for Riccati-Type Matrix Equations.

An Inversion-Free Iterative Algorithm for Riccati Matrix Equations in Discrete-Time Markov Jump Systems

An implicit iterative algorithm for generalized algebraic Riccati matrix equations

Convergence Analysis of the Doubling Algorithm for Several Nonlinear Matrix Equations in the Critical Case.

Solving Large-Scale Nonsymmetric Algebraic Riccati Equations By Doubling

On the Iterative Methods for the Solution of Three Types of Nonlinear Matrix Equations

An Implicit Iterative Algorithm for Algebraic Riccati Matrix Equations

An Inversion-Free Iterative Algorithm with a Scalar Tuning Parameter for Coupled Riccati Matrix Equations Arising in LQ Optimal Control of Markov Jump Systems

Highly Accurate Doubling Algorithms for M -Matrix Algebraic Riccati Equations

On a family of low-rank algorithms for large-scale algebraic Riccati equations

Fixed-Point Iteration Schemes to Solve Symmetric Algebraic Riccati Equation $$XBX-XA-A^{T}X-C=0$$

Gradient Based Iterative Algorithm for Solving Coupled Matrix Equations

An Iterative Algorithm for Coupled Riccati Equations in Continuous-Time Markovian Jump Linear Systems