Step 1 Corresponds to solve the non-linear system:
Step 2 is computed by the closed formula:
A continuation method over the parameter avoids to be “early” trapped by a local minimal.
This is, we start with a quadratic function ( large) and we adjust the functional by systematically decreasing
with a developed Quasi-Gauss Newton algorithm suitable for memory limited applications