Changeset 2354 for trunk/Cbc/src/CbcHeuristic.cpp

Timestamp:

Jan 2, 2018 6:34:25 AM (3 years ago)

Author:

forrest

Message:

try and allow SOS in mini BAB

File:

• trunk/Cbc/src/CbcHeuristic.cpp (modified) (1 diff)

trunk/Cbc/src/CbcHeuristic.cpp

@@ -1048 +1048 @@
                 // move seed across
                 model.randomNumberGenerator()->setSeed(model_->randomNumberGenerator()->getSeed());
+#ifdef COIN_HAS_CLP
+                // redo SOS
+                OsiClpSolverInterface * clpSolver
+                  = dynamic_cast<OsiClpSolverInterface *> (model.solver());
+                if (clpSolver && clpSolver->numberSOS()) {
+                  int numberColumns = clpSolver->getNumCols();
+                  const int * originalColumns = process.originalColumns();
+                  CoinSet * setInfo =
+                    const_cast<CoinSet *>(clpSolver->setInfo());
+                  int numberSOS = clpSolver->numberSOS();
+                  for (int iSOS = 0; iSOS < numberSOS; iSOS++) {
+                    //int type = setInfo[iSOS].setType();
+                    int n = setInfo[iSOS].numberEntries();
+                    int * which = setInfo[iSOS].modifiableWhich();
+                    double * weights = setInfo[iSOS].modifiableWeights();
+                    int n2=0;
+                    for (int j=0;j<n;j++) {
+                      int iColumn=which[j];
+                      int i;
+                      for (i = 0; i < numberColumns; i++) {
+                        if (originalColumns[i] == iColumn)
+                          break;
+                      }
+                      if (i < numberColumns) {
+                        which[n2] = i;
+                        weights[n2++] = weights[j];
+                      }
+                    }
+                    setInfo[iSOS].setNumberEntries(n2);
+                  }
+                }
+#endif
                 if (numberNodes >= 0) {
                     // normal