Changeset 2474

Timestamp:

Jan 21, 2019 4:11:14 PM (3 months ago)

Author:

unxusr

Message:

fixed modk cuts

File:

• trunk/Cbc/examples/modk.c (modified) (7 diffs)

trunk/Cbc/examples/modk.c

@@ -16 +16 @@
 #include <stdlib.h>
 #include <limits.h>
+#include <assert.h>
 #include <float.h>
 #include <ctype.h>
@@ -34 +35 @@
 void cutcallback( void *osiSolver, void *osiCuts, void *appData )
 {
-  printf("entered callback");
+  printf("entered callback\n");
   int m = Osi_getNumRows( osiSolver );
   int n = Osi_getNumCols( osiSolver );
@@ -103 +104 @@
 
     // removing those which are divisors of larger values
-    for ( int j=0 ; (j<nz-1) ; ++j )
-    {
-      for ( int jl=j+1 ; (jl<nz) ; ++jl )
-      {
+    for ( int j=0 ; (j<ndiff-1) ; ++j )
+    {
+      for ( int jl=j+1 ; (jl<ndiff) ; ++jl )
+      {
+        if (!dcoef[j])
+          continue;
         if ( dcoef[jl] % dcoef[j] == 0 )
         {
@@ -166 +169 @@
             continue;
 
-          OsiCuts_addRowCut( osiCuts, nz, cidx, coef, 'L', rhsCut );
+          OsiCuts_addRowCut( osiCuts, nzCut, cidx, coef, 'L', rhsCut );
         } // constraint in the form <= 
       } // every numerator < k
@@ -234 +237 @@
   double sumLHS = 0.0;
 
-  *nzCut = 0;
+  (*nzCut) = 0;
   for ( int j=0 ; j<nz ; ++j )
   {
     int col = rIdx[j];
+    assert( col >= 0 && col<Osi_getNumCols(osiSolver) );
 
     int cc = (icoef[j]*num)/den;
@@ -248 +252 @@
     sumLHS += ((double)cc) * x[col];
 
-    ++(*nzCut);
+    (*nzCut)++;
   }
 
@@ -268 +272 @@
   Cbc_readMps( mip, argv[1] );
 
+  Cbc_setParameter( mip, "cuts", "off" );
+  Cbc_setParameter( mip, "gomory", "ifmove" );
+  Cbc_setParameter( mip, "preprocess", "off" );
+
   Cbc_addCutCallback( mip,  cutcallback, "modk", NULL );
+
 
   Cbc_solve( mip );