3 #ifndef DUNE_ALU2DGRID_INTERSECTION_HH
4 #define DUNE_ALU2DGRID_INTERSECTION_HH
21 template<
int cd,
int dim,
class Gr
idImp>
22 class ALU2dGridEntity;
23 template<
int cd, PartitionIteratorType pitype,
class Gr
idImp >
24 class ALU2dGridLevelIterator;
25 template<
int cd,
class Gr
idImp >
26 class ALU2dGridEntityPointer;
27 template<
int mydim,
int coorddim,
class Gr
idImp>
28 class ALU2dGridGeometry;
29 template<
class Gr
idImp>
30 class ALU2dGridHierarchicIterator;
31 template<
class Gr
idImp>
32 class ALU2dGridIntersectionBase;
33 template<
class Gr
idImp>
34 class ALU2dGridLeafIntersectionIterator;
35 template<
class Gr
idImp>
36 class ALU2dGridLevelIntersectionIterator;
37 template<
int codim, PartitionIteratorType pitype,
class Gr
idImp>
38 class ALU2dGridLeafIterator;
39 template<
int dim,
int dimworld, ALU2DSPACE ElementType eltype >
47 template<
class Gr
idImp,
class LocalGeometryImpl >
53 LocalGeometryImpl geoms_[ 2 ][ 4 ][ 2 ];
55 #ifdef USE_SMP_PARALLEL
63 const LocalGeometryImpl &
localGeom (
const int aluFace,
const int twist,
const int corners )
const
65 assert( corners == 3 || corners == 4 );
66 assert( 0 <= aluFace && aluFace < corners );
67 assert( twist == 0 || twist == 1 );
68 return geoms_[ corners-3 ][ aluFace ][ twist ];
74 #ifdef USE_SMP_PARALLEL
76 static std::vector< ThisType > storage( GridObjectFactoryType :: maxThreads() );
77 return storage[ GridObjectFactoryType :: threadNumber () ];
100 template<
class Gr
idImp>
101 class ALU2dGridIntersectionBase
103 static const int dim = GridImp::dimension;
104 static const int dimworld = GridImp::dimensionworld;
107 typedef typename GridImp::Traits::template Codim< 1 >::GeometryImpl GeometryImpl;
108 typedef typename GridImp::Traits::template Codim< 1 >::LocalGeometryImpl LocalGeometryImpl;
120 typedef typename GridImp::template Codim<0>::Entity
Entity;
123 typedef typename GridImp::template Codim<1>::Geometry
Geometry;
167 return inside()->neighbour(
index_ )->thinis( ThinelementType::bndel_like );
189 nFaces_ = (inside != 0 ? inside->numfaces() : 0);
224 return reinterpret_cast< const Intersection &
>( *this );
323 template<
class Gr
idImp >
330 static const int dim = GridImp::dimension;
331 static const int dimworld = GridImp::dimensionworld;
343 typedef std::pair< HElementType *, int > IntersectionInfo;
352 typedef typename GridImp::template Codim<0>::Entity
Entity;
354 typedef typename GridImp::template Codim<1>::Geometry
Geometry;
396 void setFirstItem(
const HElementType & elem,
int wLevel);
399 template <
class EntityType>
400 void first(
const EntityType & en,
int wLevel);
402 void addNeighboursToStack();
404 static int getOppositeInFather (
int nrInChild,
int nrOfChild );
405 static int getOppositeInChild (
int nrInFather,
int nrOfChild );
407 void setupIntersection ();
416 mutable std::stack<IntersectionInfo> nbStack_;
428 template<
class Gr
idImp >
429 class ALU2dGridLeafIntersectionIterator
430 :
public ALU2dGridIntersectionBase< GridImp >
432 typedef ALU2dGridLeafIntersectionIterator< GridImp > ThisType;
433 typedef ALU2dGridIntersectionBase< GridImp > BaseType;
438 static const int dim = GridImp::dimension;
439 static const int dimworld = GridImp::dimensionworld;
455 typedef std::pair< HElementType *, int > IntersectionInfo;
458 typedef typename GridImp::template Codim<0>::Entity
Entity;
459 typedef typename GridImp::template Codim<1>::Geometry
Geometry;
499 void setFirstItem(
const HElementType & elem,
int wLevel);
502 template <
class EntityType>
503 void first(
const EntityType & en,
int wLevel);
505 void setupIntersection ();
514 std::stack<IntersectionInfo> nbStack_;
521 #if COMPILE_ALU2DGRID_INLINE
522 #include "intersection.cc"
525 #endif // #ifndef DUNE_ALU2DGRID_INTERSECTION_HH
Geometry geometry() const
Definition: intersection_imp.cc:291
void increment()
increment iterator
Definition: intersection_imp.cc:374
const LocalGeometryStorageType & localGeomStorage_
Definition: alugrid/2d/intersection.hh:309
ALU2dImplTraits< dimworld, eltype >::ThinelementType ThinelementType
Definition: alugrid/2d/intersection.hh:131
HmeshType::hbndel_t HBndElType
Definition: alu2dinclude.hh:109
ALU2dGridIntersectionBase< GridImp > ImplementationType
Definition: alugrid/2d/intersection.hh:113
ALU2dImplTraits< dimworld, eltype >::HBndElType HBndElType
Definition: alugrid/2d/intersection.hh:133
bool isBoundary() const
Definition: alugrid/2d/intersection.hh:164
FieldVector< alu2d_ctype, dim-1 > LocalCoordinate
Definition: alugrid/2d/intersection.hh:127
int opposite() const
Definition: alugrid/2d/intersection.hh:181
impl(HElementType *inside=0)
Definition: alugrid/2d/intersection.hh:147
struct Dune::ALU2dGridIntersectionBase::impl current
ElementType
Definition: alu2dinclude.hh:55
LocalGeometry geometryInOutside() const
Definition: intersection_imp.cc:268
#define ALU2DSPACE
Definition: alu2dinclude.hh:34
void assign(const ALU2dGridLevelIntersectionIterator< GridImp > &org)
The copy constructor.
Definition: intersection_imp.cc:366
GeometryType
Type representing VTK's entity geometry types.
Definition: common.hh:178
ALU2dGridEntityPointer< 0, GridImp > EntityPointer
Definition: alugrid/2d/intersection.hh:360
bool isConform() const
Definition: alugrid/2d/intersection.hh:491
size_t boundarySegmentIndex() const
return the boundary segment index
Definition: intersection_imp.cc:118
ALU2dGridGeometry< dim-1, dimworld, GridImp > GeometryImp
Definition: alugrid/2d/intersection.hh:357
int boundaryId() const
return boundary type
Definition: intersection_imp.cc:110
int walkLevel_
Definition: alugrid/2d/intersection.hh:311
const FactoryType & factory_
Definition: alugrid/2d/intersection.hh:308
HmeshType::helement_t HElementType
Definition: alu2dinclude.hh:108
GridImp::template Codim< 0 >::Entity Entity
Definition: alugrid/2d/intersection.hh:352
Wrapper class for pointers to entities.
Definition: common/entitypointer.hh:91
Definition: alugrid/2d/intersection.hh:350
GridImp::template Codim< 0 >::Entity Entity
Definition: alugrid/2d/intersection.hh:458
GeometryImpl intersectionGlobal_
Definition: alugrid/2d/intersection.hh:303
virtual bool conforming() const =0
Class that wraps IntersectionIteratorImp of a grid and gets it's internal object from a object stack ...
Definition: intersectioniteratorwrapper.hh:18
ALU2dGridGeometry< dim-1, dimworld, GridImp > GeometryImp
Definition: alugrid/2d/intersection.hh:462
void unsetUp2Date()
Definition: intersection_imp.cc:317
void done(const HElementType *inside)
Definition: intersection_imp.cc:147
HElementType & getItem() const
Definition: alugrid/2d/entity.hh:479
Definition: alugrid/2d/entity.hh:30
Wrapper class for geometries.
Definition: common/geometry.hh:101
NormalType outerNormal(const LocalCoordinate &local) const
Definition: intersection_imp.cc:203
GridImp::template Codim< 1 >::LocalGeometry LocalGeometry
Definition: alugrid/2d/intersection.hh:124
void assign(const ThisType &org)
The copy constructor.
Definition: intersection_imp.cc:50
GridImp::template Codim< 0 >::EntityPointer EntityPointer
Definition: alugrid/2d/intersection.hh:121
bool conforming() const
leaf is conforming, when conform grid version used
Definition: alugrid/2d/intersection.hh:485
virtual void setFirstItem(const HElementType &elem, int wLevel)=0
Definition: alugrid/2d/intersection.hh:145
int indexInInside() const
local index of codim 1 entity in self where intersection is contained in
Definition: intersection_imp.cc:165
Definition: objectfactory.hh:23
ALU2dGridEntity< 0, dim, GridImp > EntityImp
Definition: alugrid/2d/intersection.hh:356
FieldVector< alu2d_ctype, dimworld > NormalType
Definition: alugrid/2d/intersection.hh:359
ALU2dGridGeometry< dim-1, dim, GridImp > LocalGeometryImp
Definition: alugrid/2d/intersection.hh:463
bool isConform() const
Definition: alugrid/2d/intersection.hh:387
LocalGeometryImpl intersectionNeighborLocal_
Definition: alugrid/2d/intersection.hh:305
ALU2dGridEntity< 0, dim, GridImp > EntityImp
Definition: alugrid/2d/intersection.hh:461
HElementType * inside() const
Definition: alugrid/2d/intersection.hh:154
FieldVector< alu2d_ctype, dimworld > NormalType
Definition: alugrid/2d/intersection.hh:126
Class that wraps IntersectionIteratorImp of a grid and gets it's internal object from a object stack ...
Definition: intersectioniteratorwrapper.hh:238
Definition: alugrid/2d/entity.hh:32
NormalType unitOuterNormal(const LocalCoordinate &local) const
Definition: intersection_imp.cc:235
ALU2dGridIntersectionBase(const FactoryType &factory, int wLevel)
The default Constructor , creating an empty ALU2dGridIntersectionIterator.
Definition: intersection_imp.cc:28
MakeableInterfaceObject< Geometry > GeometryObject
Definition: alugrid/2d/intersection.hh:467
int level() const
return level of inside(entity)
Definition: intersection_imp.cc:71
Definition: alugrid/2d/grid.hh:53
const Intersection & dereference() const
Definition: alugrid/2d/intersection.hh:222
ALU2dGridLeafIntersectionIterator(const FactoryType &factory, int wLevel)
The default Constructor , createing an empty ALU2dGridIntersectionIterator.
Definition: intersection_imp.cc:420
bool equals(const ThisType &i) const
check whether entities are the same or whether iterator is done
Definition: intersection_imp.cc:63
Dune::Intersection< GridImp, Dune::ALU2dGridIntersectionBase< GridImp > > Intersection
type of the intersection
Definition: alugrid/2d/intersection.hh:115
void increment()
increment iterator
Definition: intersection_imp.cc:445
GridImp::template Codim< 1 >::Geometry Geometry
Definition: alugrid/2d/intersection.hh:123
Definition: alugrid/2d/entity.hh:26
GridImp::GridObjectFactoryType FactoryType
Definition: alugrid/2d/intersection.hh:111
bool boundary() const
return true if intersection is with boundary
Definition: intersection_imp.cc:104
Definition: objectfactory.hh:26
Class that wraps IntersectionIteratorImp of a grid and gets it's internal object from a object stack ...
Definition: intersectioniteratorwrapper.hh:333
int indexInOutside() const
local index of codim 1 entity in neighbor where intersection is contained in
Definition: intersection_imp.cc:176
double alu2d_ctype
Definition: alu2dinclude.hh:63
int index_
Definition: alugrid/2d/intersection.hh:206
MakeableInterfaceObject< Geometry > GeometryObject
Definition: alugrid/2d/intersection.hh:362
HElementType * outside() const
Definition: alugrid/2d/intersection.hh:176
GridImp::template Codim< 0 >::Entity Entity
Definition: alugrid/2d/intersection.hh:120
void checkValid()
return true if intersection is with boundary
Definition: intersection_imp.cc:92
Definition: alugrid/2d/intersection.hh:48
virtual ~ALU2dGridIntersectionBase()
Definition: alugrid/2d/intersection.hh:217
GridImp::template Codim< 1 >::LocalGeometry LocalGeometry
Definition: alugrid/2d/intersection.hh:355
Provides proxy classes for IntersectionsIterators.
organize the memory management for entitys used by the NeighborIterator
Definition: alugrid/2d/grid.hh:69
ALU2dGridLevelIntersectionIterator(const FactoryType &factory, int wLevel)
The default Constructor , creating an empty ALU2dGridIntersectionIterator.
Definition: intersection_imp.cc:349
GridImp::GridObjectFactoryType FactoryType
Definition: alugrid/2d/intersection.hh:346
void first(const EntityImp &en, int wLevel)
reset IntersectionIterator to first neighbour
Definition: intersection_imp.cc:81
bool conforming() const
Definition: alugrid/2d/intersection.hh:381
FieldVector< alu2d_ctype, dimworld > NormalType
Definition: alugrid/2d/intersection.hh:464
int nFaces() const
Definition: alugrid/2d/intersection.hh:159
const GridImp & grid() const
Definition: alugrid/2d/intersection.hh:282
Intersection of a mesh entities of codimension 0 ("elements") with a "neighboring" element or with th...
Definition: albertagrid/dgfparser.hh:26
GridImp::GridObjectFactoryType FactoryType
Definition: alugrid/2d/intersection.hh:443
ALU2dGridEntityPointer< 0, GridImp > EntityPointerImp
Definition: alugrid/2d/intersection.hh:129
Definition: alugrid/2d/intersection.hh:118
LocalGeometry geometryInInside() const
Definition: intersection_imp.cc:245
int twistInNeighbor() const
Definition: alugrid/2d/intersection.hh:266
ALU2DIntersectionGeometryStorage< GridImp, LocalGeometryImpl > LocalGeometryStorageType
Definition: alugrid/2d/intersection.hh:136
GridImp::template Codim< 1 >::Geometry Geometry
Definition: alugrid/2d/intersection.hh:459
static const ThisType & instance()
Definition: alugrid/2d/intersection.hh:72
int twistInOutside() const
Definition: intersection_imp.cc:192
void setOutside(HElementType *outside, int opposite)
Definition: alugrid/2d/intersection.hh:192
Different resources needed by all grid implementations.
ALUMemoryProvider< ThisType > StorageType
Definition: alugrid/2d/intersection.hh:347
ALU2dGridGeometry< dim-1, dim, GridImp > LocalGeometryImp
Definition: alugrid/2d/intersection.hh:358
ALU2dGridEntity< 0, dim, GridImp > EntityImp
Definition: alugrid/2d/intersection.hh:125
LocalGeometryImpl intersectionSelfLocal_
Definition: alugrid/2d/intersection.hh:304
void done(const EntityImp &en)
Definition: alugrid/2d/intersection.hh:291
EntityPointer outside() const
return EntityPointer to the Entity on the outside of this intersection.
Definition: intersection_imp.cc:157
const LocalGeometryImpl & localGeom(const int aluFace, const int twist, const int corners) const
Definition: alugrid/2d/intersection.hh:63
Definition: alugrid/2d/intersection.hh:446
bool neighbor() const
return true if intersection is with neighbor on this level
Definition: intersection_imp.cc:132
void setInside(HElementType *inside)
Definition: alugrid/2d/intersection.hh:186
Definition: alu2dinclude.hh:97
GeometryType type() const
obtain the type of reference element for this intersection
Definition: intersection_imp.cc:309
void assign(const ALU2dGridLeafIntersectionIterator< GridImp > &org)
The copy constructor.
Definition: intersection_imp.cc:437
NormalType integrationOuterNormal(const LocalCoordinate &local) const
Definition: intersection_imp.cc:228
HBndElType * boundary() const
Definition: alugrid/2d/intersection.hh:170
Definition: alugrid/2d/entity.hh:210
Definition: alugrid/2d/intersection.hh:117
Definition: alugrid/2d/intersection.hh:349
ALU2dGridIntersectionBase< GridImp > ThisType
Definition: alugrid/2d/intersection.hh:138
int twistInSelf() const
Definition: alugrid/2d/intersection.hh:261
ALU2dImplTraits< dimworld, eltype >::HElementType HElementType
Definition: alugrid/2d/intersection.hh:132
GridImp::template Codim< 1 >::Geometry Geometry
Definition: alugrid/2d/intersection.hh:354
GridImp::template Codim< 1 >::LocalGeometry LocalGeometry
Definition: alugrid/2d/intersection.hh:460
Definition: alugrid/2d/intersection.hh:447
ALU2dGridEntityPointer< 0, GridImp > EntityPointer
Definition: alugrid/2d/intersection.hh:465
EntityPointer inside() const
return EntityPointer to the Entity on the inside of this intersection.
Definition: intersection_imp.cc:140
int twistInInside() const
Definition: intersection_imp.cc:186
bool useOutside_
Definition: alugrid/2d/intersection.hh:207
ALUMemoryProvider< ThisType > StorageType
Definition: alugrid/2d/intersection.hh:444