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.

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

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

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

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

A structure-preserving doubling algorithm for nonsymmetric algebraic Riccati equation

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

Low-rank alternating direction doubling algorithm for solving large-scale continuous time algebraic Riccati equations

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

The Intrinsic Toeplitz Structure and Its Applications in Algebraic Riccati Equations

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

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

ISC: an RADI-type Method for Stochastic Continuous-Time Algebraic Riccati Equations.

An RADI-type method for stochastic continuous-time algebraic Riccati equations

Decomposition Methods for Large-Scale Semidefinite Programs with Chordal Aggregate Sparsity and Partial Orthogonality

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

Low-rank plus diagonal approximations for Riccati-like matrix differential equations

Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method

Accurate Numerical Solution for Structured M-matrix Algebraic Riccati Equations

Using $LDL^{T}$ factorizations in Newton's method for solving general large-scale algebraic Riccati equations

Improved structural methods for nonlinear differential-algebraic equations via combinatorial relaxation

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