Changeset 370


Ignore:
Timestamp:
May 18, 2004 11:26:02 AM (14 years ago)
Author:
forrest
Message:

stuff

Location:
trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/ClpFactorization.cpp

    r338 r370  
    182182          numberRowBasic = numberBasic - numberColumnBasic;
    183183        numberElements = 3 * numberBasic + 3 * numberElements + 10000;
     184#if 0
    184185        // If iteration not zero then may be compressed
    185186        getAreas ( !model->numberIterations() ? numberRows : numberBasic,
    186187                   numberRowBasic+numberColumnBasic, numberElements,
    187188                   2 * numberElements );
     189#else
     190        getAreas ( numberRows,
     191                   numberRowBasic+numberColumnBasic, numberElements,
     192                   2 * numberElements );
     193#endif
    188194        //fill
    189195        //copy
  • trunk/ClpSimplexDual.cpp

    r361 r370  
    246246        if (getStatus(i)==atLowerBound) {
    247247          if (saveDuals[i]<0.0) {
    248             if (saveDuals[i]<-1.0e-3)
    249               printf("bad dj at lb %d %g\n",i,saveDuals[i]);
     248            //if (saveDuals[i]<-1.0e-3)
     249            //printf("bad dj at lb %d %g\n",i,saveDuals[i]);
    250250            saveDuals[i]=0.0;
    251251          }
    252252        } else if (getStatus(i)==atUpperBound) {
    253253          if (saveDuals[i]>0.0) {
    254             if (saveDuals[i]>1.0e-3)
    255               printf("bad dj at ub %d %g\n",i,saveDuals[i]);
     254            //if (saveDuals[i]>1.0e-3)
     255            //printf("bad dj at ub %d %g\n",i,saveDuals[i]);
    256256            saveDuals[i]=0.0;
    257257          }
  • trunk/ClpSimplexPrimal.cpp

    r348 r370  
    24822482          int * which=columnArray->getIndices();
    24832483          for (iColumn=0;iColumn<numberRows_+numberColumns_;iColumn++) {
    2484             if (getStatus(iColumn)==superBasic) {
    2485               if (fabs(solution_[iColumn]-lower_[iColumn])<=primalTolerance_) {
    2486                 solution_[iColumn]=lower_[iColumn];
    2487                 setStatus(iColumn,atLowerBound);
    2488               } else if (fabs(solution_[iColumn]-upper_[iColumn])
    2489                          <=primalTolerance_) {
    2490                 solution_[iColumn]=upper_[iColumn];
    2491                 setStatus(iColumn,atUpperBound);
    2492               } else if (lower_[iColumn]<-1.0e20&&upper_[iColumn]>1.0e20) {
    2493                 setStatus(iColumn,isFree);
    2494                 break;
    2495               } else if (!flagged(iColumn)) {
    2496                 // put ones near bounds at end after sorting
    2497                 work[number]= - min(0.1*(solution_[iColumn]-lower_[iColumn]),
    2498                                     upper_[iColumn]-solution_[iColumn]);
    2499                 which[number++] = iColumn;
     2484            if (!flagged(iColumn)) {
     2485              if (getStatus(iColumn)==superBasic) {
     2486                if (fabs(solution_[iColumn]-lower_[iColumn])<=primalTolerance_) {
     2487                  solution_[iColumn]=lower_[iColumn];
     2488                  setStatus(iColumn,atLowerBound);
     2489                } else if (fabs(solution_[iColumn]-upper_[iColumn])
     2490                           <=primalTolerance_) {
     2491                  solution_[iColumn]=upper_[iColumn];
     2492                  setStatus(iColumn,atUpperBound);
     2493                } else if (lower_[iColumn]<-1.0e20&&upper_[iColumn]>1.0e20) {
     2494                  setStatus(iColumn,isFree);
     2495                  break;
     2496                } else if (!flagged(iColumn)) {
     2497                  // put ones near bounds at end after sorting
     2498                  work[number]= - min(0.1*(solution_[iColumn]-lower_[iColumn]),
     2499                                      upper_[iColumn]-solution_[iColumn]);
     2500                  which[number++] = iColumn;
     2501                }
    25002502              }
    25012503            }
     
    25192521      } else {
    25202522        for (;iColumn<numberRows_+numberColumns_;iColumn++) {
    2521           if (getStatus(iColumn)==superBasic) {
    2522             if (fabs(solution_[iColumn]-lower_[iColumn])<=primalTolerance_) {
    2523               solution_[iColumn]=lower_[iColumn];
    2524               setStatus(iColumn,atLowerBound);
    2525             } else if (fabs(solution_[iColumn]-upper_[iColumn])
    2526                        <=primalTolerance_) {
    2527               solution_[iColumn]=upper_[iColumn];
    2528               setStatus(iColumn,atUpperBound);
    2529             } else if (lower_[iColumn]<-1.0e20&&upper_[iColumn]>1.0e20) {
    2530               setStatus(iColumn,isFree);
    2531               break;
    2532             } else {
    2533               break;
     2523          if (!flagged(iColumn)) {
     2524            if (getStatus(iColumn)==superBasic) {
     2525              if (fabs(solution_[iColumn]-lower_[iColumn])<=primalTolerance_) {
     2526                solution_[iColumn]=lower_[iColumn];
     2527                setStatus(iColumn,atLowerBound);
     2528              } else if (fabs(solution_[iColumn]-upper_[iColumn])
     2529                         <=primalTolerance_) {
     2530                solution_[iColumn]=upper_[iColumn];
     2531                setStatus(iColumn,atUpperBound);
     2532              } else if (lower_[iColumn]<-1.0e20&&upper_[iColumn]>1.0e20) {
     2533                setStatus(iColumn,isFree);
     2534                break;
     2535              } else {
     2536                break;
     2537              }
    25342538            }
    25352539          }
Note: See TracChangeset for help on using the changeset viewer.