QsoptexLpSolver.h

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

class QsoptexLpSolver final : public LpSolver {
 public:
  explicit QsoptexLpSolver(Config config = {}, const std::string& class_name = "QsoptexLpSolver");
  ~QsoptexLpSolver() override;
 
  [[nodiscard]] int num_columns() const override;
  [[nodiscard]] int num_rows() const override;
 
  [[nodiscard]] Column column(int column_idx) const override;
  [[nodiscard]] Row row(int row_idx) const 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;

  void UpdateStats(unsigned int precision);
 
#ifndef NDEBUG
  void Dump() final;
#endif
 
 private:
  LpResult SolveCore() override;
 
  template <TypedIterable<std::pair<const Variable, mpq_class>> T>
  int AddRow(const T& addends, char sense, const mpq_class& rhs);
  template <TypedIterable<std::pair<const Variable, mpq_class>> T>
  void SetRowCoeff(int row, const T& literal_monomials);
  void SetVarCoeff(int row, const Variable& var, const mpq_class& value) const;

  void UpdateFeasible();
 
 protected:
  void EnsureSenseCore() override;
 
 private:
#if 0
  void UpdateInfeasible();
#endif
 
  mpq_QSprob qsx_;  
  QSbasis basis_;   
 
  qsopt_ex::MpqArray ray_;  
  qsopt_ex::MpqArray x_;    
};
 
}  // namespace delpi