DoFIsEquivalent_impl.hh

Go to the documentation of this file.

#include "casm/basis_set/DoFIsEquivalent.hh"
 
namespace CASM {
 
template <typename OccType>
bool OccupantDoFIsEquivalent<OccType>::operator()(
    OccDoFType const &_other) const {
  if (_other.size() != m_dof.size()) return false;
  Index j;
  for (Index i = 0; i < m_dof.size(); ++i) {
    for (j = 0; j < _other.size(); ++j) {
      if (m_dof[i].identical(_other[j], m_tol)) {
        m_P.set(i) = j;
        break;
      }
    }
    if (j == _other.size()) {
      return false;
    }
  }
  return true;
}
 
template <typename OccType>
bool OccupantDoFIsEquivalent<OccType>::operator()(
    xtal::SymOp const &_op) const {
  Index j;
  for (Index i = 0; i < m_dof.size(); ++i) {
    OccType t_occ = copy_apply(_op, m_dof[i]);
    for (j = 0; j < m_dof.size(); ++j) {
      if (t_occ.identical(m_dof[j], m_tol)) {
        m_P.set(i) = j;
        break;
      }
    }
    if (j == m_dof.size()) {
      return false;
    }
  }
  return true;
}
 
template <typename OccType>
bool OccupantDoFIsEquivalent<OccType>::operator()(
    xtal::SymOp const &_op, OccDoFType const &_other) const {
  if (_other.size() != m_dof.size()) return false;
  Index j;
  for (Index i = 0; i < m_dof.size(); ++i) {
    OccType t_occ = sym::copy_apply(_op, m_dof[i]);
    for (j = 0; j < _other.size(); ++j) {
      if (t_occ.identical(_other[j], m_tol)) {
        m_P.set(i) = j;
        break;
      }
    }
    if (j == _other.size()) {
      return false;
    }
  }
  return true;
}
 
}  // namespace CASM