<html>
<head>
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">James,<br>
<br>
In standalone Cbc bounds are added to SOS variables to make sure
they are in region -1.0e15 <= 1.0e15. This cures the problem
that the relaxed solution is unbounded.<br>
<br>
Adding<br>
<br>
HS_pos <= 1.0e15<br>
HS_P0_diff <= 1.0e15<br>
<br>
to the input gives correct result.<br>
<br>
I have added readLp. Not sure what correct course of action is
for unbounded case.<br>
<br>
John Forrest<br>
<br>
<br>
On 14/02/18 01:53, James Howey wrote:<br>
</div>
<blockquote type="cite"
cite="mid:000801d3a536$a9c9fb80$fd5df280$@gmail.com">
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
<meta name="Generator" content="Microsoft Word 15 (filtered
medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div class="WordSection1">
<p class="MsoNormal">After applying recent SOS fixes, the below
model is solved correctly with an objective value of zero.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">However, I am trying to use the cbc
interface dll. As it has no readLp, I added one.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:#6F008A">COINLIBAPI</span><span
style="font-size:9.5pt;font-family:Consolas;color:black"> </span><span
style="font-size:9.5pt;font-family:Consolas;color:blue">int</span><span
style="font-size:9.5pt;font-family:Consolas;color:black"> </span><span
style="font-size:9.5pt;font-family:Consolas;color:#6F008A">COINLINKAGE</span><span
style="font-size:9.5pt;font-family:Consolas;color:black"><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">Cbc_readLp(</span><span
style="font-size:9.5pt;font-family:Consolas;color:#2B91AF">Cbc_Model</span><span
style="font-size:9.5pt;font-family:Consolas;color:black"> *
</span><span
style="font-size:9.5pt;font-family:Consolas;color:gray">model</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">, </span><span
style="font-size:9.5pt;font-family:Consolas;color:blue">const</span><span
style="font-size:9.5pt;font-family:Consolas;color:black"> </span><span
style="font-size:9.5pt;font-family:Consolas;color:blue">char</span><span
style="font-size:9.5pt;font-family:Consolas;color:black"> *</span><span
style="font-size:9.5pt;font-family:Consolas;color:gray">filename</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">)<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">{<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">
</span><span
style="font-size:9.5pt;font-family:Consolas;color:blue">const</span><span
style="font-size:9.5pt;font-family:Consolas;color:black"> </span><span
style="font-size:9.5pt;font-family:Consolas;color:blue">char</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">
prefix[] = </span><span
style="font-size:9.5pt;font-family:Consolas;color:#A31515">"Cbc_C_Interface::Cbc_readLp():
"</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">;<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">
</span><span
style="font-size:9.5pt;font-family:Consolas;color:green">//
const int VERBOSE = 2;</span><span
style="font-size:9.5pt;font-family:Consolas;color:black"><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">
</span><span
style="font-size:9.5pt;font-family:Consolas;color:blue">if</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">
(VERBOSE > 0) printf(</span><span
style="font-size:9.5pt;font-family:Consolas;color:#A31515">"%s
begin\n"</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">,
prefix);<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">
</span><span
style="font-size:9.5pt;font-family:Consolas;color:blue">if</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">
(VERBOSE > 1) printf(</span><span
style="font-size:9.5pt;font-family:Consolas;color:#A31515">"%s
filename = '%s'\n"</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">,
prefix, </span><span
style="font-size:9.5pt;font-family:Consolas;color:gray">filename</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">);<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black"><o:p> </o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">
</span><span
style="font-size:9.5pt;font-family:Consolas;color:blue">auto</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">
clpSolver = </span><span
style="font-size:9.5pt;font-family:Consolas;color:blue">dynamic_cast</span><span
style="font-size:9.5pt;font-family:Consolas;color:black"><
</span><span
style="font-size:9.5pt;font-family:Consolas;color:#2B91AF">OsiClpSolverInterface</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">*>
(</span><span
style="font-size:9.5pt;font-family:Consolas;color:gray">model</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">->model_->solver());<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">
</span><span
style="font-size:9.5pt;font-family:Consolas;color:blue">int</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">
result = 1;<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">
</span><span
style="font-size:9.5pt;font-family:Consolas;color:green">//result
= model->model_->solver()->readLp(filename);</span><span
style="font-size:9.5pt;font-family:Consolas;color:black"><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">
result = clpSolver->readLp(</span><span
style="font-size:9.5pt;font-family:Consolas;color:gray">filename</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">);<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">
</span><span
style="font-size:9.5pt;font-family:Consolas;color:#6F008A">assert</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">(result
== 0);<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black"><o:p> </o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">
</span><span
style="font-size:9.5pt;font-family:Consolas;color:blue">if</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">
(VERBOSE > 0) printf(</span><span
style="font-size:9.5pt;font-family:Consolas;color:#A31515">"%s
return %i\n"</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">,
prefix, result);<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">
</span><span
style="font-size:9.5pt;font-family:Consolas;color:blue">return</span><span
style="font-size:9.5pt;font-family:Consolas;color:black">
result;<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black">}<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span
style="font-size:9.5pt;font-family:Consolas;color:black"><o:p> </o:p></span></p>
<p class="MsoNormal">When I route the model thru the interface,
It appears to lose the SOS sets and gives me an unbounded
result.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I would value any insights or better, a
workaround/fix to enable my scenario. Is it CbcModel vs.
Cbc_Model?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks for any assistance,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">jkh<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">MAXIMIZE<o:p></o:p></p>
<p class="MsoNormal">+HS_pos <o:p></o:p></p>
<p class="MsoNormal">SUBJECT TO<o:p></o:p></p>
<p class="MsoNormal">HS_free: -HS_@free +HS_free = 0<o:p></o:p></p>
<p class="MsoNormal">HS_P0_diff: +HS_free +HS_P0_diff -HS_pos =
0<o:p></o:p></p>
<p class="MsoNormal">HS_P0_sum: +HS_P0_y +HS_P0_z = 1<o:p></o:p></p>
<p class="MsoNormal">HS_pos: -HS_@pos +HS_pos = 0<o:p></o:p></p>
<p class="MsoNormal">BOUNDS<o:p></o:p></p>
<p class="MsoNormal">HS_@free free<o:p></o:p></p>
<p class="MsoNormal">HS_@pos free<o:p></o:p></p>
<p class="MsoNormal">HS_free = -50000<o:p></o:p></p>
<p class="MsoNormal">BINARY<o:p></o:p></p>
<p class="MsoNormal">HS_P0_y HS_P0_z <o:p></o:p></p>
<p class="MsoNormal">SOS<o:p></o:p></p>
<p class="MsoNormal">S0: S1:: HS_pos:1 HS_P0_y:2 <o:p></o:p></p>
<p class="MsoNormal">S1: S1:: HS_P0_diff:1 HS_P0_z:2 <o:p></o:p></p>
<p class="MsoNormal">END<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Cbc mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Cbc@list.coin-or.org">Cbc@list.coin-or.org</a>
<a class="moz-txt-link-freetext" href="https://list.coin-or.org/mailman/listinfo/cbc">https://list.coin-or.org/mailman/listinfo/cbc</a>
</pre>
</blockquote>
<p><br>
</p>
</body>
</html>