SoplexLpSolver.h

#pragma once
 
#ifndef DELPI_ENABLED_SOPLEX
#error SoPlex is not enabled. Please enable it by adding "--//tools:enable_soplex" to the bazel command.
#endif
 
#include <string>
#include <utility>
#include <vector>
 
#include "delpi/libs/gmp.h"
#include "delpi/libs/soplex.h"
#include "delpi/solver/LpSolver.h"
#include "delpi/symbolic/Expression.h"
#include "delpi/symbolic/Variable.h"
#include "delpi/util/concepts.h"
 
namespace delpi {

class SoplexLpSolver final : public LpSolver {
 public:
  explicit SoplexLpSolver(Config config = {}, const std::string& class_name = "SoplexLpSolver");
 
  [[nodiscard]] int num_columns() const override;
  [[nodiscard]] int num_rows() const override;
 
  [[nodiscard]] Column column(ColumnIndex column_idx) const override;
  [[nodiscard]] Row row(RowIndex row_idx) const override;
  void ReserveColumns(int num_columns) override;
  void ReserveRows(int num_rows) override;
  ColumnIndex AddColumn(const Variable& var, const mpq_class& obj, const mpq_class& lb, const mpq_class& ub) override;
  RowIndex AddRow(const std::vector<Expression::Addend>& addends, const mpq_class& lb, const mpq_class& ub) override;
  RowIndex AddRow(const Expression::Addends& lhs, FormulaKind sense, const mpq_class& rhs) override;
  void SetBound(Variable var, const mpq_class& lb, const mpq_class& ub) override;
  void SetCoefficient(RowIndex row, ColumnIndex column, const mpq_class& value) override;
  void SetObjective(int column, const mpq_class& value) override;
 
#ifndef NDEBUG
  void Dump() override;
#endif
 
 private:
  LpResult SolveCore() override;
  template <TypedIterable<std::pair<const Variable, mpq_class>> T>
  soplex::DSVectorRational ParseRowCoeff(const T& literal_monomials);
  void SetVarCoeff(soplex::DSVectorRational& coeffs, const Variable& var, const mpq_class& value) const;

  void UpdateFeasible(const soplex::Rational& max_violation);
 
 protected:
  void EnsureSenseCore() override;
 
 private:
#if 0
  void UpdateInfeasible();
#endif
 
  bool consolidated_;  
 
  soplex::SoPlex spx_;  
 
  soplex::LPColSetRational spx_cols_;  
  soplex::LPRowSetRational spx_rows_;  
 
  soplex::Rational rninfinity_;  
  soplex::Rational rinfinity_;   
};
 
}  // namespace delpi