HexahedralElement.h

//
// Class HexahedralElement
//   The HexahedralElement class. This is a class representing a hexahedron element
//   for finite element methods.
#ifndef IPPL_HEXAHEDRALELEMENT_H
#define IPPL_HEXAHEDRALELEMENT_H
 
#include "FEM/Elements/Element.h"
 
namespace ippl {
 
    // this is the basic cube 3D element with 8 vertices
    template <typename T>
    class HexahedralElement : public Element3D<T, 8> {
    public:
        static constexpr unsigned NumVertices = 8;
 
        typedef typename Element3D<T, NumVertices>::point_t point_t;
        typedef typename Element3D<T, NumVertices>::vertex_points_t vertex_points_t;
 
        KOKKOS_FUNCTION vertex_points_t getLocalVertices() const;
 
        KOKKOS_FUNCTION point_t
        getTransformationJacobian(const vertex_points_t& global_vertices) const;
 
        KOKKOS_FUNCTION point_t
        getInverseTransformationJacobian(const vertex_points_t& global_vertices) const;
 
        KOKKOS_FUNCTION point_t globalToLocal(const vertex_points_t& global_vertices,
                                              const point_t& point) const;
 
        KOKKOS_FUNCTION point_t localToGlobal(const vertex_points_t& global_vertices,
                                              const point_t& point) const;
 
        KOKKOS_FUNCTION T
        getDeterminantOfTransformationJacobian(const vertex_points_t& global_vertices) const;
 
        KOKKOS_FUNCTION point_t
        getInverseTransposeTransformationJacobian(const vertex_points_t& global_vertices) const;
 
        KOKKOS_FUNCTION bool isPointInRefElement(const Vector<T, 3>& point) const;
    };
 
}  // namespace ippl
 
#include "FEM/Elements/HexahedralElement.hpp"
 
#endif