With updated ClpSimplex.cpp, ClpSimplexPrimal.cpp and ClpSimplexDual.cpp I get the same error still.<br><br>If there is something I can look for, through gdb or printing to log, to help - let me know.<br><br><br><br><div class="gmail_quote">
On Fri, Jan 29, 2010 at 5:03 AM, John J Forrest <span dir="ltr"><<a href="mailto:jjforre@us.ibm.com">jjforre@us.ibm.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div>
<p>Matt,<br>
<br>
Can you get ClpSimplex*.cpp from svn stable/1.11 and try again. By luck I managed to get something close to your problem. It was on a branch where (luckily) I had not put in the fix I did to get round the assert you also mentioned. I had everything in debug mode apart from one function which I had copied across optimized. Interestingly this really stressed the code! I think the reason was that some computations were always 64 bit floating point, while optimized can have some at 80 bits. This seemed to set up an accuracy problem.<br>
<br>
Anyway I got the assert and tracking it down I found a serious (but unlikely to occur) bug which looks like your one.<br>
<br>
So try it and tell me if it helps.</p><div class="im"><br>
<br>
John<br>
<br>
<img src="cid:1__=0ABBFC29DFA77FDE8f9e8a93df938@us.ibm.com" alt="Inactive hide details for Matthew Galati ---01/26/2010 10:12:59 PM---This is based on CLP release 1.11.0. If I dump to MPS file" width="16" border="0" height="16"><font color="#424282">Matthew Galati ---01/26/2010 10:12:59 PM---This is based on CLP release 1.11.0. If I dump to MPS file, CLP claims the problem is infeasible - w</font><br>
<br>
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tbody><tr valign="top"><td width="1%"><img src="cid:2__=0ABBFC29DFA77FDE8f9e8a93df938@us.ibm.com" alt="" width="96" border="0" height="1"><br>
<font color="#5f5f5f" size="2">From:</font></td><td width="100%"><img src="cid:2__=0ABBFC29DFA77FDE8f9e8a93df938@us.ibm.com" alt="" width="1" border="0" height="1"><br>
<font size="2">Matthew Galati <<a href="mailto:Matthew.Galati@sas.com" target="_blank">Matthew.Galati@sas.com</a>></font></td></tr>
<tr valign="top"><td width="1%"><img src="cid:2__=0ABBFC29DFA77FDE8f9e8a93df938@us.ibm.com" alt="" width="96" border="0" height="1"><br>
<font color="#5f5f5f" size="2">To:</font></td><td width="100%"><img src="cid:2__=0ABBFC29DFA77FDE8f9e8a93df938@us.ibm.com" alt="" width="1" border="0" height="1"><br>
<font size="2">"<a href="mailto:clp@list.coin-or.org" target="_blank">clp@list.coin-or.org</a>" <<a href="mailto:clp@list.coin-or.org" target="_blank">clp@list.coin-or.org</a>></font></td></tr>
<tr valign="top"><td width="1%"><img src="cid:2__=0ABBFC29DFA77FDE8f9e8a93df938@us.ibm.com" alt="" width="96" border="0" height="1"><br>
<font color="#5f5f5f" size="2">Date:</font></td><td width="100%"><img src="cid:2__=0ABBFC29DFA77FDE8f9e8a93df938@us.ibm.com" alt="" width="1" border="0" height="1"><br>
<font size="2">01/26/2010 10:12 PM</font></td></tr>
<tr valign="top"><td width="1%"><img src="cid:2__=0ABBFC29DFA77FDE8f9e8a93df938@us.ibm.com" alt="" width="96" border="0" height="1"><br>
<font color="#5f5f5f" size="2">Subject:</font></td><td width="100%"><img src="cid:2__=0ABBFC29DFA77FDE8f9e8a93df938@us.ibm.com" alt="" width="1" border="0" height="1"><br>
<font size="2">[Clp] status wrong using OsiClp</font></td></tr>
<tr valign="top"><td width="1%"><img src="cid:2__=0ABBFC29DFA77FDE8f9e8a93df938@us.ibm.com" alt="" width="96" border="0" height="1"><br>
<font color="#5f5f5f" size="2">Sent by:</font></td><td width="100%"><img src="cid:2__=0ABBFC29DFA77FDE8f9e8a93df938@us.ibm.com" alt="" width="1" border="0" height="1"><br>
<font size="2"><a href="mailto:clp-bounces@list.coin-or.org" target="_blank">clp-bounces@list.coin-or.org</a></font></td></tr>
</tbody></table>
<hr style="color: rgb(128, 145, 165);" width="100%" align="left" noshade size="2"><br>
<br>
<br>
</div><div><div></div><div class="h5"><tt>This is based on CLP release 1.11.0.<br>
<br>
If I dump to MPS file, CLP claims the problem is infeasible - which seems correct (cplex agrees). However, if it is run through OsiClp using resolve, it returns that the problem is solved optimally. The log is below. Can you tell anything from the log? If you want to debug this directly - I can try to send instructions on how to repeat the issue - it happens when using my DIP project.<br>
<br>
Thanks,<br>
Matt<br>
<br>
<br>
<br>
Running CLP:<br>
<br>
[magala@orclus71 build-O]$ ./bin/clp Dip/examples/MILPBlock/PRICE_AND_CUT_master Prob.n0.c6.p14.lp<br>
Coin LP version 1.11.0, build Jan 26 2010<br>
command line - ./bin/clp Dip/examples/MILPBlock/PRICE_AND_CUT_masterProb.n0.c6.p 14.lp<br>
Presolve 119 (-17416) rows, 76 (-17538) columns and 1511 (-20448) elements<br>
0 Obj -1.0275e+07 Primal inf 160.178 (48) Dual inf 5.21844e+07 (71)<br>
77 Obj -1.70728e+07 Primal inf 339.779 (36)<br>
85 Obj -6.75452e+06 Primal inf 3729.11 (19)<br>
87 Obj 2.00406e+07 Primal inf 11780.8 (21)<br>
90 Obj 2.64131e+07 Primal inf 13955 (22)<br>
91 Obj 3.32559e+07 Primal inf 18130.8 (31)<br>
92 Obj 8.34552e+07 Primal inf 43072.8 (31)<br>
94 Obj 4.23954e+08 Primal inf 180013 (41)<br>
95 Obj 9.20098e+08 Primal inf 385843 (39)<br>
96 Obj 9.48195e+08 Primal inf 402516 (43)<br>
97 Obj 8.08035e+09 Primal inf 3.45018e+06 (48)<br>
98 Obj 8.08035e+09 Primal inf 3.46527e+06 (50)<br>
99 Obj 1.33977e+10 Primal inf 5.74173e+06 (51)<br>
99 Obj 1.33977e+10 Primal inf 5.74173e+06 (51)<br>
99 Obj 1.33977e+10 Primal inf 5.74173e+06 (51)<br>
99 Obj 1.33977e+10 Primal inf 5.74173e+06 (51)<br>
99 Obj 1.33977e+10 Primal inf 5.74173e+06 (51)<br>
99 Obj 1.33977e+10 Primal inf 5.74173e+06 (51)<br>
99 Obj 1.33977e+10 Primal inf 5.74173e+06 (51)<br>
99 Obj 1.33977e+10 Primal inf 5.74173e+06 (51)<br>
Primal infeasible - objective value 1.33977e+10<br>
Presolved problem not optimal, resolve after postsolve<br>
After Postsolve, objective 1.33977e+10, infeasibilities - dual -2.5429e+12 (50), primal 1.236e+08 (319)<br>
PrimalInfeasible objective 1.339766625e+10 - 198 iterations time 0.032, Presolve 0.02<br>
<br>
<br>
<br>
<br>
Running OsiClp:<br>
<br>
Coin0513I 17416 rows dropped after presolve pass 1<br>
Coin0513I 17416 rows dropped after presolve pass 2<br>
Coin0506I Presolve 119 (-17416) rows, 76 (-17538) columns and 1511 (-20459) elements<br>
Clp1001I Initial range of elements is 0.0600375 to 11508<br>
Clp1002I Range of elements is 0.0654858 to 15.2705<br>
Clp1002I Range of elements is 0.108171 to 9.24462<br>
Clp1003I Final range of elements is 0.112844 to 8.86181<br>
Clp0022I Absolute values of scaled rhs range from 0.0927259 to 22.6474, minimum gap 1e+100<br>
Clp0020I Absolute values of scaled objective range from 10250.2 to 1.61378e+07<br>
Clp0021I Absolute values of scaled bounds range from 0.0710216 to 1.80219, minimum gap 0.0710216<br>
Clp0018I 1 total structurals rejected in initial factorization<br>
Clp0006I 0 Obj -1.69954e+07 Primal inf 875.907 (53) Dual inf 2.62043e+06 (24)<br>
Clp0006I 102 Obj -7.0894e+06 Primal inf 3599.86 (21)<br>
Clp0006I 104 Obj -6.90525e+06 Primal inf 16284.6 (37)<br>
Clp0010I Flagging variable C33<br>
Clp0010I Flagging variable C61<br>
Clp0010I Flagging variable R1<br>
Clp0010I Flagging variable C34<br>
Clp0010I Flagging variable R53<br>
Clp0006I 106 Obj 8.66818e+06 Primal inf 40216.1 (34)<br>
Clp0010I Flagging variable R114<br>
Clp0010I Flagging variable R78<br>
Clp0006I 107 Obj 8.66818e+06 Primal inf 40156.6 (34)<br>
Clp0010I Flagging variable C48<br>
Clp0010I Flagging variable R102<br>
Clp0010I Flagging variable R38<br>
Clp0006I 115 Obj 2.93615e+07 Primal inf 74000.2 (44)<br>
Clp0010I Flagging variable C45<br>
Clp0010I Flagging variable C62<br>
Clp0006I 117 Obj 7.36568e+07 Primal inf 142995 (41)<br>
Clp0010I Flagging variable R50<br>
Clp0006I 119 Obj 1.7023e+08 Primal inf 293559 (38)<br>
Clp0010I Flagging variable C21<br>
Clp0006I 124 Obj 3.50591e+08 Primal inf 573525 (37)<br>
...<br>
Clp0010I Flagging variable C17470<br>
maybe forcing re-factorization - sum 1.31472e+09 1 pivots<br>
maybe forcing re-factorization - sum 1.43916e+09 2 pivots<br>
Clp0060I Primal error 1.07374e+09, dual error 3.50794e+06<br>
Clp0006I 145 Obj 1.23847e+20 Primal inf 1.98872e+24 (188) Dual inf 127.725 (1)<br>
Clp0060I Primal error 1.07374e+09, dual error 3.50794e+06<br>
Clp0060I Primal error 4.29497e+09, dual error 9.66162e+06<br>
Clp0006I 148 Obj 1.29644e+20 Primal inf 2.08182e+24 (162) Dual inf 0.0783458 (1)<br>
Clp0060I Primal error 1.07374e+09, dual error 9.66162e+06<br>
maybe forcing re-factorization - sum 6.87571e+13 1 pivots<br>
Clp0060I Primal error 16384, dual error 1.06329e+07<br>
Clp0006I 150 Obj 1.37555e+20 Primal inf 2.20884e+24 (106) Dual inf 3360 (1)<br>
Clp0060I Primal error 2.14748e+09, dual error 1.06329e+07<br>
Clp0060I Primal error 2.14748e+09, dual error 4.38545e+06<br>
...<br>
Clp0010I Flagging variable R16579<br>
Clp0006I 285 Obj 1.98074e+21 Primal inf 3.18064e+25 (196)<br>
Solution update n_cols:17614 n_rows: 17535 n_iter: 419 time: 0.871<br>
clp status = 0<br>
clp prob status = 0<br>
clp second status = 0<br>
Iteration Count : 419<br>
isAbandoned() : 0<br>
isProvenOptimal() : 1<br>
isProvenPrimalInfeasible() : 0<br>
isProvenDualInfeasible() : 0<br>
isPrimalObjectiveLimitReached : 0<br>
isDualObjectiveLimitReached : 0<br>
isIterationLimitReached : 0<br>
<br>
<br>
<br>
<br>
_______________________________________________<br>
Clp mailing list<br>
<a href="mailto:Clp@list.coin-or.org" target="_blank">Clp@list.coin-or.org</a><br>
</tt><tt><a href="http://list.coin-or.org/mailman/listinfo/clp" target="_blank">http://list.coin-or.org/mailman/listinfo/clp</a></tt><tt><br>
</tt><br>
<br>
</div></div></div>
<br>_______________________________________________<br>
Clp mailing list<br>
<a href="mailto:Clp@list.coin-or.org">Clp@list.coin-or.org</a><br>
<a href="http://list.coin-or.org/mailman/listinfo/clp" target="_blank">http://list.coin-or.org/mailman/listinfo/clp</a><br>
<br></blockquote></div><br>