Abstract:This paper is concerned with a M-matrix algebraic Riccati equation (MARE) XDX - AX - XB + C = O for which A is block-diagonal and its defining matrix W = [GRAPHICS] is a nonsingular or irreducible singular M-matrix. Such an MARE can be decomposec into many coupled algebraic Riccati equations (AREs) that can be solved by the Jacobior Gauss-Seidel-like iteration updating scheme at the outer-loop while by a doubling algorithm in the inner loop for each coupled ARE, as first proposed by Meini (2013). The goals of this paper are two-fold. One is to resolve a critical technical detail in Meini's algorithm that was not addressed. It is about whether each ARE in the inner loop has a minimal nonnegative solution. It is proved that the defining matrix of each coupled ARE during a doubling iteration is indeed a nonsingular or irreducible singular M-matrix and, as a result, they do have minimal nonnegative solutions and a doubling algorithm is an efficient way to compute them. The other goal is to design a highly accurate implementation of the doubling algorithm for the inner loop so that all entries of the minimal nonnegative solution to the original MARE are calculated with high entrywise relative accuracies, regardless of their magnitudes. This is made possible by a novel way of constructing triplet representations for the coupled ARE during doubling iterations. Numerical examples are presented to demonstrate that the resulting algorithm can indeed deliver an entrywise relatively accurate solution. (C) 2021 Elsevier B.V. All rights reserved.

Structured Condition Number And Backward Error Of Nonsymmetric Algebraic Riccati Equation

Condition Number and Backward Errors of Nonsymmetric Algebraic Riccati Equation

A Note on the Necessary Conditions for the Algebraic Riccati Equation.

A New Class of Complex Nonsymmetric Algebraic Riccati Equations

A Parameterized Class of Complex Nonsymmetric Algebraic Riccati Equations

A New Class of Complex Nonsymmetric Algebraic Riccati Equations with Its Ω-Comparison Matrix Being an Irreducible Singular M-matrix

Solving the Nonnegative Solution for a (shifted) Nonsymmetric Algebraic Riccati Equation in the Critical Case

The King-Werner Method for Solving Nonsymmetric Algebraic Riccati Equation

Normwise, Mixed and Componentwise Condition Numbers of Nonsymmetric Algebraic Riccati Equations

The relaxed Newton-like method for a nonsymmetric algebraic Riccati equation

A Class of Complex Nonsymmetric Algebraic Riccati Equations Associated with H-matrix

A structure-preserving doubling algorithm for nonsymmetric algebraic Riccati equation

Accurate Numerical Solution for Structured M-matrix Algebraic Riccati Equations

A Cone-preserving Solution to a Nonsymmetric Riccati Equation

Perturbation Analysis and Condition Numbers of Symmetric Algebraic Riccati Equations

Nonsymmetric Algebraic Riccati Equations under the Tensor Product

Alternately linearized implicit iteration methods for the minimal nonnegative solutions of the nonsymmetric algebraic Riccati equations

A modified Newton method for solving non-symmetric algebraic Riccati equations arising in transport theory

Finding an NARE whose minimal nonnegative solution represents first passage quantities in the two-dimensional Brownian motion

Comments on "upper Matrix Bounds for the Discrete Algebraic Riccati Matrix Equation" [with Author's Reply]

Accurate solutions of M -matrix algebraic Riccati equations