<html><body>
<p>Matt,<br>
<br>
I have had no luck in reproducing fault in OsiClp. I have gcc 4.4.1 and so there may be a subtle difference. I did get fathom assert.<br>
<br>
Definitely a scaling issue. Fairly sure that it will go wrong without assert. fathom is trying to do LPs fast so it normally just keeps a scaled solution around. It has completed a sub search and thinks it has a solution. However on checking it finds that an integer variable is away from integrality. Whether it allows it through or not it may cut off a good solution. I have put in a fix to abort that sub search and carry on. This is in Clp/trunk and stable/2.4 i.e. Clp/1.11<br>
<br>
Other possibility is to set<br>
<br>
-depth -999<br>
<br>
which will switch off fathom.<br>
<br>
John <br>
<br>
<br>
<img width="16" height="16" src="cid:1__=0ABBFC2BDFCD24358f9e8a93df938@us.ibm.com" border="0" alt="Inactive hide details for Matthew Galati ---01/27/2010 10:14:57 AM---Similar problem setting for last post - different error. I"><font color="#424282">Matthew Galati ---01/27/2010 10:14:57 AM---Similar problem setting for last post - different error. I ran through gdb and ended up here. #4 0x</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__=0ABBFC2BDFCD24358f9e8a93df938@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__=0ABBFC2BDFCD24358f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2">Matthew Galati <Matthew.Galati@sas.com></font></td></tr>
<tr valign="top"><td width="1%"><img width="96" height="1" src="cid:2__=0ABBFC2BDFCD24358f9e8a93df938@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__=0ABBFC2BDFCD24358f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2">"clp@list.coin-or.org" <clp@list.coin-or.org></font></td></tr>
<tr valign="top"><td width="1%"><img width="96" height="1" src="cid:2__=0ABBFC2BDFCD24358f9e8a93df938@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__=0ABBFC2BDFCD24358f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2">01/27/2010 10:14 AM</font></td></tr>
<tr valign="top"><td width="1%"><img width="96" height="1" src="cid:2__=0ABBFC2BDFCD24358f9e8a93df938@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__=0ABBFC2BDFCD24358f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2">[Clp] assertion failure in ClpSimplex::fathom</font></td></tr>
<tr valign="top"><td width="1%"><img width="96" height="1" src="cid:2__=0ABBFC2BDFCD24358f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2" color="#5F5F5F">Sent by:</font></td><td width="100%"><img width="1" height="1" src="cid:2__=0ABBFC2BDFCD24358f9e8a93df938@us.ibm.com" border="0" alt=""><br>
<font size="2">clp-bounces@list.coin-or.org</font></td></tr>
</table>
<hr width="100%" size="2" align="left" noshade style="color:#8091A5; "><br>
<br>
<br>
<tt>Similar problem setting for last post - different error. I ran through gdb and ended up here.<br>
<br>
<br>
#4 0x08308b44 in ClpSimplex::fathom (this=0xcdb7698, stuff=0xbfee59e0) at ../../../Clp/src/ClpSimplex.cpp:10079<br>
10079 assert (fabs(value-value2)<1.0e-4);<br>
(gdb) p value<br>
$1 = -40885755.234334067<br>
(gdb) p value2<br>
$2 = -40885755<br>
(gdb)<br>
<br>
<br>
Is this just a tolerance issue? Maybe that assertion should check relative value, not absolute?<br>
<br>
<br>
<br>
<br>
DW GAP = 0.00196819 isTight = 1<br>
D-ALGO : 691 [CPU: 676.8 ] --- phaseUpdate() -----------> funcT = 0<br>
D-ALGOPC : 691 [CPU: 677 ] <--- phaseDone() -------------<br>
D-ALGOPC : 691 [CPU: 677 ] <--- solutionUpdateAsIp() ----<br>
decomp_milpblock: ../../../Clp/src/ClpSimplex.cpp:10079: int ClpSimplex::fathom(void*): Assertion `fabs(value-value2)<1.0e-4' failed.<br>
<br>
Program received signal SIGABRT, Aborted.<br>
0x0092a7a2 in _dl_sysinfo_int80 () from /lib/ld-linux.so.2<br>
(gdb) where<br>
#0 0x0092a7a2 in _dl_sysinfo_int80 () from /lib/ld-linux.so.2<br>
#1 0x0096e7d5 in raise () from /lib/tls/libc.so.6<br>
#2 0x00970149 in abort () from /lib/tls/libc.so.6<br>
#3 0x00967db1 in __assert_fail () from /lib/tls/libc.so.6<br>
#4 0x08308b44 in ClpSimplex::fathom (this=0xcdb7698, stuff=0xbfee59e0) at ../../../Clp/src/ClpSimplex.cpp:10079<br>
#5 0x0820ffbf in CbcModel::doOneNode (this=0xbde8ba8, baseModel=0xbde8ba8, node=@0xbfee5e58, newNode=@0xbfee5e40)<br>
at ../../../Cbc/src/CbcModel.cpp:13531<br>
#6 0x081eaf92 in CbcModel::branchAndBound (this=0xbde8ba8, doStatistics=0) at ../../../Cbc/src/CbcModel.cpp:3686<br>
#7 0x081141d1 in CbcMain1 (argc=5, argv=0xbfeeb780, model=@0xbfeeb7d0, callBack=0x80f91f4 <dummyCallBack>) at ../../../Cbc/src/CbcSolver.cpp:8214<br>
#8 0x080f9217 in CbcMain1 (argc=5, argv=0xbfeeb780, model=@0xbfeeb7d0) at ../../../Cbc/src/CbcSolver.cpp:3104<br>
#9 0x08095d4c in DecompAlgoPC::solutionUpdateAsIP (this=0xab98eb8) at ../../../Dip/src/DecompAlgoPC.cpp:502<br>
#10 0x080956d5 in DecompAlgoPC::phaseDone (this=0xab98eb8) at ../../../Dip/src/DecompAlgoPC.cpp:428<br>
#11 0x0806da48 in DecompAlgo::processNode (this=0xab98eb8, nodeIndex=1, globalLB=-17352638.898165066, globalUB=-16685688.559999999)<br>
at ../../../Dip/src/DecompAlgo.cpp:1788<br>
#12 0x080c53d6 in AlpsDecompTreeNode::process (this=0xb9ebb48, isRoot=false, rampUp=false) at ../../../Dip/src/AlpsDecompTreeNode.cpp:244<br>
#13 0x080dfa46 in AlpsSubTree::exploreUnitWork (this=0xb937550, leaveAsIt=false, unitWork=2147483647, unitTime=1.7976931348623157e+308,<br>
exploreStatus=@0xbfeec0ac, numNodesProcessed=@0xbfeec27c, numNodesBranched=@0xbfeec280, numNodesDiscarded=@0xbfeec284, depth=@0xbfeec290,<br>
betterSolution=@0xbfeec0ab) at ../../../Alps/src/AlpsSubTree.cpp:1276<br>
#14 0x080dca45 in AlpsSubTree::exploreSubTree (this=0xb937550, root=0xb937460, nodeLimit=2147483647, timeLimit=1.7976931348623157e+308,<br>
numNodesProcessed=@0xbfeec27c, numNodesBranched=@0xbfeec280, numNodesDiscarded=@0xbfeec284, depth=@0xbfeec290)<br>
at ../../../Alps/src/AlpsSubTree.cpp:472<br>
#15 0x080e5b80 in AlpsKnowledgeBrokerSerial::rootSearch (this=0xbfeec180, root=0xb937460) at ../../../Alps/src/AlpsKnowledgeBrokerSerial.cpp:220<br>
#16 0x080c29f1 in AlpsKnowledgeBroker::search (this=0xbfeec180, model=0xbfeec3e0) at ../../../Alps/src/AlpsKnowledgeBroker.h:275<br>
#17 0x080c09f6 in AlpsDecompModel::solve (this=0xbfeec3e0) at ../../../Dip/src/AlpsDecompModel.cpp:154<br>
#18 0x0804c13e in main (argc=1, argv=0xbfeec834) at ../../../../Dip/examples/MILPBlock/MILPBlock_Main.cpp:105<br>
<br>
<br>
_______________________________________________<br>
Clp mailing list<br>
Clp@list.coin-or.org<br>
</tt><tt><a href="http://list.coin-or.org/mailman/listinfo/clp">http://list.coin-or.org/mailman/listinfo/clp</a></tt><tt><br>
</tt><br>
<br>
</body></html>