MeshStructure.h 3.47 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
// ============================================================================
// ==                                                                        ==
// == AMDiS - Adaptive multidimensional simulations                          ==
// ==                                                                        ==
// ============================================================================
// ==                                                                        ==
// ==  crystal growth group                                                  ==
// ==                                                                        ==
// ==  Stiftung caesar                                                       ==
// ==  Ludwig-Erhard-Allee 2                                                 ==
// ==  53175 Bonn                                                            ==
// ==  germany                                                               ==
// ==                                                                        ==
// ============================================================================
// ==                                                                        ==
// ==  http://www.caesar.de/cg/AMDiS                                         ==
// ==                                                                        ==
// ============================================================================

/** \file MeshStructure.h */

#ifndef AMDIS_MESHSTRUCTURE_H
#define AMDIS_MESHSTRUCTURE_H

#include <vector>
#include "MemoryManager.h"

namespace AMDiS {

  class RefinementManager;
  class TraverseStack;
  class ElInfo;

  class MeshStructure
  {
  public:
    MEMORY_MANAGED(MeshStructure);

    MeshStructure()
      : currentIndex_(0),
	currentCode_(0),
	pos_(0),
	currentElement_(0),
	numElements_(0)
    {};
  
    void clear();

    void init(Mesh *mesh);

    void init(const std::vector<unsigned long int>& code,
	      int numElements) 
    {
      code_ = code;
      numElements_ = numElements;
      reset();
    };

    void reset();

    void insertElement(bool isLeaf);

    inline void commit() {
      if(pos_ > 0) {
	code_.push_back(currentCode_);
      }
      reset();
    };

    bool skipBranch(MeshStructure *insert = NULL);

    ElInfo *skipBranch(ElInfo *elInfo, TraverseStack *stack);

    bool nextElement(MeshStructure *insert = NULL);

    inline bool isLeafElement() {
      return (currentCode_ & 1) == 0;
    };

    void merge(MeshStructure *structure) {
      MeshStructure temp(*this);
      merge(&temp, structure, this);
    };

    static void merge(MeshStructure *structure1,
		      MeshStructure *structure2,
		      MeshStructure *result);

    void fitMeshToStructure(Mesh *mesh,
			    RefinementManager *manager,
			    bool checkPartition = false);

    void print() {
      FUNCNAME("MeshStructure::print()");
      reset();
      bool cont = true;
      while(cont) {
	if(isLeafElement()) {
	  MSG("0");
	} else {
	  MSG("1");
	}
	cont = nextElement();
      }
      MSG("\n");
    };

    inline const std::vector<unsigned long int>& getCode() {
      return code_;
    };

    inline int getNumElements() { 
      return numElements_; 
    };

    inline int getCurrentElement() {
      return currentElement_;
    };

  protected:
    std::vector<unsigned long int> code_;

    int currentIndex_;

    unsigned long int currentCode_;

    int pos_;

    int currentElement_;

    int numElements_;

    static const int unsignedLongSize_;
  };

}

#endif