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.

Structure-Preserving Doubling Algorithms That Avoid Breakdowns for Algebraic Riccati-Type Matrix Equations

The Structure-preserving Doubling Algorithm for the Discrete Coupled Riccati Matrix Equations

A structure-preserving doubling algorithm for nonsymmetric algebraic Riccati equation

Solving Large-Scale Nonsymmetric Algebraic Riccati Equations By Doubling

Decoupled Structure-Preserving Doubling Algorithm with Truncation for Large-Scale Algebraic Riccati Equations

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

A decoupled form of the structure-preserving doubling algorithm with low-rank structures

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

Large-Scale Discrete-Time Algebraic Riccati Equations — Numerically Low-Ranked Solution , SDA and Error Analysis ∗

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

A STRUCTURE-PRESERVING ALGORITHM FOR SEMI-STABILIZING SOLUTIONS OF GENERALIZED ALGEBRAIC RICCATI EQUATIONS

On the generalized structure-preserving doubling algorithm for Generalized Discrete-time Algebraic Riccati equations

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

Alternating-Directional Doubling Algorithm for M-Matrix Algebraic Riccati Equations

A structure-preserving doubling algorithm for solving a class of quadratic matrix equation with $M$-matrix

On Novel Fixed-Point-Type Iterations with Structure-Preserving Doubling Algorithms for Stochastic Continuous-time Algebraic Riccati equations

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

Conversion Methods, Block Triangularization, and Structural Analysis of Differential-Algebraic Equation Systems

Two Efficient Lopsided Double-Step Methods for Solving Complex Symmetric Linear Systems

Accurate Numerical Solution for Structured M-matrix Algebraic Riccati Equations

Symbolic-numeric methods for improving structural analysis of differential-algebraic equation systems