QuadrilateralElement.h

// Class QuadriateralElement
//   The QuadrilateralElement class. This is a class representing a quadrilateral element
//   for finite element methods.
 
#ifndef IPPL_QUADRILATERALELEMENT_H
#define IPPL_QUADRILATERALELEMENT_H
 
#include "FEM/Elements/Element.h"
 
namespace ippl {
 
    // this is the basic square 2D element with 4 vertices
    template <typename T>
    class QuadrilateralElement : public Element2D<T, 4> {
    public:
        static constexpr unsigned NumVertices = 4;
 
        typedef typename Element2D<T, NumVertices>::point_t point_t;
        typedef typename Element2D<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, 2>& point) const;
    };
}  // namespace ippl
 
#include "FEM/Elements/QuadrilateralElement.hpp"
 
#endif