compute_form_data: remove do_assume_single_integral_type

Author: ksagiyam

test/test_mixed_function_space_with_mesh_sequence.py

@@ -68,7 +68,6 @@ def test_mixed_function_space_with_mesh_sequence_cell():
         do_estimate_degrees=True,
         do_replace_functions=True,
         coefficients_to_split=(f, g),
-        do_assume_single_integral_type=False,
         complex_mode=False,
     )
     (id0,) = fd.integral_data
@@ -137,7 +136,6 @@ def test_mixed_function_space_with_mesh_sequence_facet():
         do_estimate_degrees=True,
         do_replace_functions=True,
         coefficients_to_split=(f, g),
-        do_assume_single_integral_type=False,
         complex_mode=False,
     )
     (
@@ -237,7 +235,6 @@ def test_mixed_function_space_with_mesh_sequence_raise():
             do_estimate_degrees=True,
             do_replace_functions=True,
             coefficients_to_split=(f,),
-            do_assume_single_integral_type=False,
             complex_mode=False,
         )
     assert e_info.match("Found multiple domains, cannot return just one.")
@@ -254,7 +251,6 @@ def test_mixed_function_space_with_mesh_sequence_raise():
             do_estimate_degrees=True,
             do_replace_functions=True,
             coefficients_to_split=(f, g),
-            do_assume_single_integral_type=False,
             complex_mode=False,
         )
     assert e_info.match("Discontinuous type Coefficient must be restricted.")

ufl/algorithms/compute_form_data.py

@@ -262,7 +262,6 @@ def compute_form_data(
     do_append_everywhere_integrals=True,
     do_replace_functions=False,
     coefficients_to_split=None,
-    do_assume_single_integral_type=True,
     complex_mode=False,
     do_remove_component_tensors=False,
 ):
@@ -469,22 +468,11 @@ def compute_form_data(
     # Propagate restrictions to terminals
     if do_apply_restrictions:
         for itg_data in self.integral_data:
-            if do_assume_single_integral_type:  # Conventional case
-                if not itg_data.integral_type.startswith("interior_facet"):
-                    continue
             if do_apply_default_restrictions:
-                if do_assume_single_integral_type:  # Conventional case
-                    domains = set(
-                        domain
-                        for integral in itg_data.integrals
-                        for domain in extract_domains(integral)
-                    )
-                    default_restrictions = {domain: "+" for domain in domains}
-                else:
-                    default_restrictions = {
-                        domain: default_restriction_map[integral_type]
-                        for domain, integral_type in itg_data.domain_integral_type_map.items()
-                    }
+                default_restrictions = {
+                    domain: default_restriction_map[integral_type]
+                    for domain, integral_type in itg_data.domain_integral_type_map.items()
+                }
             else:
                 default_restrictions = None
             new_integrals = []