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

Adding access to H1 polynomials

Browse Source
This commit is contained in:
Stowell, Mark L
2024-06-05 16:17:15 -07:00
parent 95705c8121
commit 83a34a3e77
2 changed files with 31 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files
fem/fe/fe_h1.cpp
+28
View File
@@ -1199,6 +1199,34 @@ void H1_FuentesPyramidElement::CalcDShape(const IntegrationPoint &ip,
Ti.Mult(du, dshape);
}
void H1_FuentesPyramidElement::CalcRawShape(const IntegrationPoint &ip,
Vector &shape) const
{
const int p = order;
#ifdef MFEM_THREAD_SAFE
Vector tmp_i(p + 1);
DenseMatrix tmp1_ij(p + 1, p + 1);
#endif
calcBasis(p, ip, tmp_i, tmp1_ij, shape);
}
void H1_FuentesPyramidElement::CalcRawDShape(const IntegrationPoint &ip,
DenseMatrix &dshape) const
{
const int p = order;
#ifdef MFEM_THREAD_SAFE
Vector tmp_i(p + 1);
DenseMatrix tmp1_ij(p + 1, p + 1);
DenseMatrix tmp2_ij(p + 1, dim);
DenseTensor tmp_ijk(p + 1, p + 1, dim);
#endif
calcGradBasis(p, ip, tmp_i, tmp2_ij, tmp1_ij, tmp_ijk, dshape);
}
void H1_FuentesPyramidElement::calcBasis(const int p,
const IntegrationPoint &ip,
Vector &phi_i, DenseMatrix &phi_ij,