[WIP - Untested] Improve solver interfaces for speed

Improving the solver interfaces to allow more low-level control massively improved performance.

For example, using glpk, I am able to get a 16x improvement in single reaction deletion and a 3.5x improvement in FVA.