3 #ifndef DUNE_GRID_YASPGRIDENTITY_HH 4 #define DUNE_GRID_YASPGRIDENTITY_HH 36 static constexpr
int evaluate(
int d,
int c)
38 return _values[_offsets[d] + c];
42 static constexpr
int binomial(
int d,
int c)
45 for (
int i=d-c+1; i<=d; i++)
47 for (
long i=2; i<=c; i++)
54 BinomialTable() =
delete;
57 static constexpr
int nextValue(
int& r,
int& c)
59 const auto result = binomial(r, c);
70 template<std::size_t... I>
71 static constexpr std::array<int,
sizeof...(I)> computeValues(std::index_sequence<I...>)
74 return {{ ((void)I, nextValue(r, c))... }};
77 template<std::size_t... I>
78 static constexpr std::array<int,
sizeof...(I)> computeOffsets(std::index_sequence<I...>)
79 {
return {{ (I*(I+1)/2)... }}; }
81 static constexpr std::array<int,(n+1)*(n+2)/2> _values = computeValues(std::make_index_sequence<(n+1)*(n+2)/2>{});
82 static constexpr std::array<int,n+1> _offsets = computeOffsets(std::make_index_sequence<n+1>{});
85 #if __cplusplus < 201703L 87 constexpr std::array<int,(n+1)*(n+2)/2> BinomialTable<n>::_values;
89 constexpr std::array<int,n+1> BinomialTable<n>::_offsets;
98 template<
int dimworld>
99 constexpr
int subEnt(
int d,
int c)
101 return (d < c ? 0 : BinomialTable<dimworld>::evaluate(d,c) << c);
106 template<
typename F,
int dim>
107 struct EntityShiftTable
109 typedef std::bitset<dim> value_type;
111 static value_type evaluate(
int i,
int codim)
113 return {_values[_offsets[codim] + i]};
119 EntityShiftTable() =
delete;
122 static constexpr
int nextOffset(
int& offset,
int codim)
129 offset += subEnt<dim>(dim, codim-1);
133 template<std::size_t... I>
134 static constexpr std::array<int,
sizeof...(I)> computeOffsets(std::index_sequence<I...>)
137 return {{ (nextOffset(offset, I))... }};
141 static constexpr
unsigned char nextValue(
int& codim,
int& i)
143 const auto result = F::evaluate(i, codim);
146 if (i >= subEnt<dim>(dim, codim)) {
154 template<std::size_t... I>
155 static constexpr std::array<
unsigned char,
sizeof...(I)> computeValues(std::index_sequence<I...>)
157 int codim = 0, i = 0;
158 return {{ ((void)I, nextValue(codim, i))... }};
161 static constexpr std::array<int,dim+1> _offsets = computeOffsets(std::make_index_sequence<dim+1>{});
162 static constexpr std::array<unsigned char,StaticPower<3,dim>::power> _values = computeValues(std::make_index_sequence<StaticPower<3,dim>::power>{});
166 #if __cplusplus < 201703L 167 template<
typename F,
int dim>
168 constexpr std::array<int,dim+1> EntityShiftTable<F, dim>::_offsets;
169 template<
typename F,
int dim>
170 constexpr std::array<unsigned char,StaticPower<3,dim>::power> EntityShiftTable<F, dim>::_values;
175 struct calculate_entity_shift
177 static constexpr
unsigned long long evaluate(
int index,
int cc)
180 for (
int d = dim; d>0; d--)
184 if (index < subEnt<dim>(d-1,cc))
185 result |= 1ull << (d-1);
188 index = (index - subEnt<dim>(d-1, cc)) % subEnt<dim>(d-1,cc-1);
205 std::bitset<dim> entityShift(
int index,
int cc)
207 return EntityShiftTable<calculate_entity_shift<dim>,dim>::evaluate(index,cc);
212 struct calculate_entity_move
214 static constexpr
unsigned long long evaluate(
int index,
int cc)
217 for (
int d = dim; d>0; d--)
222 result &= ~(1ull << (d-1));
223 result |= index & (1ull << (d-1));
225 index &= ~(1<<(d-1));
227 if (index >= subEnt<dim>(d-1,cc))
229 if ((index - subEnt<dim>(d-1,cc)) / subEnt<dim>(d-1,cc-1) == 1)
231 result |= 1ull << (d-1);
233 index = (index - subEnt<dim>(d-1, cc)) % subEnt<dim>(d-1,cc-1);
250 std::bitset<dim> entityMove(
int index,
int cc)
252 return EntityShiftTable<calculate_entity_move<dim>,dim>::evaluate(index,cc);
259 template<
int codim,
int dim,
class Gr
idImp>
261 :
public EntityDefaultImplementation <codim,dim,GridImp,YaspEntity>
264 template<
int, PartitionIteratorType,
typename>
268 typedef typename GridImp::ctype
ctype;
270 typedef typename GridImp::template Codim<codim>::Geometry
Geometry;
271 typedef typename GridImp::Traits::template Codim<codim>::GeometryImpl
GeometryImpl;
273 typedef typename GridImp::template Codim<codim>::EntitySeed
EntitySeed;
292 GeometryImpl _geometry(_it.lowerleft(),_it.upperright(),_it.shift());
293 return Geometry(_geometry);
302 return Dune::Yasp::subEnt<dim>(dim-codim,cc-codim);
308 if (_g->interior[codim].inside(_it.coord(),_it.shift()))
310 if (_g->interiorborder[codim].inside(_it.coord(),_it.shift()))
312 if (_g->overlap[codim].inside(_it.coord(),_it.shift()))
314 if (_g->overlapfront[codim].inside(_it.coord(),_it.shift()))
319 typedef typename GridImp::YGridLevelIterator
YGLI;
320 typedef typename GridImp::YGrid::Iterator
I;
329 : _it(
std::move(it)), _g(
std::move(g))
335 return _it == e._it && _g == e._g;
348 std::array<int,dim> size;
350 for (
int i=0; i<dim; i++)
353 size[i] = _g->mg->levelSize(_g->level(), i);
359 PersistentIndexType
id(_it.shift().to_ulong());
363 id =
id+PersistentIndexType(_g->level());
366 for (
int i=dim-1; i>=0; i--)
369 id =
id+PersistentIndexType(_it.coord(i));
378 return _it.superindex();
386 std::bitset<dim-codim> subent_shift = Dune::Yasp::entityShift<dim-codim>(i,cc-codim);
387 std::bitset<dim-codim> subent_move = Dune::Yasp::entityMove<dim-codim>(i,cc-codim);
389 std::bitset<dim> shift = _it.shift();
390 std::array<int, dim> coord = _it.coord();
391 for (
int j=0, k=0; j<dim; j++)
396 coord[j] += subent_move[k];
397 shift[j] = subent_shift[k];
401 int which = _g->overlapfront[cc].shiftmapping(shift);
402 return _g->overlapfront[cc].superindex(coord,which);
409 const GridImp *
yaspgrid()
const {
return _g->mg; }
417 template<
int dim,
class Gr
idImp>
421 enum { dimworld = GridImp::dimensionworld };
423 typedef typename GridImp::Traits::template Codim< 0 >::GeometryImpl GeometryImpl;
425 template<
int, PartitionIteratorType,
typename>
432 typedef typename GridImp::ctype
ctype;
434 typedef typename GridImp::YGridLevelIterator
YGLI;
435 typedef typename GridImp::YGrid::Iterator
I;
437 typedef typename GridImp::template Codim< 0 >::Geometry
Geometry;
447 typedef typename GridImp::template Codim<0>::EntitySeed
EntitySeed;
457 typedef typename GridImp::YGrid::iTupel
iTupel;
468 : _it(
std::move(it)), _g(g)
472 : _it(
std::move(it)), _g(
std::move(g))
478 return _it == e._it && _g == e._g;
482 int level ()
const {
return _g->level(); }
494 if (_g->interior[0].inside(_it.coord(),_it.shift()))
496 if (_g->overlap[0].inside(_it.coord(),_it.shift()))
498 DUNE_THROW(
GridError,
"Impossible GhostEntity");
505 auto ll = _it.lowerleft();
506 auto ur = _it.upperright();
509 for (
int i=0; i<dimworld; i++) {
510 if (gridlevel()->mg->isPeriodic(i)) {
511 int coord = transformingsubiterator().coord(i);
513 auto size = _g->mg->domainSize()[i];
516 }
else if (coord + 1 > gridlevel()->mg->levelSize(gridlevel()->level(),i)) {
517 auto size = _g->mg->domainSize()[i];
524 GeometryImpl _geometry(ll,ur);
525 return Geometry( _geometry );
534 return Dune::Yasp::subEnt<dim>(dim,cc);
543 return Dune::Yasp::subEnt<dim>(dim,codim);
552 std::bitset<dim> move = Dune::Yasp::entityMove<dim>(i,cc);
555 iTupel coord = _it.coord();
556 for (
int j=0; j<dim; j++)
560 int which = _g->overlapfront[cc].shiftmapping(Dune::Yasp::entityShift<dim>(i,cc));
569 DUNE_THROW(
GridError,
"tried to call father on level 0");
576 iTupel coord = _it.coord();
579 for (
int k=0; k<dim; k++) coord[k] = coord[k]/2;
587 return (_g->level()>0);
595 FieldVector<ctype,dim> ll(0.0),ur(0.5);
597 for (
int k=0; k<dim; k++)
613 const GridImp*
yaspgrid ()
const {
return _g->mg; }
617 return (_g->level() == yaspgrid()->maxLevel());
622 bool isNew ()
const {
return yaspgrid()->adaptRefCount > 0 && yaspgrid()->maxLevel() < _g->level() + yaspgrid()->adaptRefCount; }
648 IntersectionIterator
iend ()
const 669 HierarchicIterator
hbegin (
int maxlevel)
const 675 HierarchicIterator
hend (
int maxlevel)
const 687 PersistentIndexType persistentIndex ()
const 690 PersistentIndexType
id(_it.shift().to_ulong());
694 id =
id+PersistentIndexType(_g->level());
698 for (
int i=dim-1; i>=0; i--)
701 id =
id+PersistentIndexType(_it.coord(i));
708 int compressedIndex ()
const 710 return _it.superindex();
714 PersistentIndexType subPersistentIndex (
int i,
int cc)
const 717 std::bitset<dim> shift = Dune::Yasp::entityShift<dim>(i,cc);
718 std::bitset<dim> move = Dune::Yasp::entityMove<dim>(i,cc);
720 int trailing = (cc == dim) ? 1000 : 0;
722 std::array<int,dim> size = _g->mg->levelSize(_g->level());
723 std::array<int, dim> coord = _it.coord();
724 for (
int j=0; j<dim; j++)
735 for (
int j=0; j<dim; j++)
741 for (
int k=0; k<_g->level(); k++)
742 if (coord[j] & (1<<k))
746 trailing =
std::min(trailing,zeroes);
751 PersistentIndexType
id(shift.to_ulong());
755 id =
id+PersistentIndexType(_g->level()-trailing);
758 for (
int j=dim-1; j>=0; j--)
761 id =
id+PersistentIndexType(coord[j]>>trailing);
768 int subCompressedIndex (
int i,
int cc)
const 771 std::bitset<dim> shift = Dune::Yasp::entityShift<dim>(i,cc);
772 std::bitset<dim> move = Dune::Yasp::entityMove<dim>(i,cc);
774 std::array<int, dim> coord = _it.coord();
775 for (
int j=0; j<dim; j++)
778 int which = _g->overlapfront[cc].shiftmapping(shift);
779 return _g->overlapfront[cc].superindex(coord,which);
788 template<
int dim,
class Gr
idImp>
792 enum { dimworld = GridImp::dimensionworld };
794 template<
int, PartitionIteratorType,
typename>
797 typedef typename GridImp::Traits::template Codim<dim>::GeometryImpl GeometryImpl;
800 typedef typename GridImp::ctype
ctype;
802 typedef typename GridImp::YGridLevelIterator
YGLI;
803 typedef typename GridImp::YGrid::Iterator
I;
805 typedef typename GridImp::template Codim<dim>::Geometry
Geometry;
807 typedef typename GridImp::template Codim<dim>::EntitySeed
EntitySeed;
813 typedef typename GridImp::YGrid::iTupel
iTupel;
824 : _it(
std::move(it)), _g(
std::move(g))
830 return _it == e._it && _g == e._g;
834 int level ()
const {
return _g->level();}
849 return Dune::Yasp::subEnt<dim>(dim-dim,cc-dim);
854 GeometryImpl _geometry((_it).lowerleft());
855 return Geometry( _geometry );
861 if (_g->interior[dim].inside(_it.coord(),_it.shift()))
863 if (_g->interiorborder[dim].inside(_it.coord(),_it.shift()))
865 if (_g->overlap[dim].inside(_it.coord(),_it.shift()))
867 if (_g->overlapfront[dim].inside(_it.coord(),_it.shift()))
875 return compressedIndex();
885 PersistentIndexType persistentIndex ()
const 888 iTupel size = _g->mg->levelSize(_g->level());
890 for (
int i=0; i<dim; i++)
898 for (
int i=0; i<dim; i++)
902 for (
int j=0; j<_g->level(); j++)
903 if (_it.coord(i)&(1<<j))
907 trailing =
std::min(trailing,zeros);
911 int level = _g->level()-trailing;
914 PersistentIndexType
id(0);
918 id =
id+PersistentIndexType(level);
921 for (
int i=dim-1; i>=0; i--)
924 id =
id+PersistentIndexType(_it.coord(i)>>trailing);
931 int compressedIndex ()
const {
return _it.superindex();}
939 const GridImp *
yaspgrid()
const {
return _g->mg; }
947 #endif // DUNE_GRID_YASPGRIDENTITY_HH Geometry geometry() const
geometry of this entity
Definition: yaspgridentity.hh:853
LeafIntersectionIterator ileafend() const
Reference to one past the last neighbor.
Definition: yaspgridentity.hh:654
GridImp::YGrid::Iterator I
Definition: yaspgridentity.hh:320
const GridImp * yaspgrid() const
Definition: yaspgridentity.hh:409
YaspEntity(YGLI &&g, const I &&it)
Definition: yaspgridentity.hh:328
PartitionType partitionType() const
return partition type attribute
Definition: yaspgridentity.hh:306
IdType id(const typename std::remove_const< GridImp >::type::Traits::template Codim< cd >::Entity &e) const
get id of an entity
Definition: yaspgrididset.hh:42
persistent, globally unique Ids
Definition: yaspgrid.hh:67
YGLI & gridlevel()
Definition: yaspgridentity.hh:408
PartitionType partitionType() const
return partition type attribute
Definition: yaspgridentity.hh:492
GridImp::LevelIntersectionIterator LevelIntersectionIterator
Definition: yaspgridentity.hh:449
GridImp::template Codim< 0 >::Entity Entity
Definition: yaspgridentity.hh:446
GridImp::HierarchicIterator HierarchicIterator
Definition: yaspgridentity.hh:451
Default Implementations for EntityImp.
Definition: common/entity.hh:547
GridImp::YGrid::Iterator I
Definition: yaspgridentity.hh:435
unsigned int subEntities(unsigned int codim) const
Definition: yaspgridentity.hh:541
bool mightVanish() const
Returns true, if entity might disappear during the next call to adapt()
Definition: yaspgridentity.hh:626
GridImp::YGrid::iTupel iTupel
define type used for coordinates in grid module
Definition: yaspgridentity.hh:457
YaspEntity(YGLI &&g, I &&it)
Definition: yaspgridentity.hh:823
YaspIntersectionIterator enables iteration over intersections with neighboring codim 0 entities...
Definition: yaspgrid.hh:63
const I & transformingsubiterator() const
Definition: yaspgridentity.hh:405
GridImp::YGridLevelIterator YGLI
Definition: yaspgridentity.hh:319
I & transformingsubiterator()
Definition: yaspgridentity.hh:936
bool isNew() const
Returns true, if the entity has been created during the last call to adapt()
Definition: yaspgridentity.hh:622
I & transformingsubiterator()
Definition: yaspgridentity.hh:407
int subCompressedIndex(int i, unsigned int cc) const
subentity compressed index
Definition: yaspgridentity.hh:382
PersistentIndexType persistentIndex() const
globally unique, persistent index
Definition: yaspgridentity.hh:345
GridImp::PersistentIndexType PersistentIndexType
Definition: yaspgridentity.hh:342
all interior entities
Definition: gridenums.hh:29
Implementation of Level- and LeafIndexSets for YaspGrid.
Definition: yaspgrid.hh:66
Entity father() const
Inter-level access to father element on coarser grid. Assumes that meshes are nested.
Definition: yaspgridentity.hh:565
Base class for exceptions in Dune grid modules.
Definition: exceptions.hh:16
const int yaspgrid_dim_bits
Definition: yaspgrid.hh:51
LeafIntersectionIterator ileafbegin() const
returns intersection iterator for first intersection
Definition: yaspgridentity.hh:635
GridImp::LeafIntersectionIterator LeafIntersectionIterator
Definition: yaspgridentity.hh:450
YaspEntity(const YGLI &g, const I &it)
Definition: yaspgridentity.hh:324
Geometry geometry() const
geometry of this entity
Definition: yaspgridentity.hh:290
HierarchicIterator hend(int maxlevel) const
Returns iterator to one past the last son.
Definition: yaspgridentity.hh:675
int level() const
level of this element
Definition: yaspgridentity.hh:276
GridImp::template Codim< 0 >::EntitySeed EntitySeed
Definition: yaspgridentity.hh:447
GridImp::ctype ctype
Definition: yaspgridentity.hh:268
Definition: yaspgrid.hh:60
bool isLeaf() const
Definition: yaspgridentity.hh:615
GridImp::YGridLevelIterator YGLI
Definition: yaspgridentity.hh:802
YGLI _g
Definition: yaspgridentity.hh:412
bool equals(const YaspEntity &e) const
Return true when two iterators over the same grid are equal (!).
Definition: yaspgridentity.hh:476
const GridImp * yaspgrid() const
Definition: yaspgridentity.hh:613
PartitionType partitionType() const
return partition type attribute
Definition: yaspgridentity.hh:859
Geometry geometry() const
geometry of this entity
Definition: yaspgridentity.hh:503
YGLI & gridlevel()
Definition: yaspgridentity.hh:612
bool equals(const YaspEntity &e) const
Return true when two iterators over the same grid are equal (!).
Definition: yaspgridentity.hh:828
GridImp::PersistentIndexType PersistentIndexType
define the type used for persistent indices
Definition: yaspgridentity.hh:810
int count() const
Definition: yaspgridentity.hh:532
IntersectionIterator ibegin() const
returns intersection iterator for first intersection
Definition: yaspgridentity.hh:629
GridImp::ctype ctype
Definition: yaspgridentity.hh:800
int subCompressedIndex(int, unsigned int) const
subentity compressed index simply returns compressedIndex
Definition: yaspgridentity.hh:873
unsigned int subEntities(unsigned int cc) const
Definition: yaspgridentity.hh:847
YaspHierarchicIterator enables iteration over son entities of codim 0.
Definition: yaspgrid.hh:65
const YGLI & gridlevel() const
Definition: yaspgridentity.hh:935
GridImp::template Codim< 0 >::LocalGeometry LocalGeometry
Definition: yaspgridentity.hh:438
Iterates over entities of one grid level.
Definition: yaspgrid.hh:62
The general version that handles all codimensions but 0 and dim.
Definition: yaspgrid.hh:59
I & transformingsubiterator()
Definition: yaspgridentity.hh:611
LevelIntersectionIterator ilevelbegin() const
returns intersection iterator for first intersection
Definition: yaspgridentity.hh:642
int level() const
level of this element
Definition: yaspgridentity.hh:482
YaspEntity()
Definition: yaspgridentity.hh:321
Codim< cc >::Entity subEntity(int i) const
Definition: yaspgridentity.hh:549
const YGLI & gridlevel() const
Definition: yaspgridentity.hh:406
YaspEntity()
Definition: yaspgridentity.hh:460
EntitySeed seed() const
Return the entity seed which contains sufficient information to generate the entity again and uses as...
Definition: yaspgridentity.hh:839
int compressedIndex() const
consecutive, codim-wise, level-wise index
Definition: yaspgridentity.hh:376
all entities lying in the overlap zone
Definition: gridenums.hh:31
GridImp::ctype ctype
Definition: yaspgridentity.hh:432
GridImp::LevelIntersectionIterator IntersectionIterator
Definition: yaspgridentity.hh:448
EntitySeed seed() const
Return the entity seed which contains sufficient information to generate the entity again and uses as...
Definition: yaspgridentity.hh:284
YGLI _g
Definition: yaspgridentity.hh:942
int min(const DofVectorPointer< int > &dofVector)
Definition: dofvector.hh:346
LocalGeometry geometryInFather() const
Definition: yaspgridentity.hh:592
GridImp::PersistentIndexType PersistentIndexType
define the type used for persistent indices
Definition: yaspgridentity.hh:454
YaspEntity(const YGLI &g, I &&it)
Definition: yaspgridentity.hh:467
YaspEntity(const YGLI &g, const I &it)
Definition: yaspgridentity.hh:819
Describes the minimal information necessary to create a fully functional YaspEntity.
Definition: yaspgrid.hh:61
GridImp::template Codim< dim >::Geometry Geometry
Definition: yaspgridentity.hh:805
PartitionType
Attributes used in the generic overlap model.
Definition: gridenums.hh:28
ghost entities
Definition: gridenums.hh:33
GridImp::template Codim< codim >::Geometry Geometry
Definition: yaspgridentity.hh:270
on boundary between overlap and ghost
Definition: gridenums.hh:32
const I & transformingsubiterator() const
Definition: yaspgridentity.hh:934
const int yaspgrid_level_bits
Definition: yaspgrid.hh:52
I _it
Definition: yaspgridentity.hh:411
YaspEntity(const YGLI &g, const I &it)
Definition: yaspgridentity.hh:463
bool hasFather() const
returns true if father entity exists
Definition: yaspgridentity.hh:585
HierarchicIterator hbegin(int maxlevel) const
Definition: yaspgridentity.hh:669
int level() const
level of this element
Definition: yaspgridentity.hh:834
GridImp::YGrid::Iterator I
Definition: yaspgridentity.hh:803
bool equals(const YaspEntity &e) const
Return true when two iterators over the same grid are equal (!).
Definition: yaspgridentity.hh:333
GridImp::YGrid::iTupel iTupel
define type used for coordinates in grid module
Definition: yaspgridentity.hh:813
const I & transformingsubiterator() const
Definition: yaspgridentity.hh:609
GridImp::template Codim< 0 >::Geometry Geometry
Definition: yaspgridentity.hh:437
GridImp::template Codim< cd >::Entity Entity
Definition: yaspgridentity.hh:443
LevelIntersectionIterator ilevelend() const
Reference to one past the last neighbor.
Definition: yaspgridentity.hh:660
YGLI & gridlevel()
Definition: yaspgridentity.hh:937
YaspEntity()
Definition: yaspgridentity.hh:816
const GridImp * yaspgrid() const
Definition: yaspgridentity.hh:939
GridImp::Traits::template Codim< codim >::GeometryImpl GeometryImpl
Definition: yaspgridentity.hh:271
EntitySeed seed() const
Return the entity seed which contains sufficient information to generate the entity again and uses as...
Definition: yaspgridentity.hh:487
IntersectionIterator iend() const
Reference to one past the last neighbor.
Definition: yaspgridentity.hh:648
GridImp::template Codim< dim >::EntitySeed EntitySeed
Definition: yaspgridentity.hh:807
YaspEntity(YGLI &&g, I &&it)
Definition: yaspgridentity.hh:471
GridImp::YGridLevelIterator YGLI
Definition: yaspgridentity.hh:434
const YGLI & gridlevel() const
Definition: yaspgridentity.hh:610
unsigned int subEntities(unsigned int cc) const
Definition: yaspgridentity.hh:300
on boundary between interior and overlap
Definition: gridenums.hh:30
GridImp::template Codim< codim >::EntitySeed EntitySeed
Definition: yaspgridentity.hh:273
I _it
Definition: yaspgridentity.hh:941
Include standard header files.
Definition: agrid.hh:58