Ignore:
Timestamp:
Oct 2, 2003 1:21:02 PM (17 years ago)
Author:
forrest
Message:

lots of stuff

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/pre/ClpPresolve.cpp

    r181 r212  
    894894void check_djs(CoinPostsolveMatrix *prob)
    895895{
    896   return;
     896  //return;
    897897  double *colels        = prob->colels_;
    898898  int *hrow             = prob->hrow_;
     
    935935        double coeff = colels[k];
    936936        k = link[k];
    937 
    938937        dj -= rowduals[row] * coeff;
    939938        rsol[row] += solutionValue*coeff;
    940939      }
    941       if (! (fabs(rcosts[colx] - dj) < 100*ZTOLDP))
     940      if (! (fabs(rcosts[colx] - dj) < 1.0e-4))
    942941        printf("BAD DJ:  %d %g %g\n",
    943942               colx, rcosts[colx], dj);
     
    965964        double dj = rowduals[rowx];
    966965        if (rsol[rowx]<rlo[rowx]+1.0e-6) {
    967           if (dj <-1.0e-6)
     966          if (dj <-1.0e-5)
    968967            printf("neg rDJ:  %d %g  - %g %g %g\n",
    969968                   rowx, dj, rlo[rowx], rsol[rowx], rup[rowx]);
    970969        } else if (rsol[rowx]>rup[rowx]-1.0e-6) {
    971           if (dj > 1.0e-6)
     970          if (dj > 1.0e-5)
    972971            printf("pos rDJ:  %d %g  - %g %g %g\n",
    973972                   rowx, dj, rlo[rowx], rsol[rowx], rup[rowx]);
    974973        } else {
    975           if (fabs(dj) >1.0e-6)
     974          if (fabs(dj) >1.0e-5)
    976975            printf("nonzero rDJ:  %d %g  - %g %g %g\n",
    977976                   rowx, dj, rlo[rowx], rsol[rowx], rup[rowx]);
     
    10661065  return (! (i >= 0));
    10671066}
     1067#if     DEBUG_PRESOLVE
     1068static void matrix_bounds_ok(const double *lo, const double *up, int n)
     1069{
     1070  int i;
     1071  for (i=0; i<n; i++) {
     1072    PRESOLVEASSERT(lo[i] <= up[i]);
     1073    PRESOLVEASSERT(lo[i] < PRESOLVE_INF);
     1074    PRESOLVEASSERT(-PRESOLVE_INF < up[i]);
     1075  }
     1076}
     1077#endif
    10681078CoinPresolveMatrix::CoinPresolveMatrix(int ncols0_in,
    10691079                                     double maxmin_,
Note: See TracChangeset for help on using the changeset viewer.