Ignore:
Timestamp:
Mar 26, 2013 11:23:38 AM (7 years ago)
Author:
forrest
Message:

try to fix infeasibility ray,
changes as in stable (for presolve),
stuff for Cbc parameters

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Clp/src/ClpSimplexDual.cpp

    r1924 r1926  
    15941594                    solution_[sequenceOut_] = valueOut_;
    15951595                    int whatNext = housekeeping(objectiveChange);
     1596#if 0
     1597                    for (int i=0;i<numberRows_+numberColumns_;i++) {
     1598                      if (getStatus(i)==atLowerBound) {
     1599                        assert (dj_[i]>-1.0e-5);
     1600                        assert (solution_[i]<=lower_[i]+1.0e-5);
     1601                      } else if (getStatus(i)==atUpperBound) {
     1602                        assert (dj_[i]<1.0e-5);
     1603                        assert (solution_[i]>=upper_[i]-1.0e-5);
     1604                      }
     1605                    }
     1606#endif
    15961607#ifdef CLP_REPORT_PROGRESS
    15971608                    if (ixxxxxx > ixxyyyy - 5) {
     
    17421753#endif
    17431754                    // no incoming column is valid
     1755                    spareIntArray_[3]=pivotRow_;
    17441756                    pivotRow_ = -1;
    17451757#ifdef CLP_DEBUG
     
    17531765                         if ((specialOptions_&(1024 | 4096)) == 0 || (specialOptions_ & 32) != 0) {
    17541766                              // create ray anyway
     1767#ifdef PRINT_RAY_METHOD
     1768                           printf("Dual creating infeasibility ray direction out %d - pivRow %d seqOut %d lower %g,val %g,upper %g\n",
     1769                                  directionOut_,spareIntArray_[3],sequenceOut_,lowerOut_,valueOut_,upperOut_);
     1770#endif
    17551771                              ray_ = new double [ numberRows_];
    17561772                              rowArray_[0]->expand(); // in case packed
    1757                               // swap sign NO
    17581773                              const double * array = rowArray_[0]->denseVector();
    17591774                              for (int i=0;i<numberRows_;i++)
Note: See TracChangeset for help on using the changeset viewer.