<div dir="ltr">Hi all,<div><br></div><div>I have an LP that is (almost certainly) primal infeasible, but Clp (revision 2282) seems to struggle with it when presolve is on, not solving within an hour, and perhaps cycling indefinitely. The LP is linked below both in LP format and IEEE hex MPS format.</div><div><br></div><div>LP format: <a href="https://drive.google.com/file/d/0B6wQ0pEa5jXyTnFVSnktb2JxbEU/view?usp=sharing">https://drive.google.com/file/d/0B6wQ0pEa5jXyTnFVSnktb2JxbEU/</a></div><div>IEEE hex MPS format: <a href="https://drive.google.com/file/d/0B6wQ0pEa5jXyR0txVmREcngybXM/">https://drive.google.com/file/d/0B6wQ0pEa5jXyR0txVmREcngybXM/</a></div><div><br></div><div>The commands from which I encounter the error are:</div><div>solver->readLp("solverSlow.lp");<br></div><div>solver->setHintParam(OsiDoPresolveInInitial, 1);</div><div>solver->initialSolve();</div><div><br></div><div>If I save the LP with normal or extra precision (instead of with IEEE hex), then primal infeasibility is declared within ten seconds or so. Similarly, if the line setting OsiDoPresolveInInitial to 1 is commented out, the model solves quickly. The full model with presolve on is what causes trouble, and in particular, the issue is encountered during postsolve. There is a step in which, during postsolve, the log transitions from</div><div><p style="margin:0px;font-stretch:normal;font-size:12px;line-height:normal;font-family:Helvetica"><br></p><p style="margin:0px;font-stretch:normal;font-size:12px;line-height:normal;font-family:Helvetica">Clp0006I 1  Obj 0 Primal inf 10.365599 (1)</p><p style="margin:0px;font-stretch:normal;font-size:12px;line-height:normal;font-family:Helvetica"><br></p>to<p style="margin:0px;font-stretch:normal;font-size:12px;line-height:normal;font-family:Helvetica"><br></p>
<p style="margin:0px;font-stretch:normal;font-size:12px;line-height:normal;font-family:Helvetica">Clp0006I 201  Obj 4.203557e-12 Primal inf 131157.64 (2351)</p><p style="margin:0px;font-stretch:normal;font-size:12px;line-height:normal;font-family:Helvetica"><br></p>and this latter type of entry is the one that continues indefinitely.</div><div><br></div><div>I am looking for advice on what parameter to set or cleaning of the LP to do in advance, as well as reporting the experience in case something can be changed in Clp to handle this type of case.</div><div><br></div><div><div>Here are some properties of this LP.</div><div>1. There are 4823 constraints, 668 variables, and 318711 non-zeroes in the coefficient matrix. </div><div>2. The coefficient range is 0.000178573 to 1155.46. The RHS for each row is either 0 or 1, and all constraints are >=.</div><div>3. The objective is all-zeroes. </div><div>4. There are many free variables.</div></div><div><br></div><div>I can imagine that the coefficient range may cause numerical instability, but it does not seem so terribly bad at first glance. Something I have encountered before is when coefficients are not small themselves, but the differences between the coefficients are very small, then some solvers encounter issues. The behavior I see may be related to this, as there are such "similar" coefficients in this model, and it is also suggested by the fact that saving the model with different precision so dramatically changes things.</div><div><br></div><div>Thank you in advance for your help,</div><div>Aleksandr Kazachkov</div></div>