Introduce edge integrals

AUTHORS

@@ -28,7 +28,8 @@ Contributors:
   | Corrado Maurini       <[email protected]>
   | Jack S. Hale          <[email protected]>
   | Tuomas Airaksinen     <[email protected]>
-  | Nacime Bouziani       <[email protected]>
   | Reuben W. Hill        <[email protected]>
+  | Cécile Daversin-Catty <[email protected]>
   | Nacime Bouziani       <[email protected]>
   | Matthew Scroggs       <[email protected]>
+

test/test_classcoverage.py

@@ -11,6 +11,9 @@
     Circumradius,
     Coefficient,
     Constant,
+    EdgeJacobian,
+    EdgeJacobianDeterminant,
+    EdgeJacobianInverse,
     FacetArea,
     FacetNormal,
     FunctionSpace,
@@ -292,9 +295,14 @@ def testAll(self):
     _test_object(g, (), ())
     g = FacetJacobianInverse(domain)
     _test_object(g, (dim - 1, dim), ())
-
-    g = FacetNormal(domain)
-    _test_object(g, (dim,), ())
+    g = EdgeJacobian(domain3D)
+    _test_object(g, (3, 1), ())
+    g = EdgeJacobianDeterminant(domain3D)
+    _test_object(g, (), ())
+    g = EdgeJacobianInverse(domain3D)
+    _test_object(g, (1, 3), ())
+    g = FacetNormal(domain3D)
+    _test_object(g, (3,), ())
     # g = CellNormal(domain)
     # _test_object(g, (dim,), ())
 

test/test_measures.py

@@ -16,6 +16,8 @@ def test_construct_forms_from_default_measures():
     ds = Measure("ds")
     dS = Measure("dS")
 
+    dl = Measure("dl")
+
     dP = Measure("dP")
     # dV = Measure("dV")
 
@@ -37,6 +39,8 @@ def test_construct_forms_from_default_measures():
     assert ds.integral_type() == "exterior_facet"
     assert dS.integral_type() == "interior_facet"
 
+    assert dl.integral_type() == "edge"
+
     assert dP.integral_type() == "vertex"
     # TODO: Change dP to dV:
     # assert dP.integral_type() == "point"

test/test_piecewise_checks.py

@@ -29,6 +29,9 @@
 from ufl.checks import is_cellwise_constant
 from ufl.classes import (
     CellCoordinate,
+    EdgeJacobian,
+    EdgeJacobianDeterminant,
+    EdgeJacobianInverse,
     FacetJacobian,
     FacetJacobianDeterminant,
     FacetJacobianInverse,
@@ -284,6 +287,13 @@ def mappings_are_cellwise_constant(domain, test):
         assert is_cellwise_constant(e) == test
         e = FacetJacobianInverse(domain)
         assert is_cellwise_constant(e) == test
+    if domain.topological_dimension() > 2:
+        e = EdgeJacobian(domain)
+        assert is_cellwise_constant(e) == test
+        e = EdgeJacobianDeterminant(domain)
+        assert is_cellwise_constant(e) == test
+        e = EdgeJacobianInverse(domain)
+        assert is_cellwise_constant(e) == test
 
 
 def test_mappings_are_cellwise_constant_on_linear_affine_cells(affine_domains):