Accuracy and componentwise accuracy in multilinear PageRank
Abstract
We study the stability with respect to perturbations and the accuracy of numerical algorithms for computing solutions to the multilinear PageRank problem x = (1-α)v + α P x2. 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.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.