Added adjoint Broyden method B as a limited memory variant.
This code runs, but is still buggy. More testing to come...
Bug fix in storage routine.
An improved starting condition: use Adjoint Broyden C correction of identity matrix so that resulting Jacobian estimate correctly reproduces the steepest descent direction for zeroing the constraints.
Added a version of the split quasi-Newton strategy with full-memory infeasibility Hessian.
Numerous fixes to limited memory Jacobian test.
Merge branch 'LM-Jacobian-dev' into tron-testing