Make float accuracy fix the main path. Remove alternative path.
Added utilitied collapseJacs() and vertcat().

The utilities are not optimized for speed, but make it easy
to construct a linear system from all the jacobian blocks.
Make the test run a simple bhp-controlled case.
Added well equations.

Residuals are not yet used for loop control.