
<html>

  <head>

    <meta content="text/html; charset=ISO-8859-1"

      http-equiv="Content-Type">

  </head>

  <body text="#000000" bgcolor="#FFFFFF">

    <br>

    <div class="moz-cite-prefix">On 11/27/2012 12:19 PM, William H.

      Patton wrote:<br>

    </div>

    <blockquote cite="mid:50B5042B.2090602@comcast.net" type="cite">

      <meta content="text/html; charset=ISO-8859-1"

        http-equiv="Content-Type">

      I think your issue is poor scaling of the vars.&nbsp; Try scaling the

      quadratic vars&nbsp; by 1e-3 so that the objective should be around

      -1.3 e+03 instead.<br>

      Then maybe the complementarity will down scale by 1e-3 or so and

      quit within its tolerance.<br>

      <br>

      <br>

      Any sensible person would want to quit around here at iter 18 or

      22. The relative Complementary is 40/ 10^9&nbsp; which is 8 digits of

      precision.<br>

      Pretty good when you should expect only 5 to 6 digits after a

      bunch of&nbsp; double precision scalar products on the <font size="1">911766


        elements in sparse Cholesky,<font size="1">.</font></font><font

        size="1">

        <p>16 Primal -1.3318262e+09 Dual -1.3318287e+09 Complementarity

          2519.259 - 76 fixed, rank 3924</p>

        <p>17 Primal -1.3318263e+09 Dual -1.3318263e+09 Complementarity

          40.523345 - 134 fixed, rank 3924</p>

        <p><b>18 Primal -1.3318263e+09 Dual -1.3318263e+09

            Complementarity 40.006314 - 134 fixed, rank 3924</b></p>

        <p>19 Primal -1.3318263e+09 Dual -1.3318264e+09 Complementarity

          39.542777 - 134 fixed, rank 3924</p>

        <p>20 Primal -1.3318263e+09 Dual -1.3318265e+09 Complementarity

          39.128027 - 134 fixed, rank 3924</p>

        <p>21 Primal -1.3318263e+09 Dual -1.3318266e+09 Complementarity

          38.757635 - 134 fixed, rank 3924</p>

        <p><b>22 Primal -1.3318263e+09 Dual -1.3318263e+09

            Complementarity 0.87384878 - 134 fixed, rank 3924</b></p>

        <p>23 Primal -1.3318263e+09 Dual -1.331649e+09 Complementarity

          0.0061784205 - 134 fixed, rank 3924</p>

        <p>24 Primal -1.3318263e+09 Dual -1.3047041e+09 Complementarity

          4.0066434e-05 - 134 fixed, rank 39<br>

          <br>

        </p>

      </font>Check the code again.&nbsp; <font size="1"><b>Complementarity </b></font><br>

      Searching for: Complementarity <br>

      ClpHelperFunctions.hpp(197): // Form residuals for the

      complementarity equations.<br>

      ClpHelperFunctions.hpp(200): // Cinf&nbsp; is the complementarity

      residual for X1 z1 = mu e, etc.<br>

      ClpMessage.cpp(73): {CLP_BARRIER_ITERATION, 35, 1, "%d Primal %g

      Dual %g Complementarity %g - %d fixed, rank %d"},<br>

      ClpMessage.cpp(76): {CLP_BARRIER_CLOSE_TO_OPTIMAL, 38, 2, "Close

      to optimal after %d iterations with complementarity %g"},<br>

      ClpMessage.cpp(77): {CLP_BARRIER_COMPLEMENTARITY, 39, 2,

      "Complementarity %g - %s"},<br>

      ClpMessage.cpp(84): {CLP_BARRIER_END, 46, 1, "At end primal/dual

      infeasibilities %g/%g, complementarity gap %g, objective %g"},<br>

      ClpPdco.cpp(480): // Initialize mu and complementarity residuals:<br>

      ClpPredictorCorrector.cpp(461): &lt;&lt; "- small complementarity

      gap"<br>

      ClpPredictorCorrector.cpp(468): &lt;&lt; "- objective gap and

      complementarity gap both small"<br>

      ClpPredictorCorrector.cpp(480): &lt;&lt; "- objective gap and

      complementarity gap both smallish and small steps"<br>

      ClpPredictorCorrector.cpp(952): //std::cout&lt;&lt;"Absolute

      complementarity at end

      "&lt;&lt;complementarityGap_&lt;&lt;std::endl;<br>

      ClpPredictorCorrector.cpp(954): //std::cout&lt;&lt;"maximum

      complementarity "&lt;&lt;worstComplementarity_&lt;&lt;std::endl;<br>

      ClpPredictorCorrector.cpp(1125): // minimize complementarity +

      norm*dual inf ? primal inf<br>

      ClpPredictorCorrector.cpp(1127): // Complementarity gap will be

      a*change*change + b*change +c<br>

      ClpPredictorCorrector.cpp(1265): // Complementarity gap will be

      a*change*change + b*change +c<br>

      ClpPredictorCorrector.cpp(2621): printf("good complementarity

      range %g to %g\n", minBeta, maxBeta);<br>

      ClpPredictorCorrector.cpp(3423): CoinWorkDouble complementarity =

      zVec_[iColumn] * lowerSlack_[iColumn];<br>

      ClpPredictorCorrector.cpp(3424): if (complementarity &lt; nextMu)

      {<br>

      ClpPredictorCorrector.cpp(3437): CoinWorkDouble complementarity =

      wVec_[iColumn] * upperSlack_[iColumn];<br>

      ClpPredictorCorrector.cpp(3438): if (complementarity &gt; nextMu)

      {<br>

      ClpPredictorCorrector.cpp(3441): (complementarity - nextMu) /

      upperSlack_[iColumn]);<br>

      ClpPredictorCorrector.cpp(3453): CoinWorkDouble complementarity =

      zVec_[iColumn] * lowerSlack_[iColumn];<br>

      ClpPredictorCorrector.cpp(3454): if (complementarity &gt; nextMu)

      {<br>

      ClpPredictorCorrector.cpp(3467): CoinWorkDouble complementarity =

      wVec_[iColumn] * upperSlack_[iColumn];<br>

      ClpPredictorCorrector.cpp(3468): if (complementarity &lt; nextMu)

      {<br>

      ClpPredictorCorrector.cpp(3471): (complementarity - nextMu) /

      upperSlack_[iColumn]);<br>

      ClpSimplex.hpp(771): /// Large bound value (for complementarity

      etc)<br>

      ClpSimplex.hpp(1346): /// Large bound value (for complementarity

      etc)<br>

      Found 28 occurrence(s) in 5 file(s)<br>

      <br>

      So<br>

      <b>ClpPredictorCorrector.cpp(246)</b>:

      handler_-&gt;message(CLP_BARRIER_ITERATION, messages_)<br>

      context&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <font size="1"><b>18 Primal -1.3318263e+09

          Dual -1.3318263e+09 Complementarity 40.006314 - 134 fixed,

          rank 3924</b></font><br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; complementarityGap_ =

      complementarityGap(numberComplementarityPairs_,<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; numberComplementarityItems_, 0);<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; handler_-&gt;message(CLP_BARRIER_ITERATION, messages_)<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;&lt; numberIterations_<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;&lt;

      static_cast&lt;double&gt;(primalObjective_ * objScale -

      dblParam_[ClpObjOffset])<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;&lt;

      static_cast&lt;double&gt;(dualObjective_ * objScale -

      dblParam_[ClpObjOffset])<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;&lt;

      static_cast&lt;double&gt;(complementarityGap_)<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;&lt; numberFixedTotal<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;&lt; cholesky_-&gt;rank()<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;&lt; CoinMessageEol;<br>

      <br>

      line 336<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (CoinAbs(primalObjective_ - dualObjective_) &lt;

      <b>dualTolerance</b>())

      {&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;


      &lt;&lt;&nbsp; so maybe relax the toleerance<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; gapO = 0.0;<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; }<br>

      <br>

      line

      361&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;


      &lt;&lt;&lt; figure out how to enable sloppyOptimal<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (sloppyOptimal) {<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; // vaguely optimal<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (maximumBoundInfeasibility_ &gt;

      1.0e-2 ||<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; scaledRHSError &gt; 1.0e-2 ||<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; maximumDualError_ &gt;

      objectiveNorm_ * 1.0e-2) {<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;

      handler_-&gt;message(CLP_BARRIER_EXIT2, messages_)<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;&lt; saveIteration<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;&lt; CoinMessageEol;<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; problemStatus_ = 0; // benefit of

      doubt<br>

      &nbsp; <br>

      <br>

      Feel free to continue reading this area for other stopping rules.&nbsp;

      continues to the last break at line 485&nbsp;&nbsp;&nbsp;&nbsp; ( I have source for

      1.14.7&nbsp; not.8).<br>

      Searching for: sloppyOptimal&nbsp;&nbsp; it seems to be locally detected<br>

      ClpPredictorCorrector.cpp(212): bool sloppyOptimal = false;<br>

      ClpPredictorCorrector.cpp(275): if (!sloppyOptimal) {<br>

      ClpPredictorCorrector.cpp(361): if (sloppyOptimal) {<br>

      ClpPredictorCorrector.cpp(403): if ((gapO &lt; 1.0e-6 || (gapO

      &lt; 1.0e-4 &amp;&amp; complementarityGap_ &lt; 0.1)) &amp;&amp;

      !sloppyOptimal) {<br>

      ClpPredictorCorrector.cpp(404): sloppyOptimal <b>= true;</b><br>

      ClpPredictorCorrector.cpp(417): if (saveIteration &gt;= 0

      &amp;&amp; sloppyOptimal) {<br>

      <br>

      So if you keep your current scale for the problem you need to

      tweek the block around line 403&nbsp;&nbsp; perhaps the bold blocke needs

      tweeking by relative ratio of gap to objective<br>

    </blockquote>

    The +1 is to shield small near 0 objective divisors and leave the

    test approximately as the original test with a denominator near 1.<br>

    <blockquote cite="mid:50B5042B.2090602@comcast.net" type="cite"> &nbsp;

      say the bold gap0/(primalObjective_ * objScale -

      dblParam_[ClpObjOffset]&nbsp; +1.0)&nbsp;&nbsp;&nbsp;&nbsp; and same for <b>complementarityGap_</b><br>

      <font size="1"><b>18 Primal -1.3318263e+09 Dual -1.3318263e+09

          Complementarity 40.006314 - 134 fixed, rank 3924&nbsp;&nbsp; </b></font><br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if ((gapO &lt; 1.0e-6 ||<b> (gapO &lt; 1.0e-4 &amp;&amp;

        complementarityGap_ &lt; 0.1)) &amp;&amp; !sloppyOptimal)</b> {<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sloppyOptimal = true;<br>

      <br>

      <br>

      Good luck<br>

      <br>

      William<br>

      <br>

      <div class="moz-cite-prefix">On 11/27/2012 10:50 AM, Stephanie

        Jakus wrote:<br>

      </div>

      <blockquote

cite="mid:CAKCku4B12DyuC6V_tb1Tf=W+uSkf6UGR5pBcoHMKvjGmSsO9Xw@mail.gmail.com"

        type="cite"><font size="1">

          <div><font>Dear Clp mailing list,</font></div>

          <div><font>I'm using the Clp executable with the barrier

              method, and running into problems where the

              complementarity value gets small and then diverges. I'm

              wondering if there is an option to change the convergence

              tolerance value for the complementarity when using the

              barrier method. (I see the options for tolerance on primal

              and dual feasibility... but are those enough?)</font></div>

          <div>&nbsp;</div>

          <div><font>Thanks very much!</font></div>

          <div><font>Sincerely,</font></div>

          <div><font>Stephanie Jakus</font></div>

          <div>&nbsp;</div>

          <div>Coin LP version 1.14.8, build Nov 26 2012</div>

          <p>Clp takes input from arguments ( - switches to stdin)</p>

          <p>Enter ? for list of commands or help</p>

          <p>Clp:import /cygdrive/c/Users/StephanieJ/Desktop/3892762.mps</p>

          <p>At line 8 NAME Model</p>

          <p>At line 9 ROWS</p>

          <p>At line 1364 COLUMNS</p>

          <p>At line 4998 RHS</p>

          <p>At line 5038 QUADOBJ</p>

          <p>Problem Model has 1353 rows, 1218 columns and 4830 elements</p>

          <p>At line 6257 ENDATA</p>

          <p>Model was imported from

            /cygdrive/c/Users/StephanieJ/Desktop/3892762.mps in 0.015

            seconds</p>

          <p>Clp:presolve off</p>

          <p>Clp:barr</p>

          <p>911766 elements in sparse Cholesky, flop count 1.38983e+09</p>

          <p>0 Primal 2.0281794e+14 Dual -5.1285362e+16 Complementarity

            5.7194998e+16 - 0 fixed, rank 3924</p>

          <p>1 Primal 6.4401297e+11 Dual -1.4397286e+15 Complementarity

            4.5782888e+15 - 0 fixed, rank 3924</p>

          <p>2 Primal 9.0405659e+10 Dual -6.7439781e+14 Complementarity

            1.8042711e+15 - 0 fixed, rank 3924</p>

          <p>3 Primal -4.6446979e+08 Dual -2.2353578e+14 Complementarity

            2.5748948e+14 - 0 fixed, rank 3924</p>

          <p>4 Primal -6.1493992e+08 Dual -4.9256197e+13 Complementarity

            5.1971454e+13 - 0 fixed, rank 3924</p>

          <p>5 Primal -6.2519141e+08 Dual -9.5993463e+12 Complementarity

            9.8821256e+12 - 0 fixed, rank 3924</p>

          <p>6 Primal -6.4251057e+08 Dual -6.5604007e+11 Complementarity

            6.6713932e+11 - 0 fixed, rank 3924</p>

          <p>7 Primal -7.7860822e+08 Dual -3.6508427e+10 Complementarity

            3.6354368e+10 - 0 fixed, rank 3924</p>

          <p>8 Primal -1.0644153e+09 Dual -7.4875479e+09 Complementarity

            6.5249957e+09 - 0 fixed, rank 3924</p>

          <p>9 Primal -1.2085383e+09 Dual -3.5784324e+09 Complementarity

            2.4161606e+09 - 0 fixed, rank 3924</p>

          <p>10 Primal -1.2797358e+09 Dual -2.1299084e+09

            Complementarity 8.6888912e+08 - 0 fixed, rank 3924</p>

          <p>11 Primal -1.3185777e+09 Dual -1.4373535e+09

            Complementarity 1.2056077e+08 - 0 fixed, rank 3924</p>

          <p>12 Primal -1.3295276e+09 Dual -1.3574632e+09

            Complementarity 28320570 - 0 fixed, rank 3924</p>

          <p>13 Primal -1.3314339e+09 Dual -1.3388477e+09

            Complementarity 7525528.9 - 0 fixed, rank 3924</p>

          <p>14 Primal -1.3317815e+09 Dual -1.3331553e+09

            Complementarity 1386412.8 - 13 fixed, rank 3924</p>

          <p>15 Primal -1.3318236e+09 Dual -1.3319518e+09

            Complementarity 128945.31 - 46 fixed, rank 3924</p>

          <p>16 Primal -1.3318262e+09 Dual -1.3318287e+09

            Complementarity 2519.259 - 76 fixed, rank 3924</p>

          <p>17 Primal -1.3318263e+09 Dual -1.3318263e+09

            Complementarity 40.523345 - 134 fixed, rank 3924</p>

          <p>18 Primal -1.3318263e+09 Dual -1.3318263e+09

            Complementarity 40.006314 - 134 fixed, rank 3924</p>

          <p>19 Primal -1.3318263e+09 Dual -1.3318264e+09

            Complementarity 39.542777 - 134 fixed, rank 3924</p>

          <p>20 Primal -1.3318263e+09 Dual -1.3318265e+09

            Complementarity 39.128027 - 134 fixed, rank 3924</p>

          <p>21 Primal -1.3318263e+09 Dual -1.3318266e+09

            Complementarity 38.757635 - 134 fixed, rank 3924</p>

          <p>22 Primal -1.3318263e+09 Dual -1.3318263e+09

            Complementarity 0.87384878 - 134 fixed, rank 3924</p>

          <p>23 Primal -1.3318263e+09 Dual -1.331649e+09 Complementarity

            0.0061784205 - 134 fixed, rank 3924</p>

          <p>24 Primal -1.3318263e+09 Dual -1.3047041e+09

            Complementarity 4.0066434e-05 - 134 fixed, rank 3924</p>

          <p>25 Primal -1.3318262e+09 Dual 1.3928076e+09 Complementarity

            1.3321456e-05 - 134 fixed, rank 3924</p>

          <p>26 Primal -1.3318262e+09 Dual 1.3945074e+09 Complementarity

            3.4243478e-06 - 134 fixed, rank 3924</p>

          <p>27 Primal -1.3318262e+09 Dual 2.3040717e+11 Complementarity

            7.2045846e-07 - 134 fixed, rank 3924</p>

        </font><br>

      </blockquote>

    </blockquote>

    <br>

  </body>

</html>