simplex_internal.hpp

Go to the documentation of this file.

#pragma once
 
#include <Eigen/Dense>
#include <string>
#include <vector>
 
namespace EigenOpt {
 
namespace simplex {
 
namespace internal {
 
 
struct VariableDomain {
  bool non_negative = false; 
  bool non_positive = false; 
  int idx = -1; 
};
 
 
 
template<class Scalar>
bool deduce_variables_domains(
  const Eigen::Matrix<Scalar,Eigen::Dynamic,Eigen::Dynamic>& C,
  const Eigen::Matrix<Scalar,Eigen::Dynamic,1>& d,
  const Scalar& small_number,
  std::vector<VariableDomain>& domains,
  std::string& halt_reason
);
 
 
 
template<class Scalar>
Eigen::Matrix<Scalar,Eigen::Dynamic,Eigen::Dynamic> transformation_matrix_from_domains(
  const std::vector<VariableDomain>& domains
);
 
 
 
template<class Scalar>
bool transformation_matrix(
  const Eigen::Matrix<Scalar,Eigen::Dynamic,Eigen::Dynamic>& C,
  const Eigen::Matrix<Scalar,Eigen::Dynamic,1>& d,
  const Scalar& small_number,
  Eigen::Matrix<Scalar, Eigen::Dynamic, Eigen::Dynamic>& T,
  std::string& halt_reason
);
 
 
 
template<class Scalar>
void create_tableau(
  const Eigen::Matrix<Scalar, Eigen::Dynamic, Eigen::Dynamic>& C,
  const Eigen::Matrix<Scalar, Eigen::Dynamic, 1>& d,
  Eigen::Matrix<Scalar, Eigen::Dynamic, Eigen::Dynamic>& tableau,
  std::vector<int>& basic_variables
);
 
 
 
template<class Scalar>
void eliminate_objective(
  Eigen::Matrix<Scalar, Eigen::Dynamic, Eigen::Dynamic>& tableau,
  const std::vector<int>& basic_variables
);
 
 
 
template<class Scalar>
void pivot(
  Eigen::Matrix<Scalar,Eigen::Dynamic,Eigen::Dynamic>& tableau,
  int entering_variable,
  int leaving_variable
);
 
 
 
template<class Scalar>
bool simplex(
  Eigen::Matrix<Scalar,Eigen::Dynamic,Eigen::Dynamic>& tableau,
  std::vector<int>& basic_variables,
  const Scalar& small_number,
  std::string& halt_reason
);
 
 
 
template<class Scalar>
bool two_steps_method(
  const Eigen::Matrix<Scalar,Eigen::Dynamic,1>& objective,
  Eigen::Matrix<Scalar,Eigen::Dynamic,Eigen::Dynamic>& tableau,
  std::vector<int>& basic_variables,
  unsigned int na,
  const Scalar& small_number,
  std::string& halt_reason
);
 
 
 
template<class Scalar>
bool penalty_method(
  const Eigen::Matrix<Scalar,Eigen::Dynamic,1>& objective,
  Eigen::Matrix<Scalar,Eigen::Dynamic,Eigen::Dynamic>& tableau,
  std::vector<int>& basic_variables,
  unsigned int na,
  const Scalar& small_number,
  const Scalar& large_number,
  std::string& halt_reason
);
 
} // namespace internal
 
} // namespace simplex
 
} // namespace EigenOpt
 
#include "simplex_internal.hxx"