The Gauss–Newton scheme is solved in the ME-2 by the
Conjugate Gradient algorithm (Hestenes and Stiefel, 1952),
taking advantage of the sparse structure of the Jacobian matrix
of the multilinear model. In contrast, in the PMF2, the Gauss–
Newton scheme is solved by Cholesky decomposition and
back-substitution.