Merge pull request #5528 from anne-glerum/set_face_properties_only_fo…

…r_fem Only set face properties for advection assemblers when method is fem
geodynamics · Jun 5, 2024 · a0c7101 · a0c7101
2 parents 90c3f0a + c3ae449
commit a0c7101
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;