FacetArea without restriction in dS results in unrestricted Jacobian which is an error
Issue #90
resolved
Although FacetArea is the same from either side in dS, geometry lowering turns it into an expression of the Jacobian which is not the same from either side and must be restricted.
Comments (2)
-
-
reporter - changed status to resolved
- Log in to comment
I have experienced this issue when implementing
CellVolume
andFacetArea
in TSFC. I eventually opted out for letting TSFC's non-affine implementation handle them in all cases by passingpreserve_geometry_types=(CellVolume, FacetArea)
tocompute_form_data
.