<html><body>
<p>Matt,<br>
<br>
getBasics returns an array number of rows long. The entry for row 7 would say which variable was actually pivoting on row 7. If you are looking at tableau you may need that.<br>
<br>
warmstart just gives status of each variable with no information as to which row a basic variable would pivot on.<br>
<br>
John<br>
<br>
<img width="16" height="16" src="cid:1__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt="Inactive hide details for Matthew Galati ---09/10/2009 12:39:27 PM---Thanks John. When you say "different" do you mean actually"><font color="#424282">Matthew Galati ---09/10/2009 12:39:27 PM---Thanks John. When you say "different" do you mean actually different or just the way they are stored</font><br>
<br>
<table width="100%" border="0" cellspacing="0" cellpadding="0">
<tr valign="top"><td width="1%"><img width="96" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2" color="#5F5F5F">From:</font></td><td width="100%"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2">Matthew Galati <magh@lehigh.edu></font></td></tr>
<tr valign="top"><td width="1%"><img width="96" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2" color="#5F5F5F">To:</font></td><td width="100%"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2">John J Forrest/Watson/IBM@IBMUS</font></td></tr>
<tr valign="top"><td width="1%"><img width="96" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2" color="#5F5F5F">Cc:</font></td><td width="100%" valign="middle"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2">"cbc@list.coin-or.org" <cbc@list.coin-or.org></font></td></tr>
<tr valign="top"><td width="1%"><img width="96" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2" color="#5F5F5F">Date:</font></td><td width="100%"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2">09/10/2009 12:39 PM</font></td></tr>
<tr valign="top"><td width="1%"><img width="96" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2" color="#5F5F5F">Subject:</font></td><td width="100%"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2">Re: [Cbc] failure in small IP</font></td></tr>
</table>
<hr width="100%" size="2" align="left" noshade style="color:#8091A5; "><br>
<br>
<br>
<font size="4">Thanks John. When you say "different" do you mean actually different or just the way they are stored?<br>
<br>
I.e., if I want to know which variables are basic after a solve - and I getWarmStart and check the structural status - will that information match what getBasics would return?<br>
<br>
</font><br>
<ul><font size="4">Matt,<br>
<br>
Thanks.<br>
<br>
As to getBasics - the two sets of information are different - but I will modify to point out options to user.<br>
<br>
John<br>
<br>
</font><img src="cid:1__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" width="16" height="16" alt="Inactive hide details for Matthew Galati ---09/10/2009 12:37:54 AM---Here is a fix. It was just missing an initialization in Cb"><font size="4" color="#424282">Matthew Galati ---09/10/2009 12:37:54 AM---Here is a fix. It was just missing an initialization in CbcRounding constructor. [magala@orclus71 sr</font><font size="4"><br>
</font></ul>
<br>
<table width="100%" border="0" cellspacing="0" cellpadding="0">
<tr valign="top"><td width="20%" valign="middle"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""></td><td width="43%"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""></td><td width="36%"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""></td></tr>
<tr valign="top"><td width="20%"><img src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" width="96" height="1"><font color="#5F5F5F"><br>
From:</font></td><td width="43%"><img src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" width="1" height="1"><br>
Matthew Galati <<a href="mailto:Matthew.Galati@sas.com" target="_blank"><u><font color="#0000FF">Matthew.Galati@sas.com</font></u></a>></td><td width="36%"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""></td></tr>
<tr valign="top"><td width="20%"><img src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" width="96" height="1"><font color="#5F5F5F"><br>
To:</font></td><td width="43%"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""></td><td width="36%"><img src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" width="1" height="1"><br>
Matthew Galati <<a href="mailto:magh@lehigh.edu" target="_blank"><u><font color="#0000FF">magh@lehigh.edu</font></u></a>></td></tr>
<tr valign="top"><td width="20%"><img src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" width="96" height="1"><font color="#5F5F5F"><br>
Cc:</font></td><td width="43%" valign="middle"><img src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" width="1" height="1"><br>
"<a href="mailto:cbc@list.coin-or.org" target="_blank"><u><font color="#0000FF">cbc@list.coin-or.org</font></u></a>" <<a href="mailto:cbc@list.coin-or.org" target="_blank"><u><font color="#0000FF">cbc@list.coin-or.org</font></u></a>></td><td width="36%" valign="middle"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""></td></tr>
<tr valign="top"><td width="20%"><img src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" width="96" height="1"><font color="#5F5F5F"><br>
Date:</font></td><td width="43%"><img src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" width="1" height="1"><br>
09/10/2009 12:37 AM</td><td width="36%" valign="middle"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""></td></tr>
<tr valign="top"><td width="20%" valign="middle"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""></td><td width="43%" valign="middle"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""></td><td width="36%" valign="middle"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""></td></tr>
<tr valign="top"><td width="20%"><img src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" width="96" height="1"><font color="#5F5F5F"><br>
Subject:</font></td><td width="43%"><img src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" width="1" height="1"><br>
Re: [Cbc] failure in small IP</td><td width="36%" valign="middle"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""></td></tr>
<tr valign="top"><td width="20%"><img src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" width="96" height="1"><font color="#5F5F5F"><br>
Sent by:</font></td><td width="43%"><img src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" width="1" height="1"><u><font color="#0000FF"><br>
</font></u><a href="mailto:cbc-bounces@list.coin-or.org" target="_blank"><u><font color="#0000FF">cbc-bounces@list.coin-or.org</font></u></a></td><td width="36%" valign="middle"><img width="1" height="1" src="cid:2__=0ABBFCBEDFCFB3FE8f9e8a93df938@us.ibm.com" border="0" alt=""></td></tr>
</table>
<ul><hr width="100%" size="2" align="left" noshade><br>
<font size="4"><br>
<br>
</font><font size="4" color="#1F497D"><br>
Here is a fix. It was just missing an initialization in CbcRounding constructor.</font><font size="4"><br>
<br>
<br>
</font><font size="4" color="#1F497D"><br>
[magala@orclus71 src]$ svn diff CbcHeuristic.cpp<br>
Index: CbcHeuristic.cpp<br>
===================================================================<br>
--- CbcHeuristic.cpp (revision 1223)<br>
+++ CbcHeuristic.cpp (working copy)<br>
@@ -1538,6 +1538,9 @@<br>
{<br>
// Get a copy of original matrix (and by row for rounding);<br>
assert(model.solver());<br>
+ down_=NULL;<br>
+ up_=NULL;<br>
+ equal_=NULL;<br>
if (model.solver()->getNumRows()) {<br>
matrix_ = *model.solver()->getMatrixByCol();<br>
matrixByRow_ = *model.solver()->getMatrixByRow();</font><font size="4"><br>
<br>
<br>
</font><b><font size="4" face="Tahoma"><br>
From:</font></b><font size="4" face="Tahoma"> Matthew Galati [</font><a href="mailto:magh@lehigh.edu" target="_blank"><u><font size="4" color="#0000FF" face="Tahoma">mailto:magh@lehigh.edu</font></u></a><font size="4" face="Tahoma">] </font><b><font size="4" face="Tahoma"><br>
Sent:</font></b><font size="4" face="Tahoma"> Thursday, September 10, 2009 12:19 AM</font><b><font size="4" face="Tahoma"><br>
To:</font></b><font size="4" face="Tahoma"> Matthew Galati</font><b><font size="4" face="Tahoma"><br>
Cc:</font></b><font size="4" face="Tahoma"> </font><a href="mailto:cbc@list.coin-or.org" target="_blank"><u><font size="4" color="#0000FF" face="Tahoma">cbc@list.coin-or.org</font></u></a><b><font size="4" face="Tahoma"><br>
Subject:</font></b><font size="4" face="Tahoma"> Re: [Cbc] failure in small IP</font><font size="4"><br>
</font><font size="5"><br>
The run with cbc executable that does not fail was using a May 2009 version of code. <br>
<br>
Using the current trunk it fails (same as when using CbcMain1):</font><font size="4"> </font>
<p><font size="5">[magala@orclus71 src]$ ./cbc ~/subProb.n0.c0.p1.b0.lp</font><font size="4"> </font>
<p><font size="5">Coin Cbc and Clp Solver version 1.2trunk, build Sep 10 2009</font><font size="4"> </font>
<p><font size="5">command line - ./cbc /users/magala/subProb.n0.c0.p1.b0.lp (default strategy 1)</font><font size="4"> </font>
<p><font size="5">Continuous objective value is 46.6667 - 0.00 seconds</font><font size="4"> </font>
<p><font size="5">Cgl0003I 20 fixed, 0 tightened bounds, 0 strengthened rows, 0 substitutions</font><font size="4"> </font>
<p><font size="5">Cgl0003I 5 fixed, 0 tightened bounds, 0 strengthened rows, 8 substitutions</font><font size="4"> </font>
<p><font size="5">Cgl0003I 0 fixed, 0 tightened bounds, 0 strengthened rows, 8 substitutions</font><font size="4"> </font>
<p><font size="5">Cgl0004I processed model has 0 rows, 0 columns (0 integer) and 0 elements</font><font size="4"> </font>
<p><font size="5">Segmentation fault<br>
<br>
</font>
<p><font size="5">The following small IP fails using the CbcMain1 interface. It is a unique case, because all rows and cols are removed by presolver.</font><font size="4"> </font>
<p><a href="http://coral.ie.lehigh.edu/%7Emagh/tmp/subProb.n0.c0.p1.b0.lp" target="_blank"><u><font size="5" color="#0000FF">http://coral.ie.lehigh.edu/~magh/tmp/subProb.n0.c0.p1.b0.lp</font></u></a><font size="4"> </font>
<p><a href="https://projects.coin-or.org/Cbc/ticket/79" target="_blank"><u><font size="5" color="#0000FF">https://projects.coin-or.org/Cbc/ticket/79</font></u></a><font size="4"> </font>
<p><font size="5">CbcModel cbc(*subprobSI);</font><font size="4"> </font>
<p><font size="5">CbcMain0(cbc);</font><font size="4"> </font>
<p><font size="5">const char * argv[] = {"cbc","-solve","-quit"};</font><font size="4"> </font>
<p><font size="5">CbcMain1(3, argv, cbc);</font><font size="4"> </font>
<p><font size="5">Coin Cbc and Clp Solver version 1.2trunk, build Sep 9 2009</font><font size="4"> </font>
<p><font size="5">command line - cbc -solve -quit (default strategy 1)</font><font size="4"> </font>
<p><font size="5">Continuous objective value is 46.6667 - 0.00 seconds</font><font size="4"> </font>
<p><font size="5">Cgl0003I 20 fixed, 0 tightened bounds, 0 strengthened rows, 0 substitutions</font><font size="4"> </font>
<p><font size="5">Cgl0003I 5 fixed, 0 tightened bounds, 0 strengthened rows, 8 substitutions</font><font size="4"> </font>
<p><font size="5">Cgl0003I 0 fixed, 0 tightened bounds, 0 strengthened rows, 8 substitutions</font><font size="4"> </font>
<p><font size="5">Cgl0004I processed model has 0 rows, 0 columns (0 integer) and 0 elements</font><font size="4"> </font>
<p><font size="5">gdb gives:</font><font size="4"> </font>
<p><font size="5">Program received signal SIGSEGV, Segmentation fault.</font><font size="4"> </font>
<p><font size="5">0x0000003ffa871f6b in free () from /lib64/libc.so.6</font><font size="4"> </font>
<p><font size="5">(gdb) where</font><font size="4"> </font>
<p><font size="5">#0 0x0000003ffa871f6b in free () from /lib64/libc.so.6</font><font size="4"> </font>
<p><font size="5">#1 0x000000000057bc41 in ~CbcRounding (this=0x7fffec90c080)</font><font size="4"> </font>
<p><font size="5">at ../../../Cbc/src/CbcHeuristic.cpp:1552</font><font size="4"> </font>
<p><font size="5">#2 0x00000000004ae837 in doHeuristics (model=0x1d1f3e00, type=1)</font><font size="4"> </font>
<p><font size="5">at ../../../Cbc/src/CbcSolver.cpp:3479</font><font size="4"> </font>
<p><font size="5">#3 0x00000000004c7776 in CbcMain1 (argc=3, argv=0x7fffec9150c0,</font><font size="4"> </font>
<p><font size="5">model=@0x7fffec914bb0, callBack=0x4a734c <dummyCallBack>)</font><font size="4"> </font>
<p><font size="5">at ../../../Cbc/src/CbcSolver.cpp:6680</font><font size="4"> </font>
<p><font size="5">#4 0x00000000004db968 in CbcMain1 (argc=3, argv=0x7fffec9150c0,</font><font size="4"> </font>
<p><font size="5">model=@0x7fffec914bb0) at ../../../Cbc/src/CbcSolver.cpp:3099</font><font size="4"> </font>
<p><font size="5">Note, if I run it from the command line using cbc executable, the error does not occur.</font><font size="4"> </font>
<p><font size="5">[magala@orclus71 MMKP]$ ~/COIN/coin-Cbc/build-g/bin/cbc subProb.n0.c0.p1.b0.lp</font><font size="4"> </font>
<p><font size="5">Coin Cbc and Clp Solver version 1.2trunk, build May 27 2009</font><font size="4"> </font>
<p><font size="5">command line - /users/magala/COIN/coin-Cbc/build-g/bin/cbc subProb.n0.c0.p1.b0.lp (default strategy 1)</font><font size="4"> </font>
<p><font size="5">Continuous objective value is 46.6667 - 0.00 seconds</font><font size="4"> </font>
<p><font size="5">Cgl0003I 20 fixed, 0 tightened bounds, 0 strengthened rows, 0 substitutions</font><font size="4"> </font>
<p><font size="5">Cgl0003I 5 fixed, 0 tightened bounds, 0 strengthened rows, 8 substitutions</font><font size="4"> </font>
<p><font size="5">Cgl0003I 0 fixed, 0 tightened bounds, 0 strengthened rows, 8 substitutions</font><font size="4"> </font>
<p><font size="5">Cgl0004I processed model has 0 rows, 0 columns (0 integer) and 0 elements</font><font size="4"> </font>
<p><font size="5">Cbc3007W No integer variables - nothing to do</font><font size="4"> </font>
<p><font size="5">Cuts at root node changed objective from 46.6667 to -1.79769e+308</font><font size="4"> </font>
<p><font size="5">Probing was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)</font><font size="4"> </font>
<p><font size="5">Gomory was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)</font><font size="4"> </font>
<p><font size="5">Knapsack was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)</font><font size="4"> </font>
<p><font size="5">Clique was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)</font><font size="4"> </font>
<p><font size="5">MixedIntegerRounding2 was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)</font><font size="4"> </font>
<p><font size="5">FlowCover was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)</font><font size="4"> </font>
<p><font size="5">TwoMirCuts was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)</font><font size="4"> </font>
<p><font size="5">Result - Finished objective 46.66667 after 0 nodes and 0 iterations - took 0.00 seconds (total time 0.01)</font><font size="4"> </font>
<p><font size="5">Total time 0.01<br>
<br>
_______________________________________________<br>
Cbc mailing list</font><u><font size="4" color="#0000FF"><br>
</font></u><a href="mailto:Cbc@list.coin-or.org" target="_blank"><u><font size="5" color="#0000FF">Cbc@list.coin-or.org</font></u></a><u><font size="4" color="#0000FF"><br>
</font></u><a href="http://list.coin-or.org/mailman/listinfo/cbc" target="_blank"><u><font size="5" color="#0000FF">http://list.coin-or.org/mailman/listinfo/cbc</font></u></a><tt><font size="4"><br>
_______________________________________________<br>
Cbc mailing list</font></tt><tt><u><font size="4" color="#0000FF"><br>
</font></u></tt><a href="mailto:Cbc@list.coin-or.org" target="_blank"><tt><u><font size="4" color="#0000FF">Cbc@list.coin-or.org</font></u></tt></a><tt><u><font size="4" color="#0000FF"><br>
</font></u></tt><a href="http://list.coin-or.org/mailman/listinfo/cbc" target="_blank"><tt><u><font size="4" color="#0000FF">http://list.coin-or.org/mailman/listinfo/cbc</font></u></tt></a><font size="4"><br>
<br>
</font><br>
<font size="4"><br>
_______________________________________________<br>
Cbc mailing list</font><u><font size="4" color="#0000FF"><br>
</font></u><a href="mailto:Cbc@list.coin-or.org"><u><font size="4" color="#0000FF">Cbc@list.coin-or.org</font></u></a><u><font size="4" color="#0000FF"><br>
</font></u><a href="http://list.coin-or.org/mailman/listinfo/cbc" target="_blank"><u><font size="4" color="#0000FF">http://list.coin-or.org/mailman/listinfo/cbc</font></u></a><font size="4"><br>
</font></ul>
<br>
<br>
</body></html>