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.

A tensor Alternating Anderson-Richardson method for solving multilinear systems with M-tensors

MODIFIED SPLITTING METHODS FOR SOLVING NON-HOMOGENOUS MULTI-LINEAR EQUATIONS WITH M-TENSORS

Solving Multi-linear Systems with M-Tensors.

Restarted Nonnegativity Preserving Tensor Splitting Methods via Relaxed Anderson Acceleration for Solving Multi-linear Systems

Solving Multi-linear Systems with $$\mathcal {M}$$m-Tensors

Solving Multi-linear Systems with ℳ -Tensors

Preconditioned Tensor Splitting AOR Iterative Methods for ℋ‐tensor Equations

Tensor Methods for Solving Symmetric ℳ -tensor Systems

Tensor Methods for Solving Symmetric \({\mathcal {M}}\)-tensor Systems

Efficient Alternating Least Squares Algorithms for Low Multilinear Rank Approximation of Tensors

Existence and uniqueness of positive solution for multilinear systems with generalized strong M-tensor

Randomized Kaczmarz for tensor linear systems

Solving the system of nonsingular tensor equations via randomized Kaczmarz-like method

Randomized extended average block Kaczmarz method for inconsistent tensor equations under t-product

A Multilinear Nyström Algorithm for Low-Rank Approximation of Tensors in Tucker Format

The accelerated tensor Kaczmarz algorithm with adaptive parameters for solving tensor systems

Multilinear Systems with <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" altimg="si424.gif" overflow="scroll"><mml:mi mathvariant="bold-script">ℳ</mml:mi></mml:math>-Tensors

A nested divide-and-conquer method for tensor Sylvester equations with positive definite hierarchically semiseparable coefficients

Accurate Numerical Solution for Structured M-matrix Algebraic Riccati Equations

Iterative Methods for Sparse Symmetric Multilinear Systems