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

initial working version for periodic meshes.

Browse Source
This commit is contained in:
Vladimir Z Tomov
2025-02-20 22:38:38 -08:00
parent c97450809d
commit 8f42d760fc
9 changed files with 228 additions and 111 deletions
Download Patch File Download Diff File Expand all files Collapse all files
fem/fe/fe_base.cpp
+30
View File
@@ -973,6 +973,36 @@ void NodalFiniteElement::ProjectDiv(
}
}
void NodalFiniteElement::ReorderToLexicographic(int ncomp, const Vector &dofs,
Vector &dofs_lex) const
{
MFEM_ASSERT(dofs.Size() == ncomp * dof, "Wrong input size.");
dofs_lex.SetSize(ncomp * dof);
for (int i = 0; i < dof; i++)
{
for (int c = 0; c < ncomp; c++)
{
dofs_lex(c*dof + i) = dofs(c*dof + lex_ordering[i]);
}
}
}
void NodalFiniteElement::ReorderFromLexicographic(int ncomp,
const Vector &dofs_lex,
Vector &dofs) const
{
MFEM_ASSERT(dofs_lex.Size() == ncomp * dof, "Wrong input size.");
dofs.SetSize(ncomp * dof);
for (int i = 0; i < dof; i++)
{
for (int c = 0; c < ncomp; c++)
{
dofs(c*dof + lex_ordering[i]) = dofs_lex(c*dof + i);
}
}
}
VectorFiniteElement::VectorFiniteElement(int D, Geometry::Type G,
int Do, int O, int M, int F)