FacetArea without restriction in dS results in unrestricted Jacobian which is an error

Issue #90 resolved
Martin Sandve Alnæs created an issue

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)

  1. Miklós Homolya

    I have experienced this issue when implementing CellVolume and FacetArea in TSFC. I eventually opted out for letting TSFC's non-affine implementation handle them in all cases by passing preserve_geometry_types=(CellVolume, FacetArea) to compute_form_data.

  2. Log in to comment