Accuracy and componentwise accuracy in multilinear PageRank

We study the stability with respect to perturbations and the accuracy of numerical algorithms for computing solutions to the multilinear PageRank problem $\mathbf{x} = (1-α)\mathbf{v} + α\mathcal{P} \mathbf{x}^2$. Our results reveal that the solution can be more stable with respect to perturbations and numerical errors with respect to the classical bounds for nonlinear systems of equations (based on the norm of the Jacobian). In detail, one can obtain bounds for the minimal solution which ignore the singularity of the problem for $α=1/2$, and one can show that the limiting accuracy of the Newton method depends not on the norm of the Jacobian but on a quantity that can be much smaller thanks to the nonnegativity structure of the equation. For the minimal solution, we also suggest subtraction-free modifications to the existing algorithms to achieve componentwise stability. Some of the theoretical results we obtain are interesting even outside the scope of this problem: bounds for more general quadratic vector equations, and a partial inverse for M-matrices which remains bounded when the matrix to invert approaches singularity.