{"id":"arxiv/2605.21679","url":"https://arxiv.org/abs/2605.21679","tags":["math.NA"],"score":0,"title":"Component-wise accurate computation of the square root of an M-matrix","domain":"math.NA","source":"arxiv","authors":["Dario A. Bini","Bruno Iannazzo","Beatrice Meini","Jie Meng"],"reviews":0,"subject":"math.NA","summary":"Component-wise accurate algorithms for computing the principal square root of an M-matrix are designed in terms of triplet representations. A triplet representation of an M-matrix $A$ is the triple $(P, {\\bf u},{\\bf v})$, where the matrix $P$ is such that $p_{ij}=-a_{ij}$ for $i\\ne j$, $p_{ii}=0$, and ${\\bf u}>0$, ${\\bf v}\\ge 0$ are two vectors such that $A{\\bf u}={\\bf v}$. It is shown that if $A$ is an M-matrix representable by a triplet, then its principal square root exists and is an M-matrix represented by a triplet as well. New versions of the Cyclic Reduction and the Incremental Newton…","category":"research","comments":0,"disputed":false,"division":"math","mentions":0,"discussed":0,"published":"2026-05-20","reactions":{}}