-
Notifications
You must be signed in to change notification settings - Fork 13
/
Copy pathFETri6.cpp
75 lines (52 loc) · 1.92 KB
/
FETri6.cpp
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
#include "FETri6.hpp"
namespace femus {
// =======================
FETri6::FETri6() : FEElemBase() { }
// =======================
FETri6::~FETri6() { }
// =======================
float FETri6::get_embedding_matrix(const uint a ,const uint b,const uint c ) { return _embedding_matrix[a][b][c]; }
// =======================
//STATIC data member
const float FETri6::_embedding_matrix[4][6][6] =
{ // ---------------------------------------
{// embedding matrix for child 0
// 0 1 2 3 4 5
{ 1.0, 0.0, 0.0, 0.0, 0.0, 0.0}, // 0
{ 0.0, 0.0, 0.0, 1.0, 0.0, 0.0}, // 1
{ 0.0, 0.0, 0.0, 0.0, 0.0, 1.0}, // 2
{.375, -.125, 0.0, .75, 0.0, 0.0}, // 3
{ 0.0, -.125, -.125, 0.5, .25, 0.5}, // 4
{.375, 0.0, -.125, 0.0, 0.0, .75} // 5
},
{// embedding matrix for child 1
// 0 1 2 3 4 5
{ 0.0, 0.0, 0.0, 1.0, 0.0, 0.0}, // 0
{ 0.0, 1.0, 0.0, 0.0, 0.0, 0.0}, // 1
{ 0.0, 0.0, 0.0, 0.0, 1.0, 0.0}, // 2
{-.125, .375, 0.0, .75, 0.0, 0.0}, // 3
{ 0.0, .375, -.125, 0.0, .75, 0.0}, // 4
{-.125, 0.0, -.125, 0.5, 0.5, .25} // 5
},
// embedding matrix for child 2
{
// 0 1 2 3 4 5
{ 0.0, 0.0, 0.0, 0.0, 0.0, 1.0}, // 0
{ 0.0, 0.0, 0.0, 0.0, 1.0, 0.0}, // 1
{ 0.0, 0.0, 1.0, 0.0, 0.0, 0.0}, // 2
{-.125, -.125, 0.0, .25, 0.5, 0.5}, // 3
{ 0.0, -.125, .375, 0.0, .75, 0.0}, // 4
{-.125, 0.0, .375, 0.0, 0.0, .75} // 5
},
// embedding matrix for child 3
{
// 0 1 2 3 4 5
{ 0.0, 0.0, 0.0, 1.0, 0.0, 0.0}, // 0
{ 0.0, 0.0, 0.0, 0.0, 1.0, 0.0}, // 1
{ 0.0, 0.0, 0.0, 0.0, 0.0, 1.0}, // 2
{-.125, 0.0, -.125, 0.5, 0.5, .25}, // 3
{-.125, -.125, 0.0, .25, 0.5, 0.5}, // 4
{ 0.0, -.125, -.125, 0.5, .25, 0.5} // 5
}
};
} //end namespace femus