Changeset 391 for trunk/Cbc/src/CoinSolve.cpp

Timestamp:

Jun 29, 2006 10:44:34 AM (15 years ago)

Author:

forrest

Message:

for ampl and version

File:

• trunk/Cbc/src/CoinSolve.cpp (modified) (9 diffs)

trunk/Cbc/src/CoinSolve.cpp

@@ -492 +492 @@
     double * solutionIn = NULL;
     int * prioritiesIn = NULL;
+    int numberSOS = 0;
+    int * sosStart = NULL;
+    int * sosIndices = NULL;
+    char * sosType = NULL;
+    double * sosReference = NULL;
+    int * sosPriority=NULL;
 #ifdef COIN_HAS_ASL
     ampl_info info;
@@ -1681 +1687 @@
               double timeLeft = babModel->getMaximumSeconds();
               int numberOriginalColumns = babModel->solver()->getNumCols();
-#ifdef COIN_HAS_ASL
-              if (usingAmpl&&info.numberSos&&doSOS) {
-                // SOS
-                assert (!preProcess); // do later
-                int numberSOS = info.numberSos;
-                int numberIntegers = babModel->numberIntegers();
-                int numberColumns = babModel->getNumCols();
-                /* model may not have created objects
-                   If none then create
-                */
-                if (!numberIntegers||!babModel->numberObjects()) {
-                  int type = (pseudoUp) ? 1 : 0;
-                  babModel->findIntegers(true,type);
-                  numberIntegers = babModel->numberIntegers();
-                }
-                // Do sets and priorities
-                CbcObject ** objects = new CbcObject * [numberSOS];
-                const int * starts = info.sosStart;
-                const int * which = info.sosIndices;
-                const char * type = info.sosType;
-                const double * weight = info.sosReference;
-                // see if any priorities
-                int i;
-                bool gotPriorities=false;
-                int * priorities=info.priorities;
-                if (priorities) {
-                  for (i=0;i<numberColumns;i++) {
-                    if (priorities[i]) {
-                      gotPriorities=true;
-                      break;
-                    }
-                  }
-                }
-                priorities=info.sosPriority;
-                if (priorities) {
-                  for (i=0;i<numberSOS;i++) {
-                    if (priorities[i]) {
-                      gotPriorities=true;
-                      break;
-                    }
-                  }
-                }
-                int iSOS;
-                for (iSOS =0;iSOS<numberSOS;iSOS++) {
-                  int iStart = starts[iSOS];
-                  int n=starts[iSOS+1]-iStart;
-                  objects[iSOS] = new CbcSOS(babModel,n,which+iStart,weight+iStart,
-                                             iSOS,type[iSOS]);
-                  // higher for set
-                  objects[iSOS]->setPriority(10);
-                  if (gotPriorities&&info.sosPriority&&info.sosPriority[iSOS])
-                    objects[iSOS]->setPriority(info.sosPriority[iSOS]);
-                }
-                babModel->addObjects(numberSOS,objects);
-                for (iSOS=0;iSOS<numberSOS;iSOS++)
-                  delete objects[iSOS];
-                delete [] objects;
-              }
-#endif
               if (preProcess==6) {
                 // use strategy instead
@@ -1765 +1712 @@
                   int translate[]={9999,0,0,-1,2,3};
                   process.messageHandler()->setLogLevel(babModel->logLevel());
+#ifdef COIN_HAS_ASL
+                  if (info.numberSos&&doSOS&&usingAmpl) {
+                    // SOS
+                    numberSOS = info.numberSos;
+                    sosStart = info.sosStart;
+                    sosIndices = info.sosIndices;
+                  }
+#endif
+                  if (numberSOS&&doSOS) {
+                    // SOS
+                    int numberColumns = saveSolver->getNumCols();
+                    char * prohibited = new char[numberColumns];
+                    memset(prohibited,0,numberColumns);
+                    int n=sosStart[numberSOS];
+                    for (int i=0;i<n;i++) {
+                      int iColumn = sosIndices[i];
+                      prohibited[iColumn]=1;
+                    }
+                    process.passInProhibited(prohibited,numberColumns);
+                    delete [] prohibited;
+                  }
                   solver2 = process.preProcessNonDefault(*saveSolver,translate[preProcess],10);
                   // Tell solver we are not in Branch and Cut
@@ -2045 +2013 @@
                   solutionIn=info.primalSolution;
                   prioritiesIn = info.priorities;
+                  if (info.numberSos&&doSOS) {
+                    // SOS
+                    numberSOS = info.numberSos;
+                    sosStart = info.sosStart;
+                    sosIndices = info.sosIndices;
+                    sosType = info.sosType;
+                    sosReference = info.sosReference;
+                    sosPriority = info.sosPriority;
+                  }
                 }
 #endif                
@@ -2135 +2112 @@
                     delete objects[iSOS];
                   delete [] objects;
-                } else if (priorities||branchDirection||pseudoDown||pseudoUp) {
+                } else if (priorities||branchDirection||pseudoDown||pseudoUp||numberSOS) {
                   // do anyway for priorities etc
                   int numberIntegers = babModel->numberIntegers();
@@ -2144 +2121 @@
                     int type = (pseudoUp) ? 1 : 0;
                     babModel->findIntegers(true,type);
+                  }
+                  if (numberSOS) {
+                    // Do sets and priorities
+                    CbcObject ** objects = new CbcObject * [numberSOS];
+                    int iSOS;
+                    if (originalColumns) {
+                      // redo sequence numbers
+                      int numberColumns = babModel->getNumCols();
+                      int nOld = originalColumns[numberColumns-1]+1;
+                      int * back = new int[nOld];
+                      int i;
+                      for (i=0;i<nOld;i++)
+                        back[i]=-1;
+                      for (i=0;i<numberColumns;i++)
+                        back[originalColumns[i]]=i;
+                      // Really need better checks
+                      int nMissing=0;
+                      int n=sosStart[numberSOS];
+                      for (i=0;i<n;i++) {
+                        int iColumn = sosIndices[i];
+                        int jColumn = back[iColumn];
+                        if (jColumn>=0) 
+                          sosIndices[i] = jColumn;
+                        else 
+                          nMissing++;
+                      }
+                      delete [] back;
+                      if (nMissing)
+                        printf("%d SOS variables vanished due to pre processing? - check validity?\n",nMissing);
+                    }
+                    for (iSOS =0;iSOS<numberSOS;iSOS++) {
+                      int iStart = sosStart[iSOS];
+                      int n=sosStart[iSOS+1]-iStart;
+                      objects[iSOS] = new CbcSOS(babModel,n,sosIndices+iStart,sosReference+iStart,
+                                                 iSOS,sosType[iSOS]);
+                      if (sosPriority)
+                        objects[iSOS]->setPriority(sosPriority[iSOS]);
+                      else if (!prioritiesIn)
+                        objects[iSOS]->setPriority(10);  // rather than 1000 
+                    }
+                    babModel->addObjects(numberSOS,objects);
+                    for (iSOS=0;iSOS<numberSOS;iSOS++)
+                      delete objects[iSOS];
+                    delete [] objects;
                   }
                   CbcObject ** objects = babModel->objects();
@@ -2191 +2212 @@
                   free(prioritiesIn);
                   prioritiesIn=NULL;
+                  free(sosStart);
+                  sosStart=NULL;
+                  free(sosIndices);
+                  sosIndices=NULL;
+                  free(sosType);
+                  sosType=NULL;
+                  free(sosReference);
+                  sosReference=NULL;
+                  free(sosPriority);
+                  sosPriority=NULL;
 #ifdef COIN_HAS_ASL
                 }
@@ -2362 +2393 @@
 #endif
               } else {
-                std::cout<<"Model strengthend - now has "<<clpSolver->getNumRows()
+                std::cout<<"Model strengthened - now has "<<clpSolver->getNumRows()
                          <<" rows"<<std::endl;
               }
@@ -2389 +2420 @@
                 free(prioritiesIn);
                 prioritiesIn=NULL;
+                free(sosStart);
+                sosStart=NULL;
+                free(sosIndices);
+                sosIndices=NULL;
+                free(sosType);
+                sosType=NULL;
+                free(sosReference);
+                sosReference=NULL;
+                free(sosPriority);
+                sosPriority=NULL;
 #ifdef COIN_HAS_ASL
               }