mapping.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/trafo/mapping_base.hpp>
#include <kernel/trafo/standard/evaluator.hpp>
#include <kernel/geometry/mesh_part.hpp>
#include <kernel/geometry/atlas/chart.hpp>
 
namespace FEAT
{
  namespace Trafo
  {
    namespace Standard
    {
      template<typename Mesh_>
      class Mapping :
        public MappingBase<Mesh_>
      {
      public:
        typedef MappingBase<Mesh_> BaseClass;
        typedef Mesh_ MeshType;
        typedef typename MeshType::VertexSetType::CoordType CoordType;
        typedef typename MeshType::ShapeType ShapeType;
        typedef typename Shape::FaceTraits<ShapeType, ShapeType::dimension - 1>::ShapeType FacetType;
 
      public:
        template<
          typename Shape_ = ShapeType,
          typename CoordType_ = Real>
        class Evaluator
        {
        private:
          typedef Trafo::StandardEvalPolicy<Shape_, CoordType_, MeshType::world_dim> EvalPolicy;
 
        public:
          typedef Trafo::Standard::Evaluator<Mapping, EvalPolicy> Type;
        };
 
      public:
        explicit Mapping(MeshType& mesh) :
          BaseClass(mesh)
        {
        }
 
        template<typename ShapeType_ = ShapeType, typename CoordType_ = CoordType>
        CoordType_ compute_vol(const Index cell) const
        {
          typename Evaluator<ShapeType_, CoordType_>::Type evaluator(*this);
 
          evaluator.prepare(cell);
          CoordType_ vol = evaluator.volume();
          evaluator.finish();
          return vol;
        }
 
        CoordType compute_vol()
        {
          typename Evaluator<ShapeType, CoordType>::Type evaluator(*this);
 
          CoordType vol(0);
 
          Index num_cells = this->get_mesh().get_num_entities(ShapeType::dimension);
          for(Index cell(0); cell < num_cells; ++cell)
          {
            evaluator.prepare(cell);
            vol += evaluator.volume();
            evaluator.finish();
          }
 
          return vol;
        }
 
        void add_meshpart_chart(const Geometry::MeshPart<MeshType>& DOXY(mesh_part), const Geometry::Atlas::ChartBase<MeshType>& DOXY(chart))
        {
          // nothing to do here; the only purpose of this function is to supply
          // the same interface as the Trafo::Isoparam::Mapping class
        }
      }; // class Mapping<...>
    } // namespace Standard
  } // namespace Trafo
} // namespace FEAT