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

Fixing ordering bug in RT_FuentesPyramid basis

Browse Source
This commit is contained in:
Stowell, Mark L
2024-05-07 15:52:13 -07:00
parent 15ff1a42cd
commit fb149d7d12
1 changed files with 6 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files
fem/fe/fe_rt.cpp
+6 -6
View File
@@ -1310,9 +1310,9 @@ RT_FuentesPyramidElement::RT_FuentesPyramidElement(const int p)
int o = 0;
// quadrilateral face
for (int j = 0; j <= p; j++)
for (int i = 0; i <= p; i++) // (0,1,2,3)
for (int i = 0; i <= p; i++) // (3,2,1,0)
{
Nodes.IntPoint(o).Set3(bop[i], bop[j], 0.);
Nodes.IntPoint(o).Set3(bop[i], bop[p-j], 0.);
dof2nk[o++] = 0;
}
// triangular faces
@@ -1331,17 +1331,17 @@ RT_FuentesPyramidElement::RT_FuentesPyramidElement(const int p)
dof2nk[o++] = 2;
}
for (int j = 0; j <= p; j++)
for (int i = 0; i + j <= p; i++) // (3,4,2)
for (int i = p - j; i >= 0; i--) // (2,3,4)
{
real_t w = bop[i] + bop[j] + bop[p-i-j];
Nodes.IntPoint(o).Set3(bop[j]/w, 1.0-bop[i]/w, bop[i]/w);
Nodes.IntPoint(o).Set3(bop[i]/w, 1.0-bop[j]/w, bop[j]/w);
dof2nk[o++] = 3;
}
for (int j = 0; j <= p; j++)
for (int i = 0; i + j <= p; i++) // (0,4,3)
for (int i = p - j; i >= 0; i--) // (3,0,4)
{
real_t w = bop[i] + bop[j] + bop[p-i-j];
Nodes.IntPoint(o).Set3(0., bop[j]/w, bop[i]/w);
Nodes.IntPoint(o).Set3(0., bop[i]/w, bop[j]/w);
dof2nk[o++] = 4;
}