gridview.hh

// SPDX-FileCopyrightText: Copyright © DUNE Project contributors, see file LICENSE.md in module root
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception
// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_GRID_CONCEPTS_GRIDVIEW_HH
#define DUNE_GRID_CONCEPTS_GRIDVIEW_HH
 
#include <concepts>
#include <cstddef>
#include <utility>
 
#include <dune/grid/common/capabilities.hh>
#include <dune/grid/common/gridenums.hh>
#include <dune/grid/concepts/entityiterator.hh>
#include <dune/grid/concepts/geometry.hh>
#include <dune/grid/concepts/indexidset.hh>
#include <dune/grid/concepts/intersection.hh>
#include <dune/grid/concepts/intersectioniterator.hh>
#include <dune/grid/concepts/archetypes/datahandle.hh>
 
namespace Dune::Concept {
namespace Impl {
 
  template<class GV, int codim, Dune::PartitionIteratorType partition>
  concept GridViewPartition =
    EntityIterator<typename GV::template Codim<codim>::template Partition<partition>::Iterator> &&
  requires(const GV gv)
  {
    { gv.template begin<codim,partition>() } -> std::same_as<typename GV::template Codim<codim>::template Partition<partition>::Iterator>;
    { gv.template end<codim,partition>()   } -> std::same_as<typename GV::template Codim<codim>::template Partition<partition>::Iterator>;
  };
 
  template<class GV, int codim>
  concept GridViewAllPartitions =
    GridViewPartition<GV,codim,Dune::PartitionIteratorType::InteriorBorder_Partition> &&
    GridViewPartition<GV,codim,Dune::PartitionIteratorType::Overlap_Partition> &&
    GridViewPartition<GV,codim,Dune::PartitionIteratorType::OverlapFront_Partition> &&
    GridViewPartition<GV,codim,Dune::PartitionIteratorType::All_Partition> &&
    GridViewPartition<GV,codim,Dune::PartitionIteratorType::Ghost_Partition>;
 
  template<class GV, int codim>
  concept GridViewCodim =
    Geometry<typename GV::template Codim<codim>::Geometry> &&
    Geometry<typename GV::template Codim<codim>::LocalGeometry> &&
    EntityIterator<typename GV::template Codim<codim>::Iterator> &&
  requires(const GV gv)
  {
    { gv.template begin<codim>() } -> std::same_as<typename GV::template Codim<codim>::Iterator>;
    { gv.template end<codim>()   } -> std::same_as<typename GV::template Codim<codim>::Iterator>;
 
    requires (codim != 0) || requires(const typename GV::template Codim<codim>::Entity& entity)
    {
      { gv.ibegin(entity) } -> std::same_as<typename GV::IntersectionIterator>;
      { gv.iend(entity)   } -> std::same_as<typename GV::IntersectionIterator>;
    };
  } && GridViewAllPartitions<GV,codim>;
 
  template<class GV, class Grid, int codim>
    requires Dune::Capabilities::hasEntityIterator<Grid,codim>::v
  void requireGridViewCodim()
    requires GridViewCodim<GV,codim> {}
 
  template<class GV, class Grid, int codim>
    requires (not Dune::Capabilities::hasEntityIterator<Grid,codim>::v)
  void requireGridViewCodim() {}
 
  template <class GV, std::size_t... c>
  void gridViewAllCodims(std::index_sequence<c...>)
    requires requires {
      (requireGridViewCodim<GV,typename GV::Grid,int(c)>(),...);
    };
 
} // end namespace Impl
 
template<class GV>
concept GridView = std::copyable<GV> &&
  IndexSet<typename GV::IndexSet> &&
  Intersection<typename GV::Intersection> &&
  IntersectionIterator<typename GV::IntersectionIterator> &&
requires(const GV gv, int codim, Dune::GeometryType type)
{
  typename GV::Traits;
  typename GV::ctype;
  { GV::conforming        } -> std::convertible_to<bool>;
  { GV::dimension         } -> std::convertible_to<int>;
  { GV::dimensionworld    } -> std::convertible_to<int>;
  { gv.grid()             } -> std::same_as<const typename GV::Grid&>;
  { gv.indexSet()         } -> std::same_as<const typename GV::IndexSet&>;
  { gv.size(codim)        } -> std::convertible_to<int>;
  { gv.size(type)         } -> std::convertible_to<int>;
  { gv.comm()             } -> std::convertible_to<typename GV::Communication>;
  { gv.overlapSize(codim) } -> std::convertible_to<int>;
  { gv.ghostSize(codim)   } -> std::convertible_to<int>;
 
  requires requires(Archetypes::CommDataHandle<std::byte>& handle,
                    InterfaceType iface, CommunicationDirection dir)
  {
    gv.communicate(handle, iface, dir);
  };
} &&
Impl::GridViewCodim<GV,0> &&
requires (index_constant<1> from, index_constant<GV::dimension+1> to) {
  Impl::gridViewAllCodims<GV>(range(from, to).to_integer_sequence());
};
 
}  // end namespace Dune::Concept
 
#endif // DUNE_GRID_CONCEPTS_GRIDVIEW_HH