twistutility.hh

#ifndef DUNE_MMESH_MISC_TWISTUTILITY_HH
#define DUNE_MMESH_MISC_TWISTUTILITY_HH
 
#include <dune/geometry/type.hh>
 
namespace Dune
{
 
  template< class Grid >
  class MMeshInterfaceGrid;
 
  namespace MMeshTwist
  {
    template< class LeafIntersection >
    static inline int twistInSelf(const LeafIntersection& intersection)
    {
      return 0;
    }
 
    template< class LeafIntersection >
    static inline int twistInNeighbor(const LeafIntersection& intersection )
    {
      return 0;
    }
  }
 
  namespace MMeshInterfaceTwist
  {
    template< class LeafIntersection >
    static inline int twistInSelf(const LeafIntersection& intersection)
    {
      return 0;
    }
 
    template< class LeafIntersection >
    static inline int twistInNeighbor(const LeafIntersection& intersection )
    {
      return 0;
    }
  }
 
  namespace Fem
  {
 
    template< class Grid >
    struct TwistUtility;
 
    template< class HostGrid, int dim >
    struct TwistUtility< MMesh< HostGrid, dim > >
    {
      typedef MMesh< HostGrid, dim > GridType;
      typedef typename GridType::Traits::LeafIntersectionIterator LeafIntersectionIterator;
      typedef typename LeafIntersectionIterator::Intersection LeafIntersection;
      typedef typename GridType::Traits::LevelIntersectionIterator LevelIntersectionIterator;
      typedef typename LevelIntersectionIterator::Intersection LevelIntersection;
 
    public:
      static inline int twistInSelf(const GridType& grid, const LeafIntersection& intersection)
      {
        return MMeshTwist::twistInSelf(intersection);
      }
 
      static inline int twistInNeighbor(const GridType& grid, const LeafIntersection& intersection )
      {
        return MMeshTwist::twistInNeighbor(intersection);
      }
 
      template <class Intersection>
      static inline GeometryType
      elementGeometry(const Intersection& intersection, const bool inside)
      {
        return Dune::GeometryTypes::simplex(dim);
      }
 
    private:
      TwistUtility(const TwistUtility&);
      TwistUtility& operator=(const TwistUtility&);
    };
 
 
    template< class MMesh >
    struct TwistUtility< MMeshInterfaceGrid<MMesh> >
    {
      typedef MMeshInterfaceGrid<MMesh> GridType;
      static constexpr int dim = GridType::dimension;
      typedef typename GridType::Traits::LeafIntersectionIterator LeafIntersectionIterator;
      typedef typename LeafIntersectionIterator::Intersection LeafIntersection;
      typedef typename GridType::Traits::LevelIntersectionIterator LevelIntersectionIterator;
      typedef typename LevelIntersectionIterator::Intersection LevelIntersection;
 
    public:
      static inline int twistInSelf(const GridType& grid, const LeafIntersection& intersection)
      {
        return MMeshInterfaceTwist::twistInSelf(intersection);
      }
 
      static inline int twistInNeighbor(const GridType& grid, const LeafIntersection& intersection )
      {
        return MMeshInterfaceTwist::twistInNeighbor(intersection);
      }
 
      template <class Intersection>
      static inline GeometryType
      elementGeometry(const Intersection& intersection, const bool inside)
      {
        return Dune::GeometryTypes::simplex(dim);
      }
 
    private:
      TwistUtility(const TwistUtility&);
      TwistUtility& operator=(const TwistUtility&);
    };
 
  }  // namespace Fem
 
} // namespace Dune
 
#endif // #ifndef DUNE_MMESH_MISC_TWISTUTILITY_HH