Ignore:
Timestamp:
Aug 11, 2012 3:56:05 AM (7 years ago)
Author:
forrest
Message:

cleanup some coding in outDuplicateRows

File:
1 edited

Legend:

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

    r1869 r1876  
    76887688#endif
    76897689  double inverseCleanup = (cleanUp>0.0) ? 1.0/cleanUp : 0.0;
     7690  /*
     7691    Rules for cleaning up rhs
     7692    if cleanup != 0.0 just use cleanup
     7693    otherwise if problemStatus_==0 (i.e. was optimal) use rowActivity
     7694    then go to integer value if close
     7695   */
    76907696  //#define PRINT_DUP
    76917697#ifdef PRINT_DUP
     
    78357841            break;
    78367842          } else if (fabs(rup2-rlo2)<=tolerance) {
    7837             // equal - choose closer to zero
    7838             if (fabs(rup2)<fabs(rlo2))
    7839               rlo2=rup2;
    7840             else
     7843            // equal
     7844            if (!inverseCleanup) {
     7845              if (!problemStatus_) {
     7846                // choose one closer to rowActivity
     7847                double rval2=rowActivity_[iThis];
     7848                if (fabs(rup2-rval2)<fabs(rlo2-rval2))
     7849                  rlo2=rup2;
     7850              } else {
     7851                // no solution - closer to zero
     7852                if (fabs(rup2)<fabs(rlo2))
     7853                  rlo2=rup2;
     7854              }
     7855              // see if close to integer value
     7856              double value2 = floor(rlo2+0.5);
     7857              if (fabs(rlo2-value2)<1.0e-9)
     7858                rlo2=value2;
    78417859              rup2=rlo2;
     7860            } else {
     7861              // we will be cleaning up later - choose closer to zero
     7862              if (fabs(rup2)<fabs(rlo2))
     7863                rlo2=rup2;
     7864              else
     7865                rup2=rlo2;
     7866            }
    78427867#if 0
    78437868          if (rowLength[iThis]<4)
Note: See TracChangeset for help on using the changeset viewer.