The Gauss-Newton scheme is solved in the multilinear engine by the conjugate Gradient algorithm,taking advantage of the sparse structure of the jacobian matrix of the multilinear model.in contrast, in the positive matrix factorization,the Gauss-Newton scheme is solved by Cholesky decomposition and back-substitution.