It is also worth noting that (P P − I) 2 requires each projection P •,k to be a unit vector because of P •,k P •,k = 1, which is infeasible when the nonzero values in X are large.
we multiply the matrix I by a coefficient c, which is calculated from the mean of the diagonal of P P in the current iteration.