where r(u) is as in (3.3), γ > 0 is an appropriate constant whose value depends on the selection of N and needs to be sufficiently large to bring the required stability to the new algorithm, and N = N(u) is the differential operator arising from the minimization of a convex functional such as