HermiteCounter.hh

Go to the documentation of this file.

#ifndef HERMITECOUNTER_HH
#define HERMITECOUNTER_HH
 
#include "casm/container/Counter.hh"
#include "casm/external/Eigen/Dense"
#include "casm/global/definitions.hh"
 
namespace CASM {
namespace xtal {
class HermiteCounter {
 public:
  typedef Eigen::VectorXi::Scalar value_type;
  typedef CASM::Index Index;
 
  // HermiteCounter() {};
 
  HermiteCounter(int init_determinant, int init_dim);
 
  // You probably will never need these. They're just here for testing more than
  // anything. Either way, they're safe to call.
  Index position() const;
  Eigen::VectorXi diagonal() const;
  // value_type low() const;
  // value_type high() const;
 
  Eigen::MatrixXi current() const;
 
  value_type determinant() const;
 
  Index dim() const;
 
  void reset_current();
 
  void next_determinant();
 
  void jump_to_determinant(value_type new_det);
 
  HermiteCounter &operator++();
 
  Eigen::MatrixXi operator()() const;
 
 private:
  Index m_pos;
 
  Eigen::VectorXi m_diagonal;
 
  EigenVectorXiCounter m_upper_tri;
 
  Index _increment_diagonal();
};
 
namespace HermiteCounter_impl {
HermiteCounter::Index _spill_factor(Eigen::VectorXi &diag,
                                    HermiteCounter::Index position,
                                    HermiteCounter::value_type attempt);
 
HermiteCounter::Index next_spill_position(Eigen::VectorXi &diag,
                                          HermiteCounter::Index position);
 
HermiteCounter::Index upper_size(HermiteCounter::Index init_dim);
 
EigenVectorXiCounter _upper_tri_counter(const Eigen::VectorXi &current_diag);
 
Eigen::MatrixXi _zip_matrix(const Eigen::VectorXi &current_diag,
                            const Eigen::VectorXi &current_upper_tri);
 
Eigen::MatrixXi _expand_dims_old(const Eigen::MatrixXi &hermit_mat,
                                 const Eigen::VectorXi &active_dims);
 
Eigen::MatrixXi _expand_dims(const Eigen::MatrixXi &H,
                             const Eigen::MatrixXi &G);
 
Eigen::VectorXi _canonical_unroll(const Eigen::MatrixXi &hermit_mat);
 
bool _canonical_compare(const Eigen::MatrixXi &H0, const Eigen::MatrixXi &H1);
}  // namespace HermiteCounter_impl
 
}  // namespace xtal
}  // namespace CASM
 
#endif