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
|
// Gmsh - Copyright (C) 1997-2021 C. Geuzaine, J.-F. Remacle
//
// See the LICENSE.txt file for license information. Please report all
// issues on https://gitlab.onelab.info/gmsh/gmsh/issues.
//
// Contributed by Ismail Badia.
// Reference : "Higher-Order Finite Element Methods"; Pavel Solin, Karel
// Segeth ,
// Ivo Dolezel , Chapman and Hall/CRC; Edition : Har/Cdr (2003).
#ifndef HIERARCHICAL_BASIS_HCURL_BRICK_H
#define HIERARCHICAL_BASIS_HCURL_BRICK_H
#include "HierarchicalBasisHcurl.h"
/*
* MHexahedron
*
* v
* 3----------2
* |\ ^ |\
* | \ | | \
* | \ | | \
* | 7------+---6
* | | +-- |-- | -> u
* 0---+---\--1 |
* \ | \ \ |
* \ | \ \ |
* \| w \|
* 4----------5
*
* Oriented Edges:
* e0={0, 1}, e1={0, 3}, e2={0, 4}, e3={1, 2}, e4 ={1, 5}, e5={3, 2},e6={2, 6},
* e7={3, 7},e8={4, 5}, e9= {4, 7}, e10={5, 6}, e11={7, 6}
*
* Oritented Surface:
* s0={0, 1, 3, 2}, s1={0, 1, 4, 5}, s2={0, 3, 4, 7},
* s3={1, 2, 5, 6}, s4={2, 3, 7, 6}, s5={4, 5, 7, 6}
* Local (directional) orders on mesh faces are not allowed to exceed the mini-
* mum of the (appropriate directional) orders of approximation associated with
* the interior of the adjacent elements. Local orders of approximation on mesh
* edges are limited by the minimum of all (appropriate directional) orders cor-
* responding to faces sharing that edge
*/
class HierarchicalBasisHcurlBrick : public HierarchicalBasisHcurl {
public:
HierarchicalBasisHcurlBrick(int order);
virtual ~HierarchicalBasisHcurlBrick();
virtual unsigned int getNumberOfOrientations() const;
virtual void generateBasis(double const &u, double const &v, double const &w,
std::vector<std::vector<double> > &vertexBasis,
std::vector<std::vector<double> > &edgeBasis,
std::vector<std::vector<double> > &faceBasis,
std::vector<std::vector<double> > &bubbleBasis,
std::string typeFunction)
{
if(typeFunction == "HcurlLegendre") {
generateHcurlBasis(u, v, w, edgeBasis, faceBasis, bubbleBasis);
}
else if("CurlHcurlLegendre" == typeFunction) {
generateCurlBasis(u, v, w, edgeBasis, faceBasis, bubbleBasis);
}
else {
throw std::string("unknown typeFunction");
}
};
virtual void
orientEdge(int const &flagOrientation, int const &edgeNumber,
std::vector<std::vector<double> > &edgeBasis,
const std::vector<std::vector<double> > &eTablePositiveFlag,
const std::vector<std::vector<double> > &eTableNegativeFlag);
virtual void orientEdgeFunctionsForNegativeFlag(
std::vector<std::vector<double> > &edgeFunctions);
virtual void orientFace(
int const &flag1, int const &flag2, int const &flag3, int const &faceNumber,
const std::vector<std::vector<double> > &quadFaceFunctionsAllOrientation,
const std::vector<std::vector<double> > &triFaceFunctionsAllOrientation,
std::vector<std::vector<double> > &fTableCopy);
virtual void getKeysInfo(std::vector<int> &functionTypeInfo,
std::vector<int> &orderInfo);
private:
int _pb1; // bubble function order in direction u
int _pb2; // bubble function order in direction v
int _pb3; // bubble function order in direction w
int _pOrderEdge[12]; // Edge functions order (pOrderEdge[0] matches the order
// of the edge 0)
int _pOrderFace1[6]; // Face functions order in direction u (pOrderFace1[0]
// matches the order of face 0 in direction u)
int _pOrderFace2[6]; // Face functions order in direction v (pOrderFace[0]
// matches the order of face 0 in direction v)
static double
_affineCoordinate(const int &j, const double &u, const double &v,
const double &w); // affine coordinate lambdaj j=1..6
// edgeBasis=[phie0_{0},...phie0_{pe0},phie1_{0},...phie1_{pe1}...]
// faceBasis=[phieFf1{n1,n2} (with 0<=n1<=pf1 , 2<=n2<=pf2+1), phieFf2{n1,n2}
// (with 2<=n1<=pf1+1 , 0<=n2<=pf2) ] bubbleBasis=[phieb1{n1,n2,n3} (with
// 0<=n1<=pb1 , 2<=n2<=pb2+1 , 2<=n3<=pb3+1)...]
virtual void
generateHcurlBasis(double const &u, double const &v, double const &w,
std::vector<std::vector<double> > &edgeBasis,
std::vector<std::vector<double> > &faceBasis,
std::vector<std::vector<double> > &bubbleBasis);
virtual void
generateCurlBasis(double const &u, double const &v, double const &w,
std::vector<std::vector<double> > &edgeBasis,
std::vector<std::vector<double> > &faceBasis,
std::vector<std::vector<double> > &bubbleBasis);
virtual void orientOneFace(double const &u, double const &v, double const &w,
int const &flag1, int const &flag2,
int const &flag3, int const &faceNumber,
std::vector<std::vector<double> > &faceFunctions,
std::string typeFunction);
};
#endif
|
