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
|
#ifndef FILE_MESHING3
#define FILE_MESHING3
#include "meshclass.hpp"
#include "adfront3.hpp"
#include "ruler3.hpp"
namespace netgen
{
enum MESHING3_RESULT
{
MESHING3_OK = 0,
MESHING3_GIVEUP = 1,
MESHING3_NEGVOL = 2,
MESHING3_OUTERSTEPSEXCEEDED = 3,
MESHING3_TERMINATE = 4,
MESHING3_BADSURFACEMESH = 5
};
/// 3d volume mesh generation
class Meshing3
{
/// current state of front
unique_ptr<AdFront3> adfront;
/// 3d generation rules
Array<unique_ptr<vnetrule>> rules;
/// counts how often a rule is used
Array<int> ruleused, canuse, foundmap;
/// describes, why a rule is not applied
Array<string> problems;
/// tolerance criterion
double tolfak;
public:
///
Meshing3 (const string & rulefilename);
///
Meshing3 (const char ** rulep);
///
virtual ~Meshing3 ();
///
void LoadRules (const char * filename, const char ** prules);
///
MESHING3_RESULT GenerateMesh (Mesh & mesh, const MeshingParameters & mp);
///
int ApplyRules (Array<Point3d, PointIndex> & lpoints,
Array<int, PointIndex> & allowpoint,
Array<MiniElement2d> & lfaces, INDEX lfacesplit,
INDEX_2_HASHTABLE<int> & connectedpairs,
NgArray<Element> & elements,
NgArray<INDEX> & delfaces, int tolerance,
double sloppy, int rotind1,
float & retminerr);
///
PointIndex AddPoint (const Point3d & p, PointIndex globind);
///
void AddBoundaryElement (const Element2d & elem);
///
void AddBoundaryElement (const MiniElement2d & elem);
///
int AddConnectedPair (const INDEX_2 & pair);
///
void BlockFill (Mesh & mesh, double gh);
///
void BlockFillLocalH (Mesh & mesh, const MeshingParameters & mp);
/// uses points of adfront, and puts new elements into mesh
void Delaunay (Mesh & mesh, int domainnr, const MeshingParameters & mp);
///
friend class PlotVolMesh;
///
friend void TestRules ();
};
/// status of mesh generation
class MeshingStat3d
{
public:
///
MeshingStat3d ();
///
int cntsucc;
///
int cnttrials;
///
int cntelem;
///
int nff;
///
int qualclass;
///
double vol0;
///
double vol;
///
double h;
///
int problemindex;
};
/*
template <typename POINTArray, typename FACEArray>
extern int FindInnerPoint (POINTArray & grouppoints,
FACEArray & groupfaces,
Point3d & p);
*/
} // namespace netgen
#endif
|
