Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

getobjbound for conic_to_LPQP #338

Closed
chriscoey opened this issue Mar 29, 2017 · 13 comments
Closed

getobjbound for conic_to_LPQP #338

chriscoey opened this issue Mar 29, 2017 · 13 comments

Comments

@chriscoey
Copy link
Contributor

@mlubin, in new QP tests

QP maximize: Error During Test
  Got an exception of type MethodError outside of a @test
  MethodError: no method matching getobjbound(::MathProgBase.SolverInterface.ConicToLPQPBridge)
  Closest candidates are:
    getobjbound(!Matched::GLPKMathProgInterface.GLPKInterfaceMIP.GLPKMathProgModelMIP) at /home/coey/.julia/v0.5/GLPKMathProgInterface/src/GLPKInterfaceMIP.jl:549
    getobjbound(!Matched::Cbc.CbcMathProgSolverInterface.CbcMathProgModel) at /home/coey/.julia/v0.5/Cbc/src/CbcSolverInterface.jl:190
    getobjbound(!Matched::Mosek.MosekMathProgSolverInterface.MosekLinearQuadraticModel) at /home/coey/.julia/v0.5/Mosek/src/MosekLPQCQPInterface.jl:661

@chriscoey
Copy link
Contributor Author

also one test isn't passing there because the answer is a factor of 2 off. it could be caused by maybe forgetting the 2 in the RSOC?

getobjectivevalue(m) = 8.022765808609499
QP large (cardls): Test Failed
  Expression: isapprox(getobjectivevalue(m),16.045564,atol=TOL)

@mlubin
Copy link
Member

mlubin commented Mar 29, 2017

Try out JuliaOpt/MathProgBase.jl#160 for the getobjbound issue.

With

@variable(m, w == 2)
@constraint(m, sum(t.^2) <= u*w)

You get u >= sum(t.^2)/2. Is that what you want?

@chriscoey
Copy link
Contributor Author

chriscoey commented Mar 29, 2017 via email

@mlubin
Copy link
Member

mlubin commented Mar 29, 2017

So the objective is to minimize sum(t.^2)/2 + rho*sum(x.^2)/2?

@chriscoey
Copy link
Contributor Author

chriscoey commented Mar 29, 2017 via email

@mlubin
Copy link
Member

mlubin commented Mar 29, 2017

The constraints are being enforced as written:

julia> sum(getvalue(t).^2)/2 + rho*sum(getvalue(x).^2)/2
8.022766055303913

julia> getobjectivevalue(m)
8.022766159395793

@chriscoey
Copy link
Contributor Author

chriscoey commented Mar 29, 2017 via email

@mlubin
Copy link
Member

mlubin commented Mar 29, 2017

That suggests an issue with the conic formulation, not conic to lpqp.

@chriscoey
Copy link
Contributor Author

chriscoey commented Mar 29, 2017 via email

@mlubin
Copy link
Member

mlubin commented Mar 29, 2017

What solver is reporting 16.045564?

@chriscoey
Copy link
Contributor Author

chriscoey commented Mar 29, 2017 via email

@mlubin
Copy link
Member

mlubin commented Mar 29, 2017 via email

@chriscoey
Copy link
Contributor Author

OK this seems to be fixed.

Tests passing when on MPB master.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Development

No branches or pull requests

2 participants