phynexis-deps/mfem
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fixing quad face dof order in pyramid finite elements

Browse Source
This commit is contained in:
Stowell, Mark L
2024-09-05 16:10:49 -07:00
parent 2758154175
commit f89f0eaad4
5 changed files with 18 additions and 17 deletions
Download Patch File Download Diff File Expand all files Collapse all files
fem/fe/fe_nd.cpp
+2 -2
View File
@@ -1684,7 +1684,7 @@ ND_FuentesPyramidElement::ND_FuentesPyramidElement(const int p,
for (int j = 1; j < p; j++)
for (int i = 0; i < p; i++)
{
Nodes.IntPoint(o).Set3(qop[i], qcp[j], 0.);
Nodes.IntPoint(o).Set3(qop[i], qcp[p-j], 0.);
dof2tk[o++] = 0; // (1 0 0)
}
@@ -1692,7 +1692,7 @@ ND_FuentesPyramidElement::ND_FuentesPyramidElement(const int p,
for (int j = 0; j < p; j++)
for (int i = 1; i < p; i++)
{
Nodes.IntPoint(o).Set3(qcp[i], qop[j], 0.);
Nodes.IntPoint(o).Set3(qcp[i], qop[p-1-j], 0.);
dof2tk[o++] = 7; // (0 -1 0)
}