From c3ae449bf65ee1f9ea0bed3d0c2431fe23fa3ce7 Mon Sep 17 00:00:00 2001 From: anne-glerum Date: Mon, 18 Dec 2023 21:02:23 +0100 Subject: [PATCH] Only set face properties when method is fem --- source/simulator/assembly.cc | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/source/simulator/assembly.cc b/source/simulator/assembly.cc index 0d0600748ac..fe09c592e2a 100644 --- 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::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;