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/parallel/communication.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/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;
 
  template<int dim, class HostGrid>
  struct IdentityGridFamily
  {
 
  public:
 
    typedef GridTraits<
        dim,
        HostGrid::dimensionworld,
        Dune::IdentityGrid<HostGrid>,
        IdentityGridGeometry,
        IdentityGridEntity,
        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<No_Comm>,
        DefaultLevelGridViewTraits,
        DefaultLeafGridViewTraits,
        IdentityGridEntitySeed
        > Traits;
 
  };
 
  //**********************************************************************
  //
  // --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 >
    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(seed.impl().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_->leafGridView().overlapSize(codim);
    }
 
 
    unsigned int ghostSize(int codim) const {
      return hostgrid_->leafGridView().ghostSize(codim);
    }
 
 
    unsigned int overlapSize(int level, int codim) const {
      return hostgrid_->levelGridView(level).overlapSize(codim);
    }
 
 
    unsigned int ghostSize(int level, int codim) const {
      return hostgrid_->levelGridView(level).ghostSize(codim);
    }
 
 
#if 0
    void loadBalance(int strategy, int minlevel, int depth, int maxlevel, int minelement){
      DUNE_THROW(NotImplemented, "IdentityGrid::loadBalance()");
    }
#endif
 
 
    const CollectiveCommunication<No_Comm>& comm () const
    {
      return ccobj;
    }
 
 
    // **********************************************************
    // End of Interface Methods
    // **********************************************************
 
    HostGridType& getHostGrid() const
    {
      return *hostgrid_;
    }
 
 
    template <int codim>
    const typename HostGrid::Traits::template Codim<codim>::Entity& getHostEntity(const typename Traits::template Codim<codim>::Entity& e) const
    {
      return e.impl().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<No_Comm> 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, int codim>
    struct hasEntityIterator<IdentityGrid<HostGrid>, codim>
    {
      static const bool v = hasEntityIterator<HostGrid, codim>::v;
    };
 
    template<class HostGrid, int codim>
    struct canCommunicate<IdentityGrid<HostGrid>, codim>
    {
      static const bool v = canCommunicate<HostGrid, codim>::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