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
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
|
#ifndef NETGEN_MESHING2_HPP
#define NETGEN_MESHING2_HPP
/**************************************************************************/
/* File: meshing2.hpp */
/* Author: Joachim Schoeberl */
/* Date: 01. Okt. 95 */
/**************************************************************************/
#include "adfront2.hpp"
#include "ruler2.hpp"
#include "basegeom.hpp"
namespace netgen
{
enum MESHING2_RESULT
{
MESHING2_OK = 0,
MESHING2_GIVEUP = 1
};
/*
The basis class for 2D mesh generation.
Has the method GenerateMesh
For surface mesh generation, or non-Euklidean meshing,
derive from Meshing2, and replace transformation.
*/
class Meshing2
{
/// the current advancing front
AdFront2 adfront;
/// rules for mesh generation
Array<unique_ptr<netrule>> rules;
/// statistics
NgArray<int> ruleused, canuse, foundmap;
///
Box<3> boundingbox;
///
double starttime;
///
double maxarea;
Vec3d ex, ey, ez;
Point<3> p1, p2;
const NetgenGeometry& geo;
public:
///
DLL_HEADER Meshing2 (const NetgenGeometry& geo,
const MeshingParameters & mp,
const Box<3> & aboundingbox);
///
DLL_HEADER virtual ~Meshing2 ();
/// Load rules, either from file, or compiled rules
void LoadRules (const char * filename, bool quad);
///
DLL_HEADER MESHING2_RESULT GenerateMesh (Mesh & mesh, const MeshingParameters & mp, double gh, int facenr, int layer=1);
DLL_HEADER void Delaunay (Mesh & mesh, int domainnr, const MeshingParameters & mp);
DLL_HEADER void BlockFillLocalH (Mesh & mesh, const MeshingParameters & mp);
///
DLL_HEADER int AddPoint (const Point3d & p, PointIndex globind, MultiPointGeomInfo * mgi = NULL,
bool pointonsurface = true);
DLL_HEADER PointIndex GetGlobalIndex(int pi) const;
///
DLL_HEADER void AddBoundaryElement (INDEX i1, INDEX i2,
const PointGeomInfo & gi1, const PointGeomInfo & gi2);
///
void SetStartTime (double astarttime);
///
void SetMaxArea (double amaxarea);
protected:
///
virtual void StartMesh ();
///
virtual void EndMesh ();
///
virtual double CalcLocalH (const Point<3> & p, double gh) const;
///
virtual void DefineTransformation (const Point<3> & p1, const Point<3> & p2,
const PointGeomInfo * geominfo1,
const PointGeomInfo * geominfo2);
///
virtual void TransformToPlain (const Point<3> & locpoint, const MultiPointGeomInfo & geominfo,
Point<2> & plainpoint, double h, int & zone);
/// return 0 .. ok
/// return >0 .. cannot transform point to true surface
virtual int TransformFromPlain (const Point<2>& plainpoint,
Point<3> & locpoint,
PointGeomInfo & geominfo,
double h);
/// projects to surface
/// return 0 .. ok
virtual int BelongsToActiveChart (const Point3d & p,
const PointGeomInfo & gi);
/// computes geoinfo data for line with respect to
/// selected chart
virtual int ComputePointGeomInfo (const Point3d & p,
PointGeomInfo & gi);
/// Tries to select unique geominfo on active chart
/// return 0: success
/// return 1: failed
virtual int ChooseChartPointGeomInfo (const MultiPointGeomInfo & mpgi,
PointGeomInfo & pgi);
/*
tests, whether endpoint (= 1 or 2) of line segment p1-p2
is inside of the selected chart. The endpoint must be on the
chart
*/
virtual int IsLineVertexOnChart (const Point3d & p1, const Point3d & p2,
int endpoint, const PointGeomInfo & geominfo);
/*
get (projected) boundary of current chart
*/
virtual void GetChartBoundary (NgArray<Point<2>> & points,
NgArray<Point<3>> & points3d,
NgArray<INDEX_2> & lines, double p) const;
virtual double Area () const;
/** Applies 2D rules.
Tests all 2D rules */
int ApplyRules (NgArray<Point<2>> & lpoints,
NgArray<int> & legalpoints,
int maxlegalpoint,
NgArray<INDEX_2> & llines,
int maxlegelline,
NgArray<Element2d> & elements, NgArray<INDEX> & dellines,
int tolerance,
const MeshingParameters & mp);
};
} // namespace netgen
#endif // NETGEN_MESHING2_HPP
|
