Add a residuals function prototype.
Add the residuals() function.
Add function detectNewtonOscillations.
Finish the first stablizeNewton().
Revising the function stablizeNewton().

For the SOR type relaxation, we need the dx from the previous iteration.
Adding debugging output.
Merge branch 'limited_variable_update' of github.com:GitPaean/opm-autodiff into oscillation_treatment_withlimitedupdate

Conflicts:
	opm/autodiff/FullyImplicitBlackoilSolver.hpp
Commit to save changes.
Finishing the workflow.
Merge branch 'master' into oscillation_treatment_withlimitedupdate

Conflicts:
	opm/autodiff/FullyImplicitBlackoilSolver.hpp
Removing two input paramters in stablizeNewton().

Removing the bool input parameters osicllate and stagnate from the
stablizeNewton(). Basically, the value omega will decide if the
relaxation will be applied.
Removing the bool paramters calling stablizeNewton.
Removing one blank line.
Removing the space before ().
Merge branch 'master' into oscillation_treatment_withlimitedupdate

Conflicts:
	opm/autodiff/FullyImplicitBlackoilSolver.hpp
Rewriting the detection function with for loop.

It is more automatic while it remains to change to some more flexible
form.
Correcting the determination of stagnate.
Removing some debugging output.
Removing some debugging output.
Removing some spaces.