A nonnegativity preserving algorithm for multilinear systems with nonsingular ℳ ${\mathcal M}$ -tensors