Ignore:
Timestamp:
Dec 18, 2006 9:39:29 AM (13 years ago)
Author:
forrest
Message:

mostly for rins

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/devel/Cbc/src/CbcLinked.cpp

    r501 r502  
    255255    writeMps(temp);
    256256    printf("wrote bb%d\n",xxxxxx);
    257 #if 0
    258     const double * lower = getColLower();
    259     const double * upper = getColUpper();
    260     for (int i=0;i<8;i++) {
    261       printf("%d bounds %g %g\n",i,lower[i],upper[i]);
    262     }
    263     if (lower[2]<=3.0&&upper[2]>=3) {
    264       if (lower[3]<=2.0&&upper[3]>=2.0) {
    265         if (lower[4]==0.0) {
    266           if (lower[7]==0.0) {
    267             if (lower[5]<=4.0&&upper[5]>=4.0) {
    268               if (lower[6]<=3.0&&upper[6]>=3.0) {
    269                 printf("writemps contains solution\n");
    270                 for (int i=8;i<24;i++) {
    271                   printf("%d bounds %g %g\n",i,lower[i],upper[i]);
    272                   if (false&&((!lower[i]&&upper[i]>lower[i]))) {
    273                     setColLower(i,0.0);
    274                     setColUpper(i,10.0);
    275                   }
    276                 }
    277               }
    278             }
    279           }
    280         }
    281       }
    282     }
    283 #endif
    284257  }
    285258#endif
     
    356329    }
    357330    if (isProvenOptimal()) {
     331      if (satisfied)
     332        printf("satisfied %d\n",satisfied);
    358333      if (satisfied&&(specialOptions2_&2)!=0) {
    359334        assert (quadraticModel_);
     
    586561              for ( i=rowStart[objectiveRow_];i<rowStart[objectiveRow_+1];i++)
    587562                gradient[column2[i]] = element[i];
     563              const double * columnLower = modelPtr_->columnLower();
     564              const double * columnUpper = modelPtr_->columnUpper();
    588565              for ( i =0;i<numberObjects_;i++) {
    589566                OsiBiLinear * obj = dynamic_cast<OsiBiLinear *> (object_[i]);
     
    591568                  int xColumn = obj->xColumn();
    592569                  int yColumn = obj->yColumn();
     570                  // We can do better if one fixed
     571                  assert (columnUpper[xColumn]>columnLower[xColumn]);
     572                  assert (columnUpper[yColumn]>columnLower[yColumn]);
    593573                  if (xColumn!=yColumn) {
    594574                    double coefficient = 2.0*obj->coefficient();
     
    648628            double * gradient = new double [numberColumns2];
    649629            int * column = new int[numberColumns2];
     630            const double * columnLower = modelPtr_->columnLower();
     631            const double * columnUpper = modelPtr_->columnUpper();
    650632            for (int iNon=0;iNon<numberNonLinearRows_;iNon++) {
    651633              int iRow = rowNonLinear_[iNon];
     
    666648                int xColumn = obj->xColumn();
    667649                int yColumn = obj->yColumn();
     650                // We can do better if one fixed
     651                assert (columnUpper[xColumn]>columnLower[xColumn]);
     652                assert (columnUpper[yColumn]>columnLower[yColumn]);
    668653                if (xColumn!=yColumn) {
    669654                  double coefficient = 2.0*obj->coefficient();
Note: See TracChangeset for help on using the changeset viewer.