InteriorBoundary.h 6.76 KB
Newer Older
1
2
3
4
5
6
// ============================================================================
// ==                                                                        ==
// == AMDiS - Adaptive multidimensional simulations                          ==
// ==                                                                        ==
// ============================================================================
// ==                                                                        ==
7
// ==  TU Dresden                                                            ==
8
// ==                                                                        ==
9
10
11
// ==  Institut für Wissenschaftliches Rechnen                               ==
// ==  Zellescher Weg 12-14                                                  ==
// ==  01069 Dresden                                                         ==
12
13
14
15
// ==  germany                                                               ==
// ==                                                                        ==
// ============================================================================
// ==                                                                        ==
16
// ==  https://gforge.zih.tu-dresden.de/projects/amdis/                      ==
17
18
19
20
21
22
23
24
25
// ==                                                                        ==
// ============================================================================

/** \file InteriorBoundary.h */

#ifndef AMDIS_INTERIORBOUNDARY_H
#define AMDIS_INTERIORBOUNDARY_H

#include <vector>
Thomas Witkowski's avatar
Thomas Witkowski committed
26
#include <map>
27

28
#include "AMDiS_fwd.h"
29
30
#include "MacroElement.h"
#include "Element.h"
31
#include "Boundary.h"
32
33
34

namespace AMDiS {

35
36
  typedef std::vector<std::pair<GeoIndex, int> > ExcludeList;

37
38
  /// Defines the geometrical objects that forms the boundary;
  struct BoundaryObject {
39

40
41
42
    BoundaryObject()      
      : elType(0),
	reverseMode(false),
43
	excludedSubstructures(0)
44
45
46
47
48
49
50
51
52
    {}

    BoundaryObject(Element *e, int eType, GeoIndex sObj, int iObj)
      : el(e),
	elIndex(e->getIndex()),
	elType(eType),
	subObj(sObj),
	ithObj(iObj),
	reverseMode(false),
53
	excludedSubstructures(0)
54
    {}
55

56
57
58
59
60
61
    void setReverseMode(BoundaryObject &otherBound, FiniteElemSpace *feSpace);

    bool operator==(const BoundaryObject& other) const;
    
    bool operator!=(const BoundaryObject& other) const;

62
    /// The macro element to which the boundary element corresponds to.
Thomas Witkowski's avatar
Thomas Witkowski committed
63
    Element* el;
64

65
66
67
    /// Index of the macro element.
    int elIndex;

68
69
70
    /// Element type index, only used in 3d.
    int elType;

71
72
73
74
    /** \brief
     * Defines the geometrical object at the boundary. It must be "a part" of the
     * macro element \ref el, i.e., either 1 (a vertex), 2 (an edge) or 3 (a face).
     */
75
    GeoIndex subObj;
76
77
78
79
80

    /** \brief
     * Defines which of vertix, edge or face of the macro element is part of the
     * boundary.
     *
81
82
83
     * Example: If the macro element is a triangle, than \ref subObj may be either
     * 1 (vertex) or 2 (edge). Assume its the last one. So this variable defines
     * which of the three possible edges of the triangle is at the interior
84
85
     * boundary.
     */
86
    int ithObj;
87
88
89

    bool reverseMode;

90
91
92
93
94
95
96
97
98
    /** \brief
     * In many situations it may be necessary to exclude some parts of the element to
     * be part of the boundary. In 3d, when a face is part of the boundary, an edge or
     * an vertex may be exludeded. In 2d only vertices may be exluded to be part of
     * an edge boundary. This list contains pairs of exludeded structures. The first
     * component of every pair denotes if it is a vertex or an edge, and the second
     * component denotes the local index of the structure.
     */
    ExcludeList excludedSubstructures;
99
100
101
102
103
104
105
  };

  /** \brief 
   * Defines one atomic part of the boundary, i.e., two boundary objects where the
   * boundary goes through.
   */
  struct AtomicBoundary {
106
107
108
109
110
111
    AtomicBoundary()
      : type(INTERIOR)
    {}

    bool operator==(const AtomicBoundary& other) const;

112
    /// The rank's part of the boundary.
113
    BoundaryObject rankObj;
114
115

    /// The object on the other side of the boundary.
116
    BoundaryObject neighObj;
117
118
119
120
121

    /// Integer flag that is used to distinguish between different types of 
    /// boundaries. Till now it is used only for periodic boundaries, which are also
    /// handles as interior boundaries.
    BoundaryType type;
122
123
124
125
126
127
128
  };

  /** \brief
   * Defines the interior boundary, i.e. a bound within the domain. It is used for
   * the classical domain decomposition parallelization.
   */
  class InteriorBoundary {
Thomas Witkowski's avatar
Thomas Witkowski committed
129
130
131
  public:
    typedef std::map<int, std::vector<AtomicBoundary> > RankToBoundMap;

Thomas Witkowski's avatar
Thomas Witkowski committed
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
    /// Iterator for the interior boundary object.
    class iterator {      
    public:
      iterator(InteriorBoundary &b)
	: bound(b)
      {
	reset();
      }

      /// Set the iterator to the first position.
      void reset()
      {
	mapIt = bound.boundary.begin();
	nextNonempty();

	if (mapIt != bound.boundary.end())
	  vecIt = mapIt->second.begin();	
      }

      /// Test if iterator is at the final position.
      bool end() const
      {
	return (mapIt == bound.boundary.end());
      }

      /// Move iterator to the next position.
      void operator++()
      {
	++vecIt;
	if (vecIt == mapIt->second.end()) {
	  ++mapIt;
	  nextNonempty();

	  if (mapIt != bound.boundary.end())
	    vecIt = mapIt->second.begin();	  
	}	
      }

170
171
172
173
174
175
176
177
178
179
      inline AtomicBoundary& operator*()
      {
	return *vecIt;
      }

      inline AtomicBoundary* operator->()
      {
	return &(*vecIt);
      }

Thomas Witkowski's avatar
Thomas Witkowski committed
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
      void nextRank()
      {
	++mapIt;
	nextNonempty();

	if (mapIt != bound.boundary.end())
	  vecIt = mapIt->second.begin();	
      }

      int getRank() 
      {
	return mapIt->first;
      }

    protected:

      inline void nextNonempty()
      {
Thomas Witkowski's avatar
Thomas Witkowski committed
198
199
200
	if (mapIt == bound.boundary.end())
	  return;

Thomas Witkowski's avatar
Thomas Witkowski committed
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
 	while (mapIt->second.size() == 0) {
 	  ++mapIt;
	  if (mapIt == bound.boundary.end())
	    return;
	}
      }

    protected:
      RankToBoundMap::iterator mapIt;

      std::vector<AtomicBoundary>::iterator vecIt;

      InteriorBoundary &bound;
    };

216
  public:
Thomas Witkowski's avatar
Thomas Witkowski committed
217
218
    InteriorBoundary() {}

Thomas Witkowski's avatar
Thomas Witkowski committed
219
220
221
222
223
    void clear()
    {
      boundary.clear();
    }

224
    AtomicBoundary& getNewAtomic(int rank);
Thomas Witkowski's avatar
Thomas Witkowski committed
225

226
    /// Writes this object to a file.
227
228
    void serialize(std::ostream &out);

229
230
    /// Reads the state of an interior boundary from a file.
    void deserialize(std::istream &in, std::map<int, Element*> &elIndexMap);
231

232
233
234
235
    /// Compares this interior boundaries with some other. The order of the
    /// boundary elements within the object does not play a role.
    bool operator==(const InteriorBoundary& other) const;

236
237
238
239
240
  protected:
    void serializeExcludeList(std::ostream &out, ExcludeList &list);

    void deserializeExcludeList(std::istream &in, ExcludeList &list);

Thomas Witkowski's avatar
Thomas Witkowski committed
241
  public:
Thomas Witkowski's avatar
Thomas Witkowski committed
242
    RankToBoundMap boundary;
243
244
245
246
  };
}

#endif // AMDIS_INTERIORBOUNDARY_H