rangegenerators.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_MMESH_RANGEGENERATORS_HH
#define DUNE_MMESH_RANGEGENERATORS_HH
 
#include <dune/common/iteratorrange.hh>
 
namespace Dune {
template <typename Entity>
inline auto incidentElements(const Entity& entity)
#ifndef DOXYGEN_SHOULD_SKIP_THIS
    -> IteratorRange<decltype(entity.impl().incidentBegin())>
#endif /* DOXYGEN_SHOULD_SKIP_THIS */
{
  static_assert(Entity::mydimension <= 1,
                "Incident element range iterator is only available for "
                "vertices and edges!");
  return IteratorRange<decltype(entity.impl().incidentBegin())>(
      entity.impl().incidentBegin(), entity.impl().incidentEnd());
}
 
template <typename Vertex>
inline auto incidentFacets(const Vertex& vertex)
#ifndef DOXYGEN_SHOULD_SKIP_THIS
    -> IteratorRange<decltype(vertex.impl().incidentFacetsBegin())>
#endif /* DOXYGEN_SHOULD_SKIP_THIS */
{
  static_assert(
      Vertex::mydimension == 0,
      "Incident facet range iterator is only available for vertices!");
  return IteratorRange<decltype(vertex.impl().incidentFacetsBegin())>(
      vertex.impl().incidentFacetsBegin(), vertex.impl().incidentFacetsEnd());
}
 
template <typename Vertex>
inline auto incidentVertices(const Vertex& vertex, bool includeInfinite = false)
#ifndef DOXYGEN_SHOULD_SKIP_THIS
    -> IteratorRange<
        decltype(vertex.impl().incidentVerticesBegin(includeInfinite))>
#endif /* DOXYGEN_SHOULD_SKIP_THIS */
{
  static_assert(
      Vertex::mydimension == 0,
      "Incident vertices range iterator is only available for vertices!");
  return IteratorRange<decltype(vertex.impl().incidentVerticesBegin(
      includeInfinite))>(vertex.impl().incidentVerticesBegin(includeInfinite),
                         vertex.impl().incidentVerticesEnd(includeInfinite));
}
 
template <typename GridView, int codim = 1>
inline auto interfaceElements(const GridView& gv, bool includeBoundary = false)
#ifndef DOXYGEN_SHOULD_SKIP_THIS
    -> IteratorRange<
        decltype(gv.grid().template interfaceBegin<codim>(includeBoundary))>
#endif /* DOXYGEN_SHOULD_SKIP_THIS */
{
  return IteratorRange<decltype(gv.grid().template interfaceBegin<codim>(
      includeBoundary))>(
      gv.grid().template interfaceBegin<codim>(includeBoundary),
      gv.grid().template interfaceEnd<codim>(includeBoundary));
}
 
template <typename GridView>
inline auto interfaceVertices(const GridView& gv, bool includeBoundary = false)
#ifndef DOXYGEN_SHOULD_SKIP_THIS
    -> IteratorRange<
        decltype(gv.grid().interfaceVerticesBegin(includeBoundary))>
#endif /* DOXYGEN_SHOULD_SKIP_THIS */
{
  return IteratorRange<decltype(gv.grid().interfaceVerticesBegin(
      includeBoundary))>(gv.grid().interfaceVerticesBegin(includeBoundary),
                         gv.grid().interfaceVerticesEnd(includeBoundary));
}
 
template <typename Vertex>
inline auto incidentInterfaceVertices(const Vertex& vertex)
#ifndef DOXYGEN_SHOULD_SKIP_THIS
    -> IteratorRange<decltype(vertex.impl().incidentInterfaceVerticesBegin())>
#endif /* DOXYGEN_SHOULD_SKIP_THIS */
{
  static_assert(Vertex::mydimension == 0,
                "Incident interface vertices range iterator is only available "
                "for interface vertices!");
  return IteratorRange<
      decltype(vertex.impl().incidentInterfaceVerticesBegin())>(
      vertex.impl().incidentInterfaceVerticesBegin(),
      vertex.impl().incidentInterfaceVerticesEnd());
}
 
template <typename Entity>
inline auto incidentInterfaceElements(const Entity& entity)
#ifndef DOXYGEN_SHOULD_SKIP_THIS
    -> IteratorRange<decltype(entity.impl().incidentInterfaceElementsBegin())>
#endif /* DOXYGEN_SHOULD_SKIP_THIS */
{
  static_assert(Entity::mydimension <= 1,
                "Incident interface vertices range iterator is only available "
                "for interface vertices and edges!");
  return IteratorRange<
      decltype(entity.impl().incidentInterfaceElementsBegin())>(
      entity.impl().incidentInterfaceElementsBegin(),
      entity.impl().incidentInterfaceElementsEnd());
}
 
}  // namespace Dune
 
#endif  // DUNE_MMESH_RANGEGENERATORS_HH