elem/CohElem.cpp

#include "DAPTA.hpp"

namespace DAPTA { // Define namespace DAPTA

/*
CohElem2<V>::CohElem2(ElementHandle t1, int f1, ElementHandle t2, int f2)
Create the cohesive element for a face with 3 vertices, given the face definitions
*/
template <typename E>
CohElem<E>::CohElem(ElementHandle t1, int f1, ElementHandle t2, int f2) {
   assert(t1 != NULL);
   assert(t2 != NULL);

   this->elements.at(0) = t1;
   this->elements.at(1) = t2;
   this->face_indices.at(0) = f1;
   this->face_indices.at(1) = f2;
   
   globalID = 0;
}

template <typename E>
std::vector<typename CohElem<E>::VertexHandle> CohElem<E>::vertices() {
   typename CohElem<E>::ElementType::FaceElementType face_1 = elements.at(0)->GetFace(face_indices.at(0));
   typename CohElem<E>::ElementType::FaceElementType face_2 = elements.at(1)->GetFace(face_indices.at(1));
   
   std::vector<VertexHandle> vertices;
   vertices.resize(num_vertices);
   
   for(int i=0; i<(num_vertices/2); i++) {
      vertices.at(i) = face_1.vertices.at(i);
      vertices.at((num_vertices)-i-1) = face_2.vertices.at(i);
   }
   
   return vertices;
}

} // namespace DAPTA

---