Derivation of bounds from equalities is not done by default in nonlinear solver, unless bmibnb is used
Argument missing (caught in unit test)
Performance fixes