warmstart coded, not sure which solvers it should work with
merged changes to square
allow warmstart for non-MILPs
typo
dual vector populated in solution (but not in constraints)