<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<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:"Malgun Gothic";
        panose-1:2 11 5 3 2 0 0 2 0 4;}
@font-face
        {font-family:"\@Malgun Gothic";}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
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;}
/* List Definitions */
@list l0
        {mso-list-id:738552756;
        mso-list-type:hybrid;
        mso-list-template-ids:-250719748 67698705 67698713 67698715 67698703 67698713 67698715 67698703 67698713 67698715;}
@list l0:level1
        {mso-level-text:"%1\)";
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
@list l0:level2
        {mso-level-number-format:alpha-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
@list l0:level3
        {mso-level-number-format:roman-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:right;
        text-indent:-9.0pt;}
@list l0:level4
        {mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
@list l0:level5
        {mso-level-number-format:alpha-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
@list l0:level6
        {mso-level-number-format:roman-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:right;
        text-indent:-9.0pt;}
@list l0:level7
        {mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
@list l0:level8
        {mso-level-number-format:alpha-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
@list l0:level9
        {mso-level-number-format:roman-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:right;
        text-indent:-9.0pt;}
@list l1
        {mso-list-id:1702709724;
        mso-list-type:hybrid;
        mso-list-template-ids:1867269810 67698705 67698713 67698715 67698703 67698713 67698715 67698703 67698713 67698715;}
@list l1:level1
        {mso-level-text:"%1\)";
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
@list l1:level2
        {mso-level-number-format:alpha-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
@list l1:level3
        {mso-level-number-format:roman-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:right;
        text-indent:-9.0pt;}
@list l1:level4
        {mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
@list l1:level5
        {mso-level-number-format:alpha-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
@list l1:level6
        {mso-level-number-format:roman-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:right;
        text-indent:-9.0pt;}
@list l1:level7
        {mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
@list l1:level8
        {mso-level-number-format:alpha-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
@list l1:level9
        {mso-level-number-format:roman-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:right;
        text-indent:-9.0pt;}
ol
        {margin-bottom:0in;}
ul
        {margin-bottom:0in;}
--></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]-->
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Hi ,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I have two questions, which are irrelevant to each other.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">1)   Can I set the Primal bound before solving Branch and Price algorithm?
<o:p></o:p></p>
<p class="MsoNormal">I have an initial set of columns obtained from heuristics which are added to the root node and know a good IP solution using the existing columns.<o:p></o:p></p>
<p class="MsoNormal">I would like to use that IP solution as the primal bound before it starts so that some nodes with RMP LP obj exceeding the primal bound (known IP solution) can be pruned earlier.
<o:p></o:p></p>
<p class="MsoNormal">To do this, should I create a solution using SCIPcreatesol or is there a simpler way to apply this ? ( For example, there is a lower bound pointer to store a valid lower bound. Can I do similar setting for the upper bound?)   <o:p></o:p></p>
<p class="MsoNormal"> <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">2)   I have built a MILP problem using Python-SCIP interface and wrote it out to a .lp file.<o:p></o:p></p>
<p class="MsoNormal">Then, I found that the prosolved problems are different for the same instance when I solve it with the model code (written in python) VS solve it in a SCIP interactive shell by reading the .lp file.<o:p></o:p></p>
<p class="MsoNormal">This different presolved problem results in performance difference. A model written and solved within the code has taken twice longer time to optimality than the one solved within at SCIP interactive shell.   <o:p></o:p></p>
<p class="MsoNormal">They are equivalent instances and I have check that all the parameters equally set to default values in both cases.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">What could make this difference in the presolving part? Is this a problem of using Python interface ?
<o:p></o:p></p>
<p class="MsoNormal">Is it possible that a model written in the code has more information that is not included in .lp file so that the presolved problem is different?  <o:p></o:p></p>
<p class="MsoNormal">(I don’t think this is the case since reading .lp file and solving it in a SCIP interactive shell was faster)
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I have attached the result log of the initial presolving part below, if it helps to understand my question.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">====Build a model and solve ======================================================================================================<o:p></o:p></p>
<p class="MsoNormal">presolving:<o:p></o:p></p>
<p class="MsoNormal">(round 1, fast)       9157 del vars, 9157 del conss, 0 add conss, 9157 chg bounds, 0 chg sides, 0 chg coeffs, 0 upgd conss, 0 impls, 27 clqs<o:p></o:p></p>
<p class="MsoNormal">(round 2, fast)       9157 del vars, 17473 del conss, 0 add conss, 9157 chg bounds, 0 chg sides, 0 chg coeffs, 0 upgd conss, 0 impls, 27 clqs<o:p></o:p></p>
<p class="MsoNormal">(round 3, exhaustive) 9157 del vars, 17473 del conss, 0 add conss, 9157 chg bounds, 0 chg sides, 0 chg coeffs, 6112 upgd conss, 0 impls, 27 clqs<o:p></o:p></p>
<p class="MsoNormal">   (1.0s) probing: 69/14877 (0.5%) - 0 fixings, 0 aggregations, 943 implications, 0 bound changes<o:p></o:p></p>
<p class="MsoNormal">   (1.0s) probing aborted: 50/50 successive totally useless probings<o:p></o:p></p>
<p class="MsoNormal">   Deactivated symmetry handling methods, since SCIP was built without symmetry detector (SYM=none).<o:p></o:p></p>
<p class="MsoNormal">presolving (4 rounds: 4 fast, 2 medium, 2 exhaustive):<o:p></o:p></p>
<p class="MsoNormal">9157 deleted vars, 17473 deleted constraints, 0 added constraints, 9157 tightened bounds, 0 added holes, 0 changed sides, 0 changed coefficients<o:p></o:p></p>
<p class="MsoNormal">7009 implications, 45 cliques<o:p></o:p></p>
<p class="MsoNormal"><span style="background:yellow;mso-highlight:yellow">presolved problem has 21045 variables (14877 bin, 0 int, 0 impl, 6168 cont) and 35953 constraints<o:p></o:p></span></p>
<p class="MsoNormal"><span style="background:yellow;mso-highlight:yellow">   6084 constraints of type <varbound><o:p></o:p></span></p>
<p class="MsoNormal"><span style="background:yellow;mso-highlight:yellow">      1 constraints of type <knapsack><o:p></o:p></span></p>
<p class="MsoNormal"><span style="background:yellow;mso-highlight:yellow">     27 constraints of type <setppc><o:p></o:p></span></p>
<p class="MsoNormal"><span style="background:yellow;mso-highlight:yellow">  29841 constraints of type <linear></span><o:p></o:p></p>
<p class="MsoNormal">Presolving Time: 1.00<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">====Read a model file (.lp) and solve ================================================================================================<o:p></o:p></p>
<p class="MsoNormal">SCIP> read "SCIP_MIP.lp"<o:p></o:p></p>
<p class="MsoNormal">original problem has 15802 variables (502 bin, 0 int, 0 impl, 15300 cont) and 39010 constraints<o:p></o:p></p>
<p class="MsoNormal">SCIP> optimize<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">presolving:<o:p></o:p></p>
<p class="MsoNormal">(round 1, fast)       9141 del vars, 9141 del conss, 0 add conss, 9141 chg bounds, 0 chg sides, 0 chg coeffs, 0 upgd conss, 0 impls, 43 clqs<o:p></o:p></p>
<p class="MsoNormal">(round 2, fast)       9141 del vars, 25773 del conss, 0 add conss, 9141 chg bounds, 0 chg sides, 0 chg coeffs, 0 upgd conss, 0 impls, 43 clqs<o:p></o:p></p>
<p class="MsoNormal">   (0.0s) running MILP presolver<o:p></o:p></p>
<p class="MsoNormal">   (0.0s) MILP presolver (2 rounds): 0 aggregations, 636 fixings, 0 bound changes<o:p></o:p></p>
<p class="MsoNormal">(round 3, medium)     9777 del vars, 25773 del conss, 0 add conss, 9141 chg bounds, 0 chg sides, 0 chg coeffs, 0 upgd conss, 0 impls, 43 clqs<o:p></o:p></p>
<p class="MsoNormal">(round 4, fast)       9777 del vars, 26861 del conss, 0 add conss, 9141 chg bounds, 0 chg sides, 0 chg coeffs, 0 upgd conss, 0 impls, 43 clqs<o:p></o:p></p>
<p class="MsoNormal">(round 5, exhaustive) 9777 del vars, 26861 del conss, 0 add conss, 9141 chg bounds, 0 chg sides, 0 chg coeffs, 7110 upgd conss, 0 impls, 43 clqs<o:p></o:p></p>
<p class="MsoNormal">   (1.0s) probing cycle finished: starting next cycle<o:p></o:p></p>
<p class="MsoNormal">   Deactivated symmetry handling methods, since SCIP was built without symmetry detector (SYM=none).<o:p></o:p></p>
<p class="MsoNormal">presolving (6 rounds: 6 fast, 3 medium, 2 exhaustive):<o:p></o:p></p>
<p class="MsoNormal">9777 deleted vars, 26861 deleted constraints, 0 added constraints, 9141 tightened bounds, 0 added holes, 0 changed sides, 0 changed coefficients<o:p></o:p></p>
<p class="MsoNormal">110984 implications, 70 cliques<o:p></o:p></p>
<p class="MsoNormal"><span style="background:yellow;mso-highlight:yellow">presolved problem has 6025 variables (401 bin, 0 int, 0 impl, 5624 cont) and 12149 constraints<o:p></o:p></span></p>
<p class="MsoNormal"><span style="background:yellow;mso-highlight:yellow">   7066 constraints of type <varbound><o:p></o:p></span></p>
<p class="MsoNormal"><span style="background:yellow;mso-highlight:yellow">      1 constraints of type <knapsack><o:p></o:p></span></p>
<p class="MsoNormal"><span style="background:yellow;mso-highlight:yellow">     43 constraints of type <setppc><o:p></o:p></span></p>
<p class="MsoNormal"><span style="background:yellow;mso-highlight:yellow">   5039 constraints of type <linear></span><o:p></o:p></p>
<p class="MsoNormal">transformed objective value is always integral (scale: 1)<o:p></o:p></p>
<p class="MsoNormal">Presolving Time: 1.0<o:p></o:p></p>
<p class="MsoNormal">===============================================================================================================================<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#002060">Jungeun Shin, PhD candidate<o:p></o:p></span></b></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#002060">Department of Industrial and Enterprise Systems Engineering<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#002060">College of Engineering | University of Illinois, Urbana-Champaign<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#002060"><a href="mailto:Jungeun4@illinois.edu"><span style="color:#002060">Jungeun4@illinois.edu</span></a> | 1-217-200-6772<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>