Accurate Numerical Solution for Shifted M -Matrix Algebraic Riccati Equations

Changli Liu,Jungong Xue,Ren-Cang Li
DOI: https://doi.org/10.1007/s10915-020-01263-4
2020-01-01
Journal of Scientific Computing
Abstract:An algebraic Riccati equation ( are ) is called a shifted M -matrix algebraic Riccati equation ( mare ) if it can be turned into an mare after its matrix variable is partially shifted by a diagonal matrix. Such an are can arise from computing the invariant density of a Markov modulated Brownian motion. Sufficient and necessary conditions for an are to be a shifted mare are obtained. Based on the conditions, a highly accurate implementation of the alternating directional doubling algorithm ( adda ) is established to compute the extremal solution of a shifted mare , as well as a quantity needed for computing the invariant density in the application, with high entrywise relative accuracy. Numerical examples are presented to demonstrate the theory and algorithms.
What problem does this paper attempt to address?