Add Newton: p_iteration_with_recalculated_jacobian.
Work on MUMPS.