OsiDylpSolverInterface.hpp

#ifdef COIN_USE_DYLP
#ifndef OsiDylpSolverInterface_H
#define OsiDylpSolverInterface_H

/*
  sccs: %W%     %G%
  cvs: $Id: OsiDylpSolverInterface.hpp,v 1.12 2003/10/23 17:33:47 lou-oss Exp $
*/

#include <CoinPackedMatrix.hpp>
#include <OsiSolverInterface.hpp>
#include <CoinWarmStart.hpp>
#include <CoinMessageHandler.hpp>
#include <CoinMpsIO.hpp>

#define DYLP_INTERNAL
extern "C" {
#include "dylp.h"
}

class OsiDylpSolverInterface: public OsiSolverInterface

{ friend void OsiDylpSolverInterfaceUnitTest(const std::string &mpsDir,
                                             const std::string &netLibDir) ;

/*
  Consult the COIN OSI documentation or relevant source code for details
  not covered here. Supported functions are listed first, followed by
  unsupported functions.
*/

public:


  OsiDylpSolverInterface() ;

  OsiDylpSolverInterface(const OsiDylpSolverInterface &src) ;

  OsiSolverInterface* clone(bool copyData = true) const ;

  OsiDylpSolverInterface &operator=(const OsiDylpSolverInterface &rhs) ;

  ~OsiDylpSolverInterface() ;

  void reset() ;



  int readMps(const char *filename, const char *extension = "mps") ;

  void writeMps(const char *basename,
                const char *extension = "mps",
                double objsense = 0.0) const ;

  void loadProblem(const CoinPackedMatrix &matrix,
                   const double *collb, const double *colub, const double *obj,
                   const char *rowsen, const double *rowrhs,
                   const double *rowrng) ;

  void loadProblem(const CoinPackedMatrix &matrix,
                   const double *collb, const double *colub, const double *obj,
                   const double *rowlb, const double *rowub) ;
  
  void loadProblem(const int colcnt, const int rowcnt,
               const int *start, const int *index, const double *value,
               const double *collb, const double *colub, const double *obj,
               const char *sense, const double *rhsin, const double *range) ;

  void loadProblem(const int colcnt, const int rowcnt,
                   const int *start, const int *index, const double *value,
                   const double *collb, const double *colub, const double *obj,
                   const double *row_lower, const double *row_upper) ;

  void assignProblem(CoinPackedMatrix *&matrix,
                     double *&collb, double *&colub, double *&obj, 
                     char *&rowsen, double *&rowrhs, double *&rowrng) ;

  void assignProblem(CoinPackedMatrix *&matrix,
                     double *&collb, double *&colub, double *&obj,
                     double *&rowlb, double *&rowub) ;



  int getNumCols() const ;

  int getNumRows() const ;

  int getNumElements() const ;

  const double* getColLower() const ;

  const double* getColUpper() const ;

  bool isContinuous(int colIndex) const ;

  bool isBinary(int colIndex) const ;

  bool isIntegerNonBinary (int colIndex) const ;

  bool isInteger (int colIndex) const ;

  const char* getRowSense() const ;

  const double* getRightHandSide() const ;

  const double* getRowRange() const ;

  const double* getRowLower() const ;

  const double* getRowUpper() const ;

  const double* getObjCoefficients() const ;

  double getObjSense() const ;

  const CoinPackedMatrix *getMatrixByRow() const ;

  const CoinPackedMatrix *getMatrixByCol() const ;


  void setContinuous(int index) ;

  void setContinuous(const int *indices, int len) ;

  void setInteger(int index) ;

  void setInteger(const int *indices, int len) ;

  void setColLower(int index, double value) ;

  void setColUpper(int index, double value) ;

  void setRowLower(int index, double value) ;

  void setRowUpper(int index, double value) ;

  void setRowType(int index, char rowsen, double rowrhs, double rowrng) ;

  void setObjCoeff (int index, double value) ;

  void setObjSense(double sense) ;

  void setColSolution(const double *colsol) ;

  void setRowPrice(const double*) ;

  void addCol(const CoinPackedVectorBase &vec,
              const double collb, const double colub, const double obj) ;

  void deleteCols(const int num, const int *colIndices) ;

  void addRow(const CoinPackedVectorBase &row,
              const double rowlb, const double rowub) ;

  void addRow(const CoinPackedVectorBase &row,
              const char rowsen, const double rowrhs, const double rowrng) ;

  void deleteRows(const int num, const int *rowIndices) ;

  void applyRowCut(const OsiRowCut &cut) ;

  void applyColCut(const OsiColCut &cut) ;


  void initialSolve() ;

  CoinWarmStart *getEmptyWarmStart () const ;

  CoinWarmStart *getWarmStart() const ;

  bool setWarmStart(const CoinWarmStart *warmStart) ;

  void resolve() ;

  void markHotStart() ;

  void solveFromHotStart() ;

  void unmarkHotStart() ;



  bool isAbandoned() const ;

  bool isProvenOptimal() const ;

  bool isProvenPrimalInfeasible() const ;

  bool isProvenDualInfeasible() const ;

  bool isIterationLimitReached() const ;

  int getIterationCount() const ;

  bool isPrimalObjectiveLimitReached() const ;

  bool isDualObjectiveLimitReached() const ;



  double getInfinity() const ;

  bool setIntParam(OsiIntParam key, int value) ;

  bool setDblParam(OsiDblParam key, double value) ;

  bool setStrParam(OsiStrParam key, const std::string& value) ;

  bool setHintParam(OsiHintParam key, bool sense = true,
                    OsiHintStrength strength = OsiHintTry, void *info = 0) ;

  bool getIntParam(OsiIntParam key, int &value) const ;

  bool getDblParam(OsiDblParam key, double &value) const ;

  bool getStrParam(OsiStrParam key, std::string &value) const ;

  /* For overload resolution with OSI::getHintParam functions. */

  using OsiSolverInterface::getHintParam ;

  bool getHintParam(OsiHintParam key, bool &sense,
                    OsiHintStrength &strength, void *&info) const ;

  inline void newLanguage(CoinMessages::Language language)
  { setOsiDylpMessages(language) ; } ;

  inline void setLanguage(CoinMessages::Language language)
  { setOsiDylpMessages(language) ; } ;



  const double* getColSolution() const ;

  const double* getRowPrice() const ;

  const double* getReducedCost() const ;

  const double* getRowActivity() const ;

  double getObjValue() const ;



  void activateRowCutDebugger (const char * modelName) ;



  void dylp_controlfile(const char* name, const bool silent,
                        const bool mustexist = true) ;

  void dylp_logfile(const char* name, bool echo = false) ;

  void setOsiDylpMessages(CoinMessages::Language local_language) ;




  void branchAndBound() ;

  std::vector<double *> getDualRays(int) const ;

  std::vector<double *> getPrimalRays(int) const ;

/*
  Leave them visible to the public for the nonce, until a better programmatic
  interface is available. Initialized by the constructor.
*/

  lpopts_struct *initialSolveOptions,*resolveOptions ;
  lptols_struct* tolerances ;


private:

/*
  Private implementation state and helper functions. If you're contemplating
  using any of these, you should have a look at the code.
  See OsiDylpSolverInterface.ccp for descriptions.
*/ 

  consys_struct* consys ;
  lpprob_struct* lpprob ;
  lpstats_struct* statistics ;



  static int reference_count ;
  static bool basis_ready ;
  static OsiDylpSolverInterface *dylp_owner ;




  ioid local_logchn ;

  bool initial_gtxecho ;
  bool resolve_gtxecho ;

  lpret_enum lp_retval ;

  double obj_sense ;

  const std::string solvername ;

  void *info_[OsiLastHintParam] ;

  bool mps_debug ;

  CoinWarmStart *hotstart_fallback ;

  CoinWarmStart *activeBasis ;

  bool activeIsModified ;





  mutable double _objval ;
  mutable double* _col_x ;
  mutable double* _col_obj ;
  mutable double* _col_cbar ;
  mutable double* _row_lhs ;
  mutable double* _row_lower ;
  mutable double* _row_price ;
  mutable double* _row_range ;
  mutable double* _row_rhs ;
  mutable char* _row_sense ;
  mutable double* _row_upper ;
  mutable CoinPackedMatrix* _matrix_by_row ;
  mutable CoinPackedMatrix* _matrix_by_col ;


  void construct_lpprob() ;
  void construct_options() ;
  void construct_consys(int cols, int rows) ;
  void dylp_ioinit() ;
  void gen_rowparms(int rowcnt,
                    double *rhs, double *rhslow, contyp_enum *ctyp,
                    const double *rowlb, const double *rowub) ;
  void gen_rowparms(int rowcnt,
               double *rhs, double *rhslow, contyp_enum *ctyp,
               const char *sense, const double *rhsin, const double *range) ;
  void load_problem(const CoinMpsIO &mps) ;
  void load_problem(const CoinPackedMatrix &matrix,
         const double* col_lower, const double* col_upper, const double* obj,
         const contyp_enum *ctyp, const double* rhs, const double* rhslow) ;
  void load_problem (const int colcnt, const int rowcnt,
         const int *start, const int *index, const double *value,
         const double* col_lower, const double* col_upper, const double* obj,
         const contyp_enum *ctyp, const double* rhs, const double* rhslow) ;

  void destruct_col_cache() ;
  void destruct_row_cache() ;
  void destruct_cache() ;
  void destruct_problem(bool preserve_interface) ;
  void detach_dylp() ;


/*
  There are separate groups for member and static methods so that doxygen
  won't promote the group to the top level.
*/
  
  void add_col(const CoinPackedVectorBase& coin_coli,
    vartyp_enum vtypi,double vlbi, double vubi, double obji) ;
  void add_row(const CoinPackedVectorBase& coin_rowi, 
    char clazzi, contyp_enum ctypi, double rhsi, double rhslowi) ;
  void worst_case_primal() ;
  void calc_objval() ;
  void unimp_hint(bool dylpSense, bool hintSense,
                 OsiHintStrength hintStrength, const char *msgString) ;


  static contyp_enum bound_to_type(double lower, double upper) ;
  static contyp_enum sense_to_type(char type) ;
  static char type_to_sense(contyp_enum type) ;
  static void gen_rowiparms(contyp_enum* ctypi, double* rhsi, double* rhslowi, 
                            char sensei, double rhsini, double rangei) ;
  static void gen_rowiparms(contyp_enum* ctypi, double* rhsi, double* rhslowi, 
                            double rowlbi, double rowubi) ;

  template<class T> static void copy(const T* src, T* dst, int n) ;
  template<class T> static T* copy(const T* src, int n) ;
  template<class T> static T* copy(const T* src) ;
/*
  Specializations for more complicated structures.
*/
  static basis_struct* copy_basis(const basis_struct* src, int dstsze) ;
  static void copy_basis(const basis_struct* src, basis_struct* dst) ;
  static lpprob_struct* copy_lpprob(const lpprob_struct* src) ;

#ifndef _MSC_VER

  template<class T> static void assert_same(const T& t1, const T& t2,
                                            bool exact) ;
  template<class T> static void assert_same(const T* t1, const T* t2,
                                            int n, bool exact) ;

  static void assert_same(double d1, double d2, bool exact) ;

  static void assert_same(const basis_struct& b1, 
                          const basis_struct& b2, bool exact) ;
  static void assert_same(const consys_struct& c1, const 
                          consys_struct& c2, bool exact) ;
  static void assert_same(const conbnd_struct& c1, const 
                          conbnd_struct& c2, bool exact) ;
  static void assert_same(const lpprob_struct& l1, 
                          const lpprob_struct& l2, bool exact) ;
  static void assert_same(const OsiDylpSolverInterface& o1, 
                          const OsiDylpSolverInterface& o2, bool exact) ;
#endif

  template<class T> static T* idx_vec(T* data) ;
  static int idx(int i) ;
  template<class T> static T* inv_vec(T* data) ;
  static int inv(int i) ;

  static pkvec_struct* packed_vector(
    const CoinShallowPackedVector vector, int dimension) ;
  static void packed_vector(
    const CoinShallowPackedVector vector, int dimension, pkvec_struct *dst) ;

  static std::string make_filename(const char *filename, 
                                   const char *ext1, const char *ext2) ;

} ;


/*
  OsiDylpSolverInterfaceTest.cpp
*/

void OsiDylpSolverInterfaceUnitTest(const std::string & mpsDir) ;

#endif // OsiDylpSolverInterface_H
#endif // COIN_USE_DYLP

---