Changeset 2289 for trunk


Ignore:
Timestamp:
Nov 10, 2017 1:47:57 PM (9 months ago)
Author:
forrest
Message:

zero divide and last change to ClpSimplexDual? not such a good idea - no idea why

Location:
trunk/Clp/src
Files:
2 edited

Legend:

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

    r2288 r2289  
    38883888          }
    38893889     }
     3890     if (devex_<1.001e-30) {
     3891       COIN_DETAIL_PRINT(printf("devex of incoming tiny %d %g\n",sequenceIn,devex_));
     3892       devex_=1.0e-30;
     3893     }
    38903894     double oldDevex = weights_[sequenceIn];
    38913895#ifdef CLP_DEBUG
     
    39243928     if (pivotRow >= 0) {
    39253929          // set outgoing weight here
    3926           weights_[model_->sequenceOut()] = devex_ / (model_->alpha() * model_->alpha());
     3930       double alpha = model_->alpha();
     3931       if (fabs(alpha)>1.0e15) {
     3932         COIN_DETAIL_PRINT(printf("alpha %g for %d !!\n",alpha,model_->sequenceOut()));
     3933         alpha = 1.0e15;
     3934       }
     3935       weights_[model_->sequenceOut()] = devex_ / (alpha*alpha);
    39273936     }
    39283937}
  • trunk/Clp/src/ClpSimplexDual.cpp

    r2287 r2289  
    500500              (moreSpecialOptions_&256)==0 &&
    501501              ((progress_.lastObjective(0)<-1.0e10 &&
    502                 progress_.lastObjective(1)>-1.0e5)||sumPrimalInfeasibilities_>1.0e10*smallestPrimalInfeasibility)) {
     502                progress_.lastObjective(1)>-1.0e5))) {
    503503            // problems - try primal
    504504            problemStatus_=10;
Note: See TracChangeset for help on using the changeset viewer.