The "issue" with this model when using the BCP method is that the variable x1 is free. Because we do not support extreme rays, we have to make sure that each "block" is bounded. Even though the variable x1 is not part of any of the defined blocks, internally, it is treated as a special block called the 'master-only' block. Like the other blocks, we are trying to generate extreme points of the polytope defined by that block. For your master-only block, x1 is free and for certain cost vectors can be unbounded.<br>
<br>The easiest way to ensure each polytope is bounded is to bound all your variables. You can do this manually in your MPS file, or - I have added an option to the MILPBlock application which can be used to bound all variables that have no bounds.<br>
<br>Do an svn update (at ~/coin-Dip/Dip) to get the changes. Then in the parameter file, add a ColumnLB and ColumnUB as follows:<br><br>#############################################################################<br>Settings for MILPBlock application.<br>
section=MILPBlock<br>#############################################################################<br>[MILPBlock]<br># Path to instance data.<br>DataDir = .<br><br># {Pair,List}: Block file format. See README.txt.<br>
BlockFileFormat = List<br><br># Name of block file (assumed to be in DataDir).<br>BlockFile = test_4_1.block<br><br># Name of MPS file for instance (assumed to be in DataDir).<br>Instance = test_4_1<br><br>ColumnLB = -1.0e4<br>
ColumnUB = 1.0e4<br><br><br>This will add the bounds -1000 <= x1 <= 1000 which should fix the issue related to unboundeness. And you should then be able to solve the problem using the Price and Cut algorithm.<br>
<br>Let me know if you have any issues.<br><br>In a future release, we will add extreme rays - so that the user does not have to worry about this.<br><br>Thanks,<br>Matt<br><br><br><br><br><br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Ayse,<br>
<br>
Replying to the DIP mailing list since this might be relevant to other users - although I think you are the only one so far. :)<br>
<br>
<br>
I am assuming the following happens when solving the subproblem - i.e., one of your blocks.<br>
"Coin0508I Presolve thinks problem is unbounded<br>
Error: CBC IP solver 2nd status = 7"<br>
<br>
This is because the BCP method in DIP does not consider extreme rays. Rather it assumes that the feasible region is bounded in all directions. I have yet to see any real-world example where extreme rays are needed. It should be easy to bound your feasible region by giving lower and/or upper bounds to all your variables - even if they are artificially big. Try that and see if that helps. In some future release of DIP, I will add extreme rays so this can be dealt with more naturally.<br>
<br>
<br>
<br>
"I am kind of confused why one method does not work and the others do and also I am wondering how to decide which method to use, is there any documentation that I can read about this?"<br>
<br>
There is no real documentation to help a user decide this. In fact, that is part of why we wrote DIP - so it is easy to try different methods without having to recode anything.<br>
<br>
See<br>
<a href="https://coral.ie.lehigh.edu/%7Emagh/decomp/decomp_INFORMS09.pdf" target="_blank">https://coral.ie.lehigh.edu/~magh/decomp/decomp_INFORMS09.pdf</a> and<br>
<a href="http://coral.ie.lehigh.edu/%7Eted/files/papers/DECOMP04.pdf" target="_blank">http://coral.ie.lehigh.edu/~ted/files/papers/DECOMP04.pdf</a><br>
for some background information on the different algorithms. There will be much more detail in my thesis - which is due out in one month. I will post that when it is complete.<br>
<br>
Thanks,<br>
Matt<br>
<br>
<br>
<br>
From: Ayse Durukan [mailto:<a href="mailto:adurukan@mail.uh.edu">adurukan@mail.uh.edu</a>]<br>
Sent: Thursday, October 29, 2009 2:08 PM<br>
To: Matthew Galati<br>
Subject: Re: RE: good news<br>
<br>
Hi Mr. Galati,<br>
I tried converting the mps files using cplex and it worked. Also manually removing the lines before NAME header helps too.<br>
After converting the files, I kept getting the same error until I start playing with the main settings. When I solve it with DIP's CPM method it works fine.<br>
When I solve it with DIP's BCP method it gives following error<br>
Coin0508I Presolve thinks problem is unbounded<br>
Error: CBC IP solver 2nd status = 7<br>
When I set doCut=1 and doDirect=1 it works too.<br>
I am kind of confused why one method does not work and the others do and also I am wondering how to decide which method to use, is there any documentation that I can read about this?<br>
Also, please see the attachment for the problem's mps files and block files.<br>
<br>
Thank You<br>
Ayse<br>
<br>
<br>
<br>
<br>
<br>
_______________________________________________<br>
Dip mailing list<br>
<a href="mailto:Dip@list.coin-or.org">Dip@list.coin-or.org</a><br>
<a href="http://list.coin-or.org/mailman/listinfo/dip" target="_blank">http://list.coin-or.org/mailman/listinfo/dip</a><br>
</blockquote></div><br>