dune-grid  2.6-git
yaspgridindexsets.hh
Go to the documentation of this file.
1 // -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
2 // vi: set et ts=4 sw=2 sts=2:
3 #ifndef DUNE_GRID_YASPGRIDINDEXSET_HH
4 #define DUNE_GRID_YASPGRIDINDEXSET_HH
5 
12 namespace Dune {
13 
20  template<class GridImp, bool isLeafIndexSet>
21  class YaspIndexSet
22  : public IndexSet< GridImp, YaspIndexSet< GridImp, isLeafIndexSet >, unsigned int >
23  {
24  typedef YaspIndexSet< GridImp, isLeafIndexSet > This;
25  typedef IndexSet< GridImp, This, unsigned int > Base;
26 
27  public:
28  typedef typename Base::IndexType IndexType;
29 
30  using Base::subIndex;
31 
33  YaspIndexSet ( const GridImp &g, int l )
34  : grid( g ),
35  level( l )
36  {
37  assert(not isLeafIndexSet);
38 
39  // contains a single element type;
40  for (int codim=0; codim<=GridImp::dimension; codim++)
41  mytypes[codim].push_back(GeometryTypes::cube(GridImp::dimension-codim));
42  }
43 
45  YaspIndexSet ( const GridImp &g )
46  : grid( g )
47  {
48  assert(isLeafIndexSet);
49 
50  // contains a single element type;
51  for (int codim=0; codim<=GridImp::dimension; codim++)
52  mytypes[codim].push_back(GeometryTypes::cube(GridImp::dimension-codim));
53  }
54 
56  template<int cc>
57  IndexType index (const typename std::remove_const<GridImp>::type::Traits::template Codim<cc>::Entity& e) const
58  {
59  return grid.getRealImplementation(e).compressedIndex();
60  }
61 
63  template< int cc >
64  IndexType subIndex ( const typename std::remove_const< GridImp >::type::Traits::template Codim< cc >::Entity &e,
65  int i, unsigned int codim ) const
66  {
67  return grid.getRealImplementation(e).subCompressedIndex(i, codim);
68  }
69 
71  int size (GeometryType type) const
72  {
73  return (isLeafIndexSet)
74  ? grid.size( type )
75  : grid.size( level, type );
76  }
77 
79  int size (int codim) const
80  {
81  return (isLeafIndexSet)
82  ? grid.size( codim )
83  : grid.size( level, codim );
84  }
85 
87  template<class EntityType>
88  bool contains (const EntityType& e) const
89  {
90  return (isLeafIndexSet)
91  ? e.level() == grid.maxLevel()
92  : e.level() == level;
93  }
94 
95  std::vector< GeometryType > types ( int codim ) const { return mytypes[ codim ]; }
96 
98  const std::vector<GeometryType>& geomTypes (int codim) const
99  {
100  return mytypes[codim];
101  }
102 
103  private:
104  const GridImp& grid;
105  int level;
106  std::vector<GeometryType> mytypes[std::remove_const<GridImp>::type::dimension+1];
107  };
108 
109 } // namespace Dune
110 
111 #endif // DUNE_GRID_YASPGRIDINDEXSET_HH
Include standard header files.
Definition: agrid.hh:58
YaspIndexSet(const GridImp &g)
Level grid view constructor stores reference to a grid and level.
Definition: yaspgridindexsets.hh:45
const std::vector< GeometryType > & geomTypes(int codim) const
deliver all geometry types used in this grid
Definition: yaspgridindexsets.hh:98
Base::IndexType IndexType
Definition: yaspgridindexsets.hh:28
int size(int codim) const
return size of set for a given codim
Definition: yaspgridindexsets.hh:79
std::vector< GeometryType > types(int codim) const
Definition: yaspgridindexsets.hh:95
IndexType subIndex(const typename Traits::template Codim< cc >::Entity &e, int i, unsigned int codim) const
Map a subentity to an index.
Definition: indexidset.hh:151
int size(GeometryType type) const
get number of entities of given type and level (the level is known to the object) ...
Definition: yaspgridindexsets.hh:71
bool contains(const EntityType &e) const
return true if the given entity is contained in .
Definition: yaspgridindexsets.hh:88
IndexTypeImp IndexType
The type used for the indices.
Definition: indexidset.hh:90
GeometryType
Type representing VTK&#39;s entity geometry types.
Definition: common.hh:178
IndexType subIndex(const typename std::remove_const< GridImp >::type::Traits::template Codim< cc >::Entity &e, int i, unsigned int codim) const
get index of subentity of an entity
Definition: yaspgridindexsets.hh:64
IndexType index(const typename std::remove_const< GridImp >::type::Traits::template Codim< cc >::Entity &e) const
get index of an entity
Definition: yaspgridindexsets.hh:57
YaspIndexSet(const GridImp &g, int l)
Level grid view constructor stores reference to a grid and level.
Definition: yaspgridindexsets.hh:33