// Create C++ lines to get to current state
std::string
CglClique::generateCpp( FILE * fp)
{
CglClique other;
fprintf(fp,"0#include \"CglClique.hpp\"\n");
fprintf(fp,"3 CglClique clique;\n");
std::string types[] = {"SCL_MIN_DEGREE","SCL_MAX_DEGREE",
"SCL_MAX_XJ_MAX_DEG"};
if (scl_next_node_rule!=other.scl_next_node_rule)
fprintf(fp,"3 clique.setStarCliqueNextNodeMethod(CglClique::%s);\n",
types[scl_next_node_rule].c_str());
else
fprintf(fp,"4 clique.setStarCliqueNextNodeMethod(CglClique::%s);\n",
types[scl_next_node_rule].c_str());
if (scl_candidate_length_threshold!=other.scl_candidate_length_threshold)
fprintf(fp,"3 clique.setStarCliqueCandidateLengthThreshold(%d);\n",
scl_candidate_length_threshold);
else
fprintf(fp,"4 clique.setStarCliqueCandidateLengthThreshold(%d);\n",
scl_candidate_length_threshold);
if (rcl_candidate_length_threshold!=other.rcl_candidate_length_threshold)
fprintf(fp,"3 clique.setRowCliqueCandidateLengthThreshold(%d);\n",
rcl_candidate_length_threshold);
else
fprintf(fp,"4 clique.setRowCliqueCandidateLengthThreshold(%d);\n",
rcl_candidate_length_threshold);
if (scl_report_result!=other.scl_report_result)
fprintf(fp,"3 clique.setStarCliqueReport(%s);\n",
scl_report_result ? "true" : "false");
else
fprintf(fp,"4 clique.setStarCliqueReport(%s);\n",
scl_report_result ? "true" : "false");
if (rcl_report_result!=other.rcl_report_result)
fprintf(fp,"3 clique.setRowCliqueReport(%s);\n",
rcl_report_result ? "true" : "false");
else
fprintf(fp,"4 clique.setRowCliqueReport(%s);\n",
rcl_report_result ? "true" : "false");
if (do_star_clique!=other.do_star_clique)
fprintf(fp,"3 clique.setDoStarClique(%s);\n",
do_star_clique ? "true" : "false");
else
fprintf(fp,"4 clique.setDoStarClique(%s);\n",
do_star_clique ? "true" : "false");
if (do_row_clique!=other.do_row_clique)
fprintf(fp,"3 clique.setDoRowClique(%s);\n",
do_row_clique ? "true" : "false");
else
fprintf(fp,"4 clique.setDoRowClique(%s);\n",
do_row_clique ? "true" : "false");
if (petol!=other.petol)
fprintf(fp,"3 clique.setMinViolation(%g);\n",petol);
else
fprintf(fp,"4 clique.setMinViolation(%g);\n",petol);
if (getAggressiveness()!=other.getAggressiveness())
fprintf(fp,"3 clique.setAggressiveness(%d);\n",getAggressiveness());
else
fprintf(fp,"4 clique.setAggressiveness(%d);\n",getAggressiveness());
return "clique";
}
int main( int argc, char **argv )
{
if ( argc < 2 )
{
printf("Invalid number of parameters!\n");
exit( EXIT_FAILURE );
}
char problemName[ 256 ];
getFileName( problemName, argv[1] );
clock_t start = clock();
OsiClpSolverInterface *realSolver = new OsiClpSolverInterface();
realSolver->getModelPtr()->setPerturbation(50); /* makes CLP faster for hard instances */
OsiSolverInterface *solver = (OsiSolverInterface*) realSolver;
parseParameters( argc, argv );
readLP( solver, argv[1] );
FILE *log = NULL;
if(!output.empty())
{
log = fopen(output.c_str(), "a");
if(!log)
{
printf("Could not open the file!\n");
exit(EXIT_FAILURE);
}
}
const int numCols = solver->getNumCols(), numRows = solver->getNumRows();
int pass = 0, newCuts = 0, totalCuts = 0;
double pTime, opt, cgTime;
CGraph *cgraph = NULL;
if(sepMethod == Npsep)
cgraph = build_cgraph_osi( solver );
if(!optFile.empty())
{
getOptimals();
if(optimals.find(problemName) == optimals.end())
{
fprintf(stderr, "ERROR: optimal value not found!\n");
exit(EXIT_FAILURE);
}
opt = optimals[problemName];
}
solver->initialSolve();
if (!solver->isProvenOptimal())
{
if (solver->isAbandoned())
{
fprintf( stderr, "LP solver abandoned due to numerical dificulties.\n" );
exit( EXIT_FAILURE );
}
if (solver->isProvenPrimalInfeasible())
{
fprintf( stderr, "LP solver says PRIMAL INFEASIBLE.\n" );
exit( EXIT_FAILURE );
}
if (solver->isProvenDualInfeasible())
{
fprintf( stderr, "LP solver says DUAL INFEASIBLE.\n" );
exit( EXIT_FAILURE );
}
if (solver->isPrimalObjectiveLimitReached())
{
fprintf( stderr, "LP solver says isPrimalObjectiveLimitReached.\n" );
exit( EXIT_FAILURE );
}
if (solver->isDualObjectiveLimitReached())
{
fprintf( stderr, "LP solver says isDualObjectiveLimitReached.\n" );
exit( EXIT_FAILURE );
}
if (solver->isIterationLimitReached())
{
fprintf( stderr, "LP solver says isIterationLimitReached.\n" );
exit( EXIT_FAILURE );
}
fprintf( stderr, "ERROR: Could not solve LP relaxation to optimality. Checking status...\n" );
exit( EXIT_FAILURE );
}
double initialBound = solver->getObjValue();
printf("%.2lf %d %d %.7lf", ((double)(clock()-start)) / ((double)CLOCKS_PER_SEC), pass, 0, solver->getObjValue());
if(!optFile.empty())
{
printf(" %.7lf %.7lf", opt, abs_mip_gap(solver->getObjValue(), opt));
}
printf("\n");
do
{
clock_t startSep = clock();
newCuts = 0;
switch (sepMethod)
{
case Npsep:
{
CglEClique cliqueGen;
OsiCuts cuts;
CglTreeInfo info;
info.level = 0;
info.pass = 1;
vector<string> varNames = getVarNames(solver->getColNames(), numCols);
cliqueGen.parseParameters( argc, (const char**)argv );
cliqueGen.setCGraph( cgraph );
cliqueGen.setGenOddHoles( true ); //allow (or not) inserting odd hole cuts
cliqueGen.colNames = &varNames;
cliqueGen.generateCuts( *solver, cuts, info );
newCuts = cuts.sizeCuts();
solver->applyCuts( cuts );
}
break;
case CglSepM:
{
CglClique cliqueGen;
OsiCuts cuts;
CglTreeInfo info;
info.level = 0;
info.pass = 1;
cliqueGen.setMinViolation( MIN_VIOLATION );
cliqueGen.setStarCliqueReport(false);
cliqueGen.setRowCliqueReport(false);
cliqueGen.generateCuts( *solver, cuts, info );
newCuts = cuts.sizeCuts();
solver->applyCuts( cuts );
}
break;
Default:
{
fprintf( stderr, "Separation Method does not recognized!\n" );
exit( EXIT_FAILURE );
}
}
pTime = ((double)(clock()-start)) / ((double)CLOCKS_PER_SEC);
if(pTime > MAX_TIME) break;
totalCuts += newCuts;
++pass;
if (newCuts)
{
solver->resolve();
if (!solver->isProvenOptimal())
{
if (solver->isAbandoned())
{
fprintf( stderr, "LP solver abandoned due to numerical dificulties.\n" );
exit( EXIT_FAILURE );
}
if (solver->isProvenPrimalInfeasible())
{
fprintf( stderr, "LP solver says PRIMAL INFEASIBLE.\n" );
exit( EXIT_FAILURE );
}
if (solver->isProvenDualInfeasible())
{
fprintf( stderr, "LP solver says DUAL INFEASIBLE.\n" );
exit( EXIT_FAILURE );
}
if (solver->isPrimalObjectiveLimitReached())
{
fprintf( stderr, "LP solver says isPrimalObjectiveLimitReached.\n" );
exit( EXIT_FAILURE );
}
if (solver->isDualObjectiveLimitReached())
{
fprintf( stderr, "LP solver says isDualObjectiveLimitReached.\n" );
exit( EXIT_FAILURE );
}
if (solver->isIterationLimitReached())
{
fprintf( stderr, "LP solver says isIterationLimitReached.\n" );
exit( EXIT_FAILURE );
}
fprintf( stderr, "ERROR: Could not solve LP relaxation. Exiting.\n" );
exit( EXIT_FAILURE );
}
pTime = ((double)(clock()-start)) / ((double)CLOCKS_PER_SEC);
if(pTime > MAX_TIME) break;
double sepTime = ((double)(clock()-startSep)) / ((double)CLOCKS_PER_SEC);
printf("%.2lf %d %d %.7lf", sepTime, pass, newCuts, solver->getObjValue());
if(!optFile.empty())
printf(" %.7lf %.7lf", opt, abs_mip_gap(solver->getObjValue(), opt));
printf("\n");
}
}
while ( (newCuts>0) && (pass<MAX_PASSES) ) ;
if(log)
{
double totalTime = ((double)(clock()-start)) / ((double)CLOCKS_PER_SEC);
fprintf(log, "%s %.2lf %d %d %.7lf", problemName, totalTime, pass - 1, totalCuts, solver->getObjValue());
if(!optFile.empty())
fprintf(log, " %.7lf", abs_mip_gap(solver->getObjValue(), opt));
fprintf(log, "\n");
}
if(cgraph)
cgraph_free( &cgraph );
delete realSolver;
return EXIT_SUCCESS;
}
void
CglCliqueUnitTest(const OsiSolverInterface *baseSiP,
const std::string mpsDir)
{
// Test default constructor
{
CglClique aGenerator;
}
// Test copy & assignment
{
CglClique rhs;
{
CglClique bGenerator;
CglClique cGenerator(bGenerator);
//rhs=bGenerator;
}
}
// Test get/set methods
{
CglClique getset;
// None to test
}
// Test generateCuts
{
CglClique gct;
OsiSolverInterface *siP = baseSiP->clone();
std::string fn = mpsDir+"l152lav";
std::string fn2 = mpsDir+"l152lav.mps";
FILE *in_f = fopen(fn2.c_str(), "r");
if(in_f == NULL) {
std::cout<<"Can not open file "<<fn2<<std::endl<<"Skip test of CglClique::generateCuts()"<<std::endl;
}
else {
fclose(in_f);
siP->readMps(fn.c_str(),"mps");
siP->initialSolve();
double lpRelax = siP->getObjValue();
OsiCuts cs;
gct.generateCuts(*siP, cs);
int nRowCuts = cs.sizeRowCuts();
std::cout<<"There are "<<nRowCuts<<" Clique cuts"<<std::endl;
assert(cs.sizeRowCuts() > 0);
OsiSolverInterface::ApplyCutsReturnCode rc = siP->applyCuts(cs);
siP->resolve();
double lpRelaxAfter= siP->getObjValue();
std::cout<<"Initial LP value: "<<lpRelax<<std::endl;
std::cout<<"LP value with cuts: "<<lpRelaxAfter<<std::endl;
assert( lpRelax < lpRelaxAfter );
assert(lpRelaxAfter < 4722.1);
}
delete siP;
}
}
void LinearCutsGenerator::initialize(BabSetupBase & s){
assert(dynamic_cast<TMINLP2TNLPQuadCuts *> (s.nonlinearSolver()->problem()));
int freq;
s.options()->GetIntegerValue("Gomory_cuts", freq,"bonmin.");
if (freq) {
Coin::SmartPtr<CuttingMethod> cg = new CuttingMethod;
cg->frequency = freq;
CglGomory * gom = new CglGomory;
cg->cgl = gom;
gom->setLimitAtRoot(512);
gom->setLimit(50);
cg->id = "Mixed Integer Gomory";
methods_.push_back(cg);
}
s.options()->GetIntegerValue("mir_cuts",freq,"bonmin.");
if (freq) {
Coin::SmartPtr<CuttingMethod> cg = new CuttingMethod;
cg->frequency = freq;
CglMixedIntegerRounding2 * mir = new CglMixedIntegerRounding2;
cg->cgl = mir;
cg->id = "Mixed Integer Rounding";
methods_.push_back(cg);
}
s.options()->GetIntegerValue("2mir_cuts",freq,"bonmin.");
if (freq) {
Coin::SmartPtr<CuttingMethod> cg = new CuttingMethod;
cg->frequency = freq;
CglTwomir * mir2 = new CglTwomir;
cg->cgl = mir2;
cg->id = "2-MIR";
methods_.push_back(cg);
}
s.options()->GetIntegerValue("cover_cuts",freq,"bonmin.");
if (freq) {
Coin::SmartPtr<CuttingMethod> cg = new CuttingMethod;
cg->frequency = freq;
CglKnapsackCover * cover = new CglKnapsackCover;
cg->cgl = cover;
cg->id = "Cover";
methods_.push_back(cg);
}
s.options()->GetIntegerValue("clique_cuts",freq,"bonmin.");
if (freq) {
Coin::SmartPtr<CuttingMethod> cg = new CuttingMethod;
cg->frequency = freq;
CglClique * clique = new CglClique;
clique->setStarCliqueReport(false);
clique->setRowCliqueReport(false);
clique->setMinViolation(0.1);
cg->cgl = clique;
cg->id = "Clique";
methods_.push_back(cg);
}
s.options()->GetIntegerValue("flow_cover_cuts",freq,"bonmin.");
if (freq) {
Coin::SmartPtr<CuttingMethod> cg = new CuttingMethod;
cg->frequency = freq;
CglFlowCover * flow = new CglFlowCover;
cg->cgl = flow;
cg->id = "Flow Covers";
methods_.push_back(cg);
}
s.options()->GetIntegerValue("lift_and_project_cuts",freq,"bonmin.");
if (freq) {
Coin::SmartPtr<CuttingMethod> cg = new CuttingMethod;
cg->frequency = freq;
CglLandP * landp = new CglLandP;
cg->cgl = landp;
cg->id = "Lift-and-Project";
methods_.push_back(cg);
}
s.options()->GetIntegerValue("reduce_and_split_cuts",freq,"bonmin.");
if (freq) {
Coin::SmartPtr<CuttingMethod> cg = new CuttingMethod;
cg->frequency = freq;
CglRedSplit * rands = new CglRedSplit;
cg->cgl = rands;
cg->id = "Reduce-and-Split";
methods_.push_back(cg);
}
}
