RhodeCode

/* -*- mode: C++; indent-tabs-mode: nil; -*-

 *

 * This file is a part of LEMON, a generic C++ optimization library.

 *

 * Copyright (C) 2003-2008

 * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport

 * (Egervary Research Group on Combinatorial Optimization, EGRES).

 *

 * Permission to use, modify and distribute this software is granted

 * provided that this copyright notice appears in all copies. For

 * precise terms see the accompanying LICENSE file.

 *

 * This software is provided "AS IS" with no warranty of any kind,

 * express or implied, and with no claim as to its suitability for any

 * purpose.

 *

 */

#ifndef LEMON_GLPK_H

#define LEMON_GLPK_H

///\file

///\brief Header of the LEMON-GLPK lp solver interface.

///\ingroup lp_group

#include <lemon/lp_base.h>

// forward declaration

#ifndef _GLP_PROB

#define _GLP_PROB

typedef struct { double _prob; } glp_prob;

/* LP/MIP problem object */

#endif

namespace lemon {

  /// \brief Base interface for the GLPK LP and MIP solver

  ///

  /// This class implements the common interface of the GLPK LP and MIP solver.

  /// \ingroup lp_group

  class GlpkBase : virtual public LpBase {

  protected:

    typedef glp_prob LPX;

    glp_prob* lp;

    GlpkBase();

    GlpkBase(const GlpkBase&);

    virtual ~GlpkBase();

  protected:

    virtual int _addCol();

    virtual int _addRow();

    virtual void _eraseCol(int i);

    virtual void _eraseRow(int i);

    virtual void _eraseColId(int i);

    virtual void _eraseRowId(int i);

    virtual void _getColName(int col, std::string& name) const;

    virtual void _setColName(int col, const std::string& name);

    virtual int _colByName(const std::string& name) const;

    virtual void _getRowName(int row, std::string& name) const;

    virtual void _setRowName(int row, const std::string& name);

    virtual int _rowByName(const std::string& name) const;

    virtual void _setRowCoeffs(int i, ExprIterator b, ExprIterator e);

    virtual void _getRowCoeffs(int i, InsertIterator b) const;

    virtual void _setColCoeffs(int i, ExprIterator b, ExprIterator e);

    virtual void _getColCoeffs(int i, InsertIterator b) const;

    virtual void _setCoeff(int row, int col, Value value);

    virtual Value _getCoeff(int row, int col) const;

    virtual void _setColLowerBound(int i, Value value);

    virtual Value _getColLowerBound(int i) const;

    virtual void _setColUpperBound(int i, Value value);

    virtual Value _getColUpperBound(int i) const;

    virtual void _setRowLowerBound(int i, Value value);

    virtual Value _getRowLowerBound(int i) const;

    virtual void _setRowUpperBound(int i, Value value);

    virtual Value _getRowUpperBound(int i) const;

    virtual void _setObjCoeffs(ExprIterator b, ExprIterator e);

    virtual void _getObjCoeffs(InsertIterator b) const;

    virtual void _setObjCoeff(int i, Value obj_coef);

    virtual Value _getObjCoeff(int i) const;

    virtual void _setSense(Sense);

    virtual Sense _getSense() const;

    virtual void _clear();

  public:

    ///Pointer to the underlying GLPK data structure.

    LPX *lpx() {return lp;}

    ///Const pointer to the underlying GLPK data structure.

    const LPX *lpx() const {return lp;}

    ///Returns the constraint identifier understood by GLPK.

    int lpxRow(Row r) const { return rows(id(r)); }

    ///Returns the variable identifier understood by GLPK.

    int lpxCol(Col c) const { return cols(id(c)); }

  };

  /// \brief Interface for the GLPK LP solver

  ///

  /// This class implements an interface for the GLPK LP solver.

  ///\ingroup lp_group

  class GlpkLp : public LpSolver, public GlpkBase {

  public:

    ///\e

    GlpkLp();

    ///\e

    GlpkLp(const GlpkLp&);

    ///\e

    virtual GlpkLp* cloneSolver() const;

    ///\e

    virtual GlpkLp* newSolver() const;

  private:

    mutable std::vector<double> _primal_ray;

    mutable std::vector<double> _dual_ray;

    void _clear_temporals();

  protected:

    virtual const char* _solverName() const;

    virtual SolveExitStatus _solve();

    virtual Value _getPrimal(int i) const;

    virtual Value _getDual(int i) const;

    virtual Value _getPrimalValue() const;

    virtual VarStatus _getColStatus(int i) const;

    virtual VarStatus _getRowStatus(int i) const;

    virtual Value _getPrimalRay(int i) const;

    virtual Value _getDualRay(int i) const;

    virtual ProblemType _getPrimalType() const;

    virtual ProblemType _getDualType() const;

  public:

    ///Solve with primal simplex

    SolveExitStatus solvePrimal();

    ///Solve with dual simplex

    SolveExitStatus solveDual();

    ///Turns on or off the presolver

    ///Turns on (\c b is \c true) or off (\c b is \c false) the presolver

    ///

    ///The presolver is off by default.

    void presolver(bool b);

    ///Enum for \c messageLevel() parameter

    enum MessageLevel {

      /// no output (default value)

      MESSAGE_NO_OUTPUT = 0,

      /// error messages only

      MESSAGE_ERROR_MESSAGE = 1,

      /// normal output

      MESSAGE_NORMAL_OUTPUT = 2,

      /// full output (includes informational messages)

      MESSAGE_FULL_OUTPUT = 3

    };

  private:

    MessageLevel _message_level;

  public:

    ///Set the verbosity of the messages

    ///Set the verbosity of the messages

    ///

    ///\param m is the level of the messages output by the solver routines.

    void messageLevel(MessageLevel m);

  };

  /// \brief Interface for the GLPK MIP solver

  ///

  /// This class implements an interface for the GLPK MIP solver.

  ///\ingroup lp_group

  class GlpkMip : public MipSolver, public GlpkBase {

  public:

    ///\e

    GlpkMip();

    ///\e

    GlpkMip(const GlpkMip&);

    virtual GlpkMip* cloneSolver() const;

    virtual GlpkMip* newSolver() const;

  protected:

    virtual const char* _solverName() const;

    virtual ColTypes _getColType(int col) const;

    virtual void _setColType(int col, ColTypes col_type);

    virtual SolveExitStatus _solve();

    virtual ProblemType _getType() const;

    virtual Value _getSol(int i) const;

    virtual Value _getSolValue() const;

    ///Enum for \c messageLevel() parameter

    enum MessageLevel {

      /// no output (default value)

      MESSAGE_NO_OUTPUT = 0,

      /// error messages only

      MESSAGE_ERROR_MESSAGE = 1,

      /// normal output

      MESSAGE_NORMAL_OUTPUT = 2,

      /// full output (includes informational messages)

      MESSAGE_FULL_OUTPUT = 3

    };

  private:

    MessageLevel _message_level;

  public:

    ///Set the verbosity of the messages

    ///Set the verbosity of the messages

    ///

    ///\param m is the level of the messages output by the solver routines.

    void messageLevel(MessageLevel m);

  };

} //END OF NAMESPACE LEMON

#endif //LEMON_GLPK_H

/* -*- mode: C++; indent-tabs-mode: nil; -*-

 *

 * This file is a part of LEMON, a generic C++ optimization library.

 *

 * Copyright (C) 2003-2008

 * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport

 * (Egervary Research Group on Combinatorial Optimization, EGRES).

 *

 * Permission to use, modify and distribute this software is granted

 * provided that this copyright notice appears in all copies. For

 * precise terms see the accompanying LICENSE file.

 *

 * This software is provided "AS IS" with no warranty of any kind,

 * express or implied, and with no claim as to its suitability for any

 * purpose.

 *

 */

#ifndef LEMON_LP_SKELETON_H

#define LEMON_LP_SKELETON_H

#include <lemon/lp_base.h>

///\file

///\brief Skeleton file to implement LP/MIP solver interfaces

///  

///The classes in this file do nothing, but they can serve as skeletons when

///implementing an interface to new solvers.

namespace lemon {

  ///A skeleton class to implement LP/MIP solver base interface

  ///This class does nothing, but it can serve as a skeleton when

  ///implementing an interface to new solvers.

  class SkeletonSolverBase : public virtual LpBase {

    int col_num,row_num;

  protected:

    SkeletonSolverBase()

      : col_num(-1), row_num(-1) {}

    /// \e

    virtual int _addCol();

    /// \e

    virtual int _addRow();

    /// \e

    virtual void _eraseCol(int i);

    /// \e

    virtual void _eraseRow(int i);

    /// \e

    virtual void _getColName(int col, std::string& name) const;

    /// \e

    virtual void _setColName(int col, const std::string& name);

    /// \e

    virtual int _colByName(const std::string& name) const;

    /// \e

    virtual void _getRowName(int row, std::string& name) const;

    /// \e

    virtual void _setRowName(int row, const std::string& name);

    /// \e

    virtual int _rowByName(const std::string& name) const;

    /// \e

    virtual void _setRowCoeffs(int i, ExprIterator b, ExprIterator e);

    /// \e

    virtual void _getRowCoeffs(int i, InsertIterator b) const;

    /// \e

    virtual void _setColCoeffs(int i, ExprIterator b, ExprIterator e);

    /// \e

    virtual void _getColCoeffs(int i, InsertIterator b) const;

    /// Set one element of the coefficient matrix

    virtual void _setCoeff(int row, int col, Value value);

    /// Get one element of the coefficient matrix

    virtual Value _getCoeff(int row, int col) const;

    /// The lower bound of a variable (column) have to be given by an

    /// extended number of type Value, i.e. a finite number of type

    /// Value or -\ref INF.

    virtual void _setColLowerBound(int i, Value value);

    /// \e

    /// The lower bound of a variable (column) is an

    /// extended number of type Value, i.e. a finite number of type

    /// Value or -\ref INF.

    virtual Value _getColLowerBound(int i) const;

    /// The upper bound of a variable (column) have to be given by an

    /// extended number of type Value, i.e. a finite number of type

    /// Value or \ref INF.

    virtual void _setColUpperBound(int i, Value value);

    /// \e

    /// The upper bound of a variable (column) is an

    /// extended number of type Value, i.e. a finite number of type

    /// Value or \ref INF.

    virtual Value _getColUpperBound(int i) const;

    /// The lower bound of a constraint (row) have to be given by an

    /// extended number of type Value, i.e. a finite number of type

    /// Value or -\ref INF.

    virtual void _setRowLowerBound(int i, Value value);

    /// \e

    /// The lower bound of a constraint (row) is an

    /// extended number of type Value, i.e. a finite number of type

    /// Value or -\ref INF.

    virtual Value _getRowLowerBound(int i) const;

    /// The upper bound of a constraint (row) have to be given by an

    /// extended number of type Value, i.e. a finite number of type

    /// Value or \ref INF.

    virtual void _setRowUpperBound(int i, Value value);

    /// \e

    /// The upper bound of a constraint (row) is an

    /// extended number of type Value, i.e. a finite number of type

    /// Value or \ref INF.

    virtual Value _getRowUpperBound(int i) const;

    /// \e

    virtual void _setObjCoeffs(ExprIterator b, ExprIterator e);

    /// \e

    virtual void _getObjCoeffs(InsertIterator b) const;

    /// \e

    virtual void _setObjCoeff(int i, Value obj_coef);

    /// \e

    virtual Value _getObjCoeff(int i) const;

    ///\e

    virtual void _setSense(Sense);

    ///\e

    virtual Sense _getSense() const;

    ///\e

    virtual void _clear();

  };

  /// \brief Skeleton class for an LP solver interface

  ///

  ///This class does nothing, but it can serve as a skeleton when

  ///implementing an interface to new solvers.

  ///\ingroup lp_group

  class LpSkeleton : public LpSolver, public SkeletonSolverBase {

  public:

    ///\e

    LpSkeleton() : LpSolver(), SkeletonSolverBase() {}

    ///\e

    virtual LpSkeleton* newSolver() const;

    ///\e

    virtual LpSkeleton* cloneSolver() const;

  protected:

    ///\e

    virtual SolveExitStatus _solve();

    ///\e

    virtual Value _getPrimal(int i) const;

    ///\e

    virtual Value _getDual(int i) const;

    ///\e

    virtual Value _getPrimalValue() const;

    ///\e

    virtual Value _getPrimalRay(int i) const;

    ///\e

    virtual Value _getDualRay(int i) const;

    ///\e

    virtual ProblemType _getPrimalType() const;

    ///\e

    virtual ProblemType _getDualType() const;

    ///\e

    virtual VarStatus _getColStatus(int i) const;

    ///\e

    virtual VarStatus _getRowStatus(int i) const;

    ///\e

    virtual const char* _solverName() const;

  };

  /// \brief Skeleton class for a MIP solver interface

  ///

  ///This class does nothing, but it can serve as a skeleton when

  ///implementing an interface to new solvers.

  ///\ingroup lp_group

  class MipSkeleton : public MipSolver, public SkeletonSolverBase {

  public:

    ///\e

    MipSkeleton() : MipSolver(), SkeletonSolverBase() {}

    ///\e

    virtual MipSkeleton* newSolver() const;

    ///\e

    virtual MipSkeleton* cloneSolver() const;

  protected:

    ///\e

    virtual SolveExitStatus _solve();

    ///\e

    virtual Value _getSol(int i) const;

    ///\e

    virtual Value _getSolValue() const;

    ///\e

    virtual ProblemType _getType() const;

    ///\e

    virtual const char* _solverName() const;

  };

} //namespace lemon

#endif