
<HTML><BODY>

Hello,<div><br></div><div>I am trying to solve a relatively simple quadratic problem with the barrier method, but I seem to be hitting some really strange behavior. The problem itself consists of diagonal matrix for the quadratic part and has a number of linear constraints that may vary between different instances of the problem.</div><div>I have tried the following 2 ways of calling the solver:</div><div><br></div><div>1. Variant 1:</div><div>ClpSolve options;</div><div>options.setSolveType(ClpSolve::SolveType::useBarrier);</div><div>options.setPresolveType(ClpSolve::PresolveType::presolveOn);</div><div><br></div><div>ClpSimplex model;</div><div>// ... do some loading of the model, from file or from externally provided arrays</div><div><br></div><div><div>model.setPrimalTolerance(1e-6);</div><div>model.setDualTolerance(1e-6);</div></div><div>model.initialSolve(options);</div><div><br></div><div><br></div><div>2. Variant 2:</div><div> ClpSimplex model;</div><div>// ... do some loading of the model, from file or from externally provided arrays</div><div><br></div><div><div>model.setPrimalTolerance(1e-6);</div><div>model.setDualTolerance(1e-6);</div></div><div>model.barrier(true);</div><div><br></div><div><br></div><div>I encountered an instance of the problem for which variant 1 returns a solution that violates the constraints (with status 0) while variant 2 correctly identifies it as infeasible (status 1). What is interesting here is that in case of variant 1 there is a line in the log that states the problem is infeasible, however the solver keeps working on it and eventually returns status 0 and a solution.</div><div><br></div><div>I have also encountered another instance of the problem for which variant 1 returns the correct solution while variant 2 returns a solution that violates the constraints (with status 0).</div><div><br></div><div>At this point I am pretty confident I am missing something in the way I am invoking the solver, but I do not seem to be able to figure out what.</div><div>Thanks in advance for any help.</div><div><br></div></BODY></HTML>