edge_element.h

Go to the documentation of this file.

#ifndef EDGE_ELEMENT_H
#define EDGE_ELEMENT_H
 
#include "bitrl/utils/geometry/mesh/element.h"
#include "bitrl/utils/geometry/mesh/face_element.h"
 
namespace bitrl
{
namespace utils
{
namespace geom
{
template <int dim> class EdgeElem;
 
template <int dim> class EdgeElem : public FaceElement<dim, 1>
{
  public:
    typedef FaceElement<2, 1>::node_ptr_t node_ptr_t;
 
    EdgeElem(uint_t id, uint_t proc_id = 0);
 
    EdgeElem(uint_t id, uint_t n_nodes, uint_t proc_id);
};
 
template <> class EdgeElem<1> : public Element<1>
{
  public:
    typedef Element<1>::neighbor_ptr_t neighbor_ptr_t;
    typedef Element<1>::neighbor_ref_t neighbor_ref_t;
    typedef Element<1>::node_ptr_t node_ptr_t;
    typedef Element<1>::edge_ptr_t edge_ptr_t;
    typedef Element<1>::cedge_ptr_t cedge_ptr_t;
    typedef Element<1>::edge_ref_t edge_ref_t;
    typedef Element<1>::cedge_ref_t cedge_ref_t;
    typedef Element<1>::face_ptr_t face_ptr_t;
    typedef Element<1>::cface_ptr_t cface_ptr_t;
    typedef Element<1>::face_ref_t face_ref_t;
    typedef Element<1>::cface_ref_t cface_ref_t;
 
    EdgeElem(uint_t id, uint_t proc_id = 0);
 
    virtual uint_t n_nodes() const override { return 2; }
 
    virtual void set_node(uint_t i, node_ptr_t node) override final;
 
    virtual node_ptr_t get_node(uint_t n) override final;
 
    virtual uint_t n_vertices() const override final { return 2; }
 
    virtual uint_t n_edges() const override final { return 2; }
 
    virtual uint_t n_faces() const override final { return 2; }
 
    virtual void resize_faces() override final {}
 
    virtual void set_face(uint_t, face_ptr_t) override final {}
 
    virtual cface_ref_t get_face(uint_t f) const override;
 
    virtual face_ref_t get_face(uint_t f) override;
 
    virtual uint_t n_neighbors() const { return 2; }
 
    virtual void face_vertices(uint_t f, std::vector<uint_t> &ids) const override final;
 
    virtual real_t volume() const override final;
 
    virtual uint_t which_face_am_i(cface_ref_t face) const override final
    {
        throw std::logic_error("Not implemented");
    }
 
    virtual const DynVec<real_t> face_normal_vector(uint_t f) const override final
    {
        throw std::logic_error("Not implemented");
    }
};
 
} // namespace geom
} // namespace utils
} // namespace bitrl
#endif // EDGE_ELEMENT_H