Doxygen/OS/_o_s_cgl_cuts_8cpp_source.html

 /* $Id: OSAddCuts.cpp 2710 2009-06-10 21:13:43Z kmartin $ */

 #include <cppad/cppad.hpp>

 #include "OSConfig.h"

 #include "OSCoinSolver.h"

 #include "OSResult.h"

 #include "OSiLReader.h"

 #include "OSiLWriter.h"

 #include "OSrLReader.h"

 #include "OSrLWriter.h"

 #include "OSInstance.h"

 #include "OSFileUtil.h"

 #include "OSDefaultSolver.h"

 #include "OSWSUtil.h"

 #include "OSSolverAgent.h"

 #include "OShL.h"

 #include "OSErrorClass.h"

 #include "OSmps2OS.h"

 #include "OSBase64.h"

 #include "OSErrorClass.h"

 #include "CglGomory.hpp"

 #include "CglSimpleRounding.hpp"

 #include "CglKnapsackCover.hpp"

 #include "OSMathUtil.h"

 #include "CbcModel.hpp"


 #include <iostream>

 using std::cout;

 using std::endl;


 int main( ){

         WindowsErrorPopupBlocker();

 // test OS code samples here

         FileUtil *fileUtil = NULL;

         fileUtil = new FileUtil();

         const char dirsep =  CoinFindDirSeparator();

         // Set directory containing mps data files.

         std::string dataDir;

    dataDir = dirsep == '/' ? "../data/" : "..\\data\\";

         cout << "Start Building the Model" << endl;

         try{

                 // get the p0033 problem

                 std::string osilFileName;

                 osilFileName =  dataDir  +  "osilFiles" + dirsep +  "p0033.osil";

                 std::cout << "Try to read a sample file" << std::endl;

                 std::cout << "The file is: " ;

                 std::cout <<  osilFileName << std::endl;

                 std::string osil = fileUtil->getFileAsString( osilFileName.c_str() );

                 OSiLReader *osilreader = NULL;

                 osilreader = new OSiLReader();

                 OSInstance *osinstance;

                 osinstance = osilreader->readOSiL( osil);

                 // done writing the model

                 cout << "Done writing the Model" << endl;

                 // now solve the model

                 CoinSolver *solver  = NULL;

                 solver = new CoinSolver();

                 solver->sSolverName ="cbc";

                 solver->osinstance = osinstance;

                 solver->buildSolverInstance();

                 solver->osiSolver->setHintParam(OsiDoReducePrint, true, OsiHintTry);

                 solver->osiSolver->initialSolve();

                 cout << "Here is the initial objective value "  << solver->osiSolver->getObjValue()  << endl;


                 CglKnapsackCover cover;

             CglSimpleRounding round;

                 CglGomory gomory;

                 CbcModel *model = new CbcModel( *solver->osiSolver);


                 //model->setBestObjectiveValue(4000);

                 model->setMaximumNodes(100000);

                 //

                 model->addCutGenerator(&gomory, 1, "Gomory");

                 model->addCutGenerator(&cover, 1, "Cover");

                 model->addCutGenerator(&round, 1, "Round");

                 model->branchAndBound();

                 // now create a result object

                 OSResult *osresult = new OSResult();

                 // if we are throw an exception if the problem is nonlinear

                 double *x = NULL;

                 double *y = NULL;

                 double *z = NULL;

                 //int i = 0;

                 std::string *rcost = NULL;

                 // resultHeader infomration

                 if(osresult->setServiceName("Solved with Coin Solver: " + solver->sSolverName) != true)

                         throw ErrorClass("OSResult error: setServiceName");

                 if(osresult->setInstanceName(  solver->osinstance->getInstanceName()) != true)

                         throw ErrorClass("OSResult error: setInstanceName");

                 if(osresult->setVariableNumber( solver->osinstance->getVariableNumber()) != true)

                         throw ErrorClass("OSResult error: setVariableNumer");

                 if(osresult->setObjectiveNumber( 1) != true)

                         throw ErrorClass("OSResult error: setObjectiveNumber");

                 if(osresult->setConstraintNumber( solver->osinstance->getConstraintNumber()) != true)

                         throw ErrorClass("OSResult error: setConstraintNumber");

                 if(osresult->setSolutionNumber(  1) != true)

                         throw ErrorClass("OSResult error: setSolutionNumer");

                 int solIdx = 0;

                 std::string description = "";

                 osresult->setGeneralStatusType("success");

                 std::cout << "PROVEN OPTIMAL " << model->isProvenOptimal() << std::endl;

                 int i;

                 if (model->isProvenOptimal() == 1){

                         osresult->setSolutionStatus(solIdx, "optimal", description);

                         /* Retrieve the solution */

                         x = new double[solver->osinstance->getVariableNumber() ];

                         y = new double[solver->osinstance->getConstraintNumber() ];

                         z = new double[1];

                         rcost = new std::string[ solver->osinstance->getVariableNumber()];

                         //

                         *(z + 0)  =  model->getObjValue();

                         osresult->setObjectiveValuesDense(solIdx, z);

                         for(i=0; i < solver->osinstance->getVariableNumber(); i++){

                                 *(x + i) = model->getColSolution()[i];

                         }

                         osresult->setPrimalVariableValuesDense(solIdx, x );

                         //if( solver->sSolverName.find( "symphony") == std::string::npos){

                         for(i=0; i <  solver->osinstance->getConstraintNumber(); i++){

                                 *(y + i) = model->getRowPrice()[ i];

                         }

                         osresult->setDualVariableValuesDense(solIdx, y);

                         //

                         // now put the reduced costs into the osrl

                         int numberOfOtherVariableResult = 1;

                         int otherIdx = 0;

                         // first set the number of Other Variable Results

                         osresult->setNumberOfOtherVariableResults(solIdx, numberOfOtherVariableResult);


                         std::ostringstream outStr;

                         int numberOfVar =  solver->osinstance->getVariableNumber();

                         for(i=0; i < numberOfVar; i++){

                                 outStr << model->getReducedCost()[ i];

                                 rcost[ i] = outStr.str();

                                 outStr.str("");

                         }

                         osresult->setAnOtherVariableResultDense(solIdx, otherIdx, "reduced costs", "", "the variable reduced costs",

                                 rcost);

                         // end of settiing reduced costs

                 }

                 else{

                         if(solver->osiSolver->isProvenPrimalInfeasible() == true)

                                 osresult->setSolutionStatus(solIdx, "infeasible", description);

                         else

                                 if(solver->osiSolver->isProvenDualInfeasible() == true)

                                         osresult->setSolutionStatus(solIdx, "dualinfeasible", description);

                                 else

                                         osresult->setSolutionStatus(solIdx, "other", description);

                 }

                 OSrLWriter *osrlwriter = new OSrLWriter();

                 std::cout <<  osrlwriter->writeOSrL( osresult) << std::endl;

                 if(solver->osinstance->getVariableNumber() > 0){

                         delete[] x;

                         x = NULL;

                 }

                 if(solver->osinstance->getConstraintNumber()) delete[] y;

                 y = NULL;

                 delete[] z;

                 z = NULL;

                 if(solver->osinstance->getVariableNumber() > 0){

                         delete[] rcost;

                         rcost = NULL;

                 }

         // do garbage collection

                 delete osresult;

                 osresult = NULL;

                 delete osrlwriter;

                 osrlwriter = NULL;

                 delete solver;

                 solver = NULL;

                 delete osilreader;

                 osilreader = NULL;

                 delete fileUtil;

                 fileUtil  = NULL;

                 delete model;

                 model = NULL;

                 cout << "Done with garbage collection" << endl;

                 cout << "Program terminates normally" << endl;

                 return 0;

         }

         catch(const ErrorClass& eclass){

                 delete fileUtil;

                 std::cout << eclass.errormsg <<  std::endl;

                 return 0;

         }

 }// end main


OSWSUtil.h

OSErrorClass.h

OSResult::setSolutionStatus
bool setSolutionStatus(int solIdx, std::string type, std::string description)
Set the [i]th optimization solution status, where i equals the given solution index.

OSResult::setPrimalVariableValuesDense
bool setPrimalVariableValuesDense(int solIdx, double *x)
Set the [i]th optimization solution&#39;s primal variable values, where i equals the given solution index...
Definition: OSResult.cpp:5001

OSInstance::getVariableNumber
int getVariableNumber()
Get number of variables.
Definition: OSInstance.cpp:2018

OSResult::setServiceName
bool setServiceName(std::string serviceName)
Set service name.

OSResult::setVariableNumber
bool setVariableNumber(int variableNumber)
Set the variable number.
Definition: OSResult.cpp:4712

ErrorClass::errormsg
std::string errormsg
errormsg is the error that is causing the exception to be thrown
Definition: OSErrorClass.h:42

main
int main(int argc, char *argv[])
Definition: BB_tm.cpp:32

CoinSolver::osiSolver
OsiSolverInterface * osiSolver
osiSolver is the osi solver object – in this case clp, glpk, cbc, cplex, symphony or dylp ...
Definition: OSCoinSolver.h:93

OSCoinSolver.h

OSResult
The Result Class.
Definition: OSResult.h:2548

OSrLWriter
Take an OSResult object and write a string that validates against OSrL.
Definition: OSrLWriter.h:30

OSiLReader.h

OSFileUtil.h

OSResult::setDualVariableValuesDense
bool setDualVariableValuesDense(int solIdx, double *y)
Set the [i]th optimization solution&#39;s dual variable values, where i equals the given solution index...
Definition: OSResult.cpp:6291

OSResult::setAnOtherVariableResultDense
bool setAnOtherVariableResultDense(int solIdx, int otherIdx, std::string name, std::string value, std::string description, std::string *s)
Set the [i]th optimization solution&#39;s other (non-standard/solver specific)variable-related results...

OSResult::setObjectiveNumber
bool setObjectiveNumber(int objectiveNumber)
Set the objective number.
Definition: OSResult.cpp:4721

OSResult::setInstanceName
bool setInstanceName(std::string instanceName)
Set instance name.

OSiLReader::readOSiL
OSInstance * readOSiL(const std::string &osil)
parse the OSiL model instance.
Definition: OSiLReader.cpp:53

OSrLWriter::writeOSrL
std::string writeOSrL(OSResult *theosresult)
create an osrl string from an OSResult object
Definition: OSrLWriter.cpp:45

CoinSolver::buildSolverInstance
virtual void buildSolverInstance()
The implementation of the corresponding virtual function.
Definition: OSCoinSolver.cpp:145

OSResult::setSolutionNumber
bool setSolutionNumber(int number)
set the number of solutions.
Definition: OSResult.cpp:4740

OSiLReader
Used to read an OSiL string.
Definition: OSiLReader.h:37

OSResult::setNumberOfOtherVariableResults
bool setNumberOfOtherVariableResults(int solIdx, int numberOfOtherVariableResults)
Set the [i]th optimization solution&#39;s other (non-standard/solver specific) variable-related results...
Definition: OSResult.cpp:5236

OSrLReader.h

osresult
OSResult * osresult
Definition: OSParseosrl.tab.cpp:5025

OSDefaultSolver.h

DefaultSolver::osinstance
OSInstance * osinstance
osinstance holds the problem instance in-memory as an OSInstance object
Definition: OSDefaultSolver.h:54

OSResult::setGeneralStatusType
bool setGeneralStatusType(std::string type)
Set the general status type, which can be: success, error, warning.

OSInstance::getInstanceName
std::string getInstanceName()
Get instance name.
Definition: OSInstance.cpp:1963

OSConfig.h

OSInstance::getConstraintNumber
int getConstraintNumber()
Get number of constraints.
Definition: OSInstance.cpp:2327

OSResult::setConstraintNumber
bool setConstraintNumber(int constraintNumber)
Set the constraint number.
Definition: OSResult.cpp:4731

DefaultSolver::sSolverName
std::string sSolverName
sSolverName is the name of the Coin solver used, e.g.
Definition: OSDefaultSolver.h:68

OSSolverAgent.h

CoinSolver
Implements a solve method for the Coin solvers.
Definition: OSCoinSolver.h:37

FileUtil::getFileAsString
std::string getFileAsString(const char *fname)
read a file and return contents as a string.
Definition: OSFileUtil.cpp:35

OSrLWriter.h

OSmps2OS.h

OSMathUtil.h

OSResult::setObjectiveValuesDense
bool setObjectiveValuesDense(int solIdx, double *objectiveValues)
Set the [i]th optimization solution&#39;s objective values, where i equals the given solution index...
Definition: OSResult.cpp:5824

OSInstance
The in-memory representation of an OSiL instance..
Definition: OSInstance.h:2262

osinstance
OSInstance * osinstance
Definition: OSParseosil.tab.cpp:3860

FileUtil
class used to make it easy to read and write files.
Definition: OSFileUtil.h:37

OSiLWriter.h

ErrorClass
used for throwing exceptions.
Definition: OSErrorClass.h:31

OShL.h

OSBase64.h

x
void fint fint fint real fint real * x
Definition: BonBqpdSolver.cpp:27