standard_refinement_traits.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/shape.hpp>
 
namespace FEAT
{
  namespace Geometry
  {
    namespace Intern
    {
      template<
        typename Shape_,
        int face_dim_>
  #ifndef DOXYGEN
      struct StandardRefinementTraits;
  #else
      struct StandardRefinementTraits
      {
        typedef Shape_ ShapeType;
 
        static constexpr int cell_dim = ShapeType::dimension;
 
        static constexpr int face_dim = face_dim_;
 
        static constexpr int count = ...;
      };
  #endif // DOXYGEN
 
      template<int face_dim_>
      struct StandardRefinementTraits<Shape::Vertex, face_dim_>
      {
        // validate dimensions
        static_assert(face_dim_ == 0, "invalid face dimension");
 
        typedef Shape::Vertex ShapeType;
 
        static constexpr int cell_dim = ShapeType::dimension; // = 0
 
        static constexpr int face_dim = face_dim_; // = 0
 
        static constexpr int count = 1;
 
        static inline String name()
        {
          return "StandardRefinementTraits<Vertex," + stringify(face_dim_) + ">";
        }
      }; // struct StandardRefinementTraits<...,Vertex,...>
 
      template<
        int cell_dim_,
        int face_dim_>
      struct StandardRefinementTraits<Shape::Hypercube<cell_dim_>, face_dim_>
      {
        // validate dimensions
        static_assert(face_dim_ >= 0, "invalid face dimension");
        static_assert(cell_dim_ >= face_dim_, "invalid cell dimension");
 
        typedef Shape::Hypercube<cell_dim_> ShapeType;
 
        typedef typename Shape::FaceTraits<ShapeType, face_dim_>::ShapeType FaceType;
 
        static constexpr int cell_dim = cell_dim_;
 
        static constexpr int face_dim = face_dim_;
 
        static constexpr int count = Shape::FaceTraits<ShapeType, cell_dim - face_dim>::count;
 
        static inline String name()
        {
          return "StandardRefinementTraits<" + ShapeType::name() + "," + stringify(face_dim_) + ">";
        }
      }; // struct StandardRefinementTraits<Hypercube<...>,...>
 
      template<
        int face_dim_>
      struct StandardRefinementTraits<Shape::Simplex<1>, face_dim_>
      {
        // validate dimensions
        static_assert(face_dim_ >= 0, "invalid face dimension");
        static_assert(1 >= face_dim_, "invalid cell dimension");
 
        typedef Shape::Simplex<1> ShapeType;
 
        typedef typename Shape::FaceTraits<ShapeType, face_dim_>::ShapeType FaceType;
 
        static constexpr int cell_dim = 1;
 
        static constexpr int face_dim = face_dim_;
 
        static constexpr int count = Shape::FaceTraits<ShapeType, cell_dim - face_dim>::count;
 
        static inline String name()
        {
          return "StandardRefinementTraits<" + ShapeType::name() + "," + stringify(face_dim_) + ">";
        }
      }; // struct StandardRefinementTraits<Simplex<1>,...>
 
      template<>
      struct StandardRefinementTraits<Shape::Simplex<2>, 0>
      {
        typedef Shape::Simplex<2> ShapeType;
 
        static constexpr int cell_dim = 2;
 
        static constexpr int face_dim = 0;
 
        static constexpr int count = 0;
 
        static inline String name()
        {
          return "StandardRefinementTraits<" + ShapeType::name() + ",0>";
        }
      }; // struct StandardRefinementTraits<Simplex<2>,0>
 
      template<>
      struct StandardRefinementTraits<Shape::Simplex<2>, 1>
      {
        typedef Shape::Simplex<2> ShapeType;
 
        static constexpr int cell_dim = 2;
 
        static constexpr int face_dim = 1;
 
        static constexpr int count = 3;
 
        static inline String name()
        {
          return "StandardRefinementTraits<" + ShapeType::name() + ",1>";
        }
      }; // struct StandardRefinementTraits<Simplex<2>,1>
 
      template<>
      struct StandardRefinementTraits<Shape::Simplex<2>, 2>
      {
        typedef Shape::Simplex<2> ShapeType;
 
        static constexpr int cell_dim = 2;
 
        static constexpr int face_dim = 2;
 
        static constexpr int count = 4;
 
        static inline String name()
        {
          return "StandardRefinementTraits<" + ShapeType::name() + ",2>";
        }
      }; // struct StandardRefinementTraits<Simplex<2>,2>
 
      template<>
      struct StandardRefinementTraits<Shape::Simplex<3>, 0>
      {
        typedef Shape::Simplex<3> ShapeType;
 
        static constexpr int cell_dim = 3;
 
        static constexpr int face_dim = 0;
 
        static constexpr int count = 1;
 
        static inline String name()
        {
          return "StandardRefinementTraits<" + ShapeType::name() + ",0>";
        }
      }; // struct StandardRefinementTraits<Simplex<3>,0>
 
      template<>
      struct StandardRefinementTraits<Shape::Simplex<3>, 1>
      {
        typedef Shape::Simplex<3> ShapeType;
 
        static constexpr int cell_dim = 3;
 
        static constexpr int face_dim = 1;
 
        static constexpr int count = 6;
 
        static inline String name()
        {
          return "StandardRefinementTraits<" + ShapeType::name() + ",1>";
        }
      }; // struct StandardRefinementTraits<Simplex<3>,1>
 
      template<>
      struct StandardRefinementTraits<Shape::Simplex<3>, 2>
      {
        typedef Shape::Simplex<3> ShapeType;
 
        static constexpr int cell_dim = 3;
 
        static constexpr int face_dim = 2;
 
        static constexpr int count = 16;
 
        static inline String name()
        {
          return "StandardRefinementTraits<" + ShapeType::name() + ",2>";
        }
      }; // struct StandardRefinementTraits<Simplex<3>,2>
 
      template<>
      struct StandardRefinementTraits<Shape::Simplex<3>, 3>
      {
        typedef Shape::Simplex<3> ShapeType;
 
        static constexpr int cell_dim = 3;
 
        static constexpr int face_dim = 3;
 
        static constexpr int count = 12;
 
        static inline String name()
        {
          return "StandardRefinementTraits<" + ShapeType::name() + ",3>";
        }
      }; // struct StandardRefinementTraits<Simplex<3>,3>
    } // namespace Intern
  } // namespace Geometry
} // namespace FEAT