identitygrid.hh

Go to the documentation of this file.

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_GRID_IDENTITYGRID_HH
#define DUNE_GRID_IDENTITYGRID_HH
 
#include <string>
#include <map>
 
#include <dune/common/deprecated.hh>
#include <dune/common/parallel/collectivecommunication.hh>
#include <dune/grid/common/capabilities.hh>
#include <dune/grid/common/grid.hh>
 
// The components of the IdentityGrid interface
#include "identitygrid/identitygridgeometry.hh"
#include "identitygrid/identitygridentity.hh"
#include "identitygrid/identitygridentitypointer.hh"
#include "identitygrid/identitygridentityseed.hh"
#include "identitygrid/identitygridintersectioniterator.hh"
#include "identitygrid/identitygridleveliterator.hh"
#include "identitygrid/identitygridleafiterator.hh"
#include "identitygrid/identitygridhierarchiciterator.hh"
#include "identitygrid/identitygridindexsets.hh"
 
namespace Dune
{
  // Forward declaration
  template <class HostGrid>
  class IdentityGrid;
 
  // External forward declarations
  template< class Grid >
  struct HostGridAccess;
 
  namespace {
 
    // This is just a placeholder to stuff into the standard
    // GridTraits, as those expect an EntityPointer implementation
    // with two template parameters, but our EntityPointer uses three.
    template<int,typename>
    struct DummyIdentityGridEntityPointer
    {};
 
  }
 
 
  template<int dim, class HostGrid>
  struct IdentityGridFamily
  {
 
  private:
 
    // The type signature of our EntityPointer implementation does not fit into
    // the scheme expected by GridTraits, so we have to provide our own Traits
    // On the other hand, we don't want to retype everything, so we just inherit
    // from the default Traits and just tweak the EntityPointer type in the derived
    // Traits class
 
    typedef GridTraits<
        dim,
        HostGrid::dimensionworld,
        Dune::IdentityGrid<HostGrid>,
        IdentityGridGeometry,
        IdentityGridEntity,
        DummyIdentityGridEntityPointer, // placeholder
        IdentityGridLevelIterator,
        IdentityGridLeafIntersection,
        IdentityGridLevelIntersection,
        IdentityGridLeafIntersectionIterator,
        IdentityGridLevelIntersectionIterator,
        IdentityGridHierarchicIterator,
        IdentityGridLeafIterator,
        IdentityGridLevelIndexSet< const IdentityGrid<HostGrid> >,
        IdentityGridLeafIndexSet< const IdentityGrid<HostGrid> >,
        IdentityGridGlobalIdSet< const IdentityGrid<HostGrid> >,
        typename HostGrid::Traits::GlobalIdSet::IdType,
        IdentityGridLocalIdSet< const IdentityGrid<HostGrid> >,
        typename HostGrid::Traits::LocalIdSet::IdType,
        CollectiveCommunication<IdentityGrid<HostGrid> >,
        DefaultLevelGridViewTraits,
        DefaultLeafGridViewTraits,
        IdentityGridEntitySeed
        > BaseTraits;
 
  public:
 
    // inherit everyting from the default Traits...
    struct Traits
      : public BaseTraits
    {
 
      // Except for the Codim struct, which we reimplement, but which
      // again inherits everything from the base class...
      template<int codim>
      struct Codim
        : public BaseTraits::template Codim<codim>
      {
 
        // ... except for the EntityPointer, which we override with our
        // actual implementation
        typedef Dune::EntityPointer<
          const typename BaseTraits::Grid,
          IdentityGridEntityPointer<
            codim,
            const typename BaseTraits::Grid,
            typename HostGrid::template Codim<codim>::EntityPointer
            >
          > EntityPointer;
 
      };
 
    };
 
  };
 
  //**********************************************************************
  //
  // --IdentityGrid
  //
  //************************************************************************
  template <class HostGrid>
  class IdentityGrid
  : public GridDefaultImplementation<HostGrid::dimension, HostGrid::dimensionworld,
                                     typename HostGrid::ctype, IdentityGridFamily<HostGrid::dimension, HostGrid> >
  {
    friend class IdentityGridLevelIndexSet<const IdentityGrid<HostGrid> >;
    friend class IdentityGridLeafIndexSet<const IdentityGrid<HostGrid> >;
    friend class IdentityGridGlobalIdSet<const IdentityGrid<HostGrid> >;
    friend class IdentityGridLocalIdSet<const IdentityGrid<HostGrid> >;
    friend class IdentityGridHierarchicIterator<const IdentityGrid<HostGrid> >;
    friend class IdentityGridLevelIntersectionIterator<const IdentityGrid<HostGrid> >;
    friend class IdentityGridLeafIntersectionIterator<const IdentityGrid<HostGrid> >;
 
    template<int codim, PartitionIteratorType pitype, class GridImp_>
    friend class IdentityGridLevelIterator;
 
    template<int codim, PartitionIteratorType pitype, class GridImp_>
    friend class IdentityGridLeafIterator;
 
 
    template<int codim_, int dim_, class GridImp_>
    friend class IdentityGridEntity;
 
    friend struct HostGridAccess< IdentityGrid< HostGrid > >;
 
  public:
 
    typedef HostGrid HostGridType;
 
    //**********************************************************
    // The Interface Methods
    //**********************************************************
 
    typedef IdentityGridFamily<HostGrid::dimension,HostGrid>  GridFamily;
 
    typedef typename IdentityGridFamily<HostGrid::dimension,HostGrid>::Traits Traits;
 
    typedef typename HostGrid::ctype ctype;
 
 
    explicit IdentityGrid(HostGrid& hostgrid) :
      hostgrid_(&hostgrid),
      leafIndexSet_(*this),
      globalIdSet_(*this),
      localIdSet_(*this)
    {
      setIndices();
    }
 
    ~IdentityGrid()
    {
      // Delete level index sets
      for (size_t i=0; i<levelIndexSets_.size(); i++)
        if (levelIndexSets_[i])
          delete (levelIndexSets_[i]);
    }
 
 
    int maxLevel() const {
      return hostgrid_->maxLevel();
    }
 
    template<int codim>
    typename Traits::template Codim<codim>::LevelIterator lbegin (int level) const {
      return IdentityGridLevelIterator<codim,All_Partition, const IdentityGrid<HostGrid> >(this, level);
    }
 
 
    template<int codim>
    typename Traits::template Codim<codim>::LevelIterator lend (int level) const {
      return IdentityGridLevelIterator<codim,All_Partition, const IdentityGrid<HostGrid> >(this, level, true);
    }
 
 
    template<int codim, PartitionIteratorType PiType>
    typename Traits::template Codim<codim>::template Partition<PiType>::LevelIterator lbegin (int level) const {
      return IdentityGridLevelIterator<codim,PiType, const IdentityGrid<HostGrid> >(this, level);
    }
 
 
    template<int codim, PartitionIteratorType PiType>
    typename Traits::template Codim<codim>::template Partition<PiType>::LevelIterator lend (int level) const {
      return IdentityGridLevelIterator<codim,PiType, const IdentityGrid<HostGrid> >(this, level, true);
    }
 
 
    template<int codim>
    typename Traits::template Codim<codim>::LeafIterator leafbegin() const {
      return IdentityGridLeafIterator<codim,All_Partition, const IdentityGrid<HostGrid> >(this);
    }
 
 
    template<int codim>
    typename Traits::template Codim<codim>::LeafIterator leafend() const {
      return IdentityGridLeafIterator<codim,All_Partition, const IdentityGrid<HostGrid> >(this, true);
    }
 
 
    template<int codim, PartitionIteratorType PiType>
    typename Traits::template Codim<codim>::template Partition<PiType>::LeafIterator leafbegin() const {
      return IdentityGridLeafIterator<codim,PiType, const IdentityGrid<HostGrid> >(this);
    }
 
 
    template<int codim, PartitionIteratorType PiType>
    typename Traits::template Codim<codim>::template Partition<PiType>::LeafIterator leafend() const {
      return IdentityGridLeafIterator<codim,PiType, const IdentityGrid<HostGrid> >(this, true);
    }
 
 
    int size (int level, int codim) const {
      return hostgrid_->size(level,codim);
    }
 
    size_t numBoundarySegments () const {
      return hostgrid_->numBoundarySegments();
    }
 
    int size (int codim) const {
      return leafIndexSet().size(codim);
    }
 
 
    int size (int level, GeometryType type) const {
      return levelIndexSets_[level]->size(type);
    }
 
 
    int size (GeometryType type) const
    {
      return leafIndexSet().size(type);
    }
 
 
    const typename Traits::GlobalIdSet& globalIdSet() const {
      return globalIdSet_;
    }
 
 
    const typename Traits::LocalIdSet& localIdSet() const {
      return localIdSet_;
    }
 
 
    const typename Traits::LevelIndexSet& levelIndexSet(int level) const
    {
      if (level < 0 || level > maxLevel())
      {
        DUNE_THROW(GridError, "levelIndexSet of nonexisting level " << level << " requested!");
      }
      return *levelIndexSets_[level];
    }
 
 
    const typename Traits::LeafIndexSet& leafIndexSet() const
    {
      return leafIndexSet_;
    }
 
 
    template < class EntitySeed >
    DUNE_DEPRECATED_MSG("entityPointer() is deprecated and will be removed after the release of dune-grid 2.4. Use entity() instead to directly obtain an Entity object.")
    typename Traits::template Codim<EntitySeed::codimension>::EntityPointer
    entityPointer(const EntitySeed& seed) const
    {
      typedef typename Traits::template Codim<EntitySeed::codimension>::EntityPointer EntityPointer;
      typedef IdentityGridEntityPointer<
        EntitySeed::codimension,
        const typename Traits::Grid,
        typename HostGrid::template Codim<EntitySeed::codimension>::EntityPointer
        > EntityPointerImp;
 
      return EntityPointer(EntityPointerImp(this, hostgrid_->entity(this->getRealImplementation(seed).hostEntitySeed())));
    }
 
    template < class EntitySeed >
    typename Traits::template Codim<EntitySeed::codimension>::Entity
    entity(const EntitySeed& seed) const
    {
      typedef IdentityGridEntity<
        EntitySeed::codimension,
        HostGrid::dimension,
        const typename Traits::Grid
        > EntityImp;
 
      return EntityImp(this, hostgrid_->entity(this->getRealImplementation(seed).hostEntitySeed()));
    }
 
 
 
 
    void globalRefine (int refCount)
    {
      hostgrid_->globalRefine(refCount);
    }
 
    bool mark(int refCount, const typename Traits::template Codim<0>::Entity & e)
    {
      return hostgrid_->mark(refCount, getHostEntity<0>(e));
    }
 
    int getMark(const typename Traits::template Codim<0>::Entity & e) const
    {
      return hostgrid_->getMark(getHostEntity<0>(e));
    }
 
    bool preAdapt() {
      return hostgrid_->preAdapt();
    }
 
 
    bool adapt()
    {
      return hostgrid_->adapt();
    }
 
    void postAdapt() {
      return hostgrid_->postAdapt();
    }
 
    unsigned int overlapSize(int codim) const {
      return hostgrid_->overlapSize(codim);
    }
 
 
    unsigned int ghostSize(int codim) const {
      return hostgrid_->ghostSize(codim);
    }
 
 
    unsigned int overlapSize(int level, int codim) const {
      return hostgrid_->overlapSize(level,codim);
    }
 
 
    unsigned int ghostSize(int level, int codim) const {
      return hostgrid_->ghostSize(level,codim);
    }
 
 
#if 0
    void loadBalance(int strategy, int minlevel, int depth, int maxlevel, int minelement){
      DUNE_THROW(NotImplemented, "IdentityGrid::loadBalance()");
    }
 
    template<class T, template<class> class P, int codim>
    void communicate (T& t, InterfaceType iftype, CommunicationDirection dir, int level);
 
    template<class DataHandle>
    void communicate (DataHandle& data, InterfaceType iftype, CommunicationDirection dir, int level) const
    {}
 
    template<class DataHandle>
    void communicate (DataHandle& data, InterfaceType iftype, CommunicationDirection dir) const
    {}
#endif
 
 
    const CollectiveCommunication<IdentityGrid>& comm () const
    {
      return ccobj;
    }
 
 
    // **********************************************************
    // End of Interface Methods
    // **********************************************************
 
    HostGridType& getHostGrid() const
    {
      return *hostgrid_;
    }
 
 
    template <int codim>
    DUNE_DEPRECATED_MSG("getHostEntityPointer() is deprecated and will be removed after the release of dune-grid 2.4. Use getHostEntity() instead to obtain a reference to the host entity object.")
    typename HostGrid::Traits::template Codim<codim>::EntityPointer getHostEntityPointer(const typename Traits::template Codim<codim>::Entity& e) const
    {
      return this->getRealImplementation(e).hostEntity_;
    }
 
    template <int codim>
    const typename HostGrid::Traits::template Codim<codim>::Entity& getHostEntity(const typename Traits::template Codim<codim>::Entity& e) const
    {
      return this->getRealImplementation(e).hostEntity_;
    }
 
  protected:
 
    HostGrid* hostgrid_;
 
  private:
 
    void setIndices()
    {
      localIdSet_.update();
 
      globalIdSet_.update();
 
      // //////////////////////////////////////////
      //   Create the index sets
      // //////////////////////////////////////////
      for (int i=levelIndexSets_.size(); i<=maxLevel(); i++) {
        IdentityGridLevelIndexSet<const IdentityGrid<HostGrid> >* p
          = new IdentityGridLevelIndexSet<const IdentityGrid<HostGrid> >();
        levelIndexSets_.push_back(p);
      }
 
      for (int i=0; i<=maxLevel(); i++)
        if (levelIndexSets_[i])
          levelIndexSets_[i]->update(*this, i);
 
      leafIndexSet_.update(*this);
 
    }
 
    CollectiveCommunication<IdentityGrid> ccobj;
 
    std::vector<IdentityGridLevelIndexSet<const IdentityGrid<HostGrid> >*> levelIndexSets_;
 
    IdentityGridLeafIndexSet<const IdentityGrid<HostGrid> > leafIndexSet_;
 
    IdentityGridGlobalIdSet<const IdentityGrid<HostGrid> > globalIdSet_;
 
    IdentityGridLocalIdSet<const IdentityGrid<HostGrid> > localIdSet_;
 
  }; // end Class IdentityGrid
 
 
 
 
  namespace Capabilities
  {
    template<class HostGrid, int codim>
    struct hasEntity<IdentityGrid<HostGrid>, codim>
    {
      static const bool v = hasEntity<HostGrid,codim>::v;
    };
 
    template<class HostGrid>
    struct DUNE_DEPRECATED_MSG("Capabilities::isParallel will be removed after dune-grid-2.4.") isParallel<IdentityGrid<HostGrid> >
    {
      static const bool DUNE_DEPRECATED_MSG("Capabilities::isParallel will be removed after dune-grid-2.4.") v = isParallel<HostGrid>::v;
    };
 
    template<class HostGrid>
    struct isLevelwiseConforming<IdentityGrid<HostGrid> >
    {
      static const bool v = isLevelwiseConforming<HostGrid>::v;
    };
 
    template<class HostGrid>
    struct isLeafwiseConforming<IdentityGrid<HostGrid> >
    {
      static const bool v = isLeafwiseConforming<HostGrid>::v;
    };
  } // end namespace Capabilities
 
} // namespace Dune
 
#endif // DUNE_GRID_IDENTITYGRID_HH