evaluator.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/parametric_evaluator.hpp>
#include <kernel/space/dof_mapping_common.hpp>
#include <kernel/space/lagrange2/details.hpp>
 
namespace FEAT
{
  namespace Space
  {
    namespace Lagrange2
    {
      static constexpr SpaceTags ref_caps = SpaceTags::ref_value | SpaceTags::ref_grad | SpaceTags::ref_hess;
 
      template<
        typename Space_,
        typename TrafoEvaluator_,
        typename SpaceEvalTraits_,
        typename Shape_ = typename Space_::ShapeType>
      class Evaluator :
        public ParametricEvaluator<
          Evaluator<
            Space_,
            TrafoEvaluator_,
            SpaceEvalTraits_,
            Shape_ >,
          TrafoEvaluator_,
          SpaceEvalTraits_,
          ref_caps>
      {
      public:
        typedef ParametricEvaluator<Evaluator, TrafoEvaluator_, SpaceEvalTraits_, ref_caps> BaseClass;
 
        typedef Space_ SpaceType;
 
        typedef Shape_ ShapeType;
 
        typedef SpaceEvalTraits_ SpaceEvalTraits;
 
        typedef typename SpaceEvalTraits::EvalPolicy EvalPolicy;
 
        typedef typename EvalPolicy::DomainPointType DomainPointType;
 
        typedef typename SpaceEvalTraits::DataType DataType;
 
        typedef EvalHelper<DomainPointType, DataType, ShapeType> EvalHp;
 
      public:
        explicit Evaluator(const SpaceType& DOXY(space))
        {
        }
 
        int get_num_local_dofs() const
        {
          return EvalHp::get_num_local_dofs();
        }
 
        template<typename EvalData_>
        void eval_ref_values(
          EvalData_& data,
          const DomainPointType& point) const
        {
          EvalHp::template eval_ref_values<EvalData_>(data, point);
        }
 
        template<typename EvalData_>
        void eval_ref_gradients(
          EvalData_& data,
          const DomainPointType& point) const
        {
          EvalHp::template eval_ref_gradients<EvalData_>(data, point);
        }
 
 
        template<typename EvalData_>
        void eval_ref_hessians(
          EvalData_& data,
          const DomainPointType& point) const
        {
          EvalHp::template eval_ref_hessians<EvalData_>(data, point);
        }
      };
    } // namespace Lagrange2
  } // namespace Space
} // namespace FEAT