This is for Simplex stuff which is neither dual nor primal. More...
#include <ClpSimplexOther.hpp>
Methods | |
void | dualRanging (int numberCheck, const int *which, double *costIncrease, int *sequenceIncrease, double *costDecrease, int *sequenceDecrease, double *valueIncrease=NULL, double *valueDecrease=NULL) |
Dual ranging. | |
void | primalRanging (int numberCheck, const int *which, double *valueIncrease, int *sequenceIncrease, double *valueDecrease, int *sequenceDecrease) |
Primal ranging. | |
int | parametrics (double startingTheta, double &endingTheta, double reportIncrement, const double *changeLowerBound, const double *changeUpperBound, const double *changeLowerRhs, const double *changeUpperRhs, const double *changeObjective) |
Parametrics This is an initial slow version. | |
int | writeBasis (const char *filename, bool writeValues=false, int formatType=0) const |
Write the basis in MPS format to the specified file. | |
int | readBasis (const char *filename) |
Read a basis from the given filename. | |
ClpSimplex * | dualOfModel (double fractionRowRanges=1.0, double fractionColumnRanges=1.0) const |
Creates dual of a problem if looks plausible (defaults will always create model) fractionRowRanges is fraction of rows allowed to have ranges fractionColumnRanges is fraction of columns allowed to have ranges. | |
int | restoreFromDual (const ClpSimplex *dualProblem) |
Restores solution from dualized problem non-zero return code indicates minor problems. | |
ClpSimplex * | crunch (double *rhs, int *whichRows, int *whichColumns, int &nBound, bool moreBounds=false, bool tightenBounds=false) |
Does very cursory presolve. | |
void | afterCrunch (const ClpSimplex &small, const int *whichRows, const int *whichColumns, int nBound) |
After very cursory presolve. | |
void | cleanupAfterPostsolve () |
Quick try at cleaning up duals if postsolve gets wrong. | |
int | tightenIntegerBounds (double *rhsSpace) |
Tightens integer bounds - returns number tightened or -1 if infeasible. | |
int | expandKnapsack (int knapsackRow, int &numberOutput, double *buildObj, CoinBigIndex *buildStart, int *buildRow, double *buildElement, int reConstruct=-1) const |
Expands out all possible combinations for a knapsack If buildObj NULL then just computes space needed - returns number elements On entry numberOutput is maximum allowed, on exit it is number needed or -1 (as will be number elements) if maximum exceeded. | |
int | parametricsLoop (double startingTheta, double &endingTheta, double reportIncrement, const double *changeLower, const double *changeUpper, const double *changeObjective, ClpDataSave &data, bool canTryQuick) |
Parametrics - inner loop This first attempt is when reportIncrement non zero and may not report endingTheta correctly If it can not reach input endingTheta return code will be 1 for infeasible, 2 for unbounded, otherwise 0. | |
void | statusOfProblemInParametrics (int type, ClpDataSave &saveData) |
Refactorizes if necessary Checks if finished. | |
int | whileIterating (double startingTheta, double &endingTheta, double reportIncrement, const double *changeLower, const double *changeUpper, const double *changeObjective) |
This has the flow between re-factorizations. | |
int | nextTheta (int type, double maxTheta, double *primalChange, double *dualChange, const double *changeLower, const double *changeUpper, const double *changeObjective) |
Computes next theta and says if objective or bounds (0= bounds, 1 objective, -1 none). | |
void | checkDualRatios (CoinIndexedVector *rowArray, CoinIndexedVector *columnArray, double &costIncrease, int &sequenceIncrease, double &alphaIncrease, double &costDecrease, int &sequenceDecrease, double &alphaDecrease) |
Row array has row part of pivot row Column array has column part. | |
void | checkPrimalRatios (CoinIndexedVector *rowArray, int direction) |
Row array has pivot column This is used in primal ranging. | |
double | primalRanging1 (int whichIn, int whichOther) |
Returns new value of whichOther when whichIn enters basis. |
Additional Inherited Members | |
![]() | |
enum | Status { isFree = 0x00, basic = 0x01, atUpperBound = 0x02, atLowerBound = 0x03, superBasic = 0x04, isFixed = 0x05 } |
enums for status of various sorts. More... | |
enum | FakeBound { noFake = 0x00, lowerFake = 0x01, upperFake = 0x02, bothFake = 0x03 } |
![]() | |
ClpSimplex (bool emptyMessages=false) | |
Default constructor. | |
ClpSimplex (const ClpSimplex &rhs, int scalingMode=-1) | |
Copy constructor. | |
ClpSimplex (const ClpModel &rhs, int scalingMode=-1) | |
Copy constructor from model. | |
ClpSimplex (const ClpModel *wholeModel, int numberRows, const int *whichRows, int numberColumns, const int *whichColumns, bool dropNames=true, bool dropIntegers=true, bool fixOthers=false) | |
Subproblem constructor. | |
ClpSimplex (const ClpSimplex *wholeModel, int numberRows, const int *whichRows, int numberColumns, const int *whichColumns, bool dropNames=true, bool dropIntegers=true, bool fixOthers=false) | |
Subproblem constructor. | |
ClpSimplex (ClpSimplex *wholeModel, int numberColumns, const int *whichColumns) | |
This constructor modifies original ClpSimplex and stores original stuff in created ClpSimplex. | |
void | originalModel (ClpSimplex *miniModel) |
This copies back stuff from miniModel and then deletes miniModel. | |
void | setPersistenceFlag (int value) |
Array persistence flag If 0 then as now (delete/new) 1 then only do arrays if bigger needed 2 as 1 but give a bit extra if bigger needed. | |
void | makeBaseModel () |
Save a copy of model with certain state - normally without cuts. | |
void | deleteBaseModel () |
Switch off base model. | |
ClpSimplex * | baseModel () const |
See if we have base model. | |
void | setToBaseModel (ClpSimplex *model=NULL) |
Reset to base model (just size and arrays needed) If model NULL use internal copy. | |
ClpSimplex & | operator= (const ClpSimplex &rhs) |
Assignment operator. This copies the data. | |
~ClpSimplex () | |
Destructor. | |
void | loadProblem (const ClpMatrixBase &matrix, const double *collb, const double *colub, const double *obj, const double *rowlb, const double *rowub, const double *rowObjective=NULL) |
Loads a problem (the constraints on the rows are given by lower and upper bounds). | |
void | loadProblem (const CoinPackedMatrix &matrix, const double *collb, const double *colub, const double *obj, const double *rowlb, const double *rowub, const double *rowObjective=NULL) |
Default constructor. | |
void | loadProblem (const int numcols, const int numrows, const CoinBigIndex *start, const int *index, const double *value, const double *collb, const double *colub, const double *obj, const double *rowlb, const double *rowub, const double *rowObjective=NULL) |
Just like the other loadProblem() method except that the matrix is given in a standard column major ordered format (without gaps). | |
void | loadProblem (const int numcols, const int numrows, const CoinBigIndex *start, const int *index, const double *value, const int *length, const double *collb, const double *colub, const double *obj, const double *rowlb, const double *rowub, const double *rowObjective=NULL) |
This one is for after presolve to save memory. | |
int | loadProblem (CoinModel &modelObject, bool keepSolution=false) |
This loads a model from a coinModel object - returns number of errors. | |
int | readMps (const char *filename, bool keepNames=false, bool ignoreErrors=false) |
Read an mps file from the given filename. | |
int | readGMPL (const char *filename, const char *dataName, bool keepNames=false) |
Read GMPL files from the given filenames. | |
int | readLp (const char *filename, const double epsilon=1e-5) |
Read file in LP format from file with name filename. | |
void | borrowModel (ClpModel &otherModel) |
Borrow model. | |
void | borrowModel (ClpSimplex &otherModel) |
Default constructor. | |
void | passInEventHandler (const ClpEventHandler *eventHandler) |
Pass in Event handler (cloned and deleted at end) | |
void | getbackSolution (const ClpSimplex &smallModel, const int *whichRow, const int *whichColumn) |
Puts solution back into small model. | |
int | loadNonLinear (void *info, int &numberConstraints, ClpConstraint **&constraints) |
Load nonlinear part of problem from AMPL info Returns 0 if linear 1 if quadratic objective 2 if quadratic constraints 3 if nonlinear objective 4 if nonlinear constraints -1 on failure. | |
int | initialSolve (ClpSolve &options) |
General solve algorithm which can do presolve. | |
int | initialSolve () |
Default initial solve. | |
int | initialDualSolve () |
Dual initial solve. | |
int | initialPrimalSolve () |
Primal initial solve. | |
int | initialBarrierSolve () |
Barrier initial solve. | |
int | initialBarrierNoCrossSolve () |
Barrier initial solve, not to be followed by crossover. | |
int | dual (int ifValuesPass=0, int startFinishOptions=0) |
Dual algorithm - see ClpSimplexDual.hpp for method. | |
int | dualDebug (int ifValuesPass=0, int startFinishOptions=0) |
General solve algorithm which can do presolve. | |
int | primal (int ifValuesPass=0, int startFinishOptions=0) |
Primal algorithm - see ClpSimplexPrimal.hpp for method. | |
int | nonlinearSLP (int numberPasses, double deltaTolerance) |
Solves nonlinear problem using SLP - may be used as crash for other algorithms when number of iterations small. | |
int | nonlinearSLP (int numberConstraints, ClpConstraint **constraints, int numberPasses, double deltaTolerance) |
Solves problem with nonlinear constraints using SLP - may be used as crash for other algorithms when number of iterations small. | |
int | barrier (bool crossover=true) |
Solves using barrier (assumes you have good cholesky factor code). | |
int | reducedGradient (int phase=0) |
Solves non-linear using reduced gradient. | |
int | solve (CoinStructuredModel *model) |
Solve using structure of model and maybe in parallel. | |
int | loadProblem (CoinStructuredModel &modelObject, bool originalOrder=true, bool keepSolution=false) |
This loads a model from a CoinStructuredModel object - returns number of errors. | |
int | cleanup (int cleanupScaling) |
When scaling is on it is possible that the scaled problem is feasible but the unscaled is not. | |
CoinWarmStartBasis * | getBasis () const |
Returns a basis (to be deleted by user) | |
void | setFactorization (ClpFactorization &factorization) |
Passes in factorization. | |
ClpFactorization * | swapFactorization (ClpFactorization *factorization) |
General solve algorithm which can do presolve. | |
void | copyFactorization (ClpFactorization &factorization) |
Copies in factorization to existing one. | |
int | tightenPrimalBounds (double factor=0.0, int doTight=0, bool tightIntegers=false) |
Tightens primal bounds to make dual faster. | |
int | crash (double gap, int pivot) |
Crash - at present just aimed at dual, returns -2 if dual preferred and crash basis created -1 if dual preferred and all slack basis preferred 0 if basis going in was not all slack 1 if primal preferred and all slack basis preferred 2 if primal preferred and crash basis created. | |
void | setDualRowPivotAlgorithm (ClpDualRowPivot &choice) |
Sets row pivot choice algorithm in dual. | |
void | setPrimalColumnPivotAlgorithm (ClpPrimalColumnPivot &choice) |
Sets column pivot choice algorithm in primal. | |
int | strongBranching (int numberVariables, const int *variables, double *newLower, double *newUpper, double **outputSolution, int *outputStatus, int *outputIterations, bool stopOnFirstInfeasible=true, bool alwaysFinish=false, int startFinishOptions=0) |
For strong branching. | |
int | fathom (void *stuff) |
Fathom - 1 if solution. | |
int | fathomMany (void *stuff) |
Do up to N deep - returns -1 - no solution nNodes_ valid nodes >= if solution and that node gives solution ClpNode array is 2**N long. | |
double | doubleCheck () |
Double checks OK. | |
int | startFastDual2 (ClpNodeStuff *stuff) |
Starts Fast dual2. | |
int | fastDual2 (ClpNodeStuff *stuff) |
Like Fast dual. | |
void | stopFastDual2 (ClpNodeStuff *stuff) |
Stops Fast dual2. | |
ClpSimplex * | fastCrunch (ClpNodeStuff *stuff, int mode) |
Deals with crunch aspects mode 0 - in 1 - out with solution 2 - out without solution returns small model or NULL. | |
int | pivot () |
Pivot in a variable and out a variable. | |
int | primalPivotResult () |
Pivot in a variable and choose an outgoing one. | |
int | dualPivotResult () |
Pivot out a variable and choose an incoing one. | |
int | startup (int ifValuesPass, int startFinishOptions=0) |
Common bits of coding for dual and primal. | |
void | finish (int startFinishOptions=0) |
Pivot in a variable and out a variable. | |
bool | statusOfProblem (bool initial=false) |
Factorizes and returns true if optimal. | |
void | defaultFactorizationFrequency () |
If user left factorization frequency then compute. | |
bool | primalFeasible () const |
If problem is primal feasible. | |
bool | dualFeasible () const |
If problem is dual feasible. | |
ClpFactorization * | factorization () const |
factorization | |
bool | sparseFactorization () const |
Sparsity on or off. | |
void | setSparseFactorization (bool value) |
If problem is primal feasible. | |
int | factorizationFrequency () const |
Factorization frequency. | |
void | setFactorizationFrequency (int value) |
If problem is primal feasible. | |
double | dualBound () const |
Dual bound. | |
void | setDualBound (double value) |
If problem is primal feasible. | |
double | infeasibilityCost () const |
Infeasibility cost. | |
void | setInfeasibilityCost (double value) |
If problem is primal feasible. | |
int | perturbation () const |
Amount of print out: 0 - none 1 - just final 2 - just factorizations 3 - as 2 plus a bit more 4 - verbose above that 8,16,32 etc just for selective debug. | |
void | setPerturbation (int value) |
If problem is primal feasible. | |
int | algorithm () const |
Current (or last) algorithm. | |
void | setAlgorithm (int value) |
Set algorithm. | |
bool | isObjectiveLimitTestValid () const |
Return true if the objective limit test can be relied upon. | |
double | sumDualInfeasibilities () const |
Sum of dual infeasibilities. | |
void | setSumDualInfeasibilities (double value) |
If problem is primal feasible. | |
double | sumOfRelaxedDualInfeasibilities () const |
Sum of relaxed dual infeasibilities. | |
void | setSumOfRelaxedDualInfeasibilities (double value) |
If problem is primal feasible. | |
int | numberDualInfeasibilities () const |
Number of dual infeasibilities. | |
void | setNumberDualInfeasibilities (int value) |
If problem is primal feasible. | |
int | numberDualInfeasibilitiesWithoutFree () const |
Number of dual infeasibilities (without free) | |
double | sumPrimalInfeasibilities () const |
Sum of primal infeasibilities. | |
void | setSumPrimalInfeasibilities (double value) |
If problem is primal feasible. | |
double | sumOfRelaxedPrimalInfeasibilities () const |
Sum of relaxed primal infeasibilities. | |
void | setSumOfRelaxedPrimalInfeasibilities (double value) |
If problem is primal feasible. | |
int | numberPrimalInfeasibilities () const |
Number of primal infeasibilities. | |
void | setNumberPrimalInfeasibilities (int value) |
If problem is primal feasible. | |
int | saveModel (const char *fileName) |
Save model to file, returns 0 if success. | |
int | restoreModel (const char *fileName) |
Restore model from file, returns 0 if success, deletes current model. | |
void | checkSolution (int setToBounds=0) |
Just check solution (for external use) - sets sum of infeasibilities etc. | |
void | checkSolutionInternal () |
Just check solution (for internal use) - sets sum of infeasibilities etc. | |
CoinIndexedVector * | rowArray (int index) const |
Useful row length arrays (0,1,2,3,4,5) | |
CoinIndexedVector * | columnArray (int index) const |
Useful column length arrays (0,1,2,3,4,5) | |
double | alphaAccuracy () const |
Initial value for alpha accuracy calculation (-1.0 off) | |
void | setAlphaAccuracy (double value) |
If problem is primal feasible. | |
void | setDisasterHandler (ClpDisasterHandler *handler) |
Objective value. | |
ClpDisasterHandler * | disasterHandler () const |
Get disaster handler. | |
double | largeValue () const |
Large bound value (for complementarity etc) | |
void | setLargeValue (double value) |
If problem is primal feasible. | |
double | largestPrimalError () const |
Largest error on Ax-b. | |
double | largestDualError () const |
Largest error on basic duals. | |
void | setLargestPrimalError (double value) |
Largest error on Ax-b. | |
void | setLargestDualError (double value) |
Largest error on basic duals. | |
double | zeroTolerance () const |
Get zero tolerance. | |
void | setZeroTolerance (double value) |
Set zero tolerance. | |
int * | pivotVariable () const |
Basic variables pivoting on which rows. | |
bool | automaticScaling () const |
If automatic scaling on. | |
void | setAutomaticScaling (bool onOff) |
If problem is primal feasible. | |
double | currentDualTolerance () const |
Current dual tolerance. | |
void | setCurrentDualTolerance (double value) |
If problem is primal feasible. | |
double | currentPrimalTolerance () const |
Current primal tolerance. | |
void | setCurrentPrimalTolerance (double value) |
If problem is primal feasible. | |
int | numberRefinements () const |
How many iterative refinements to do. | |
void | setNumberRefinements (int value) |
If problem is primal feasible. | |
double | alpha () const |
Alpha (pivot element) for use by classes e.g. steepestedge. | |
void | setAlpha (double value) |
If problem is primal feasible. | |
double | dualIn () const |
Reduced cost of last incoming for use by classes e.g. steepestedge. | |
int | pivotRow () const |
Pivot Row for use by classes e.g. steepestedge. | |
void | setPivotRow (int value) |
If problem is primal feasible. | |
double | valueIncomingDual () const |
value of incoming variable (in Dual) | |
double * | solutionRegion (int section) const |
Return row or column sections - not as much needed as it once was. | |
double * | djRegion (int section) const |
Return row or column sections - not as much needed as it once was. | |
double * | lowerRegion (int section) const |
Return row or column sections - not as much needed as it once was. | |
double * | upperRegion (int section) const |
Return row or column sections - not as much needed as it once was. | |
double * | costRegion (int section) const |
Return row or column sections - not as much needed as it once was. | |
double * | solutionRegion () const |
Return region as single array. | |
double * | djRegion () const |
Return row or column sections - not as much needed as it once was. | |
double * | lowerRegion () const |
Return row or column sections - not as much needed as it once was. | |
double * | upperRegion () const |
Return row or column sections - not as much needed as it once was. | |
double * | costRegion () const |
Return row or column sections - not as much needed as it once was. | |
Status | getStatus (int sequence) const |
Return row or column sections - not as much needed as it once was. | |
void | setStatus (int sequence, Status newstatus) |
Return row or column sections - not as much needed as it once was. | |
bool | startPermanentArrays () |
Start or reset using maximumRows_ and Columns_ - true if change. | |
void | setInitialDenseFactorization (bool onOff) |
Normally the first factorization does sparse coding because the factorization could be singular. | |
bool | initialDenseFactorization () const |
Return row or column sections - not as much needed as it once was. | |
int | sequenceIn () const |
Return sequence In or Out. | |
int | sequenceOut () const |
Return row or column sections - not as much needed as it once was. | |
void | setSequenceIn (int sequence) |
Set sequenceIn or Out. | |
void | setSequenceOut (int sequence) |
Return row or column sections - not as much needed as it once was. | |
int | directionIn () const |
Return direction In or Out. | |
int | directionOut () const |
Return row or column sections - not as much needed as it once was. | |
void | setDirectionIn (int direction) |
Set directionIn or Out. | |
void | setDirectionOut (int direction) |
Return row or column sections - not as much needed as it once was. | |
double | valueOut () const |
Value of Out variable. | |
void | setValueOut (double value) |
Set value of out variable. | |
void | setLowerOut (double value) |
Set lower of out variable. | |
void | setUpperOut (double value) |
Set upper of out variable. | |
void | setTheta (double value) |
Set theta of out variable. | |
int | isColumn (int sequence) const |
Returns 1 if sequence indicates column. | |
int | sequenceWithin (int sequence) const |
Returns sequence number within section. | |
double | solution (int sequence) |
Return row or column values. | |
double & | solutionAddress (int sequence) |
Return address of row or column values. | |
double | reducedCost (int sequence) |
Return row or column sections - not as much needed as it once was. | |
double & | reducedCostAddress (int sequence) |
Return row or column sections - not as much needed as it once was. | |
double | lower (int sequence) |
Return row or column sections - not as much needed as it once was. | |
double & | lowerAddress (int sequence) |
Return address of row or column lower bound. | |
double | upper (int sequence) |
Return row or column sections - not as much needed as it once was. | |
double & | upperAddress (int sequence) |
Return address of row or column upper bound. | |
double | cost (int sequence) |
Return row or column sections - not as much needed as it once was. | |
double & | costAddress (int sequence) |
Return address of row or column cost. | |
double | originalLower (int iSequence) const |
Return original lower bound. | |
double | originalUpper (int iSequence) const |
Return original lower bound. | |
double | theta () const |
Theta (pivot change) | |
double | bestPossibleImprovement () const |
Best possible improvement using djs (primal) or obj change by flipping bounds to make dual feasible (dual) | |
ClpNonLinearCost * | nonLinearCost () const |
Return pointer to details of costs. | |
int | moreSpecialOptions () const |
Return more special options 1 bit - if presolve says infeasible in ClpSolve return 2 bit - if presolved problem infeasible return 4 bit - keep arrays like upper_ around 8 bit - if factorization kept can still declare optimal at once 16 bit - if checking replaceColumn accuracy before updating 32 bit - say optimal if primal feasible! 64 bit - give up easily in dual (and say infeasible) 128 bit - no objective, 0-1 and in B&B 256 bit - in primal from dual or vice versa 512 bit - alternative use of solveType_. | |
void | setMoreSpecialOptions (int value) |
Set more special options 1 bit - if presolve says infeasible in ClpSolve return 2 bit - if presolved problem infeasible return 4 bit - keep arrays like upper_ around 8 bit - no free or superBasic variables 16 bit - if checking replaceColumn accuracy before updating 32 bit - say optimal if primal feasible! 64 bit - give up easily in dual (and say infeasible) 128 bit - no objective, 0-1 and in B&B 256 bit - in primal from dual or vice versa 512 bit - alternative use of solveType_. | |
void | setFakeBound (int sequence, FakeBound fakeBound) |
To flag a variable (not inline to allow for column generation) | |
FakeBound | getFakeBound (int sequence) const |
To flag a variable (not inline to allow for column generation) | |
void | setRowStatus (int sequence, Status newstatus) |
To flag a variable (not inline to allow for column generation) | |
Status | getRowStatus (int sequence) const |
To flag a variable (not inline to allow for column generation) | |
void | setColumnStatus (int sequence, Status newstatus) |
To flag a variable (not inline to allow for column generation) | |
Status | getColumnStatus (int sequence) const |
To flag a variable (not inline to allow for column generation) | |
void | setPivoted (int sequence) |
To flag a variable (not inline to allow for column generation) | |
void | clearPivoted (int sequence) |
To flag a variable (not inline to allow for column generation) | |
bool | pivoted (int sequence) const |
To flag a variable (not inline to allow for column generation) | |
void | setFlagged (int sequence) |
To flag a variable (not inline to allow for column generation) | |
void | clearFlagged (int sequence) |
To flag a variable (not inline to allow for column generation) | |
bool | flagged (int sequence) const |
To flag a variable (not inline to allow for column generation) | |
void | setActive (int iRow) |
To say row active in primal pivot row choice. | |
void | clearActive (int iRow) |
To flag a variable (not inline to allow for column generation) | |
bool | active (int iRow) const |
To flag a variable (not inline to allow for column generation) | |
void | createStatus () |
Set up status array (can be used by OsiClp). | |
void | allSlackBasis (bool resetSolution=false) |
Sets up all slack basis and resets solution to as it was after initial load or readMps. | |
int | lastBadIteration () const |
So we know when to be cautious. | |
int | progressFlag () const |
Progress flag - at present 0 bit says artificials out. | |
void | forceFactorization (int value) |
Force re-factorization early. | |
double | rawObjectiveValue () const |
Raw objective value (so always minimize in primal) | |
void | computeObjectiveValue (bool useWorkingSolution=false) |
Compute objective value from solution and put in objectiveValue_. | |
double | computeInternalObjectiveValue () |
Compute minimization objective value from internal solution without perturbation. | |
int | numberExtraRows () const |
Number of extra rows. | |
int | maximumBasic () const |
Maximum number of basic variables - can be more than number of rows if GUB. | |
int | baseIteration () const |
Iteration when we entered dual or primal. | |
void | generateCpp (FILE *fp, bool defaultFactor=false) |
Create C++ lines to get to current state. | |
ClpFactorization * | getEmptyFactorization () |
Gets clean and emptyish factorization. | |
void | setEmptyFactorization () |
May delete or may make clean and emptyish factorization. | |
void | moveInfo (const ClpSimplex &rhs, bool justStatus=false) |
Move status and solution across. | |
void | getBInvARow (int row, double *z, double *slack=NULL) |
Get a row of the tableau (slack part in slack if not NULL) | |
void | getBInvRow (int row, double *z) |
Get a row of the basis inverse. | |
void | getBInvACol (int col, double *vec) |
Get a column of the tableau. | |
void | getBInvCol (int col, double *vec) |
Get a column of the basis inverse. | |
void | getBasics (int *index) |
Get basic indices (order of indices corresponds to the order of elements in a vector retured by getBInvACol() and getBInvCol()). | |
void | setObjectiveCoefficient (int elementIndex, double elementValue) |
Set an objective function coefficient. | |
void | setObjCoeff (int elementIndex, double elementValue) |
Set an objective function coefficient. | |
void | setColumnLower (int elementIndex, double elementValue) |
Set a single column lower bound Use -DBL_MAX for -infinity. | |
void | setColumnUpper (int elementIndex, double elementValue) |
Set a single column upper bound Use DBL_MAX for infinity. | |
void | setColumnBounds (int elementIndex, double lower, double upper) |
Set a single column lower and upper bound. | |
void | setColumnSetBounds (const int *indexFirst, const int *indexLast, const double *boundList) |
Set the bounds on a number of columns simultaneously The default implementation just invokes setColLower() and setColUpper() over and over again. | |
void | setColLower (int elementIndex, double elementValue) |
Set a single column lower bound Use -DBL_MAX for -infinity. | |
void | setColUpper (int elementIndex, double elementValue) |
Set a single column upper bound Use DBL_MAX for infinity. | |
void | setColBounds (int elementIndex, double newlower, double newupper) |
Set a single column lower and upper bound. | |
void | setColSetBounds (const int *indexFirst, const int *indexLast, const double *boundList) |
Set the bounds on a number of columns simultaneously | |
void | setRowLower (int elementIndex, double elementValue) |
Set a single row lower bound Use -DBL_MAX for -infinity. | |
void | setRowUpper (int elementIndex, double elementValue) |
Set a single row upper bound Use DBL_MAX for infinity. | |
void | setRowBounds (int elementIndex, double lower, double upper) |
Set a single row lower and upper bound. | |
void | setRowSetBounds (const int *indexFirst, const int *indexLast, const double *boundList) |
Set the bounds on a number of rows simultaneously | |
void | resize (int newNumberRows, int newNumberColumns) |
Resizes rim part of model. | |
int | getSolution (const double *rowActivities, const double *columnActivities) |
Given an existing factorization computes and checks primal and dual solutions. | |
int | getSolution () |
Given an existing factorization computes and checks primal and dual solutions. | |
int | createPiecewiseLinearCosts (const int *starts, const double *lower, const double *gradient) |
Constructs a non linear cost from list of non-linearities (columns only) First lower of each column is taken as real lower Last lower is taken as real upper and cost ignored. | |
ClpDualRowPivot * | dualRowPivot () const |
dual row pivot choice | |
ClpPrimalColumnPivot * | primalColumnPivot () const |
primal column pivot choice | |
bool | goodAccuracy () const |
Returns true if model looks OK. | |
void | returnModel (ClpSimplex &otherModel) |
Return model - updates any scalars. | |
int | internalFactorize (int solveType) |
Factorizes using current basis. | |
ClpDataSave | saveData () |
Save data. | |
void | restoreData (ClpDataSave saved) |
Restore data. | |
void | cleanStatus () |
Clean up status. | |
int | factorize () |
Factorizes using current basis. For external use. | |
void | computeDuals (double *givenDjs) |
Computes duals from scratch. | |
void | computePrimals (const double *rowActivities, const double *columnActivities) |
Computes primals from scratch. | |
void | add (double *array, int column, double multiplier) const |
Adds multiple of a column into an array. | |
void | unpack (CoinIndexedVector *rowArray) const |
Unpacks one column of the matrix into indexed array Uses sequenceIn_ Also applies scaling if needed. | |
void | unpack (CoinIndexedVector *rowArray, int sequence) const |
Unpacks one column of the matrix into indexed array Slack if sequence>= numberColumns Also applies scaling if needed. | |
void | unpackPacked (CoinIndexedVector *rowArray) |
Unpacks one column of the matrix into indexed array as packed vector Uses sequenceIn_ Also applies scaling if needed. | |
void | unpackPacked (CoinIndexedVector *rowArray, int sequence) |
Unpacks one column of the matrix into indexed array as packed vector Slack if sequence>= numberColumns Also applies scaling if needed. | |
void | setValuesPassAction (double incomingInfeasibility, double allowedInfeasibility) |
For advanced use. | |
![]() | |
int | spareIntArray_ [4] |
Spare int array for passing information [0]!=0 switches on. | |
double | spareDoubleArray_ [4] |
Spare double array for passing information [0]!=0 switches on. | |
![]() | |
int | gutsOfSolution (double *givenDuals, const double *givenPrimals, bool valuesPass=false) |
May change basis and then returns number changed. | |
void | gutsOfDelete (int type) |
Does most of deletion (0 = all, 1 = most, 2 most + factorization) | |
void | gutsOfCopy (const ClpSimplex &rhs) |
Does most of copying. | |
bool | createRim (int what, bool makeRowCopy=false, int startFinishOptions=0) |
puts in format I like (rowLower,rowUpper) also see StandardMatrix 1 bit does rows (now and columns), (2 bit does column bounds), 4 bit does objective(s). | |
void | createRim1 (bool initial) |
Does rows and columns. | |
void | createRim4 (bool initial) |
Does objective. | |
void | createRim5 (bool initial) |
Does rows and columns and objective. | |
void | deleteRim (int getRidOfFactorizationData=2) |
releases above arrays and does solution scaling out. | |
bool | sanityCheck () |
Sanity check on input rim data (after scaling) - returns true if okay. | |
int | housekeeping (double objectiveChange) |
This does basis housekeeping and does values for in/out variables. | |
void | checkPrimalSolution (const double *rowActivities=NULL, const double *columnActivies=NULL) |
This sets largest infeasibility and most infeasible and sum and number of infeasibilities (Primal) | |
void | checkDualSolution () |
This sets largest infeasibility and most infeasible and sum and number of infeasibilities (Dual) | |
void | checkBothSolutions () |
This sets sum and number of infeasibilities (Dual and Primal) | |
double | scaleObjective (double value) |
If input negative scales objective so maximum <= -value and returns scale factor used. | |
int | solveDW (CoinStructuredModel *model) |
Solve using Dantzig-Wolfe decomposition and maybe in parallel. | |
int | solveBenders (CoinStructuredModel *model) |
Solve using Benders decomposition and maybe in parallel. | |
![]() | |
double | bestPossibleImprovement_ |
Best possible improvement using djs (primal) or obj change by flipping bounds to make dual feasible (dual) | |
double | zeroTolerance_ |
Zero tolerance. | |
int | columnPrimalSequence_ |
Sequence of worst (-1 if feasible) | |
int | rowPrimalSequence_ |
Sequence of worst (-1 if feasible) | |
double | bestObjectiveValue_ |
"Best" objective value | |
int | moreSpecialOptions_ |
More special options - see set for details. | |
int | baseIteration_ |
Iteration when we entered dual or primal. | |
double | primalToleranceToGetOptimal_ |
Primal tolerance needed to make dual feasible (<largeTolerance) | |
double | largeValue_ |
Large bound value (for complementarity etc) | |
double | largestPrimalError_ |
Largest error on Ax-b. | |
double | largestDualError_ |
Largest error on basic duals. | |
double | alphaAccuracy_ |
For computing whether to re-factorize. | |
double | dualBound_ |
Dual bound. | |
double | alpha_ |
Alpha (pivot element) | |
double | theta_ |
Theta (pivot change) | |
double | lowerIn_ |
Lower Bound on In variable. | |
double | valueIn_ |
Value of In variable. | |
double | upperIn_ |
Upper Bound on In variable. | |
double | dualIn_ |
Reduced cost of In variable. | |
double | lowerOut_ |
Lower Bound on Out variable. | |
double | valueOut_ |
Value of Out variable. | |
double | upperOut_ |
Upper Bound on Out variable. | |
double | dualOut_ |
Infeasibility (dual) or ? (primal) of Out variable. | |
double | dualTolerance_ |
Current dual tolerance for algorithm. | |
double | primalTolerance_ |
Current primal tolerance for algorithm. | |
double | sumDualInfeasibilities_ |
Sum of dual infeasibilities. | |
double | sumPrimalInfeasibilities_ |
Sum of primal infeasibilities. | |
double | infeasibilityCost_ |
Weight assigned to being infeasible in primal. | |
double | sumOfRelaxedDualInfeasibilities_ |
Sum of Dual infeasibilities using tolerance based on error in duals. | |
double | sumOfRelaxedPrimalInfeasibilities_ |
Sum of Primal infeasibilities using tolerance based on error in primals. | |
double | acceptablePivot_ |
Acceptable pivot value just after factorization. | |
double * | lower_ |
Working copy of lower bounds (Owner of arrays below) | |
double * | rowLowerWork_ |
Row lower bounds - working copy. | |
double * | columnLowerWork_ |
Column lower bounds - working copy. | |
double * | upper_ |
Working copy of upper bounds (Owner of arrays below) | |
double * | rowUpperWork_ |
Row upper bounds - working copy. | |
double * | columnUpperWork_ |
Column upper bounds - working copy. | |
double * | cost_ |
Working copy of objective (Owner of arrays below) | |
double * | rowObjectiveWork_ |
Row objective - working copy. | |
double * | objectiveWork_ |
Column objective - working copy. | |
CoinIndexedVector * | rowArray_ [6] |
Useful row length arrays. | |
CoinIndexedVector * | columnArray_ [6] |
Useful column length arrays. | |
int | sequenceIn_ |
Sequence of In variable. | |
int | directionIn_ |
Direction of In, 1 going up, -1 going down, 0 not a clude. | |
int | sequenceOut_ |
Sequence of Out variable. | |
int | directionOut_ |
Direction of Out, 1 to upper bound, -1 to lower bound, 0 - superbasic. | |
int | pivotRow_ |
Pivot Row. | |
int | lastGoodIteration_ |
Last good iteration (immediately after a re-factorization) | |
double * | dj_ |
Working copy of reduced costs (Owner of arrays below) | |
double * | rowReducedCost_ |
Reduced costs of slacks not same as duals (or - duals) | |
double * | reducedCostWork_ |
Possible scaled reduced costs. | |
double * | solution_ |
Working copy of primal solution (Owner of arrays below) | |
double * | rowActivityWork_ |
Row activities - working copy. | |
double * | columnActivityWork_ |
Column activities - working copy. | |
int | numberDualInfeasibilities_ |
Number of dual infeasibilities. | |
int | numberDualInfeasibilitiesWithoutFree_ |
Number of dual infeasibilities (without free) | |
int | numberPrimalInfeasibilities_ |
Number of primal infeasibilities. | |
int | numberRefinements_ |
How many iterative refinements to do. | |
ClpDualRowPivot * | dualRowPivot_ |
dual row pivot choice | |
ClpPrimalColumnPivot * | primalColumnPivot_ |
primal column pivot choice | |
int * | pivotVariable_ |
Basic variables pivoting on which rows. | |
ClpFactorization * | factorization_ |
factorization | |
double * | savedSolution_ |
Saved version of solution. | |
int | numberTimesOptimal_ |
Number of times code has tentatively thought optimal. | |
ClpDisasterHandler * | disasterArea_ |
Disaster handler. | |
int | changeMade_ |
If change has been made (first attempt at stopping looping) | |
int | algorithm_ |
Algorithm >0 == Primal, <0 == Dual. | |
int | forceFactorization_ |
Now for some reliability aids This forces re-factorization early. | |
int | perturbation_ |
Perturbation: -50 to +50 - perturb by this power of ten (-6 sounds good) 100 - auto perturb if takes too long (1.0e-6 largest nonzero) 101 - we are perturbed 102 - don't try perturbing again default is 100. | |
unsigned char * | saveStatus_ |
Saved status regions. | |
ClpNonLinearCost * | nonLinearCost_ |
Very wasteful way of dealing with infeasibilities in primal. | |
int | lastBadIteration_ |
So we know when to be cautious. | |
int | lastFlaggedIteration_ |
So we know when to open up again. | |
int | numberFake_ |
Can be used for count of fake bounds (dual) or fake costs (primal) | |
int | numberChanged_ |
Can be used for count of changed costs (dual) or changed bounds (primal) | |
int | progressFlag_ |
Progress flag - at present 0 bit says artificials out, 1 free in. | |
int | firstFree_ |
First free/super-basic variable (-1 if none) | |
int | numberExtraRows_ |
Number of extra rows. | |
int | maximumBasic_ |
Maximum number of basic variables - can be more than number of rows if GUB. | |
int | dontFactorizePivots_ |
If may skip final factorize then allow up to this pivots (default 20) | |
double | incomingInfeasibility_ |
For advanced use. | |
double | allowedInfeasibility_ |
Best possible improvement using djs (primal) or obj change by flipping bounds to make dual feasible (dual) | |
int | automaticScale_ |
Automatic scaling of objective and rhs and bounds. | |
int | maximumPerturbationSize_ |
Maximum perturbation array size (take out when code rewritten) | |
double * | perturbationArray_ |
Perturbation array (maximumPerturbationSize_) | |
ClpSimplex * | baseModel_ |
A copy of model with certain state - normally without cuts. | |
ClpSimplexProgress | progress_ |
For dealing with all issues of cycling etc. | |
![]() | |
double | optimizationDirection_ |
Direction of optimization (1 - minimize, -1 - maximize, 0 - ignore. | |
double | dblParam_ [ClpLastDblParam] |
Array of double parameters. | |
double | objectiveValue_ |
Objective value. | |
double | smallElement_ |
Small element value. | |
double | objectiveScale_ |
Scaling of objective. | |
double | rhsScale_ |
Scaling of rhs and bounds. | |
int | numberRows_ |
Number of rows. | |
int | numberColumns_ |
Number of columns. | |
double * | rowActivity_ |
Row activities. | |
double * | columnActivity_ |
Column activities. | |
double * | dual_ |
Duals. | |
double * | reducedCost_ |
Reduced costs. | |
double * | rowLower_ |
Row lower. | |
double * | rowUpper_ |
Row upper. | |
ClpObjective * | objective_ |
Objective. | |
double * | rowObjective_ |
Row Objective (? sign) - may be NULL. | |
double * | columnLower_ |
Column Lower. | |
double * | columnUpper_ |
Column Upper. | |
ClpMatrixBase * | matrix_ |
Packed matrix. | |
ClpMatrixBase * | rowCopy_ |
Row copy if wanted. | |
ClpPackedMatrix * | scaledMatrix_ |
Scaled packed matrix. | |
double * | ray_ |
Infeasible/unbounded ray. | |
double * | rowScale_ |
Row scale factors for matrix. | |
double * | columnScale_ |
Column scale factors. | |
double * | inverseRowScale_ |
Inverse row scale factors for matrix (end of rowScale_) | |
double * | inverseColumnScale_ |
Inverse column scale factors for matrix (end of columnScale_) | |
int | scalingFlag_ |
Scale flag, 0 none, 1 equilibrium, 2 geometric, 3, auto, 4 dynamic, 5 geometric on rows. | |
unsigned char * | status_ |
Status (i.e. | |
char * | integerType_ |
Integer information. | |
void * | userPointer_ |
User pointer for whatever reason. | |
ClpTrustedData * | trustedUserPointer_ |
Trusted user pointer e.g. for heuristics. | |
int | intParam_ [ClpLastIntParam] |
Array of integer parameters. | |
int | numberIterations_ |
Number of iterations. | |
int | solveType_ |
Solve type - 1 simplex, 2 simplex interface, 3 Interior. | |
unsigned int | whatsChanged_ |
Whats changed since last solve. | |
int | problemStatus_ |
Status of problem. | |
int | secondaryStatus_ |
Secondary status of problem. | |
int | lengthNames_ |
length of names (0 means no names) | |
int | numberThreads_ |
Number of threads (not very operational) | |
unsigned int | specialOptions_ |
For advanced options See get and set for meaning. | |
CoinMessageHandler * | handler_ |
Message handler. | |
bool | defaultHandler_ |
Flag to say if default handler (so delete) | |
CoinThreadRandom | randomNumberGenerator_ |
Thread specific random number generator. | |
ClpEventHandler * | eventHandler_ |
Event handler. | |
std::vector< std::string > | rowNames_ |
Row names. | |
std::vector< std::string > | columnNames_ |
Column names. | |
CoinMessages | messages_ |
Messages. | |
CoinMessages | coinMessages_ |
Coin messages. | |
int | maximumColumns_ |
Maximum number of columns in model. | |
int | maximumRows_ |
Maximum number of rows in model. | |
int | maximumInternalColumns_ |
Maximum number of columns (internal arrays) in model. | |
int | maximumInternalRows_ |
Maximum number of rows (internal arrays) in model. | |
CoinPackedMatrix | baseMatrix_ |
Base packed matrix. | |
CoinPackedMatrix | baseRowCopy_ |
Base row copy. | |
double * | savedRowScale_ |
Saved row scale factors for matrix. | |
double * | savedColumnScale_ |
Saved column scale factors. | |
std::string | strParam_ [ClpLastStrParam] |
Array of string parameters. | |
![]() |
This is for Simplex stuff which is neither dual nor primal.
It inherits from ClpSimplex. It has no data of its own and is never created - only cast from a ClpSimplex object at algorithm time.
Definition at line 23 of file ClpSimplexOther.hpp.
void ClpSimplexOther::dualRanging | ( | int | numberCheck, |
const int * | which, | ||
double * | costIncrease, | ||
int * | sequenceIncrease, | ||
double * | costDecrease, | ||
int * | sequenceDecrease, | ||
double * | valueIncrease = NULL , |
||
double * | valueDecrease = NULL |
||
) |
Dual ranging.
This computes increase/decrease in cost for each given variable and corresponding sequence numbers which would change basis. Sequence numbers are 0..numberColumns and numberColumns.. for artificials/slacks. For non-basic variables the information is trivial to compute and the change in cost is just minus the reduced cost and the sequence number will be that of the non-basic variables. For basic variables a ratio test is between the reduced costs for non-basic variables and the row of the tableau corresponding to the basic variable. The increase/decrease value is always >= 0.0
Up to user to provide correct length arrays where each array is of length numberCheck. which contains list of variables for which information is desired. All other arrays will be filled in by function. If fifth entry in which is variable 7 then fifth entry in output arrays will be information for variable 7.
If valueIncrease/Decrease not NULL (both must be NULL or both non NULL) then these are filled with the value of variable if such a change in cost were made (the existing bounds are ignored)
When here - guaranteed optimal
Reimplemented from ClpSimplex.
void ClpSimplexOther::primalRanging | ( | int | numberCheck, |
const int * | which, | ||
double * | valueIncrease, | ||
int * | sequenceIncrease, | ||
double * | valueDecrease, | ||
int * | sequenceDecrease | ||
) |
Primal ranging.
This computes increase/decrease in value for each given variable and corresponding sequence numbers which would change basis. Sequence numbers are 0..numberColumns and numberColumns.. for artificials/slacks. This should only be used for non-basic variabls as otherwise information is pretty useless For basic variables the sequence number will be that of the basic variables.
Up to user to provide correct length arrays where each array is of length numberCheck. which contains list of variables for which information is desired. All other arrays will be filled in by function. If fifth entry in which is variable 7 then fifth entry in output arrays will be information for variable 7.
When here - guaranteed optimal
Reimplemented from ClpSimplex.
int ClpSimplexOther::parametrics | ( | double | startingTheta, |
double & | endingTheta, | ||
double | reportIncrement, | ||
const double * | changeLowerBound, | ||
const double * | changeUpperBound, | ||
const double * | changeLowerRhs, | ||
const double * | changeUpperRhs, | ||
const double * | changeObjective | ||
) |
Parametrics This is an initial slow version.
The code uses current bounds + theta * change (if change array not NULL) and similarly for objective. It starts at startingTheta and returns ending theta in endingTheta. If reportIncrement 0.0 it will report on any movement If reportIncrement >0.0 it will report at startingTheta+k*reportIncrement. If it can not reach input endingTheta return code will be 1 for infeasible, 2 for unbounded, if error on ranges -1, otherwise 0. Normal report is just theta and objective but if event handler exists it may do more On exit endingTheta is maximum reached (can be used for next startingTheta)
|
private |
Parametrics - inner loop This first attempt is when reportIncrement non zero and may not report endingTheta correctly If it can not reach input endingTheta return code will be 1 for infeasible, 2 for unbounded, otherwise 0.
Normal report is just theta and objective but if event handler exists it may do more
|
private |
Refactorizes if necessary Checks if finished.
Updates status.
type - 0 initial so set up save arrays etc
2 restoring from saved
|
private |
This has the flow between re-factorizations.
Reasons to come out: -1 iterations etc -2 inaccuracy -3 slight inaccuracy (and done iterations) +0 looks optimal (might be unbounded - but we will investigate) +1 looks infeasible +3 max iterations
|
private |
Computes next theta and says if objective or bounds (0= bounds, 1 objective, -1 none).
theta is in theta_. type 1 bounds, 2 objective, 3 both.
|
private |
Row array has row part of pivot row Column array has column part.
This is used in dual ranging
|
private |
Row array has pivot column This is used in primal ranging.
|
private |
Returns new value of whichOther when whichIn enters basis.
int ClpSimplexOther::writeBasis | ( | const char * | filename, |
bool | writeValues = false , |
||
int | formatType = 0 |
||
) | const |
Write the basis in MPS format to the specified file.
If writeValues true writes values of structurals (and adds VALUES to end of NAME card)
Row and column names may be null. formatType is
Returns non-zero on I/O error
Reimplemented from ClpSimplex.
int ClpSimplexOther::readBasis | ( | const char * | filename | ) |
Read a basis from the given filename.
Reimplemented from ClpSimplex.
ClpSimplex* ClpSimplexOther::dualOfModel | ( | double | fractionRowRanges = 1.0 , |
double | fractionColumnRanges = 1.0 |
||
) | const |
Creates dual of a problem if looks plausible (defaults will always create model) fractionRowRanges is fraction of rows allowed to have ranges fractionColumnRanges is fraction of columns allowed to have ranges.
int ClpSimplexOther::restoreFromDual | ( | const ClpSimplex * | dualProblem | ) |
Restores solution from dualized problem non-zero return code indicates minor problems.
ClpSimplex* ClpSimplexOther::crunch | ( | double * | rhs, |
int * | whichRows, | ||
int * | whichColumns, | ||
int & | nBound, | ||
bool | moreBounds = false , |
||
bool | tightenBounds = false |
||
) |
Does very cursory presolve.
rhs is numberRows, whichRows is 3*numberRows and whichColumns is 2*numberColumns.
void ClpSimplexOther::afterCrunch | ( | const ClpSimplex & | small, |
const int * | whichRows, | ||
const int * | whichColumns, | ||
int | nBound | ||
) |
After very cursory presolve.
rhs is numberRows, whichRows is 3*numberRows and whichColumns is 2*numberColumns.
void ClpSimplexOther::cleanupAfterPostsolve | ( | ) |
Quick try at cleaning up duals if postsolve gets wrong.
int ClpSimplexOther::tightenIntegerBounds | ( | double * | rhsSpace | ) |
Tightens integer bounds - returns number tightened or -1 if infeasible.
int ClpSimplexOther::expandKnapsack | ( | int | knapsackRow, |
int & | numberOutput, | ||
double * | buildObj, | ||
CoinBigIndex * | buildStart, | ||
int * | buildRow, | ||
double * | buildElement, | ||
int | reConstruct = -1 |
||
) | const |
Expands out all possible combinations for a knapsack If buildObj NULL then just computes space needed - returns number elements On entry numberOutput is maximum allowed, on exit it is number needed or -1 (as will be number elements) if maximum exceeded.
numberOutput will have at least space to return values which reconstruct input. Rows returned will be original rows but no entries will be returned for any rows all of whose entries are in knapsack. So up to user to allow for this. If reConstruct >=0 then returns number of entrie which make up item "reConstruct" in expanded knapsack. Values in buildRow and buildElement;