ClpModel.hpp

Go to the documentation of this file.
00001 /* $Id: ClpModel.hpp 1941 2013-04-10 16:52:27Z stefan $ */
00002 // Copyright (C) 2002, International Business Machines
00003 // Corporation and others.  All Rights Reserved.
00004 // This code is licensed under the terms of the Eclipse Public License (EPL).
00005 
00006 #ifndef ClpModel_H
00007 #define ClpModel_H
00008 
00009 #include "ClpConfig.h"
00010 
00011 #include <iostream>
00012 #include <cassert>
00013 #include <cmath>
00014 #include <vector>
00015 #include <string>
00016 //#ifndef COIN_USE_CLP
00017 //#define COIN_USE_CLP
00018 //#endif
00019 #include "ClpPackedMatrix.hpp"
00020 #include "CoinMessageHandler.hpp"
00021 #include "CoinHelperFunctions.hpp"
00022 #include "CoinTypes.hpp"
00023 #include "CoinFinite.hpp"
00024 #include "ClpParameters.hpp"
00025 #include "ClpObjective.hpp"
00026 class ClpEventHandler;
00036 class CoinBuild;
00037 class CoinModel;
00038 class ClpModel {
00039 
00040 public:
00041 
00047 
00048      ClpModel (bool emptyMessages = false  );
00049 
00054      ClpModel(const ClpModel & rhs, int scalingMode = -1);
00056      ClpModel & operator=(const ClpModel & rhs);
00061      ClpModel (const ClpModel * wholeModel,
00062                int numberRows, const int * whichRows,
00063                int numberColumns, const int * whichColumns,
00064                bool dropNames = true, bool dropIntegers = true);
00066      ~ClpModel (  );
00068 
00082      void loadProblem (  const ClpMatrixBase& matrix,
00083                          const double* collb, const double* colub,
00084                          const double* obj,
00085                          const double* rowlb, const double* rowub,
00086                          const double * rowObjective = NULL);
00087      void loadProblem (  const CoinPackedMatrix& matrix,
00088                          const double* collb, const double* colub,
00089                          const double* obj,
00090                          const double* rowlb, const double* rowub,
00091                          const double * rowObjective = NULL);
00092 
00095      void loadProblem (  const int numcols, const int numrows,
00096                          const CoinBigIndex* start, const int* index,
00097                          const double* value,
00098                          const double* collb, const double* colub,
00099                          const double* obj,
00100                          const double* rowlb, const double* rowub,
00101                          const double * rowObjective = NULL);
00107      int loadProblem (  CoinModel & modelObject, bool tryPlusMinusOne = false);
00109      void loadProblem (  const int numcols, const int numrows,
00110                          const CoinBigIndex* start, const int* index,
00111                          const double* value, const int * length,
00112                          const double* collb, const double* colub,
00113                          const double* obj,
00114                          const double* rowlb, const double* rowub,
00115                          const double * rowObjective = NULL);
00117      void loadQuadraticObjective(const int numberColumns,
00118                                  const CoinBigIndex * start,
00119                                  const int * column, const double * element);
00120      void loadQuadraticObjective (  const CoinPackedMatrix& matrix);
00122      void deleteQuadraticObjective();
00124      void setRowObjective(const double * rowObjective);
00126      int readMps(const char *filename,
00127                  bool keepNames = false,
00128                  bool ignoreErrors = false);
00130      int readGMPL(const char *filename, const char * dataName,
00131                   bool keepNames = false);
00133      void copyInIntegerInformation(const char * information);
00135      void deleteIntegerInformation();
00137      void setContinuous(int index);
00139      void setInteger(int index);
00141      bool isInteger(int index) const;
00143      void resize (int newNumberRows, int newNumberColumns);
00145      void deleteRows(int number, const int * which);
00147      void addRow(int numberInRow, const int * columns,
00148                  const double * elements, double rowLower = -COIN_DBL_MAX,
00149                  double rowUpper = COIN_DBL_MAX);
00151      void addRows(int number, const double * rowLower,
00152                   const double * rowUpper,
00153                   const CoinBigIndex * rowStarts, const int * columns,
00154                   const double * elements);
00156      void addRows(int number, const double * rowLower,
00157                   const double * rowUpper,
00158                   const CoinBigIndex * rowStarts, const int * rowLengths,
00159                   const int * columns,
00160                   const double * elements);
00161 #ifndef CLP_NO_VECTOR
00162      void addRows(int number, const double * rowLower,
00163                   const double * rowUpper,
00164                   const CoinPackedVectorBase * const * rows);
00165 #endif
00166 
00171      int addRows(const CoinBuild & buildObject, bool tryPlusMinusOne = false,
00172                  bool checkDuplicates = true);
00181      int addRows(CoinModel & modelObject, bool tryPlusMinusOne = false,
00182                  bool checkDuplicates = true);
00183 
00185      void deleteColumns(int number, const int * which);
00187      void deleteRowsAndColumns(int numberRows, const int * whichRows,
00188                                int numberColumns, const int * whichColumns);
00190      void addColumn(int numberInColumn,
00191                     const int * rows,
00192                     const double * elements,
00193                     double columnLower = 0.0,
00194                     double  columnUpper = COIN_DBL_MAX,
00195                     double  objective = 0.0);
00197      void addColumns(int number, const double * columnLower,
00198                      const double * columnUpper,
00199                      const double * objective,
00200                      const CoinBigIndex * columnStarts, const int * rows,
00201                      const double * elements);
00202      void addColumns(int number, const double * columnLower,
00203                      const double * columnUpper,
00204                      const double * objective,
00205                      const CoinBigIndex * columnStarts, const int * columnLengths,
00206                      const int * rows,
00207                      const double * elements);
00208 #ifndef CLP_NO_VECTOR
00209      void addColumns(int number, const double * columnLower,
00210                      const double * columnUpper,
00211                      const double * objective,
00212                      const CoinPackedVectorBase * const * columns);
00213 #endif
00214 
00219      int addColumns(const CoinBuild & buildObject, bool tryPlusMinusOne = false,
00220                     bool checkDuplicates = true);
00228      int addColumns(CoinModel & modelObject, bool tryPlusMinusOne = false,
00229                     bool checkDuplicates = true);
00231      inline void modifyCoefficient(int row, int column, double newElement,
00232                                    bool keepZero = false) {
00233           matrix_->modifyCoefficient(row, column, newElement, keepZero);
00234      }
00236      void chgRowLower(const double * rowLower);
00238      void chgRowUpper(const double * rowUpper);
00240      void chgColumnLower(const double * columnLower);
00242      void chgColumnUpper(const double * columnUpper);
00244      void chgObjCoefficients(const double * objIn);
00248      void borrowModel(ClpModel & otherModel);
00251      void returnModel(ClpModel & otherModel);
00252 
00254      void createEmptyMatrix();
00262      int cleanMatrix(double threshold = 1.0e-20);
00264      void copy(const ClpMatrixBase * from, ClpMatrixBase * & to);
00265 #ifndef CLP_NO_STD
00267      void dropNames();
00269      void copyNames(const std::vector<std::string> & rowNames,
00270                     const std::vector<std::string> & columnNames);
00272      void copyRowNames(const std::vector<std::string> & rowNames, int first, int last);
00274      void copyColumnNames(const std::vector<std::string> & columnNames, int first, int last);
00276      void copyRowNames(const char * const * rowNames, int first, int last);
00278      void copyColumnNames(const char * const * columnNames, int first, int last);
00280      void setRowName(int rowIndex, std::string & name) ;
00282      void setColumnName(int colIndex, std::string & name) ;
00283 #endif
00284 
00291      int findNetwork(char * rotate, double fractionNeeded = 0.75);
00294      CoinModel * createCoinModel() const;
00295 
00308      int writeMps(const char *filename,
00309                   int formatType = 0, int numberAcross = 2,
00310                   double objSense = 0.0) const ;
00312 
00314 
00315      inline int numberRows() const {
00316           return numberRows_;
00317      }
00318      inline int getNumRows() const {
00319           return numberRows_;
00320      }
00322      inline int getNumCols() const {
00323           return numberColumns_;
00324      }
00325      inline int numberColumns() const {
00326           return numberColumns_;
00327      }
00329      inline double primalTolerance() const {
00330           return dblParam_[ClpPrimalTolerance];
00331      }
00332      void setPrimalTolerance( double value) ;
00334      inline double dualTolerance() const  {
00335           return dblParam_[ClpDualTolerance];
00336      }
00337      void setDualTolerance( double value) ;
00339      inline double primalObjectiveLimit() const {
00340           return dblParam_[ClpPrimalObjectiveLimit];
00341      }
00342      void setPrimalObjectiveLimit(double value);
00344      inline double dualObjectiveLimit() const {
00345           return dblParam_[ClpDualObjectiveLimit];
00346      }
00347      void setDualObjectiveLimit(double value);
00349      inline double objectiveOffset() const {
00350           return dblParam_[ClpObjOffset];
00351      }
00352      void setObjectiveOffset(double value);
00354      inline double presolveTolerance() const {
00355           return dblParam_[ClpPresolveTolerance];
00356      }
00357 #ifndef CLP_NO_STD
00358      inline const std::string & problemName() const {
00359           return strParam_[ClpProbName];
00360      }
00361 #endif
00363      inline int numberIterations() const  {
00364           return numberIterations_;
00365      }
00366      inline int getIterationCount() const {
00367           return numberIterations_;
00368      }
00369      inline void setNumberIterations(int numberIterationsNew) {
00370           numberIterations_ = numberIterationsNew;
00371      }
00373      inline int solveType() const {
00374           return solveType_;
00375      }
00376      inline void setSolveType(int type) {
00377           solveType_ = type;
00378      }
00380      inline int maximumIterations() const {
00381           return intParam_[ClpMaxNumIteration];
00382      }
00383      void setMaximumIterations(int value);
00385      inline double maximumSeconds() const {
00386           return dblParam_[ClpMaxSeconds];
00387      }
00388      void setMaximumSeconds(double value);
00390      bool hitMaximumIterations() const;
00400      inline int status() const            {
00401           return problemStatus_;
00402      }
00403      inline int problemStatus() const            {
00404           return problemStatus_;
00405      }
00407      inline void setProblemStatus(int problemStatusNew) {
00408           problemStatus_ = problemStatusNew;
00409      }
00424      inline int secondaryStatus() const            {
00425           return secondaryStatus_;
00426      }
00427      inline void setSecondaryStatus(int newstatus) {
00428           secondaryStatus_ = newstatus;
00429      }
00431      inline bool isAbandoned() const             {
00432           return problemStatus_ == 4;
00433      }
00435      inline bool isProvenOptimal() const         {
00436           return problemStatus_ == 0;
00437      }
00439      inline bool isProvenPrimalInfeasible() const {
00440           return problemStatus_ == 1;
00441      }
00443      inline bool isProvenDualInfeasible() const  {
00444           return problemStatus_ == 2;
00445      }
00447      bool isPrimalObjectiveLimitReached() const ;
00449      bool isDualObjectiveLimitReached() const ;
00451      inline bool isIterationLimitReached() const {
00452           return problemStatus_ == 3;
00453      }
00455      inline double optimizationDirection() const {
00456           return  optimizationDirection_;
00457      }
00458      inline double getObjSense() const    {
00459           return optimizationDirection_;
00460      }
00461      void setOptimizationDirection(double value);
00463      inline double * primalRowSolution() const    {
00464           return rowActivity_;
00465      }
00466      inline const double * getRowActivity() const {
00467           return rowActivity_;
00468      }
00470      inline double * primalColumnSolution() const {
00471           return columnActivity_;
00472      }
00473      inline const double * getColSolution() const {
00474           return columnActivity_;
00475      }
00476      inline void setColSolution(const double * input) {
00477           memcpy(columnActivity_, input, numberColumns_ * sizeof(double));
00478      }
00480      inline double * dualRowSolution() const      {
00481           return dual_;
00482      }
00483      inline const double * getRowPrice() const    {
00484           return dual_;
00485      }
00487      inline double * dualColumnSolution() const   {
00488           return reducedCost_;
00489      }
00490      inline const double * getReducedCost() const {
00491           return reducedCost_;
00492      }
00494      inline double* rowLower() const              {
00495           return rowLower_;
00496      }
00497      inline const double* getRowLower() const     {
00498           return rowLower_;
00499      }
00501      inline double* rowUpper() const              {
00502           return rowUpper_;
00503      }
00504      inline const double* getRowUpper() const     {
00505           return rowUpper_;
00506      }
00507      //-------------------------------------------------------------------------
00511      void setObjectiveCoefficient( int elementIndex, double elementValue );
00513      inline void setObjCoeff( int elementIndex, double elementValue ) {
00514           setObjectiveCoefficient( elementIndex, elementValue);
00515      }
00516 
00519      void setColumnLower( int elementIndex, double elementValue );
00520 
00523      void setColumnUpper( int elementIndex, double elementValue );
00524 
00526      void setColumnBounds( int elementIndex,
00527                            double lower, double upper );
00528 
00537      void setColumnSetBounds(const int* indexFirst,
00538                              const int* indexLast,
00539                              const double* boundList);
00540 
00543      inline void setColLower( int elementIndex, double elementValue ) {
00544           setColumnLower(elementIndex, elementValue);
00545      }
00548      inline void setColUpper( int elementIndex, double elementValue ) {
00549           setColumnUpper(elementIndex, elementValue);
00550      }
00551 
00553      inline void setColBounds( int elementIndex,
00554                                double lower, double upper ) {
00555           setColumnBounds(elementIndex, lower, upper);
00556      }
00557 
00564      inline void setColSetBounds(const int* indexFirst,
00565                                  const int* indexLast,
00566                                  const double* boundList) {
00567           setColumnSetBounds(indexFirst, indexLast, boundList);
00568      }
00569 
00572      void setRowLower( int elementIndex, double elementValue );
00573 
00576      void setRowUpper( int elementIndex, double elementValue ) ;
00577 
00579      void setRowBounds( int elementIndex,
00580                         double lower, double upper ) ;
00581 
00588      void setRowSetBounds(const int* indexFirst,
00589                           const int* indexLast,
00590                           const double* boundList);
00591 
00593 
00594      inline const double * rowScale() const {
00595           return rowScale_;
00596      }
00597      inline const double * columnScale() const {
00598           return columnScale_;
00599      }
00600      inline const double * inverseRowScale() const {
00601           return inverseRowScale_;
00602      }
00603      inline const double * inverseColumnScale() const {
00604           return inverseColumnScale_;
00605      }
00606      inline double * mutableRowScale() const {
00607           return rowScale_;
00608      }
00609      inline double * mutableColumnScale() const {
00610           return columnScale_;
00611      }
00612      inline double * mutableInverseRowScale() const {
00613           return inverseRowScale_;
00614      }
00615      inline double * mutableInverseColumnScale() const {
00616           return inverseColumnScale_;
00617      }
00618      inline double * swapRowScale(double * newScale) {
00619           double * oldScale = rowScale_;
00620           rowScale_ = newScale;
00621           return oldScale;
00622      }
00623      void setRowScale(double * scale) ;
00624      void setColumnScale(double * scale);
00626      inline double objectiveScale() const {
00627           return objectiveScale_;
00628      }
00629      inline void setObjectiveScale(double value) {
00630           objectiveScale_ = value;
00631      }
00633      inline double rhsScale() const {
00634           return rhsScale_;
00635      }
00636      inline void setRhsScale(double value) {
00637           rhsScale_ = value;
00638      }
00640      void scaling(int mode = 1);
00643      void unscale();
00645      inline int scalingFlag() const {
00646           return scalingFlag_;
00647      }
00649      inline double * objective() const {
00650           if (objective_) {
00651                double offset;
00652                return objective_->gradient(NULL, NULL, offset, false);
00653           } else {
00654                return NULL;
00655           }
00656      }
00657      inline double * objective(const double * solution, double & offset, bool refresh = true) const {
00658           offset = 0.0;
00659           if (objective_) {
00660                return objective_->gradient(NULL, solution, offset, refresh);
00661           } else {
00662                return NULL;
00663           }
00664      }
00665      inline const double * getObjCoefficients() const {
00666           if (objective_) {
00667                double offset;
00668                return objective_->gradient(NULL, NULL, offset, false);
00669           } else {
00670                return NULL;
00671           }
00672      }
00674      inline double * rowObjective() const         {
00675           return rowObjective_;
00676      }
00677      inline const double * getRowObjCoefficients() const {
00678           return rowObjective_;
00679      }
00681      inline double * columnLower() const          {
00682           return columnLower_;
00683      }
00684      inline const double * getColLower() const    {
00685           return columnLower_;
00686      }
00688      inline double * columnUpper() const          {
00689           return columnUpper_;
00690      }
00691      inline const double * getColUpper() const    {
00692           return columnUpper_;
00693      }
00695      inline CoinPackedMatrix * matrix() const {
00696           if ( matrix_ == NULL ) return NULL;
00697           else return matrix_->getPackedMatrix();
00698      }
00700      inline int getNumElements() const {
00701           return matrix_->getNumElements();
00702      }
00705      inline double getSmallElementValue() const {
00706           return smallElement_;
00707      }
00708      inline void setSmallElementValue(double value) {
00709           smallElement_ = value;
00710      }
00712      inline ClpMatrixBase * rowCopy() const       {
00713           return rowCopy_;
00714      }
00716      void setNewRowCopy(ClpMatrixBase * newCopy);
00718      inline ClpMatrixBase * clpMatrix() const     {
00719           return matrix_;
00720      }
00722      inline ClpPackedMatrix * clpScaledMatrix() const     {
00723           return scaledMatrix_;
00724      }
00726      inline void setClpScaledMatrix(ClpPackedMatrix * scaledMatrix) {
00727           delete scaledMatrix_;
00728           scaledMatrix_ = scaledMatrix;
00729      }
00731      inline ClpPackedMatrix * swapScaledMatrix(ClpPackedMatrix * scaledMatrix) {
00732           ClpPackedMatrix * oldMatrix = scaledMatrix_;
00733           scaledMatrix_ = scaledMatrix;
00734           return oldMatrix;
00735      }
00741      void replaceMatrix(ClpMatrixBase * matrix, bool deleteCurrent = false);
00747      inline void replaceMatrix(CoinPackedMatrix * newmatrix,
00748                                bool deleteCurrent = false) {
00749           replaceMatrix(new ClpPackedMatrix(newmatrix), deleteCurrent);
00750      }
00752      inline double objectiveValue() const {
00753           return objectiveValue_ * optimizationDirection_ - dblParam_[ClpObjOffset];
00754      }
00755      inline void setObjectiveValue(double value) {
00756           objectiveValue_ = (value + dblParam_[ClpObjOffset]) / optimizationDirection_;
00757      }
00758      inline double getObjValue() const {
00759           return objectiveValue_ * optimizationDirection_ - dblParam_[ClpObjOffset];
00760      }
00762      inline char * integerInformation() const     {
00763           return integerType_;
00764      }
00767      double * infeasibilityRay(bool fullRay=false) const;
00768      double * unboundedRay() const;
00770      inline double * ray() const
00771      { return ray_;}
00773      inline bool rayExists() const {
00774          return (ray_!=NULL);
00775      }
00777      inline void deleteRay() {
00778          delete [] ray_;
00779          ray_=NULL;
00780      }
00782          inline const double * internalRay() const {
00783                  return ray_;
00784          }
00786      inline bool statusExists() const {
00787           return (status_ != NULL);
00788      }
00790      inline unsigned char *  statusArray() const {
00791           return status_;
00792      }
00795      unsigned char *  statusCopy() const;
00797      void copyinStatus(const unsigned char * statusArray);
00798 
00800      inline void setUserPointer (void * pointer) {
00801           userPointer_ = pointer;
00802      }
00803      inline void * getUserPointer () const {
00804           return userPointer_;
00805      }
00807      inline void setTrustedUserPointer (ClpTrustedData * pointer) {
00808           trustedUserPointer_ = pointer;
00809      }
00810      inline ClpTrustedData * getTrustedUserPointer () const {
00811           return trustedUserPointer_;
00812      }
00814      inline int whatsChanged() const {
00815           return whatsChanged_;
00816      }
00817      inline void setWhatsChanged(int value) {
00818           whatsChanged_ = value;
00819      }
00821      inline int numberThreads() const {
00822           return numberThreads_;
00823      }
00824      inline void setNumberThreads(int value) {
00825           numberThreads_ = value;
00826      }
00828 
00830 
00831      void passInMessageHandler(CoinMessageHandler * handler);
00833      CoinMessageHandler * pushMessageHandler(CoinMessageHandler * handler,
00834                                              bool & oldDefault);
00836      void popMessageHandler(CoinMessageHandler * oldHandler, bool oldDefault);
00838      void newLanguage(CoinMessages::Language language);
00839      inline void setLanguage(CoinMessages::Language language) {
00840           newLanguage(language);
00841      }
00843      void setDefaultMessageHandler();
00845      inline CoinMessageHandler * messageHandler() const       {
00846           return handler_;
00847      }
00849      inline CoinMessages messages() const                     {
00850           return messages_;
00851      }
00853      inline CoinMessages * messagesPointer()                  {
00854           return & messages_;
00855      }
00857      inline CoinMessages coinMessages() const                  {
00858           return coinMessages_;
00859      }
00861      inline CoinMessages * coinMessagesPointer()                  {
00862           return & coinMessages_;
00863      }
00872      inline void setLogLevel(int value)    {
00873           handler_->setLogLevel(value);
00874      }
00875      inline int logLevel() const           {
00876           return handler_->logLevel();
00877      }
00879      inline bool defaultHandler() const {
00880           return defaultHandler_;
00881      }
00883      void passInEventHandler(const ClpEventHandler * eventHandler);
00885      inline ClpEventHandler * eventHandler() const {
00886           return eventHandler_;
00887      }
00889      inline CoinThreadRandom * randomNumberGenerator() {
00890           return &randomNumberGenerator_;
00891      }
00893      inline CoinThreadRandom & mutableRandomNumberGenerator() {
00894           return randomNumberGenerator_;
00895      }
00897      inline void setRandomSeed(int value) {
00898           randomNumberGenerator_.setSeed(value);
00899      }
00901      inline int lengthNames() const {
00902           return lengthNames_;
00903      }
00904 #ifndef CLP_NO_STD
00906      inline void setLengthNames(int value) {
00907           lengthNames_ = value;
00908      }
00910      inline const std::vector<std::string> * rowNames() const {
00911           return &rowNames_;
00912      }
00913      inline const std::string& rowName(int iRow) const {
00914           return rowNames_[iRow];
00915      }
00917      std::string getRowName(int iRow) const;
00919      inline const std::vector<std::string> * columnNames() const {
00920           return &columnNames_;
00921      }
00922      inline const std::string& columnName(int iColumn) const {
00923           return columnNames_[iColumn];
00924      }
00926      std::string getColumnName(int iColumn) const;
00927 #endif
00929      inline ClpObjective * objectiveAsObject() const {
00930           return objective_;
00931      }
00932      void setObjective(ClpObjective * objective);
00933      inline void setObjectivePointer(ClpObjective * newobjective) {
00934           objective_ = newobjective;
00935      }
00938      int emptyProblem(int * infeasNumber = NULL, double * infeasSum = NULL, bool printMessage = true);
00939 
00941 
00950      void times(double scalar,
00951                 const double * x, double * y) const;
00955      void transposeTimes(double scalar,
00956                          const double * x, double * y) const ;
00958 
00959 
00960      //---------------------------------------------------------------------------
00978 
00979      bool setIntParam(ClpIntParam key, int value) ;
00981      bool setDblParam(ClpDblParam key, double value) ;
00982 #ifndef CLP_NO_STD
00984      bool setStrParam(ClpStrParam key, const std::string & value);
00985 #endif
00986      // Get an integer parameter
00987      inline bool getIntParam(ClpIntParam key, int& value) const {
00988           if (key < ClpLastIntParam) {
00989                value = intParam_[key];
00990                return true;
00991           } else {
00992                return false;
00993           }
00994      }
00995      // Get an double parameter
00996      inline bool getDblParam(ClpDblParam key, double& value) const {
00997           if (key < ClpLastDblParam) {
00998                value = dblParam_[key];
00999                return true;
01000           } else {
01001                return false;
01002           }
01003      }
01004 #ifndef CLP_NO_STD
01005      // Get a string parameter
01006      inline bool getStrParam(ClpStrParam key, std::string& value) const {
01007           if (key < ClpLastStrParam) {
01008                value = strParam_[key];
01009                return true;
01010           } else {
01011                return false;
01012           }
01013      }
01014 #endif
01016      void generateCpp( FILE * fp);
01017 
01050      inline unsigned int specialOptions() const {
01051           return specialOptions_;
01052      }
01053      void setSpecialOptions(unsigned int value);
01054 #define COIN_CBC_USING_CLP 0x01000000
01055      inline bool inCbcBranchAndBound() const {
01056           return (specialOptions_ & COIN_CBC_USING_CLP) != 0;
01057      }
01059 
01062 protected:
01064      void gutsOfDelete(int type);
01068      void gutsOfCopy(const ClpModel & rhs, int trueCopy = 1);
01070      void getRowBound(int iRow, double& lower, double& upper) const;
01072      void gutsOfLoadModel ( int numberRows, int numberColumns,
01073                             const double* collb, const double* colub,
01074                             const double* obj,
01075                             const double* rowlb, const double* rowub,
01076                             const double * rowObjective = NULL);
01078      void gutsOfScaling();
01080      inline double rawObjectiveValue() const {
01081           return objectiveValue_;
01082      }
01084      inline bool permanentArrays() const {
01085           return (specialOptions_ & 65536) != 0;
01086      }
01088      void startPermanentArrays();
01090      void stopPermanentArrays();
01092      const char * const * rowNamesAsChar() const;
01094      const char * const * columnNamesAsChar() const;
01096      void deleteNamesAsChar(const char * const * names, int number) const;
01098      void onStopped();
01100 
01101 
01103 protected:
01104 
01107 
01108      double optimizationDirection_;
01110      double dblParam_[ClpLastDblParam];
01112      double objectiveValue_;
01114      double smallElement_;
01116      double objectiveScale_;
01118      double rhsScale_;
01120      int numberRows_;
01122      int numberColumns_;
01124      double * rowActivity_;
01126      double * columnActivity_;
01128      double * dual_;
01130      double * reducedCost_;
01132      double* rowLower_;
01134      double* rowUpper_;
01136      ClpObjective * objective_;
01138      double * rowObjective_;
01140      double * columnLower_;
01142      double * columnUpper_;
01144      ClpMatrixBase * matrix_;
01146      ClpMatrixBase * rowCopy_;
01148      ClpPackedMatrix * scaledMatrix_;
01150      double * ray_;
01152      double * rowScale_;
01154      double * columnScale_;
01156      double * inverseRowScale_;
01158      double * inverseColumnScale_;
01161      int scalingFlag_;
01169      unsigned char * status_;
01171      char * integerType_;
01173      void * userPointer_;
01175      ClpTrustedData * trustedUserPointer_;
01177      int intParam_[ClpLastIntParam];
01179      int numberIterations_;
01181      int solveType_;
01198 #define ROW_COLUMN_COUNTS_SAME 1
01199 #define MATRIX_SAME 2
01200 #define MATRIX_JUST_ROWS_ADDED 4
01201 #define MATRIX_JUST_COLUMNS_ADDED 8
01202 #define ROW_LOWER_SAME 16
01203 #define ROW_UPPER_SAME 32
01204 #define OBJECTIVE_SAME 64 
01205 #define COLUMN_LOWER_SAME 128
01206 #define COLUMN_UPPER_SAME 256
01207 #define BASIS_SAME 512
01208 #define ALL_SAME 65339
01209 #define ALL_SAME_EXCEPT_COLUMN_BOUNDS 65337
01210      unsigned int whatsChanged_;
01212      int problemStatus_;
01214      int secondaryStatus_;
01216      int lengthNames_;
01218      int numberThreads_;
01222      unsigned int specialOptions_;
01224      CoinMessageHandler * handler_;
01226      bool defaultHandler_;
01228      CoinThreadRandom randomNumberGenerator_;
01230      ClpEventHandler * eventHandler_;
01231 #ifndef CLP_NO_STD
01233      std::vector<std::string> rowNames_;
01235      std::vector<std::string> columnNames_;
01236 #endif
01238      CoinMessages messages_;
01240      CoinMessages coinMessages_;
01242      int maximumColumns_;
01244      int maximumRows_;
01246      int maximumInternalColumns_;
01248      int maximumInternalRows_;
01250      CoinPackedMatrix baseMatrix_;
01252      CoinPackedMatrix baseRowCopy_;
01254      double * savedRowScale_;
01256      double * savedColumnScale_;
01257 #ifndef CLP_NO_STD
01259      std::string strParam_[ClpLastStrParam];
01260 #endif
01261 
01262 };
01265 class ClpDataSave {
01266 
01267 public:
01271 
01272      ClpDataSave (  );
01273 
01275      ClpDataSave(const ClpDataSave &);
01277      ClpDataSave & operator=(const ClpDataSave & rhs);
01279      ~ClpDataSave (  );
01280 
01282 
01284 public:
01285 
01288      double dualBound_;
01289      double infeasibilityCost_;
01290      double pivotTolerance_;
01291      double zeroFactorizationTolerance_;
01292      double zeroSimplexTolerance_;
01293      double acceptablePivot_;
01294      double objectiveScale_;
01295      int sparseThreshold_;
01296      int perturbation_;
01297      int forceFactorization_;
01298      int scalingFlag_;
01299      unsigned int specialOptions_;
01301 };
01302 
01303 #endif
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines

Generated on 11 Apr 2013 for Clp by  doxygen 1.6.1