Ignore:
Timestamp:
Oct 14, 2002 9:12:24 AM (17 years ago)
Author:
forrest
Message:

Testing presolve on nasty burlington problem

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/devel-1/Test/ClpMain.cpp

    r34 r41  
    10821082#ifdef USE_PRESOLVE
    10831083              Presolve pinfo;
    1084               if (preSolve)
     1084              if (preSolve) {
    10851085                model2 = pinfo.presolvedModel(models[iModel],1.0e-8,preSolve);
     1086                model2->checkSolution();
     1087#ifdef CLP_DEBUG
     1088                printf("%g %g (%d) %g (%d)\n"
     1089                       ,model2->objectiveValue()
     1090                       ,model2->sumDualInfeasibilities()
     1091                       ,model2->numberDualInfeasibilities()
     1092                       ,model2->sumPrimalInfeasibilities()
     1093                       ,model2->numberPrimalInfeasibilities());
     1094#endif
     1095                if (type==DUALSIMPLEX) {
     1096                  int numberInfeasibilities = model2->tightenPrimalBounds();
     1097                  if (numberInfeasibilities)
     1098                    std::cout<<"** Analysis indicates model infeasible"
     1099                             <<std::endl;
     1100                }
     1101              }
    10861102#endif
    10871103#ifdef READLINE     
     
    11031119#endif
    11041120                models[iModel].primal(1);
     1121#ifdef CLP_DEBUG
     1122                models[iModel].checkSolution();
     1123                printf("%g dual %g(%d) Primal %g(%d)\n",
     1124                       models[iModel].objectiveValue(),
     1125                       models[iModel].sumDualInfeasibilities(),
     1126                       models[iModel].numberDualInfeasibilities(),
     1127                       models[iModel].sumPrimalInfeasibilities(),
     1128                       models[iModel].numberPrimalInfeasibilities());
     1129                {
     1130                  Presolve pinfoA;
     1131                  model2 = pinfoA.presolvedModel(models[iModel],1.0e-8);
     1132                 
     1133                  printf("Resolving from presolved optimal solution\n");
     1134                  model2->primal(1);
     1135                 
     1136                  delete model2;
     1137                }
     1138#endif
    11051139              }
    11061140#endif
     
    11221156          case TIGHTEN:
    11231157            if (goodModels[iModel]) {
    1124               int numberInfeasibilities = models[iModel].tightenPrimalBounds();
     1158              int numberInfeasibilities = models[iModel].tightenPrimalBounds();
    11251159              if (numberInfeasibilities)
    11261160                std::cout<<"** Analysis indicates model infeasible"<<std::endl;
Note: See TracChangeset for help on using the changeset viewer.