element.hpp

// FEAT3: Finite Element Analysis Toolbox, Version 3
// Copyright (C) 2010 by Stefan Turek & the FEAT group
// FEAT3 is released under the GNU General Public License version 3,
// see the file 'copyright.txt' in the top level directory for details.
 
#pragma once
 
// includes, FEAT
#include <kernel/space/element_base.hpp>
#include <kernel/space/dof_mapping_common.hpp>
#include <kernel/space/dof_assignment_common.hpp>
#include <kernel/space/argyris/dof_traits.hpp>
#include <kernel/space/argyris/evaluator.hpp>
#include <kernel/space/argyris/node_functional.hpp>
 
namespace FEAT
{
  namespace Space
  {
    namespace Argyris
    {
      template<typename Trafo_>
      class Element :
        public ElementBase<Trafo_>
      {
      public:
        typedef ElementBase<Trafo_> BaseClass;
        typedef Trafo_ TrafoType;
        typedef typename TrafoType::MeshType MeshType;
        typedef typename TrafoType::ShapeType ShapeType;
 
        static constexpr int num_vert_dofs = 6;
        static constexpr int num_edge_dofs = 1;
 
        static constexpr bool have_node_func = true;
 
        static constexpr int local_degree = 5;
 
        template<
          typename TrafoEvaluator_,
          typename DataType_ = typename TrafoEvaluator_::DataType>
        class Evaluator
        {
        private:
          typedef typename TrafoEvaluator_::EvalPolicy EvalPolicy;
 
          static constexpr int num_loc_dofs = 21;
 
          typedef StandardScalarEvalTraits<EvalPolicy, num_loc_dofs, DataType_> Traits;
 
        public:
          typedef Argyris::Evaluator<Element, TrafoEvaluator_, Traits> Type;
        };
 
        typedef DofMappingUniform<Element, Argyris::DofTraits, ShapeType> DofMappingType;
 
        template<
          int shape_dim_,
          typename DataType_ = Real>
        class DofAssignment
        {
        public:
          typedef DofAssignmentUniform<Element, shape_dim_, DataType_, Argyris::DofTraits, ShapeType> Type;
        };
 
        template<
          int shape_dim_,
          typename DataType_ = Real>
        class NodeFunctional
        {
        public:
          typedef Argyris::NodeFunctional<Element, shape_dim_, DataType_> Type;
        };
 
      public:
        explicit Element(TrafoType& trafo)
          : BaseClass(trafo)
        {
        }
 
        virtual ~Element()
        {
        }
 
        Index get_num_dofs() const
        {
          return
            // number of vertex dofs
            Index(num_vert_dofs) * this->get_mesh().get_num_entities(0) +
            // number of edge dofs
            Index(num_edge_dofs) * this->get_mesh().get_num_entities(1);
        }
 
        static String name()
        {
          return "Argyris";
        }
      }; // class Element<...>
    } // namespace Argyris
  } // namespace Space
} // namespace FEAT