Hello
I'm implementing a simple Benders decomposition problem and I would like to access the Farkas infeasibility proof. According to the Gurobi API, there is an attribute of the constraint called FarkasDual. I was looking at the plugins folder of pyomo and I noticed the gurobi_direct.py
I was thinking I could add something like
if (soln.status == SolutionStatus.infeasible):
if extract_duals:
for con in cons:
# FarkasDual attributes in Gurobi is the Farkas infeasibility proof
soln_constraints[ con.ConstrName ]["Dual"] = con.FarkasDual
to get the unbounded ray in the dual suffix.
Now I have two questions:
1. How do I use the gurobi_direct plugin?
opt = SolverFactory('gurobi_direct') ?
thanks!
Braulio