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
|
#ifndef PSURFACE_FACTORY_H
#define PSURFACE_FACTORY_H
// Check for VC9 / VS2008 with installed feature pack.
#if defined(_MSC_VER) && (_MSC_VER>=1500)
// Dummy-include to define _CPPLIB_VER.
#include <vector>
#if defined(_CPPLIB_VER) && _CPPLIB_VER>=505
#include <array>
#else
#error Please install the Visual Studio 2008 SP1 for TR1 support.
#endif
#else
#include <tr1/array>
#endif
#ifdef PSURFACE_STANDALONE
namespace psurface { class Surface; }
#else
class Surface;
#endif
namespace psurface {
template <int dim, class ctype> class PSurface;
template <class ctype, int dim> class StaticVector;
template <int dim, class ctype>
class PSurfaceFactory
{
public:
PSurfaceFactory(PSurface<dim,ctype>* psurface)
: psurface_(psurface)
{}
void setTargetSurface(Surface* surface);
void insertVertex(const StaticVector<ctype,dim+1>& position);
/** \brief Insert a new domain triangle
\return The index of the newly inserted triangle
*/
unsigned int insertSimplex(const std::tr1::array<unsigned int, dim+1>& v);
/**
\param domainVertex if the normal projection hits a base grid vertex, this is the vertex
*/
void insertTargetVertexMapping(unsigned int targetVertex,
unsigned int domainTriangle,
const StaticVector<ctype,dim>& domainLocalPosition,
NodeBundle& projectedTo,
int& domainVertex);
void insertGhostNode(unsigned int domainVertex,
unsigned int targetTriangle,
const StaticVector<ctype,dim>& targetLocalPosition);
void insertEdge();
NodeIdx addInteriorNode(int tri, const StaticVector<ctype,2>& dom, int nodeNumber);
NodeIdx addGhostNode(int tri, int corner, int targetTri, const StaticVector<ctype,2>& localTargetCoords);
NodeIdx addCornerNode(int tri, int corner, int nodeNumber);
NodeBundle addIntersectionNodePair(int tri1, int tri2,
const StaticVector<ctype,2>& dP1, const StaticVector<ctype,2>& dP2,
int edge1, int edge2, const StaticVector<ctype,3>& range);
NodeBundle addBoundaryNode(int tri, const StaticVector<ctype,2>& dP, int edge,
const StaticVector<ctype,3>& range, int targetVert);
NodeIdx addTouchingNode(int tri, const StaticVector<ctype,2>& dP, int edge, int nodeNumber);
NodeBundle addTouchingNodePair(int tri1, int tri2,
const StaticVector<ctype,2>& dP1, const StaticVector<ctype,2>& dP2,
int edge1, int edge2, int nodeNumber);
void addParTriangle(int tri, const std::tr1::array<int,3>& p);
protected:
void addCornerNodeBundle(int domainVertex, int targetVertex);
private:
PSurface<dim,ctype>* psurface_;
};
} // namespace psurface;
#endif
|
