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>
 
namespace FEAT
{
  namespace Space
  {
    namespace CaiDouSanSheYe
    {
      static constexpr SpaceTags ref_caps = SpaceTags::ref_value | SpaceTags::ref_grad;
 
      template<
        typename Space_,
        typename TrafoEvaluator_,
        typename SpaceEvalTraits_,
        typename Shape_ = typename Space_::ShapeType>
      class Evaluator DOXY({});
 
      template<
        typename Space_,
        typename TrafoEvaluator_,
        typename SpaceEvalTraits_>
      class Evaluator<Space_, TrafoEvaluator_, SpaceEvalTraits_, Shape::Quadrilateral> :
        public ParametricEvaluator<
          Evaluator< Space_, TrafoEvaluator_, SpaceEvalTraits_, Shape::Quadrilateral>,
          TrafoEvaluator_,
          SpaceEvalTraits_,
          ref_caps>
      {
      public:
        typedef ParametricEvaluator<Evaluator, TrafoEvaluator_, SpaceEvalTraits_, ref_caps> BaseClass;
 
        typedef Space_ SpaceType;
 
        typedef SpaceEvalTraits_ SpaceEvalTraits;
 
        typedef typename SpaceEvalTraits::EvalPolicy EvalPolicy;
 
        typedef typename EvalPolicy::DomainPointType DomainPointType;
 
        typedef typename SpaceEvalTraits::DataType DataType;
 
      public:
        explicit Evaluator(const SpaceType& DOXY(space))
        {
        }
 
        virtual ~Evaluator()
        {
        }
 
        int get_num_local_dofs() const
        {
          return 5;
        }
 
        template<typename EvalData_>
        void eval_ref_values(
          EvalData_& data,
          const DomainPointType& point) const
        {
          const DataType x = point[0];
          const DataType y = point[1];
          const DataType x2 = x*x;
          const DataType y2 = y*y;
          // edge midpoints dofs
          data.phi[0].ref_value = DataType(0.25) - DataType(0.5)*y + ( DataType(0.375) - DataType(0.625)*x2)*x2 + (-DataType(0.375) + DataType(0.625)*y2)*y2;
          data.phi[1].ref_value = DataType(0.25) + DataType(0.5)*y + ( DataType(0.375) - DataType(0.625)*x2)*x2 + (-DataType(0.375) + DataType(0.625)*y2)*y2;
          data.phi[2].ref_value = DataType(0.25) - DataType(0.5)*x + (-DataType(0.375) + DataType(0.625)*x2)*x2 + ( DataType(0.375) - DataType(0.625)*y2)*y2;
          data.phi[3].ref_value = DataType(0.25) + DataType(0.5)*x + (-DataType(0.375) + DataType(0.625)*x2)*x2 + ( DataType(0.375) - DataType(0.625)*y2)*y2;
          // bubble dof
          data.phi[4].ref_value = DataType(2.25)*x*y;
        }
 
        template<typename EvalData_>
        void eval_ref_gradients(
          EvalData_& data,
          const DomainPointType& point) const
        {
          const DataType x = point[0];
          const DataType y = point[1];
          const DataType x2 = x*x;
          const DataType y2 = y*y;
          // edge midpoint dofs
          data.phi[0].ref_grad[0] = ( DataType(0.75) - DataType(2.5)*x2)*x;
          data.phi[0].ref_grad[1] = (-DataType(0.75) + DataType(2.5)*y2)*y - DataType(0.5);
          data.phi[1].ref_grad[0] = ( DataType(0.75) - DataType(2.5)*x2)*x;
          data.phi[1].ref_grad[1] = (-DataType(0.75) + DataType(2.5)*y2)*y + DataType(0.5);
          data.phi[2].ref_grad[0] = (-DataType(0.75) + DataType(2.5)*x2)*x - DataType(0.5);
          data.phi[2].ref_grad[1] = ( DataType(0.75) - DataType(2.5)*y2)*y;
          data.phi[3].ref_grad[0] = (-DataType(0.75) + DataType(2.5)*x2)*x + DataType(0.5);
          data.phi[3].ref_grad[1] = ( DataType(0.75) - DataType(2.5)*y2)*y;
          // bubble dof
          data.phi[4].ref_grad[0] = DataType(2.25)*y;
          data.phi[4].ref_grad[1] = DataType(2.25)*x;
        }
      }; // Evaluator<..., Shape::Quadrilateral>
    } // namespace CaiDouSanSheYe
  } // namespace Space
} // namespace FEAT