Fix range parsing
Comment out safety checks in variable bounds
Prepare for normalisation to standard form
Standard form
Initial affine scaling
Reuse some allocation in affine scaline
Be smarter about residual
Factor out scaled constraint matrix, put explicit LS solve in repair
More robust conversion to affine scaling format

 * convert equality to tight (1d-6) ranges.
 * avoid issues with infty ranges
Explicit projection with normal matrix instead of DGGLSES
Try and be robust in face of numerical instability
Cholesky solve
Punt to LS solve on degenerate symmetric solve
Parameter tweaks
Cholesky for repair iterations as well
Cull unused work space slots
initial commit of SuiteSparse bindings
Detect factorisation issues in sparse cholesky
Use SuiteSparse as dense solver
More utility