LAUSR

Since the eigenvalue problem of real supersymmetric tensors was proposed, there have been many results regarding the numerical algorithms for computing the spectral radius of nonnegative tensors, among which the NQZ method is the most studied. However, the NQZ method is only suitable for calculating the spectral radius of a special weakly primitive tensor, or a weakly irreducible primitive tensor that satisfies certain conditions. In this paper, by means of diagonal similarrity transformation of tensors, we construct a numerical algorithm for computing the spectral radius of nonnegative tensors with the aid of power functions. This algorithm is suitable for the calculation of the spectral radius of all weakly irreducible nonnegative tensors. Furthermore, it is efficient and can be widely applied.