intersectiontype.hh

#ifndef DUNE_PDELAB_COMMON_INTERSECTIONTYPE_HH
#define DUNE_PDELAB_COMMON_INTERSECTIONTYPE_HH
 
#include <utility>
#include <tuple>
 
#include <dune/common/version.hh>
#include <dune/grid/common/partitionset.hh>
 
namespace Dune {
  namespace PDELab {
 
    enum class IntersectionType
    {
 
      processor = 0, 
      skeleton = 1,  
      boundary = 2,  
      periodic = 3   
 
    };
 
 
    template<typename EntitySet, typename Intersection>
    std::tuple<IntersectionType,typename EntitySet::Element> classifyIntersection(const EntitySet& entity_set, const Intersection& is)
    {
      auto type = static_cast<IntersectionType>(1* is.neighbor() + 2*is.boundary());
      if (type == IntersectionType::skeleton || type == IntersectionType::periodic)
        if (entity_set.partitions() == Partitions::all)
          return std::make_tuple(type,is.outside());
        else
          {
            auto outside_entity = is.outside();
            if (entity_set.partitions().contains(outside_entity.partitionType()))
              return std::make_tuple(type,outside_entity);
            else
              return std::make_tuple(IntersectionType::processor,std::move(outside_entity));
          }
      else
        return std::make_tuple(type,decltype(is.outside()){});
    }
 
 
  } // namespace PDELab
} // namespace Dune
 
#endif // DUNE_PDELAB_COMMON_INTERSECTIONTYPE_HH