[Clp] status wrong using OsiClp

Wed Jan 27 09:26:04 EST 2010

Sorry about all the steps...

1.  svn checkout https://projects.coin-or.org/svn/Dip/trunk coin-Dip

2.  cd coin-Dip

3.  mkdir build-g

4.  cd build-g

5.  ../configure --enable-debug

6.  make

7.  make install

8.  cd Dip/examples/MILPBlock/

9.  make

10. wget http://coral.ie.lehigh.edu/~magh/tmp/milpblockR.parm

11. wget http://coral.ie.lehigh.edu/~magh/tmp/forestry_Random_40.mps

12. wget http://coral.ie.lehigh.edu/~magh/tmp/forestry_Random_40_blocks_stand.txt

13. ./decomp_milpblock --param milpblockR.parm

There is a lot of output. You can reduce that by changing LogDebugLevel in milpblockR.parm. It took about 70 seconds to get to issue. It will throw an exception and the last LP it solves is the issue.

The linux machine I am using has:
[magala at orclus71 MILPBlock]$ uname -a
Linux orclus71 2.6.18-8.el5 #1 SMP Fri Jan 26 14:15:14 EST 2007 x86_64 x86_64 x86_64 GNU/Linux
[magala at orclus71 MILPBlock]$ g++ -v
Using built-in specs.
Target: x86_64-redhat-linux
Configured with: ../configure --prefix=/usr --mandir=/usr/share/man --infodir=/usr/share/info --enable-shared --enable-threads=posix --enable-checking=release --with-system-zlib --enable-__cxa_atexit --disable-libunwind-exceptions --enable-libgcj-multifile --enable-languages=c,c++,objc,obj-c++,java,fortran,ada --enable-java-awt=gtk --disable-dssi --enable-plugin --with-java-home=/usr/lib/jvm/java-1.4.2-gcj-1.4.2.0/jre --with-cpu=generic --host=x86_64-redhat-linux
Thread model: posix
gcc version 4.1.1 20070105 (Red Hat 4.1.1-52)

Note - I could NOT repeat the issue on a linux machine with:
[magala at ordlnx01 MILPBlock]$ uname -a
Linux ordlnx01.unx.sas.com 2.6.9-11.ELsmp #1 SMP Fri May 20 18:26:27 EDT 2005 i686 i686 i386 GNU/Linux
[magala at ordlnx01 MILPBlock]$ g++ -v
Reading specs from /usr/lib/gcc/i386-redhat-linux/3.4.3/specs
Configured with: ../configure --prefix=/usr --mandir=/usr/share/man --infodir=/usr/share/info --enable-shared --enable-threads=posix --disable-checking --with-system-zlib --enable-__cxa_atexit --disable-libunwind-exceptions --enable-java-awt=gtk --host=i386-redhat-linux
Thread model: posix
gcc version 3.4.3 20050227 (Red Hat 3.4.3-22.1)

From: clp-bounces at list.coin-or.org [mailto:clp-bounces at list.coin-or.org] On Behalf Of John J Forrest
Sent: Wednesday, January 27, 2010 4:27 AM
Cc: clp at list.coin-or.org
Subject: Re: [Clp] status wrong using OsiClp

Matt,

Looking at output, OsiClp should have said problem was infeasible. After a quick look at code, I can't see anything obvious.

So you will have to give me something so I can reproduce error.

John

[cid:image001.gif at 01CA9F2C.85005790]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

From:

Matthew Galati <Matthew.Galati at sas.com>

To:

"clp at list.coin-or.org" <clp at list.coin-or.org>

Date:

01/26/2010 10:12 PM

Subject:

[Clp] status wrong using OsiClp

Sent by:

clp-bounces at list.coin-or.org

________________________________

This is based on CLP release 1.11.0.

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.

Thanks,
Matt

Running CLP:

[magala at orclus71 build-O]$ ./bin/clp Dip/examples/MILPBlock/PRICE_AND_CUT_master                                                               Prob.n0.c6.p14.lp
Coin LP version 1.11.0, build Jan 26 2010
command line - ./bin/clp Dip/examples/MILPBlock/PRICE_AND_CUT_masterProb.n0.c6.p                                                               14.lp
Presolve 119 (-17416) rows, 76 (-17538) columns and 1511 (-20448) elements
0  Obj -1.0275e+07 Primal inf 160.178 (48) Dual inf 5.21844e+07 (71)
77  Obj -1.70728e+07 Primal inf 339.779 (36)
85  Obj -6.75452e+06 Primal inf 3729.11 (19)
87  Obj 2.00406e+07 Primal inf 11780.8 (21)
90  Obj 2.64131e+07 Primal inf 13955 (22)
91  Obj 3.32559e+07 Primal inf 18130.8 (31)
92  Obj 8.34552e+07 Primal inf 43072.8 (31)
94  Obj 4.23954e+08 Primal inf 180013 (41)
95  Obj 9.20098e+08 Primal inf 385843 (39)
96  Obj 9.48195e+08 Primal inf 402516 (43)
97  Obj 8.08035e+09 Primal inf 3.45018e+06 (48)
98  Obj 8.08035e+09 Primal inf 3.46527e+06 (50)
99  Obj 1.33977e+10 Primal inf 5.74173e+06 (51)
99  Obj 1.33977e+10 Primal inf 5.74173e+06 (51)
99  Obj 1.33977e+10 Primal inf 5.74173e+06 (51)
99  Obj 1.33977e+10 Primal inf 5.74173e+06 (51)
99  Obj 1.33977e+10 Primal inf 5.74173e+06 (51)
99  Obj 1.33977e+10 Primal inf 5.74173e+06 (51)
99  Obj 1.33977e+10 Primal inf 5.74173e+06 (51)
99  Obj 1.33977e+10 Primal inf 5.74173e+06 (51)
Primal infeasible - objective value 1.33977e+10
Presolved problem not optimal, resolve after postsolve
After Postsolve, objective 1.33977e+10, infeasibilities - dual -2.5429e+12 (50),                                                                primal 1.236e+08 (319)
PrimalInfeasible objective 1.339766625e+10 - 198 iterations time 0.032, Presolve                                                                0.02

Running OsiClp:

Coin0513I 17416 rows dropped after presolve pass 1
Coin0513I 17416 rows dropped after presolve pass 2
Coin0506I Presolve 119 (-17416) rows, 76 (-17538) columns and 1511 (-20459) elements
Clp1001I Initial range of elements is 0.0600375 to 11508
Clp1002I Range of elements is 0.0654858 to 15.2705
Clp1002I Range of elements is 0.108171 to 9.24462
Clp1003I Final range of elements is 0.112844 to 8.86181
Clp0022I Absolute values of scaled rhs range from 0.0927259 to 22.6474, minimum gap 1e+100
Clp0020I Absolute values of scaled objective range from 10250.2 to 1.61378e+07
Clp0021I Absolute values of scaled bounds range from 0.0710216 to 1.80219, minimum gap 0.0710216
Clp0018I 1 total structurals rejected in initial factorization
Clp0006I 0  Obj -1.69954e+07 Primal inf 875.907 (53) Dual inf 2.62043e+06 (24)
Clp0006I 102  Obj -7.0894e+06 Primal inf 3599.86 (21)
Clp0006I 104  Obj -6.90525e+06 Primal inf 16284.6 (37)
Clp0010I Flagging variable C33
Clp0010I Flagging variable C61
Clp0010I Flagging variable R1
Clp0010I Flagging variable C34
Clp0010I Flagging variable R53
Clp0006I 106  Obj 8.66818e+06 Primal inf 40216.1 (34)
Clp0010I Flagging variable R114
Clp0010I Flagging variable R78
Clp0006I 107  Obj 8.66818e+06 Primal inf 40156.6 (34)
Clp0010I Flagging variable C48
Clp0010I Flagging variable R102
Clp0010I Flagging variable R38
Clp0006I 115  Obj 2.93615e+07 Primal inf 74000.2 (44)
Clp0010I Flagging variable C45
Clp0010I Flagging variable C62
Clp0006I 117  Obj 7.36568e+07 Primal inf 142995 (41)
Clp0010I Flagging variable R50
Clp0006I 119  Obj 1.7023e+08 Primal inf 293559 (38)
Clp0010I Flagging variable C21
Clp0006I 124  Obj 3.50591e+08 Primal inf 573525 (37)
...
Clp0010I Flagging variable C17470
maybe forcing re-factorization - sum 1.31472e+09  1 pivots
maybe forcing re-factorization - sum 1.43916e+09  2 pivots
Clp0060I Primal error 1.07374e+09, dual error 3.50794e+06
Clp0006I 145  Obj 1.23847e+20 Primal inf 1.98872e+24 (188) Dual inf 127.725 (1)
Clp0060I Primal error 1.07374e+09, dual error 3.50794e+06
Clp0060I Primal error 4.29497e+09, dual error 9.66162e+06
Clp0006I 148  Obj 1.29644e+20 Primal inf 2.08182e+24 (162) Dual inf 0.0783458 (1)
Clp0060I Primal error 1.07374e+09, dual error 9.66162e+06
maybe forcing re-factorization - sum 6.87571e+13  1 pivots
Clp0060I Primal error 16384, dual error 1.06329e+07
Clp0006I 150  Obj 1.37555e+20 Primal inf 2.20884e+24 (106) Dual inf 3360 (1)
Clp0060I Primal error 2.14748e+09, dual error 1.06329e+07
Clp0060I Primal error 2.14748e+09, dual error 4.38545e+06
...
Clp0010I Flagging variable R16579
Clp0006I 285  Obj 1.98074e+21 Primal inf 3.18064e+25 (196)
Solution update n_cols:17614      n_rows: 17535      n_iter: 419        time: 0.871
clp status        = 0
clp prob status   = 0
clp second status = 0
Iteration Count               : 419
isAbandoned()                 : 0
isProvenOptimal()             : 1
isProvenPrimalInfeasible()    : 0
isProvenDualInfeasible()      : 0
isPrimalObjectiveLimitReached : 0
isDualObjectiveLimitReached   : 0
isIterationLimitReached       : 0

_______________________________________________
Clp mailing list
Clp at list.coin-or.org
http://list.coin-or.org/mailman/listinfo/clp

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://list.coin-or.org/pipermail/clp/attachments/20100127/5378a03e/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.gif
Type: image/gif
Size: 105 bytes
Desc: image001.gif
URL: <http://list.coin-or.org/pipermail/clp/attachments/20100127/5378a03e/attachment.gif>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image003.png
Type: image/png
Size: 168 bytes
Desc: image003.png
URL: <http://list.coin-or.org/pipermail/clp/attachments/20100127/5378a03e/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image004.png
Type: image/png
Size: 166 bytes
Desc: image004.png
URL: <http://list.coin-or.org/pipermail/clp/attachments/20100127/5378a03e/attachment-0001.png>