v0.14.0
Loading...
Searching...
No Matches
Public Member Functions | Private Member Functions | Private Attributes | List of all members
MoFEM::FlatPrismPolynomialBase Struct Reference
MoFEM » Base functions

Calculate base functions on tetrahedralFIXME: Need moab and mofem finite element structure to work (that not perfect) More...

#include <src/approximation/FlatPrismPolynomialBase.hpp>

Inheritance diagram for MoFEM::FlatPrismPolynomialBase:
[legend]
Collaboration diagram for MoFEM::FlatPrismPolynomialBase:
[legend]

Public Member Functions

MoFEMErrorCode query_interface (boost::typeindex::type_index type_index, UnknownInterface **iface) const
 
 FlatPrismPolynomialBase ()
 
 ~FlatPrismPolynomialBase ()
 
MoFEMErrorCode getValue (MatrixDouble &pts, boost::shared_ptr< BaseFunctionCtx > ctx_ptr)
 
- Public Member Functions inherited from MoFEM::BaseFunction
MoFEMErrorCode query_interface (boost::typeindex::type_index type_index, MoFEM::UnknownInterface **iface) const
 
virtual MoFEMErrorCode getValue (MatrixDouble &pts, boost::shared_ptr< BaseFunctionCtx > ctx_ptr)
 
virtual MoFEMErrorCode getValue (MatrixDouble &pts_x, MatrixDouble &pts_t, boost::shared_ptr< BaseFunctionCtx > ctx_ptr)
 
- Public Member Functions inherited from MoFEM::BaseFunctionUnknownInterface
virtual MoFEMErrorCode query_interface (boost::typeindex::type_index type_index, UnknownInterface **iface) const =0
 
virtual ~BaseFunctionUnknownInterface ()=default
 
- Public Member Functions inherited from MoFEM::UnknownInterface
virtual MoFEMErrorCode query_interface (boost::typeindex::type_index type_index, UnknownInterface **iface) const =0
 
template<class IFACE >
MoFEMErrorCode registerInterface (bool error_if_registration_failed=true)
 Register interface. More...
 
template<class IFACE >
MoFEMErrorCode getInterface (IFACE *&iface) const
 Get interface refernce to pointer of interface. More...
 
template<class IFACE >
MoFEMErrorCode getInterface (IFACE **const iface) const
 Get interface pointer to pointer of interface. More...
 
template<class IFACE , typename boost::enable_if< boost::is_pointer< IFACE >, int >::type = 0>
IFACE getInterface () const
 Get interface pointer to pointer of interface. More...
 
template<class IFACE , typename boost::enable_if< boost::is_reference< IFACE >, int >::type = 0>
IFACE getInterface () const
 Get reference to interface. More...
 
template<class IFACE >
IFACE * getInterface () const
 Function returning pointer to interface. More...
 
virtual ~UnknownInterface ()=default
 

Private Member Functions

MoFEMErrorCode getValueH1 (MatrixDouble &pts)
 
MoFEMErrorCode getValueL2 (MatrixDouble &pts)
 
MoFEMErrorCode getValueHdiv (MatrixDouble &pts)
 
MoFEMErrorCode getValueHcurl (MatrixDouble &pts)
 

Private Attributes

FlatPrismPolynomialBaseCtx * cTx
 
int numNodes
 
const EntityHandle * connPrism
 
const EntityHandle * connFace3
 
const EntityHandle * connFace4
 
int faceNodes [2][3]
 
MatrixDouble N
 
MatrixDouble diffN
 

Additional Inherited Members

- Static Public Member Functions inherited from MoFEM::UnknownInterface
static MoFEMErrorCode getLibVersion (Version &version)
 Get library version. More...
 
static MoFEMErrorCode getFileVersion (moab::Interface &moab, Version &version)
 Get database major version. More...
 
static MoFEMErrorCode setFileVersion (moab::Interface &moab, Version version=Version(MoFEM_VERSION_MAJOR, MoFEM_VERSION_MINOR, MoFEM_VERSION_BUILD))
 Get database major version. More...
 
static MoFEMErrorCode getInterfaceVersion (Version &version)
 Get database major version. More...
 

Detailed Description

Calculate base functions on tetrahedral

FIXME: Need moab and mofem finite element structure to work (that not perfect)

Definition at line 44 of file FlatPrismPolynomialBase.hpp.

Constructor & Destructor Documentation

◆ FlatPrismPolynomialBase()

FlatPrismPolynomialBase::FlatPrismPolynomialBase ( )

Definition at line 34 of file FlatPrismPolynomialBase.cpp.

34{}

◆ ~FlatPrismPolynomialBase()

FlatPrismPolynomialBase::~FlatPrismPolynomialBase ( )

Definition at line 33 of file FlatPrismPolynomialBase.cpp.

33{}

Member Function Documentation

◆ getValue()

MoFEMErrorCode FlatPrismPolynomialBase::getValue ( MatrixDouble &  pts,
boost::shared_ptr< BaseFunctionCtx >  ctx_ptr 
)
virtual

Reimplemented from MoFEM::BaseFunction.

Definition at line 37 of file FlatPrismPolynomialBase.cpp.

38 {
39 MoFEMFunctionBeginHot;
40
41 cTx = ctx_ptr->getInterface<FlatPrismPolynomialBaseCtx>();
42
43 if (!cTx->fePtr)
44 SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
45 "Pointer to element should be given "
46 "when EntPolynomialBaseCtx is constructed "
47 "(use different constructor)");
48
49 int nb_gauss_pts = pts.size2();
50 if (!nb_gauss_pts)
51 MoFEMFunctionReturnHot(0);
52
53 if (pts.size1() < 1)
54 SETERRQ(
55 PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
56 "Wrong dimension of pts, should be at least 3 rows with coordinates");
57
58 const FieldApproximationBase base = cTx->bAse;
59 EntitiesFieldData &data = cTx->dAta;
60
61 if (cTx->copyNodeBase == LASTBASE)
62 SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "Not implemented");
63 else
64 data.dataOnEntities[MBVERTEX][0].getNSharedPtr(base) =
65 data.dataOnEntities[MBVERTEX][0].getNSharedPtr(cTx->copyNodeBase);
66
67 data.dataOnEntities[MBVERTEX][0].getN(base).resize(nb_gauss_pts, 6, false);
68 data.dataOnEntities[MBVERTEX][0].getDiffN(base).resize(nb_gauss_pts, 12,
69 false);
70 if (data.dataOnEntities[MBVERTEX][0].getN(base).size1() !=
71 (unsigned int)nb_gauss_pts) {
72 SETERRQ1(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
73 "Base functions or nodes has wrong number of integration points "
74 "for base %s",
75 ApproximationBaseNames[base]);
76 }
77
78 data.dataOnEntities[MBVERTEX][0].getN(base).resize(nb_gauss_pts, 6, false);
79 data.dataOnEntities[MBVERTEX][0].getDiffN(base).resize(nb_gauss_pts, 12,
80 false);
81 N.resize(nb_gauss_pts, 3, false);
82 diffN.resize(3, 2, false);
83 CHKERR ShapeMBTRI(&*N.data().begin(), &pts(0, 0), &pts(1, 0), nb_gauss_pts);
84 std::copy(Tools::diffShapeFunMBTRI.begin(), Tools::diffShapeFunMBTRI.end(),
85 &*diffN.data().begin());
86
87 // This is needed to have proper order of nodes on faces
88 CHKERR cTx->mOab.get_connectivity(cTx->fePtr->getEnt(), connPrism, numNodes,
89 true);
90 SideNumber_multiIndex &side_table =
91 const_cast<SideNumber_multiIndex &>(cTx->fePtr->getSideNumberTable());
92 SideNumber_multiIndex::nth_index<1>::type::iterator siit3 =
93 side_table.get<1>().find(boost::make_tuple(MBTRI, 3));
94 SideNumber_multiIndex::nth_index<1>::type::iterator siit4 =
95 side_table.get<1>().find(boost::make_tuple(MBTRI, 4));
96 if (siit3 == side_table.get<1>().end())
97 SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
98 if (siit4 == side_table.get<1>().end())
99 SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
100 CHKERR cTx->mOab.get_connectivity(siit3->get()->ent, connFace3, numNodes,
101 true);
102 CHKERR cTx->mOab.get_connectivity(siit4->get()->ent, connFace4, numNodes,
103 true);
104
105 for (int nn = 0; nn < 3; nn++) {
106 faceNodes[0][nn] = std::distance(
107 connPrism, std::find(connPrism, connPrism + 3, connFace3[nn]));
108 faceNodes[1][nn] = std::distance(
109 connPrism + 3, std::find(connPrism + 3, connPrism + 6, connFace4[nn]));
110 for (int gg = 0; gg < nb_gauss_pts; gg++) {
111 double val = N(gg, nn);
112 double val_x = diffN(nn, 0);
113 double val_y = diffN(nn, 1);
114 data.dataOnEntities[MBVERTEX][0].getN(base)(gg, nn) = val;
115 data.dataOnEntities[MBVERTEX][0].getDiffN(base)(gg, 2 * nn + 0) = val_x;
116 data.dataOnEntities[MBVERTEX][0].getDiffN(base)(gg, 2 * nn + 1) = val_y;
117 data.dataOnEntities[MBVERTEX][0].getN(base)(gg, 3 + nn) = val;
118 data.dataOnEntities[MBVERTEX][0].getDiffN(base)(gg, 6 + 2 * nn + 0) =
119 val_x;
120 data.dataOnEntities[MBVERTEX][0].getDiffN(base)(gg, 6 + 2 * nn + 1) =
121 val_y;
122 }
123 }
124
125 switch (cTx->sPace) {
126 case H1:
127 CHKERR getValueH1(pts);
128 break;
129 case HDIV:
130 CHKERR getValueHdiv(pts);
131 break;
132 case HCURL:
133 CHKERR getValueHcurl(pts);
134 break;
135 case L2:
136 CHKERR getValueL2(pts);
137 break;
138 default:
139 SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "Not yet implemented");
140 }
141
142 MoFEMFunctionReturnHot(0);
143}
FieldApproximationBase
FieldApproximationBase
approximation base
Definition: definitions.h:58
LASTBASE
@ LASTBASE
Definition: definitions.h:69
MoFEMFunctionReturnHot
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:447
L2
@ L2
field with C-1 continuity
Definition: definitions.h:88
H1
@ H1
continuous field
Definition: definitions.h:85
HCURL
@ HCURL
field with continuous tangents
Definition: definitions.h:86
HDIV
@ HDIV
field with continuous normal traction
Definition: definitions.h:87
MOFEM_DATA_INCONSISTENCY
@ MOFEM_DATA_INCONSISTENCY
Definition: definitions.h:31
MOFEM_NOT_IMPLEMENTED
@ MOFEM_NOT_IMPLEMENTED
Definition: definitions.h:32
ApproximationBaseNames
static const char *const ApproximationBaseNames[]
Definition: definitions.h:72
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:535
MoFEMFunctionBeginHot
#define MoFEMFunctionBeginHot
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:440
ShapeMBTRI
PetscErrorCode ShapeMBTRI(double *N, const double *X, const double *Y, const int G_DIM)
calculate shape functions on triangle
Definition: fem_tools.c:182
SideNumber_multiIndex
multi_index_container< boost::shared_ptr< SideNumber >, indexed_by< ordered_unique< member< SideNumber, EntityHandle, &SideNumber::ent > >, ordered_non_unique< composite_key< SideNumber, const_mem_fun< SideNumber, EntityType, &SideNumber::getEntType >, member< SideNumber, signed char, &SideNumber::side_number > > > > > SideNumber_multiIndex
SideNumber_multiIndex for SideNumber.
Definition: RefEntsMultiIndices.hpp:101
MoFEM::EntPolynomialBaseCtx::copyNodeBase
const FieldApproximationBase copyNodeBase
Definition: EntPolynomialBaseCtx.hpp:40
MoFEM::EntPolynomialBaseCtx::bAse
const FieldApproximationBase bAse
Definition: EntPolynomialBaseCtx.hpp:38
MoFEM::EntitiesFieldData
data structure for finite element entity
Definition: EntitiesFieldData.hpp:40
MoFEM::EntitiesFieldData::dataOnEntities
std::array< boost::ptr_vector< EntData >, MBMAXTYPE > dataOnEntities
Definition: EntitiesFieldData.hpp:56
MoFEM::FlatPrismPolynomialBaseCtx
Class used to pass element data to calculate base functions on flat prism.
Definition: FlatPrismPolynomialBase.hpp:21
MoFEM::FlatPrismPolynomialBaseCtx::fePtr
const NumeredEntFiniteElement * fePtr
Definition: FlatPrismPolynomialBase.hpp:27
MoFEM::FlatPrismPolynomialBase::getValueHcurl
MoFEMErrorCode getValueHcurl(MatrixDouble &pts)
Definition: FlatPrismPolynomialBase.cpp:228
MoFEM::FlatPrismPolynomialBase::cTx
FlatPrismPolynomialBaseCtx * cTx
Definition: FlatPrismPolynomialBase.hpp:56
MoFEM::FlatPrismPolynomialBase::getValueHdiv
MoFEMErrorCode getValueHdiv(MatrixDouble &pts)
Definition: FlatPrismPolynomialBase.cpp:224
MoFEM::FlatPrismPolynomialBase::getValueH1
MoFEMErrorCode getValueH1(MatrixDouble &pts)
Definition: FlatPrismPolynomialBase.cpp:145
MoFEM::FlatPrismPolynomialBase::getValueL2
MoFEMErrorCode getValueL2(MatrixDouble &pts)
Definition: FlatPrismPolynomialBase.cpp:218
MoFEM::Tools::diffShapeFunMBTRI
static constexpr std::array< double, 6 > diffShapeFunMBTRI
Definition: Tools.hpp:104
MoFEM::UnknownInterface::getInterface
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface refernce to pointer of interface.
Definition: UnknownInterface.hpp:93
MoFEM::interface_EntFiniteElement::getSideNumberTable
SideNumber_multiIndex & getSideNumberTable() const
Definition: FEMultiIndices.hpp:705
MoFEM::interface_RefEntity::getEnt
EntityHandle getEnt() const
Definition: RefEntsMultiIndices.hpp:603

◆ getValueH1()

MoFEMErrorCode FlatPrismPolynomialBase::getValueH1 ( MatrixDouble &  pts)
private

Definition at line 145 of file FlatPrismPolynomialBase.cpp.

145 {
146 MoFEMFunctionBeginHot;
147
148 EntitiesFieldData &data = cTx->dAta;
149 const FieldApproximationBase base = cTx->bAse;
150 PetscErrorCode (*base_polynomials)(int p, double s, double *diff_s, double *L,
151 double *diffL, const int dim) =
152 cTx->basePolynomialsType0;
153
154 int nb_gauss_pts = pts.size2();
155
156 // edges
157 if (data.dataOnEntities[MBEDGE].size() != 9)
158 SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
159
160 int sense[9], order[9];
161 double *H1edgeN[9], *diffH1edgeN[9];
162
163 auto set_edge_base_data = [&](const int ee) {
164 MoFEMFunctionBegin;
165 auto &ent_data = data.dataOnEntities[MBEDGE][ee];
166 if (ent_data.getSense() == 0)
167 SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
168 sense[ee] = ent_data.getSense();
169 order[ee] = ent_data.getOrder();
170 int nb_dofs = NBEDGE_H1(ent_data.getOrder());
171 ent_data.getN(base).resize(nb_gauss_pts, nb_dofs, false);
172 ent_data.getDiffN(base).resize(nb_gauss_pts, 2 * nb_dofs, false);
173 H1edgeN[ee] = &*ent_data.getN(base).data().begin();
174 diffH1edgeN[ee] = &*ent_data.getDiffN(base).data().begin();
175 MoFEMFunctionReturn(0);
176 };
177
178 if ((data.spacesOnEntities[MBEDGE]).test(H1)) {
179
180 for (int ee = 0; ee != 3; ++ee)
181 CHKERR set_edge_base_data(ee);
182 for (int ee = 6; ee != 9; ++ee)
183 CHKERR set_edge_base_data(ee);
184
185 // shape functions on face 3
186 CHKERR H1_EdgeShapeFunctions_MBTRI(
187 &sense[0], &order[0], &*N.data().begin(), &*diffN.data().begin(),
188 &H1edgeN[0], &diffH1edgeN[0], nb_gauss_pts, base_polynomials);
189 // shape functions on face 4
190 CHKERR H1_EdgeShapeFunctions_MBTRI(
191 &sense[6], &order[6], &*N.data().begin(), &*diffN.data().begin(),
192 &H1edgeN[6], &diffH1edgeN[6], nb_gauss_pts, base_polynomials);
193 }
194
195 // face
196 if (data.dataOnEntities[MBTRI].size() != 5)
197 SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
198
199 if ((data.spacesOnEntities[MBTRI]).test(H1)) {
200 for (int ff = 3; ff <= 4; ff++) {
201 int nb_dofs = NBFACETRI_H1(data.dataOnEntities[MBTRI][ff].getOrder());
202 data.dataOnEntities[MBTRI][ff].getN(base).resize(nb_gauss_pts, nb_dofs,
203 false);
204 data.dataOnEntities[MBTRI][ff].getDiffN(base).resize(nb_gauss_pts,
205 2 * nb_dofs, false);
206 CHKERR H1_FaceShapeFunctions_MBTRI(
207 faceNodes[ff - 3], data.dataOnEntities[MBTRI][ff].getOrder(),
208 &*N.data().begin(), &*diffN.data().begin(),
209 &*data.dataOnEntities[MBTRI][ff].getN(base).data().begin(),
210 &*data.dataOnEntities[MBTRI][ff].getDiffN(base).data().begin(),
211 nb_gauss_pts, base_polynomials);
212 }
213 }
214
215 MoFEMFunctionReturnHot(0);
216}
p
static Index< 'p', 3 > p
Definition: BasicFeTools.hpp:90
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:346
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:416
dim
const int dim
Definition: elec_phys_2D.cpp:12
order
int order
Definition: fluid_structure_eigenproblem.cpp:84
NBEDGE_H1
#define NBEDGE_H1(P)
Numer of base function on edge for H1 space.
Definition: h1_hdiv_hcurl_l2.h:55
H1_EdgeShapeFunctions_MBTRI
PetscErrorCode H1_EdgeShapeFunctions_MBTRI(int *sense, int *p, double *N, double *diffN, double *edgeN[3], double *diff_edgeN[3], int GDIM, PetscErrorCode(*base_polynomials)(int p, double s, double *diff_s, double *L, double *diffL, const int dim))
H1_EdgeShapeFunctions_MBTRI.
Definition: h1.c:17
H1_FaceShapeFunctions_MBTRI
PetscErrorCode H1_FaceShapeFunctions_MBTRI(const int *face_nodes, int p, double *N, double *diffN, double *faceN, double *diff_faceN, int GDIM, PetscErrorCode(*base_polynomials)(int p, double s, double *diff_s, double *L, double *diffL, const int dim))
Definition: h1.c:191
NBFACETRI_H1
#define NBFACETRI_H1(P)
Number of base function on triangle for H1 space.
Definition: h1_hdiv_hcurl_l2.h:60
MoFEM::EntPolynomialBaseCtx::basePolynomialsType0
PetscErrorCode(* basePolynomialsType0)(int p, double s, double *diff_s, double *L, double *diffL, const int dim)
Definition: EntPolynomialBaseCtx.hpp:27
MoFEM::EntitiesFieldData::spacesOnEntities
std::array< std::bitset< LASTSPACE >, MBMAXTYPE > spacesOnEntities
spaces on entity types
Definition: EntitiesFieldData.hpp:50

◆ getValueHcurl()

MoFEMErrorCode FlatPrismPolynomialBase::getValueHcurl ( MatrixDouble &  pts)
private

Definition at line 228 of file FlatPrismPolynomialBase.cpp.

228 {
229 SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "Not yet implemented");
230}

◆ getValueHdiv()

MoFEMErrorCode FlatPrismPolynomialBase::getValueHdiv ( MatrixDouble &  pts)
private

Definition at line 224 of file FlatPrismPolynomialBase.cpp.

224 {
225 SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "Not yet implemented");
226}

◆ getValueL2()

MoFEMErrorCode FlatPrismPolynomialBase::getValueL2 ( MatrixDouble &  pts)
private

Definition at line 218 of file FlatPrismPolynomialBase.cpp.

218 {
219 MoFEMFunctionBeginHot;
220 SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "Not yet implemented");
221 MoFEMFunctionReturnHot(0);
222}

◆ query_interface()

MoFEMErrorCode FlatPrismPolynomialBase::query_interface ( boost::typeindex::type_index  type_index,
UnknownInterface **  iface 
) const
virtual

Reimplemented from MoFEM::BaseFunction.

Definition at line 27 of file FlatPrismPolynomialBase.cpp.

28 {
29 *iface = const_cast<FlatPrismPolynomialBase *>(this);
30 return 0;
31}
MoFEM::FlatPrismPolynomialBase
Calculate base functions on tetrahedralFIXME: Need moab and mofem finite element structure to work (t...
Definition: FlatPrismPolynomialBase.hpp:44

Member Data Documentation

◆ connFace3

const EntityHandle* MoFEM::FlatPrismPolynomialBase::connFace3
private

Definition at line 68 of file FlatPrismPolynomialBase.hpp.

◆ connFace4

const EntityHandle* MoFEM::FlatPrismPolynomialBase::connFace4
private

Definition at line 69 of file FlatPrismPolynomialBase.hpp.

◆ connPrism

const EntityHandle* MoFEM::FlatPrismPolynomialBase::connPrism
private

Definition at line 67 of file FlatPrismPolynomialBase.hpp.

◆ cTx

FlatPrismPolynomialBaseCtx* MoFEM::FlatPrismPolynomialBase::cTx
private

Definition at line 56 of file FlatPrismPolynomialBase.hpp.

◆ diffN

MatrixDouble MoFEM::FlatPrismPolynomialBase::diffN
private

Definition at line 72 of file FlatPrismPolynomialBase.hpp.

◆ faceNodes

int MoFEM::FlatPrismPolynomialBase::faceNodes[2][3]
private

Definition at line 70 of file FlatPrismPolynomialBase.hpp.

◆ N

MatrixDouble MoFEM::FlatPrismPolynomialBase::N
private

Definition at line 71 of file FlatPrismPolynomialBase.hpp.

◆ numNodes

int MoFEM::FlatPrismPolynomialBase::numNodes
private

Definition at line 66 of file FlatPrismPolynomialBase.hpp.

The documentation for this struct was generated from the following files:
Generated by Doxygen 1.9.5 and hosted at University of Glasgow