Dune Core Modules (2.9.0)

alu3dinclude.hh
1 #ifndef DUNE_ALUGRID_ALU3DINCLUDE_HH
2 #define DUNE_ALUGRID_ALU3DINCLUDE_HH
3 
4 #include "aluinline.hh"
5 
6 // all methods and classes of the ALUGrid are defined in the namespace
7 #define ALU3DSPACE ::ALUGrid::
8 
9 #include <dune/alugrid/common/declaration.hh>
10 
11 #include <dune/alugrid/impl/serial/gatherscatter.hh>
12 #include <dune/alugrid/impl/serial/key.h>
13 #include <dune/alugrid/impl/serial/myalloc.h>
14 #include <dune/alugrid/impl/serial/serialize.h>
15 
16 #include <dune/alugrid/impl/parallel/mpAccess.h>
17 #include <dune/alugrid/impl/parallel/gitter_pll_ldb.h>
18 
19 #include <dune/alugrid/impl/serial/gitter_sti.h>
20 
21 #include <dune/alugrid/impl/serial/gitter_hexa_top.h>
22 #include <dune/alugrid/impl/serial/mapp_tetra_3d_ext.h>
23 #include <dune/alugrid/impl/serial/gitter_tetra_top.h>
24 #include <dune/alugrid/impl/serial/walk.h>
25 #include <dune/alugrid/impl/serial/gitter_impl.h>
26 #include <dune/alugrid/impl/serial/gitter_mgb.h>
27 #include <dune/alugrid/impl/serial/key.h>
28 #include <dune/alugrid/impl/serial/lock.h>
29 
30 #include <dune/alugrid/impl/duneinterface/gitter_dune_impl.h>
31 
32 namespace ALUGrid
33 {
34 
35  typedef Gitter::AdaptRestrictProlong AdaptRestrictProlongType;
36 
37  static const int ProcessorBoundary_t = Gitter::hbndseg_STI::closure;
38 
39  // general GatherScatter type
40  typedef GatherScatter GatherScatterType;
41 
42 } // namespace ALUGrid
43 
44 
45 // headers for parallel grid structures
46 #include <dune/alugrid/impl/parallel/gitter_pll_sti.h>
47 #include <dune/alugrid/impl/parallel/gitter_pll_impl.h>
48 #include <dune/alugrid/impl/parallel/gitter_pll_ldb.h>
49 #include <dune/alugrid/impl/parallel/gitter_tetra_top_pll.h>
50 #include <dune/alugrid/impl/parallel/gitter_hexa_top_pll.h>
51 #include <dune/alugrid/impl/parallel/gitter_pll_mgb.h>
52 #include <dune/alugrid/impl/duneinterface/gitter_dune_pll_impl.h>
53 
54 #if ALU3DGRID_PARALLEL
55 // if MPI was found include MPI communications
56 #include <dune/alugrid/impl/parallel/mpAccess_MPI.h>
57 #endif // #if ALU3DGRID_PARALLEL
58 
59 //- local includes
60 #include <dune/alugrid/3d/topology.hh>
61 
62 namespace Dune
63 {
64 
65  // typedef of ALU3dGridElementType see topology.hh
66 
67  // i.e. double or float
68  typedef double alu3d_ctype;
69 
70 
71  // ALU3dBasicImplTraits
72  // --------------------
73 
74  template< class Comm >
75  struct ALU3dBasicImplTraits;
76 
77  template<>
78  struct ALU3dBasicImplTraits< ALUGridNoComm >
79  {
80  typedef ALU3DSPACE Gitter GitterType;
81  typedef ALU3DSPACE GitterDuneImpl GitterImplType;
82 
83  typedef GitterType::helement_STI HElementType; // Interface Element
84  typedef GitterType::hface_STI HFaceType; // Interface Face
85  typedef GitterType::hedge_STI HEdgeType; // Interface Edge
86  typedef GitterType::vertex_STI VertexType; // Interface Vertex
87  typedef GitterType::hbndseg_STI HBndSegType;
88  typedef GitterType::ghostpair_STI GhostPairType;
89 
90  typedef HElementType PllElementType;
91 
92  typedef GitterType::Geometric::hedge1_GEO GEOEdgeType;
93 
95  template <class BndFaceType>
96  static bool isGhost( const BndFaceType* ghost )
97  {
98  return false ;
99  }
100  };
101 
102  template<>
103  struct ALU3dBasicImplTraits< ALUGridMPIComm >
104  {
105  typedef ALU3DSPACE GitterDunePll GitterType;
106  typedef ALU3DSPACE GitterDunePll GitterImplType;
107 
108  typedef GitterType::helement_STI HElementType; // Interface Element
109  typedef GitterType::hface_STI HFaceType; // Interface Face
110  typedef GitterType::hedge_STI HEdgeType; // Interface Edge
111  typedef GitterType::vertex_STI VertexType; // Interface Vertex
112  typedef GitterType::hbndseg_STI HBndSegType;
113  typedef GitterType::ghostpair_STI GhostPairType;
114 
115  typedef ALU3DSPACE ElementPllXIF_t PllElementType;
116 
117  typedef GitterType::Geometric::hedge1_GEO GEOEdgeType;
118 
119  // method for ghost check
120  template <class BndFaceType>
121  static bool isGhost( const BndFaceType* ghost )
122  {
123  return ( ghost != 0 );
124  }
125  };
126 
127 
128  // ALU3dCodimImplTraits
129  // --------------------
130 
131  template< ALU3dGridElementType elType, class Comm, int dim, int codim >
132  struct ALU3dCodimImplTraits;
133 
134  template< class Comm, int dim>
135  struct ALU3dCodimImplTraits< tetra, Comm, dim, 0 >
136  {
137  typedef typename ALU3dBasicImplTraits< Comm >::GitterType GitterType;
138  typedef typename ALU3dBasicImplTraits< Comm >::GitterImplType GitterImplType;
139 
140  typedef typename GitterType::helement_STI InterfaceType;
141  typedef typename GitterType::Geometric::hasFace3 EntitySeedType;
142  typedef typename GitterImplType::Objects::tetra_IMPL ImplementationType;
143  typedef typename GitterType::hbndseg_STI GhostInterfaceType;
144  typedef typename GitterImplType::Objects::Hbnd3Default GhostImplementationType;
145  };
146 
147  template< class Comm, int dim >
148  struct ALU3dCodimImplTraits< hexa, Comm, dim, 0 >
149  {
150  typedef typename ALU3dBasicImplTraits< Comm >::GitterType GitterType;
151  typedef typename ALU3dBasicImplTraits< Comm >::GitterImplType GitterImplType;
152 
153  typedef typename GitterType::helement_STI InterfaceType;
154  typedef typename GitterType::Geometric::hasFace4 EntitySeedType;
155  typedef typename GitterImplType::Objects::hexa_IMPL ImplementationType;
156  typedef typename GitterType::hbndseg_STI GhostInterfaceType;
157  typedef typename GitterImplType::Objects::Hbnd4Default GhostImplementationType;
158  };
159 
160  template< class Comm, int dim >
161  struct ALU3dCodimImplTraits< tetra, Comm, dim, 1 >
162  {
163  typedef typename ALU3dBasicImplTraits< Comm >::GitterType GitterType;
164 
165  typedef typename GitterType::hface_STI InterfaceType;
166  typedef InterfaceType EntitySeedType;
167  typedef typename GitterType::Geometric::hface3_GEO ImplementationType;
168  };
169 
170  template< class Comm, int dim >
171  struct ALU3dCodimImplTraits< hexa, Comm, dim, 1 >
172  {
173  typedef typename ALU3dBasicImplTraits< Comm >::GitterType GitterType;
174 
175  typedef typename GitterType::hface_STI InterfaceType;
176  typedef InterfaceType EntitySeedType;
177  typedef typename GitterType::Geometric::hface4_GEO ImplementationType;
178  };
179 
180  template< ALU3dGridElementType elType, class Comm >
181  struct ALU3dCodimImplTraits< elType, Comm, 3, 2 >
182  {
183  typedef typename ALU3dBasicImplTraits< Comm >::GitterType GitterType;
184 
185  typedef typename GitterType::hedge_STI InterfaceType;
186  typedef InterfaceType EntitySeedType;
187  typedef typename GitterType::Geometric::hedge1_GEO ImplementationType;
188  };
189 
190 
191  template< ALU3dGridElementType elType, class Comm >
192  struct ALU3dCodimImplTraits< elType, Comm, 2, 2 >
193  {
194  typedef typename ALU3dBasicImplTraits< Comm >::GitterType GitterType;
195 
196  typedef typename GitterType::vertex_STI InterfaceType;
197  typedef InterfaceType EntitySeedType;
198  typedef typename GitterType::Geometric::VertexGeo ImplementationType;
199  };
200 
201  template< ALU3dGridElementType elType, class Comm >
202  struct ALU3dCodimImplTraits< elType, Comm, 3, 3 >
203  {
204  typedef typename ALU3dBasicImplTraits< Comm >::GitterType GitterType;
205 
206  typedef typename GitterType::vertex_STI InterfaceType;
207  typedef InterfaceType EntitySeedType;
208  typedef typename GitterType::Geometric::VertexGeo ImplementationType;
209  };
210 
211  // Refinement rules in general
212  template< class MarkRuleType, ALU3dGridElementType elType >
213  struct ALU3dRefinementTraits {};
214 
215  // Refinement rules for simplices
216  template< class MarkRuleType >
217  struct ALU3dRefinementTraits < MarkRuleType, tetra >
218  {
219  // refinement and coarsening enum
220  enum { bisect_element_t = MarkRuleType::bisect };
221  enum { refine_element_t = MarkRuleType::regular };
222  enum { coarse_element_t = MarkRuleType::crs };
223  enum { nosplit_element_t = MarkRuleType::nosplit };
224  };
225 
226  // Refinement rules for cubes
227  template< class MarkRuleType >
228  struct ALU3dRefinementTraits < MarkRuleType, hexa >
229  {
230  // refinement and coarsening enum
231  enum { bisect_element_t = MarkRuleType::regular };
232  enum { refine_element_t = MarkRuleType::regular };
233  enum { coarse_element_t = MarkRuleType::crs };
234  enum { nosplit_element_t = MarkRuleType::nosplit };
235  };
236 
237 
238  // ALU3dImplTraits
239  // ---------------
240 
241  template< ALU3dGridElementType elType, class Comm >
242  struct ALU3dImplTraits;
243 
244  template< class Comm >
245  struct ALU3dImplTraits< tetra, Comm >
246  : public ALU3dBasicImplTraits< Comm >
247  {
248  typedef typename ALU3dBasicImplTraits< Comm >::GitterType GitterType;
249  typedef typename ALU3dBasicImplTraits< Comm >::GitterImplType GitterImplType;
250 
251  typedef typename GitterType::Geometric::hface3_GEO GEOFaceType;
252  typedef typename GitterType::Geometric::VertexGeo GEOVertexType;
253  typedef typename GitterImplType::Objects::tetra_IMPL IMPLElementType;
254  typedef typename GitterType::Geometric::tetra_GEO GEOElementType;
255  typedef typename GitterType::Geometric::periodic3_GEO GEOPeriodicType;
256  typedef typename GitterType::Geometric::hasFace3 HasFaceType;
257  typedef typename GitterType::Geometric::Hface3Rule HfaceRuleType;
258  typedef typename GitterImplType::Objects::Hbnd3Default BNDFaceType;
259  typedef typename GitterImplType::Objects::hbndseg3_IMPL ImplBndFaceType;
260 
261  typedef typename GitterType::Geometric::TetraRule MarkRuleType;
262 
263  struct RefinementRules
264  : public ALU3dRefinementTraits<MarkRuleType, tetra>
265  {};
266 
267  typedef std::pair< GEOFaceType *, int > NeighbourFaceType;
268  typedef std::pair< HasFaceType *, int > NeighbourPairType;
269 
270  template< int dim, int codim >
271  struct Codim
272  : public ALU3dCodimImplTraits< tetra, Comm, dim, codim >
273  {};
274 
275  // access of faces
276  template <class Elem>
277  static const GEOFaceType* getFace( const Elem& elem, const int aluFace )
278  {
279  return elem.myhface( aluFace );
280  }
281  };
282 
283  template< class Comm >
284  struct ALU3dImplTraits< hexa, Comm >
285  : public ALU3dBasicImplTraits< Comm >
286  {
287  typedef typename ALU3dBasicImplTraits< Comm >::GitterType GitterType;
288  typedef typename ALU3dBasicImplTraits< Comm >::GitterImplType GitterImplType;
289 
290  typedef typename GitterType::Geometric::hface4_GEO GEOFaceType;
291  typedef typename GitterType::Geometric::VertexGeo GEOVertexType;
292  typedef typename GitterImplType::Objects::hexa_IMPL IMPLElementType;
293  typedef typename GitterType::Geometric::hexa_GEO GEOElementType;
294  typedef typename GitterType::Geometric::periodic4_GEO GEOPeriodicType;
295  typedef typename GitterType::Geometric::hasFace4 HasFaceType;
296  typedef typename GitterType::Geometric::Hface4Rule HfaceRuleType;
297  typedef typename GitterImplType::Objects::Hbnd4Default BNDFaceType;
298  typedef typename GitterImplType::Objects::hbndseg4_IMPL ImplBndFaceType;
299 
300  typedef typename GitterType::Geometric::HexaRule MarkRuleType;
301 
302  struct RefinementRules
303  : public ALU3dRefinementTraits<MarkRuleType, hexa>
304  {};
305 
306  typedef std::pair< GEOFaceType *, int > NeighbourFaceType;
307  typedef std::pair< HasFaceType *, int > NeighbourPairType;
308 
309  template< int dim, int codim >
310  struct Codim
311  : public ALU3dCodimImplTraits< hexa, Comm, dim, codim >
312  {};
313 
314  // access of faces
315  template <class Elem>
316  static const GEOFaceType* getFace( const Elem& elem, const int aluFace )
317  {
318  return elem.myhface( aluFace );
319  }
320  };
321 
322 
323 
326  template< class Comm >
327  struct ALU3dGridVertexList
328  {
329  // level vertex iterator list
330  typedef typename ALU3dBasicImplTraits< Comm >::VertexType VertexType;
331  typedef std::vector< VertexType * > VertexListType;
332  typedef typename VertexListType::iterator IteratorType;
333 
334  ALU3dGridVertexList ()
335  : up2Date_( false )
336  {}
337 
338  size_t size () const { return vertexList_.size(); }
339 
340  bool up2Date () const { return up2Date_; }
341  void unsetUp2Date () { up2Date_ = false; }
342 
343  // make grid walkthrough and calc global size
344  template <class GridType>
345  void setupVxList (const GridType & grid, int level);
346 
347  IteratorType begin () { return vertexList_.begin(); }
348  IteratorType end () { return vertexList_.end(); }
349 
350  VertexListType & getItemList() { return vertexList_; }
351 
352  private:
353  bool up2Date_;
354  //careful: due to the setupVxList structure the ordering of vertexList_ and validateList_ differ in the level Case
355  //for validateList_ we want the ALUGrid Index as ordering, as we want to use it for faces
356  VertexListType vertexList_;
357  };
358 
359 
362  template< class Comm >
363  struct ALU3dGridLeafVertexList
364  {
365  // level vertex iterator list
366  typedef typename ALU3dBasicImplTraits< Comm >::VertexType VertexType;
367  typedef std::pair< VertexType *, int > ItemType;
368  typedef std::vector< ItemType > VertexListType;
369  typedef typename VertexListType::iterator IteratorType;
370 
371  ALU3dGridLeafVertexList ()
372  : up2Date_( false )
373  {}
374 
375  size_t size () const { return vertexList_.size(); }
376 
377  bool up2Date () const { return up2Date_; }
378  void unsetUp2Date () { up2Date_ = false; }
379 
380  // make grid walkthrough and calc global size
381  template <class GridType>
382  void setupVxList (const GridType & grid);
383 
384  IteratorType begin () { return vertexList_.begin(); }
385  IteratorType end () { return vertexList_.end(); }
386 
387  VertexListType & getItemList() { return vertexList_; }
388 
389  int getLevel ( const VertexType &vertex ) const
390  {
391  const int idx = vertex.getIndex();
392  alugrid_assert ( idx >= 0 );
393  alugrid_assert ( idx < (int)size());
394  const ItemType & p = vertexList_[idx];
395  if( p.first == 0 )
396  return vertex.level();
397  else
398  return p.second;
399  }
400 
401  private:
402  bool up2Date_;
403  VertexListType vertexList_;
404  };
405 
406 
407 
408  class ALU3dGridItemList
409  {
410  public:
411  // level vertex iterator list
412  typedef std::vector < void * > ItemListType;
413  typedef ItemListType :: iterator IteratorType;
414 
415  ALU3dGridItemList () : up2Date_(false) {}
416 
417  size_t size () const { return itemList_.size(); }
418 
419  bool up2Date () const { return up2Date_; }
420  void unsetUp2Date () { up2Date_ = false; }
421 
422  void markAsUp2Date() { up2Date_ = true; }
423 
424  IteratorType begin () { return itemList_.begin(); }
425  IteratorType end () { return itemList_.end(); }
426 
427  ItemListType & getItemList() { return itemList_; }
428 
429  private:
430  bool up2Date_;
431  ItemListType itemList_;
432  };
433 
434  typedef ALU3dGridItemList ALU3dGridItemListType;
435 
437  // some helper functions
439 
440  template< class Comm >
441  struct ALU3dGridFaceGetter
442  {
443  static const typename ALU3dImplTraits< tetra, Comm >::GEOFaceType *
444  getFace( const typename ALU3dImplTraits< tetra, Comm >::GEOElementType& elem, int index)
445  {
446  alugrid_assert (index >= 0 && index < 4);
447  return elem.myhface( ElementTopologyMapping< tetra >::dune2aluFace(index) );
448  }
449 
450  static const typename ALU3dImplTraits< hexa, Comm >::GEOFaceType*
451  getFace( const typename ALU3dImplTraits< hexa, Comm >::GEOElementType &elem, int index )
452  {
453  alugrid_assert (index >= 0 && index < 6);
454  return elem.myhface( ElementTopologyMapping< hexa >::dune2aluFace(index) );
455  }
456  };
457 
458 } // end namespace Dune
459 
460 #endif // #ifndef DUNE_ALUGRID_ALU3DINCLUDE_HH
Dune::ElementTopologyMapping::dune2aluFace
static int dune2aluFace(int index)
Maps face index from Dune onto ALU3dGrid reference element.
Definition: topology.cc:387
Dune::InterfaceType
InterfaceType
Parameter to be used for the communication functions.
Definition: gridenums.hh:86
Dune::GeometryTypes::vertex
constexpr GeometryType vertex
GeometryType representing a vertex.
Definition: type.hh:506
Dune
Dune namespace.
Definition: alignedallocator.hh:13
Dune::ALU3dGridLeafVertexList
Definition: alu3dinclude.hh:364
Dune::ALU3dGridVertexList
Definition: alu3dinclude.hh:328
Creative Commons License   |  Legal Statements / Impressum  |  Hosted by TU Dresden  |  generated with Hugo v0.80.0 (May 5, 22:29, 2024)