BonSubMipSolver.hpp

Go to the documentation of this file.

// (C) Copyright International Business Machines (IBM) 2006
// All Rights Reserved.
// This code is published under the Eclipse Public License.
//
// Authors :
// P. Bonami, International Business Machines
//
// Date :  12/07/2006


// Code separated from BonOaDecBase to try to clarify OAs
#ifndef BonSubMipSolver_HPP
#define BonSubMipSolver_HPP
#include "IpSmartPtr.hpp"
#include <string>
/* forward declarations.*/
class OsiSolverInterface;
class OsiClpSolverInterface;
class OsiCpxSolverInterface;
class CbcStrategy;
class CbcStrategyDefault;

#include "OsiCuts.hpp"

namespace Bonmin {
    class RegisteredOptions;
    class BabSetupBase; 
    class SubMipSolver
    {
    public:
      enum MILP_solve_strategy{
         FindGoodSolution,
         GetOptimum};
      SubMipSolver(BabSetupBase &b, const std::string &prefix);

      SubMipSolver(const SubMipSolver &copy);

      ~SubMipSolver();

      void setLpSolver(OsiSolverInterface * lp);

      void setStrategy(CbcStrategyDefault * strategy);

      const double * getLastSolution()
      {
        return integerSolution_;
      }

      double getLowerBound()
      {
        return lowBound_;
      }

      void solve(double cutoff,
          int loglevel,
          double maxTime){
         if(milp_strat_ == FindGoodSolution){
            find_good_sol(cutoff, loglevel, maxTime);
         }
         else
            optimize(cutoff, loglevel, maxTime);
      }
 

      void find_good_sol(double cutoff,
          int loglevel,
          double maxTime);

      void optimize(double cutoff,
          int loglevel,
          double maxTime);

      void optimize_with_lazy_constraints(double cutoff,
          int loglevel,
          double maxTime, const OsiCuts & cs);

      inline double lowBound()
      {
        return lowBound_;
      }

      inline bool optimal()
      {
        return optimal_;
      }

      inline int nodeCount()
      {
        return nodeCount_;
      }

      inline int iterationCount()
      {
        return iterationCount_;
      }


      OsiSolverInterface * solver();

     static void registerOptions(Ipopt::SmartPtr<Bonmin::RegisteredOptions> roptions);
    private:
      OsiClpSolverInterface *clp_;
      OsiCpxSolverInterface * cpx_;
      double lowBound_;
      bool optimal_;
      double * integerSolution_;
      CbcStrategyDefault * strategy_;
      int nodeCount_;
      int iterationCount_;
      MILP_solve_strategy milp_strat_;
      double gap_tol_;
      bool ownClp_;
    };

}

#endif