Only set face properties when method is fem

geodynamics · Jun 4, 2024 · c3ae449 · c3ae449
1 parent 6c5e7d4
commit c3ae449
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/source/simulator/assembly.cc b/source/simulator/assembly.cc
@@ -192,7 +192,8 @@ namespace aspect
             assemblers->advection_system_assembler_on_face_properties[0].need_face_finite_element_evaluation = true;
           }
 
-        if (i > 0 && parameters.use_discontinuous_composition_discretization)
+        if (i > 0 && parameters.use_discontinuous_composition_discretization
+            && parameters.compositional_field_methods[i-1] == Parameters<dim>::AdvectionFieldMethod::fem_field)
           {
             // TODO should these only be set when method is fem_field?
             assemblers->advection_system_assembler_on_face_properties[i].need_face_material_model_data = true;