Farkas duals for variable bounds

[odow]

Currently we only get the farkas dual for constraints:

LinQuadOptInterface.jl/src/solve.jl

Line 58 in d6957f3

get_farkas_dual!(m, m.constraint_dual_solution)

However, linear12 needs it for variable bounds.

[joaquimg]

I talked about this with @blegat once. Part of the problem is due to the fact that xpress gurobi and cplex only return data for constraints. They ignore bounds.
We can have a method for computing that.

[blegat]

Yes, this is the dual analogous to computing the ConstraintPrimals using the VariablePrimals and the ConstraintFunctions

[odow]

I don't really know what I'm talking about, but for the linear case, it seems we just need a matrix multiplication: https://orinanobworld.blogspot.com/2011/07/farkas-certificates-in-cplex.html

get_farkas_dual!(m, m.constraint_dual_solution)
A = get_constraint_matrix(m)  # doesn't exist
m.variable_dual_solution .= m.constraint_dual_solution' * A
# plus some +/- corrections

It seems we have two options

1. add a get_constraint_matrix(model) method
2. add a get_farkas_dual_variables!(model, m.variable_dual_solution) method
   2a. with a fallback that loops through the linear constraints querying the constraint functions.

[joaquimg]

basic infra structure is on #59 should we close?