Lagrange multipliers can be obtained using the “dual” Suffix. Depending on which Pyomo interface you use, some solvers also return a “slack” Suffix. However, you can also compute the slack after loading the solution using the lslack() and uslack() methods available on constraints objects.
If you use the NL file interface to Gurobi (SolverFactory(“gurobi”, solver_io=“nl”) or SolverFactory(“gurobi_ampl”)), you can obtain the unbounded rays using the “dunbdd” Suffix (explained in the output from "gurobi_ampl -=“). We also have a user that has submitted a patch to Pyomo’s Gurobi-Python interface to obtain this information, but that hasn’t been incorporated into trunk yet.
Gabe