lemon/clp.h
author Balazs Dezso <deba@inf.elte.hu>
Thu, 24 Jun 2010 09:27:53 +0200
changeset 894 bb70ad62c95f
parent 540 9db62975c32b
child 746 e4554cd6b2bf
permissions -rw-r--r--
Fix critical bug in preflow (#372)

The wrong transition between the bound decrease and highest active
heuristics caused the bug. The last node chosen in bound decrease mode
is used in the first iteration in highest active mode.
alpar@461
     1
/* -*- mode: C++; indent-tabs-mode: nil; -*-
alpar@461
     2
 *
alpar@461
     3
 * This file is a part of LEMON, a generic C++ optimization library.
alpar@461
     4
 *
alpar@461
     5
 * Copyright (C) 2003-2008
alpar@461
     6
 * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
alpar@461
     7
 * (Egervary Research Group on Combinatorial Optimization, EGRES).
alpar@461
     8
 *
alpar@461
     9
 * Permission to use, modify and distribute this software is granted
alpar@461
    10
 * provided that this copyright notice appears in all copies. For
alpar@461
    11
 * precise terms see the accompanying LICENSE file.
alpar@461
    12
 *
alpar@461
    13
 * This software is provided "AS IS" with no warranty of any kind,
alpar@461
    14
 * express or implied, and with no claim as to its suitability for any
alpar@461
    15
 * purpose.
alpar@461
    16
 *
alpar@461
    17
 */
alpar@461
    18
alpar@461
    19
#ifndef LEMON_CLP_H
alpar@461
    20
#define LEMON_CLP_H
alpar@461
    21
alpar@461
    22
///\file
alpar@461
    23
///\brief Header of the LEMON-CLP lp solver interface.
alpar@461
    24
alpar@461
    25
#include <vector>
alpar@461
    26
#include <string>
alpar@461
    27
alpar@461
    28
#include <lemon/lp_base.h>
alpar@461
    29
alpar@461
    30
class ClpSimplex;
alpar@461
    31
alpar@461
    32
namespace lemon {
alpar@461
    33
alpar@461
    34
  /// \ingroup lp_group
alpar@461
    35
  ///
alpar@461
    36
  /// \brief Interface for the CLP solver
alpar@461
    37
  ///
alpar@461
    38
  /// This class implements an interface for the Clp LP solver.  The
alpar@461
    39
  /// Clp library is an object oriented lp solver library developed at
alpar@461
    40
  /// the IBM. The CLP is part of the COIN-OR package and it can be
alpar@461
    41
  /// used with Common Public License.
alpar@462
    42
  class ClpLp : public LpSolver {
alpar@461
    43
  protected:
alpar@461
    44
alpar@461
    45
    ClpSimplex* _prob;
alpar@461
    46
alpar@461
    47
    std::map<std::string, int> _col_names_ref;
alpar@461
    48
    std::map<std::string, int> _row_names_ref;
alpar@461
    49
alpar@461
    50
  public:
alpar@461
    51
alpar@461
    52
    /// \e
alpar@462
    53
    ClpLp();
alpar@461
    54
    /// \e
alpar@462
    55
    ClpLp(const ClpLp&);
alpar@461
    56
    /// \e
alpar@462
    57
    ~ClpLp();
alpar@461
    58
alpar@540
    59
    /// \e
alpar@540
    60
    virtual ClpLp* newSolver() const;
alpar@540
    61
    /// \e
alpar@540
    62
    virtual ClpLp* cloneSolver() const;
alpar@540
    63
alpar@461
    64
  protected:
alpar@461
    65
alpar@461
    66
    mutable double* _primal_ray;
alpar@461
    67
    mutable double* _dual_ray;
alpar@461
    68
alpar@461
    69
    void _init_temporals();
alpar@461
    70
    void _clear_temporals();
alpar@461
    71
alpar@461
    72
  protected:
alpar@461
    73
alpar@461
    74
    virtual const char* _solverName() const;
alpar@461
    75
alpar@461
    76
    virtual int _addCol();
alpar@461
    77
    virtual int _addRow();
alpar@461
    78
alpar@461
    79
    virtual void _eraseCol(int i);
alpar@461
    80
    virtual void _eraseRow(int i);
alpar@461
    81
alpar@461
    82
    virtual void _eraseColId(int i);
alpar@461
    83
    virtual void _eraseRowId(int i);
alpar@461
    84
alpar@461
    85
    virtual void _getColName(int col, std::string& name) const;
alpar@461
    86
    virtual void _setColName(int col, const std::string& name);
alpar@461
    87
    virtual int _colByName(const std::string& name) const;
alpar@461
    88
alpar@461
    89
    virtual void _getRowName(int row, std::string& name) const;
alpar@461
    90
    virtual void _setRowName(int row, const std::string& name);
alpar@461
    91
    virtual int _rowByName(const std::string& name) const;
alpar@461
    92
alpar@461
    93
    virtual void _setRowCoeffs(int i, ExprIterator b, ExprIterator e);
alpar@461
    94
    virtual void _getRowCoeffs(int i, InsertIterator b) const;
alpar@461
    95
alpar@461
    96
    virtual void _setColCoeffs(int i, ExprIterator b, ExprIterator e);
alpar@461
    97
    virtual void _getColCoeffs(int i, InsertIterator b) const;
alpar@461
    98
alpar@461
    99
    virtual void _setCoeff(int row, int col, Value value);
alpar@461
   100
    virtual Value _getCoeff(int row, int col) const;
alpar@461
   101
alpar@461
   102
    virtual void _setColLowerBound(int i, Value value);
alpar@461
   103
    virtual Value _getColLowerBound(int i) const;
alpar@461
   104
    virtual void _setColUpperBound(int i, Value value);
alpar@461
   105
    virtual Value _getColUpperBound(int i) const;
alpar@461
   106
alpar@461
   107
    virtual void _setRowLowerBound(int i, Value value);
alpar@461
   108
    virtual Value _getRowLowerBound(int i) const;
alpar@461
   109
    virtual void _setRowUpperBound(int i, Value value);
alpar@461
   110
    virtual Value _getRowUpperBound(int i) const;
alpar@461
   111
alpar@461
   112
    virtual void _setObjCoeffs(ExprIterator, ExprIterator);
alpar@461
   113
    virtual void _getObjCoeffs(InsertIterator) const;
alpar@461
   114
alpar@461
   115
    virtual void _setObjCoeff(int i, Value obj_coef);
alpar@461
   116
    virtual Value _getObjCoeff(int i) const;
alpar@461
   117
alpar@461
   118
    virtual void _setSense(Sense sense);
alpar@461
   119
    virtual Sense _getSense() const;
alpar@461
   120
alpar@461
   121
    virtual SolveExitStatus _solve();
alpar@461
   122
alpar@461
   123
    virtual Value _getPrimal(int i) const;
alpar@461
   124
    virtual Value _getDual(int i) const;
alpar@461
   125
alpar@461
   126
    virtual Value _getPrimalValue() const;
alpar@461
   127
alpar@461
   128
    virtual Value _getPrimalRay(int i) const;
alpar@461
   129
    virtual Value _getDualRay(int i) const;
alpar@461
   130
alpar@461
   131
    virtual VarStatus _getColStatus(int i) const;
alpar@461
   132
    virtual VarStatus _getRowStatus(int i) const;
alpar@461
   133
alpar@461
   134
    virtual ProblemType _getPrimalType() const;
alpar@461
   135
    virtual ProblemType _getDualType() const;
alpar@461
   136
alpar@461
   137
    virtual void _clear();
alpar@461
   138
deba@576
   139
    virtual void _messageLevel(MessageLevel);
deba@576
   140
    
alpar@461
   141
  public:
alpar@461
   142
alpar@461
   143
    ///Solves LP with primal simplex method.
alpar@461
   144
    SolveExitStatus solvePrimal();
alpar@461
   145
alpar@461
   146
    ///Solves LP with dual simplex method.
alpar@461
   147
    SolveExitStatus solveDual();
alpar@461
   148
alpar@461
   149
    ///Solves LP with barrier method.
alpar@461
   150
    SolveExitStatus solveBarrier();
alpar@461
   151
alpar@461
   152
    ///Returns the constraint identifier understood by CLP.
alpar@461
   153
    int clpRow(Row r) const { return rows(id(r)); }
alpar@461
   154
alpar@461
   155
    ///Returns the variable identifier understood by CLP.
alpar@461
   156
    int clpCol(Col c) const { return cols(id(c)); }
alpar@461
   157
alpar@461
   158
  };
alpar@461
   159
alpar@461
   160
} //END OF NAMESPACE LEMON
alpar@461
   161
alpar@461
   162
#endif //LEMON_CLP_H
alpar@461
   163