-
Notifications
You must be signed in to change notification settings - Fork 0
/
MeshRegion.h
92 lines (88 loc) · 3.37 KB
/
MeshRegion.h
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
#ifndef MESHREGION_H
#define MESHREGION_H
#include <map>
#include <set>
#include <string>
#include <vector>
#define NFACEMAX 6
// only topology structures
// no physics
class MeshRegion {
public:
MeshRegion(std::string name, double tolerance = 1.E-6);
int PointIsExist(std::vector<double> p, int &pId);
void RebuildEdgesIndex();
void RebuildFacesIndex();
void ExtractBndPts();
void OutCompAsGeo3D(std::string name, std::vector<std::vector<double>> center,
std::vector<double> radius);
void AddMeshRegion(MeshRegion &doc);
void GetCellCenter(int i, std::vector<double> &c);
void GetPts(int index, std::vector<int> &pts, char &type);
void GetEdgePts(int index, std::vector<int> &pts, char &type);
void GetFacePts(int index, std::vector<int> &pts, char &type);
void GetFaceEdges(int index, std::vector<int> &edges, char &type);
void GetFaceUnitNorm(int index, std::vector<double> &norm, char &type);
void GetCellPts(int index, std::vector<int> &pts, char &type);
void GetCellEdges(int index, std::vector<int> &edges, char &type);
void GetElementPts(int index, std::vector<int> &pts, char &type);
void GetBndElementPts(int index, std::vector<int> &pts, char &type);
void ReorgDomain(std::vector<void *> condition, bool detectSingular);
void OutPutSU2(std::string name);
void GetElements(std::map<int, std::vector<int>> &elements);
bool ReorgBoundary(double angle = 75. / 180. * 3.1415926);
void CheckMesh(double angle = 75. / 180. * 3.1415926);
double m_tolerance;
std::string m_name;
int m_dim;
std::map<int, std::vector<double>> m_pts;
std::map<int, std::vector<int>> m_edges;
std::map<int, std::vector<int>> m_faces;
std::map<int, std::vector<int>> m_cells;
std::map<int, char> m_cellsType;
int m_ptsIndexMax;
int m_edgeIndexMax;
int m_faceIndexMax;
int m_cellIndexMax;
std::map<int, std::set<int>> m_bndComposite;
std::map<int, std::set<int>> m_domain;
std::map<int, char> m_domainType;
std::set<int> m_singularElements;
std::set<int> m_singularEdge;
std::set<int> m_singularPts;
std::set<int> m_singularFace;
enum ElementType {
eHexahedron,
ePrism,
eTetrahedron,
ePyramid,
eQuadrilateral,
eTriangle,
eSegment,
ePoint
};
std::vector<char> ElementTag = {'H', 'R', 'A', 'P', 'Q', 'T', 'S', 'V'};
std::vector<int> ElementSU2 = {12, 13, 10, 14, 9, 5, 3, 0};
std::map<char, int> ElementTypeMap;
protected:
void MergePts(MeshRegion &doc, std::map<int, int> &ptsMap);
void MergeEdge(MeshRegion &doc, std::map<int, int> &ptsMap,
std::map<int, int> &edgeMap);
void MergeFace(MeshRegion &doc, std::map<int, int> &edgeMap,
std::map<int, int> &faceMap);
void AddCell(MeshRegion &doc, std::map<int, int> &faceMap,
std::map<int, int> &cellMap);
void MergeComposite(MeshRegion &doc, std::map<int, int> &faceMap,
std::map<int, int> &cellMap);
int PointHash(std::vector<double> p);
void InsertBndPts(int index);
int FindSingularElements(double angle = 75. / 180. * 3.1415926);
int FindSharedSingluarEdges(std::set<int> &bnd1, std::set<int> &bnd2,
int &shared, int &singular);
std::vector<double> m_minRange;
std::vector<double> m_maxRange;
std::map<std::set<int>, int> m_edgesIndex;
std::map<std::set<int>, int> m_facesIndex;
std::map<int, std::set<int>> m_bndPts;
};
#endif // MESHREGION_H