4 #include <dune/common/stdstreams.hh> 19 template< ALU3dGr
idElementType elType,
class Comm >
33 , localIdSet_( *this )
34 , levelIndexVec_(MAXL,0) , leafIndexSet_(0)
35 , referenceElement_( elType ==
tetra 39 #ifdef USE_SMP_PARALLEL
44 , lockPostAdapt_( false )
49 , refinementType_( refinementType )
51 assert( elType ==
tetra || elType ==
hexa );
53 geomTypes_.resize( dimension+1 );
55 for(
int codim = 0; codim <= dimension; ++codim )
58 tmpType.makeSimplex( dimension - codim );
60 tmpType.makeCube( dimension - codim );
62 geomTypes_[ codim ].push_back( tmpType );
66 checkMacroGridFile( macroTriangFilename );
68 mygrid_ = createALUGrid( macroTriangFilename );
71 dverb <<
"************************************************" << std::endl;
72 dverb <<
"Created grid on p=" << comm().rank() << std::endl;
73 dverb <<
"************************************************" << std::endl;
81 template< ALU3dGr
idElementType elType,
class Comm >
87 return myGrid().indexManager(codim).getMaxIndex();
91 template< ALU3dGr
idElementType elType,
class Comm >
95 return myGrid().indexManager(codim).getMaxIndex();
99 template< ALU3dGr
idElementType elType,
class Comm >
106 template< ALU3dGr
idElementType elType,
class Comm >
116 template< ALU3dGr
idElementType elType,
class Comm >
117 template<
int cd, PartitionIteratorType pitype >
121 assert( level >= 0 );
123 if( level > maxlevel_ )
124 return this->
template lend<cd,pitype> (level);
130 template< ALU3dGr
idElementType elType,
class Comm >
131 template<
int cd, PartitionIteratorType pitype >
135 assert( level >= 0 );
141 template< ALU3dGr
idElementType elType,
class Comm >
146 return this->
template lbegin<cd,All_Partition>( level );
150 template< ALU3dGr
idElementType elType,
class Comm >
155 assert( level >= 0 );
156 return this->
template lend<cd,All_Partition>( level );
165 template< ALU3dGr
idElementType elType,
class Comm >
166 template<
int cd, PartitionIteratorType pitype >
170 assert( level >= 0 );
175 template< ALU3dGr
idElementType elType,
class Comm >
176 template<
int cd, PartitionIteratorType pitype >
180 return createLeafIteratorBegin<cd, pitype> (level) ;
184 template< ALU3dGr
idElementType elType,
class Comm >
189 return createLeafIteratorBegin<cd, All_Partition> (level) ;
193 template< ALU3dGr
idElementType elType,
class Comm >
194 template<
int cd, PartitionIteratorType pitype >
198 return createLeafIteratorBegin< cd, pitype > (maxlevel_) ;
202 template< ALU3dGr
idElementType elType,
class Comm >
207 return createLeafIteratorBegin< cd, All_Partition> (maxlevel_) ;
211 template< ALU3dGr
idElementType elType,
class Comm >
215 return createLeafIteratorBegin<0, All_Partition> (level) ;
219 template< ALU3dGr
idElementType elType,
class Comm >
223 return createLeafIteratorBegin<0, All_Partition> (maxlevel_) ;
232 template< ALU3dGr
idElementType elType,
class Comm >
233 template<
int cd, PartitionIteratorType pitype >
237 assert( level >= 0 );
242 template< ALU3dGr
idElementType elType,
class Comm >
243 template<
int cd, PartitionIteratorType pitype >
247 return createLeafIteratorEnd < cd, pitype> (level);
251 template< ALU3dGr
idElementType elType,
class Comm >
256 return createLeafIteratorEnd < cd, All_Partition> (level);
260 template< ALU3dGr
idElementType elType,
class Comm >
261 template<
int cd, PartitionIteratorType pitype >
265 return createLeafIteratorEnd < cd, pitype> (maxlevel_);
269 template< ALU3dGr
idElementType elType,
class Comm >
274 return createLeafIteratorEnd < cd, All_Partition> (maxlevel_);
278 template< ALU3dGr
idElementType elType,
class Comm >
282 return createLeafIteratorEnd <0, All_Partition> (level);
286 template< ALU3dGr
idElementType elType,
class Comm >
290 return createLeafIteratorEnd <0,All_Partition> (maxlevel_);
297 template< ALU3dGr
idElementType elType,
class Comm >
301 bool marked = (this->getRealImplementation( entity )).mark(ref);
304 if(ref > 0) ++refineMarked_;
305 if(ref < 0) ++coarsenMarked_;
312 template< ALU3dGr
idElementType elType,
class Comm >
316 return this->getRealImplementation( entity ).getMark();
321 template< ALU3dGr
idElementType elType,
class Comm >
322 template<
class Gr
idImp,
class DataHandle >
327 assert( (refCount + maxLevel()) < MAXL );
329 for(
int count = refCount; count > 0; --count )
341 template< ALU3dGr
idElementType elType,
class Comm >
342 template<
class Gr
idImp,
class DataHandle >
350 EntityObject father( EntityImp( factory(), this->maxLevel() ) );
351 EntityObject son ( EntityImp( factory(), this->maxLevel() ) );
354 bool mightCoarse = preAdapt();
356 bool refined = false ;
361 ALU3DSPACE AdaptRestrictProlongGlSet< MyType, AdaptDataHandle, GlobalIdSetImp >
363 father,this->getRealImplementation(father),
364 son, this->getRealImplementation(son),
368 refined = myGrid().duneAdapt(rp);
372 ALU3DSPACE AdaptRestrictProlongImpl< MyType, AdaptDataHandle >
374 father,this->getRealImplementation(father),
375 son, this->getRealImplementation(son),
378 refined = myGrid().duneAdapt(rp);
381 if(refined || mightCoarse)
398 template< ALU3dGr
idElementType elType,
class Comm >
401 template< ALU3dGr
idElementType elType >
408 template<
class DataHandle >
409 static bool loadBalance ( Grid &grid, DataHandle &data ) {
return false; }
411 template<
class DataHandle,
class DataType >
419 template<
class DataHandle,
class DataType >
427 #if ALU3DGRID_PARALLEL 428 template< ALU3dGr
idElementType elType >
434 static bool loadBalance ( Grid &grid )
436 if( grid.
comm().size() <= 1 )
439 const bool changed = grid.
myGrid().duneLoadBalance();
443 grid.
myGrid().duneExchangeDynamicState();
462 template<
class DataHandle >
463 static bool loadBalance ( Grid &grid, DataHandle &data )
465 if( grid.
comm().size() <= 1 )
469 typedef typename EntityObject::ImplementationType EntityImp;
474 typedef ALU3DSPACE LoadBalanceElementCount< Grid, DataHandle > LDBElCountType;
478 LDBElCountType elCount( grid,
479 father, Grid::getRealImplementation( father ),
480 son, Grid::getRealImplementation( son ),
483 ALU3DSPACE GatherScatterLoadBalance< Grid, DataHandle, LDBElCountType >
484 gs( grid, en, Grid::getRealImplementation( en ), data, elCount );
487 const bool changed = grid.
myGrid().duneLoadBalance( gs, elCount );
492 grid.
myGrid().duneExchangeDynamicState();
512 template<
class DataHandle,
class DataType >
513 static void communicate (
const Grid &grid,
521 typedef typename VertexObject::ImplementationType VertexImp;
523 typedef typename EdgeObject::ImplementationType EdgeImp;
525 typedef typename FaceObject::ImplementationType FaceImp;
527 typedef typename ElementObject::ImplementationType ElementImp;
529 if( grid.
comm().size() > 1 )
537 grid.template lbegin<0>( level ),
538 grid.template lend<0>( level ), level );
542 VertexObject vx( VertexImp( grid.
factory(), level ) );
543 ALU3DSPACE GatherScatterLevelData< Grid, DataHandleType, 3 >
544 vertexData( grid, vx, Grid::getRealImplementation( vx ), data, *levelISet, level );
546 EdgeObject edge( EdgeImp( grid.
factory(), level ) );
547 ALU3DSPACE GatherScatterLevelData< Grid, DataHandleType, 2 >
548 edgeData( grid, edge, Grid::getRealImplementation( edge ), data, *levelISet, level );
550 FaceObject face( FaceImp( grid.
factory(), level ) );
551 ALU3DSPACE GatherScatterLevelData< Grid, DataHandleType, 1 >
552 faceData( grid, face, Grid::getRealImplementation( face ), data, *levelISet, level );
554 ElementObject element( ElementImp( grid.
factory(), level ) );
555 ALU3DSPACE GatherScatterLevelData< Grid, DataHandleType, 0 >
556 elementData( grid, element, Grid::getRealImplementation( element ), data, *levelISet, level );
558 doCommunication( grid, vertexData, edgeData, faceData, elementData, iftype, dir );
565 template<
class DataHandle,
class DataType >
566 static void communicate (
const Grid &grid,
573 typedef typename VertexObject::ImplementationType VertexImp;
575 typedef typename EdgeObject::ImplementationType EdgeImp;
577 typedef typename FaceObject::ImplementationType FaceImp;
579 typedef typename ElementObject::ImplementationType ElementImp;
581 if( grid.
comm().size() > 1 )
584 ALU3DSPACE GatherScatterLeafData< Grid, DataHandleType, 3 >
585 vertexData( grid, vx, Grid::getRealImplementation( vx ), data );
588 ALU3DSPACE GatherScatterLeafData< Grid, DataHandleType, 2 >
589 edgeData( grid, edge, Grid::getRealImplementation( edge ), data );
592 ALU3DSPACE GatherScatterLeafData< Grid, DataHandleType, 1 >
593 faceData( grid, face, Grid::getRealImplementation( face ), data );
595 ElementObject element( ElementImp( grid.
factory(), grid.
maxLevel() ) );
596 ALU3DSPACE GatherScatterLeafData< Grid, DataHandleType, 0 >
597 elementData( grid, element, Grid::getRealImplementation( element ), data );
599 doCommunication( grid, vertexData, edgeData, faceData, elementData, iftype, dir );
604 doCommunication (
const Grid &grid,
605 GatherScatterType &vertexData, GatherScatterType &edgeData,
606 GatherScatterType &faceData, GatherScatterType &elementData,
612 dverb <<
"ALUGrid contains no overlap, therefore no communication for" << std::endl;
613 dverb <<
"Overlap_OverlapFront_Interface or Overlap_All_Interface interfaces!" << std::endl;
617 grid.
myGrid().borderBorderCommunication(vertexData,edgeData,faceData,elementData);
622 grid.
myGrid().interiorGhostCommunication(vertexData,edgeData,faceData,elementData);
625 grid.
myGrid().ghostInteriorCommunication(vertexData,edgeData,faceData,elementData);
629 grid.
myGrid().allAllCommunication(vertexData,edgeData,faceData,elementData);
631 DUNE_THROW(
GridError,
"Wrong set of parameters in ALUGridCommHelper::doCommunication" );
634 #endif // #if ALU3DGRID_PARALLEL 638 template< ALU3dGr
idElementType elType,
class Comm >
646 template< ALU3dGr
idElementType elType,
class Comm >
647 template<
class DataHandle >
655 template< ALU3dGr
idElementType elType,
class Comm >
656 template <
class DataHandleImp,
class DataType>
666 template< ALU3dGr
idElementType elType,
class Comm >
667 template <
class DataHandleImp,
class DataType>
677 template< ALU3dGr
idElementType elType,
class Comm >
681 return "ALUCubeGrid";
683 return "ALUSimplexGrid";
void postAdapt()
clear all entity new markers
int global_size(int cd) const
number of grid entities on all levels for given codim
Definition: grid_inline.hh:82
Definition: objectfactory.hh:23
Interface class for the Grid's adapt method where the parameter is a AdaptDataHandleInterface.
Definition: adaptcallback.hh:30
double alu3d_ctype
Definition: alu3dinclude.hh:59
bool checkMacroGrid(ALU3dGridElementType elType, const std::string filename)
static void communicate(const Grid &grid, const CommDataHandleIF< DataHandle, DataType > &data, const InterfaceType iftype, const CommunicationDirection dir, const int level)
Definition: grid_inline.hh:412
DefaultIndexSet creates an index set by using the grids persistent container an a given pair of itera...
Definition: defaultindexsets.hh:66
[ provides Dune::Grid ]
Definition: alugrid/3d/entity.hh:36
send/receive interior and border entities
Definition: gridenums.hh:85
ALU3dGrid< elType, No_Comm > Grid
Definition: grid_inline.hh:404
send overlap, receive all entities
Definition: gridenums.hh:88
Base class for exceptions in Dune grid modules.
Definition: exceptions.hh:16
Interface class for vertex projection at the boundary.
Definition: boundaryprojection.hh:23
send overlap, receive overlap and front entities
Definition: gridenums.hh:87
void communicate(CommDataHandleIF< DataHandleImp, DataTypeImp > &data, InterfaceType iftype, CommunicationDirection dir, int level) const
Communicate information on distributed entities on a given level Template parameter is a model of Dun...
GeometryType
Type representing VTK's entity geometry types.
Definition: common.hh:178
Definition: alugrid/common/declaration.hh:18
CommunicationDirection
Define a type for communication direction parameter.
Definition: gridenums.hh:168
GitterImplType & myGrid() const
Definition: grid_inline.hh:108
ALU3dGrid(const std::string ¯oTriangFilename, const MPICommunicatorType mpiComm, const DuneBoundaryProjectionType *bndPrj, const DuneBoundaryProjectionVector *bndVec, const ALUGridRefinementType refinementType)
Definition: grid_inline.hh:21
bool loadBalance()
Calculates load of each process and repartition the grid if neccessary. For parameters of the load ba...
Definition: grid_inline.hh:639
Traits::template Codim< cd >::template Partition< pitype >::LevelIterator lbegin(int level) const
Iterator to first entity of given codim on level.
Definition: grid_inline.hh:119
Definition: topology.hh:13
Traits::DuneBoundaryProjectionVector DuneBoundaryProjectionVector
boundary projection type
Definition: alugrid/3d/grid.hh:496
send interior and border, receive all entities
Definition: gridenums.hh:86
int getMark(const typename Traits::template Codim< 0 >::Entity &e) const
returns adaptation mark for given entity
Definition: grid_inline.hh:314
Definition: alugrid/common/declaration.hh:18
const CollectiveCommunication & comm() const
return const reference to a collective communication object. The return type is a model of Dune::Coll...
Definition: alugrid/3d/grid.hh:785
static std::string name()
for grid identification
Definition: grid_inline.hh:678
std::vector< LevelIndexSetImp * > levelIndexVec_
Definition: alugrid/3d/grid.hh:1063
Traits::template Codim< cd >::template Partition< pitype >::LevelIterator lend(int level) const
one past the end on this level
Definition: grid_inline.hh:133
reverse communication direction
Definition: gridenums.hh:170
Definition: adaptcallback.hh:19
CommDataHandleIF describes the features of a data handle for communication in parallel runs using the...
Definition: datahandleif.hh:72
communicate as given in InterfaceType
Definition: gridenums.hh:169
int hierSetSize(int cd) const
Definition: grid_inline.hh:92
send all and receive all entities
Definition: gridenums.hh:89
ALU3dImplTraits< elType, Comm >::GitterImplType GitterImplType
Definition: alugrid/3d/grid.hh:536
Definition: topology.hh:13
Leaf iterator.
Definition: alugrid/3d/entity.hh:34
GatherScatter GatherScatterType
Definition: alu3dinclude.hh:49
static void communicate(const Grid &grid, const CommDataHandleIF< DataHandle, DataType > &data, const InterfaceType iftype, const CommunicationDirection dir)
Definition: grid_inline.hh:420
ALUGrid boundary projection implementation DuneBndProjection has to fulfil the DuneBoundaryProjection...
Definition: common/bndprojection.hh:11
Types for GridView.
Definition: common/grid.hh:420
InterfaceType
Parameter to be used for the communication functions.
Definition: gridenums.hh:84
Include standard header files.
Definition: agrid.hh:59
#define ALU3DSPACE
Definition: alu3dinclude.hh:26
InterfaceType::Implementation ImplementationType
Definition: common/grid.hh:1468
Comm MPICommunicatorType
Definition: alugrid/3d/grid.hh:559
Traits::template Codim< 0 >::LeafIterator LeafIteratorType
Definition: alugrid/3d/grid.hh:531
static bool loadBalance(Grid &grid, DataHandle &data)
Definition: grid_inline.hh:409
bool adapt()
Refine all positive marked leaf entities, coarsen all negative marked entities if possible...
Definition: alugrid/3d/grid.hh:137
static bool loadBalance(Grid &grid)
Definition: grid_inline.hh:406
Definition: alugrid/3d/grid.hh:88
int maxLevel() const
Return maximum level defined in this grid. Levels are numbered maxLevel with 0 the coarsest level...
Definition: grid_inline.hh:100
Definition: alugrid/3d/entity.hh:24
ALUGridRefinementType
available refinement types for ALUGrid
Definition: alugrid/common/declaration.hh:20
const GridObjectFactoryType & factory() const
Definition: alugrid/3d/grid.hh:1007
GlobalIdSetImp * globalIdSet_
Definition: alugrid/3d/grid.hh:1057
bool mark(int refCount, const typename Traits::template Codim< 0 >::Entity &e)
Marks an entity to be refined/coarsened in a subsequent adapt.
Definition: grid_inline.hh:299
A Traits struct that collects all associated types of one implementation.
Definition: common/grid.hh:437
void globalRefine(int refCount)
uses the interface, mark on entity and refineLocal