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

add specialized L2 restriction from fine to coarse for discontinuous space elements

Browse Source
This commit is contained in:
Robert W. Anderson
2022-05-18 21:51:20 -07:00
parent 86332c87ee
commit 66959f111e
5 changed files with 35 additions and 13 deletions
Download Patch File Download Diff File Expand all files Collapse all files
fem/fe/fe_base.cpp
+3 -2
View File
@@ -456,7 +456,7 @@ void ScalarFiniteElement::ScalarLocalInterpolation(
}
}
void ScalarFiniteElement::ScalarLocalRestriction(
void ScalarFiniteElement::ScalarLocalL2Restriction(
ElementTransformation &Trans, DenseMatrix &R,
const ScalarFiniteElement &coarse_fe) const
{
@@ -489,7 +489,8 @@ void ScalarFiniteElement::ScalarLocalRestriction(
Trans.Transform(f_ip, vv);
c_ip.Set(v, dim);
coarse_fe.CalcShape(c_ip, coarse_shape);
AddMult_a_VWt(f_ip.weight*Trans.Weight(), coarse_shape, fine_shape, coarse_fine_mass);
AddMult_a_VWt(f_ip.weight*Trans.Weight(), coarse_shape, fine_shape,
coarse_fine_mass);
}
DenseMatrixInverse coarse_mass_inv(coarse_mass);