Changeset 1435 for stable


Ignore:
Timestamp:
Sep 25, 2009 11:22:06 AM (10 years ago)
Author:
forrest
Message:

don't allow zero pivot in dualColumn

File:
1 edited

Legend:

Unmodified
Added
Removed
  • stable/1.10/Clp/src/ClpSimplexDual.cpp

    r1406 r1435  
    36553655  }
    36563656
    3657   if (sequenceIn_>=0) {
    36583657#ifdef MORE_CAREFUL
    36593658    // If we have done pivots and things look bad set alpha_ 0.0 to force factorization
    3660     if (badSumPivots) {
     3659    if (badSumPivots||
     3660        fabs(theta_*badFree)>10.0*dualTolerance_&&factorization_->pivots()) {
    36613661      if (handler_->logLevel()>1)
    36623662        printf("forcing re-factorization\n");
    3663       alpha_=0.0;
    3664     }
    3665     if (fabs(theta_*badFree)>10.0*dualTolerance_&&factorization_->pivots()) {
    3666       if (handler_->logLevel()>1)
    3667         printf("forcing re-factorizationon free\n");
    3668       alpha_=0.0;
    3669     }
    3670 #endif
     3663      sequenceIn_=-1;
     3664    }
     3665#endif
     3666  if (sequenceIn_>=0) {
    36713667    lowerIn_ = lower_[sequenceIn_];
    36723668    upperIn_ = upper_[sequenceIn_];
     
    37333729      lowerIn_=valueIn_;
    37343730    }
     3731  } else {
     3732    // no pivot
     3733    bestPossible=0.0;
     3734    alpha_=0.0;
    37353735  }
    37363736  //if (thisIncrease)
Note: See TracChangeset for help on using the changeset viewer.