Ignore:
Timestamp:
Dec 2, 2005 6:46:13 AM (15 years ago)
Author:
forrest
Message:

for quadratic barrier

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/ClpQuadraticObjective.cpp

    r468 r700  
    337337                if (iColumn!=jColumn) {
    338338                  offset += valueI*valueJ*elementValue;
     339                  //if (fabs(valueI*valueJ*elementValue)>1.0e-12)
     340                  //printf("%d %d %g %g %g -> %g\n",
     341                  //       iColumn,jColumn,valueI,valueJ,elementValue,
     342                  //       valueI*valueJ*elementValue);
    339343                  double gradientI = valueJ*elementValue;
    340344                  double gradientJ = valueI*elementValue;
     
    343347                } else {
    344348                  offset += 0.5*valueI*valueI*elementValue;
     349                  //if (fabs(valueI*valueI*elementValue)>1.0e-12)
     350                  //printf("XX %d %g %g -> %g\n",
     351                  //       iColumn,valueI,elementValue,
     352                  //       0.5*valueI*valueI*elementValue);
    345353                  double gradientI = valueI*elementValue;
    346354                  gradient_[iColumn] += gradientI;
     
    369377        }
    370378      }
    371       offset *= model->optimizationDirection()*model->objectiveScale();
     379      if (model)
     380        offset *= model->optimizationDirection()*model->objectiveScale();
    372381      return gradient_;
    373382    }
     
    462471      }
    463472    }
    464     offset *= model->optimizationDirection();
     473    if (model)
     474      offset *= model->optimizationDirection();
    465475    return gradient_;
    466476  }
Note: See TracChangeset for help on using the changeset viewer.