Opened 4 years ago
Last modified 4 years ago
#18733 new enhancement
MixedIntegerLinearProgram should support tableau-query methods getBInvARow, getBInvRow, getBInvACol, getBInvCol
Reported by: | mkoeppe | Owned by: | |
---|---|---|---|
Priority: | major | Milestone: | sage-6.8 |
Component: | numerical | Keywords: | lp |
Cc: | ncohen, yzh | Merged in: | |
Authors: | Reviewers: | ||
Report Upstream: | N/A | Work issues: | |
Branch: | Commit: | ||
Dependencies: | #18732 | Stopgaps: |
Description (last modified by )
If one uses Sage's MixedIntegerLinearProgram
class for solving an LP, it is useful to access dictionary (tableau) data corresponding to the current (often optimal) basis. For example, one needs this to compute tableau cutting planes such as Gomory's fractional cut.
The Sage interface should, of course, be designed to work consistently across all solvers. There are some subtleties here. Best to be compatible with a library that has already sorted most of it out: In the COIN open solver interface, there are the functions getBInvARow, getBInvRow, getBInvACol, getBInvCol. Unfortunately these functions do not seem to be implemented in Osi's GLPK backend; so one needs to use one of Osi more complete backends such as CPLEX or CLP to find the precise definition.
#18732 adds the necessary backend functions for the GLPK backend. #18763 adds the necessary backend functions for the COIN (CBC/CLP) backend.
This kind of information is available in every solver.
Change History (5)
comment:1 Changed 4 years ago by
- Type changed from PLEASE CHANGE to enhancement
comment:2 Changed 4 years ago by
- Description modified (diff)
comment:3 Changed 4 years ago by
The methods would raise NotImplementedError
or other errors when no tableau is available. Even if the backend uses the simplex method, these methods may raise an error when solve
has not been called yet.
comment:4 Changed 4 years ago by
- Description modified (diff)
comment:5 Changed 4 years ago by
- Description modified (diff)
Note that not all Sage's LP solvers use a simplex method. E.g. cvxopt does not. In this case is makes no sense talking about tableaux...