From 516844109b83aae6738828a283910d3e96210184 Mon Sep 17 00:00:00 2001
From: Florian Zill <florian.zill@ufz.de>
Date: Mon, 3 Feb 2020 16:53:25 +0100
Subject: [PATCH] [HM] CheckRequiredProperties on all elements
---
.../HydroMechanics/CreateHydroMechanicsProcess.cpp | 13 ++++++++-----
1 file changed, 8 insertions(+), 5 deletions(-)
diff --git a/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp b/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp
index 40b1264aec3..c381578582e 100644
--- a/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp
+++ b/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp
@@ -143,13 +143,16 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
MaterialPropertyLib::porosity, MaterialPropertyLib::biot_coefficient,
MaterialPropertyLib::density};
- for (auto const& m : media)
+ for (auto const& element : mesh.getElements())
{
- checkRequiredProperties(*m.second, requiredMediumProperties);
- checkRequiredProperties(m.second->phase("Gas"), requiredGasProperties);
- checkRequiredProperties(m.second->phase("Solid"),
- requiredSolidProperties);
+ auto const element_id = element->getID();
+ media_map->checkElementHasMedium(element_id);
+ auto const& medium = *media_map->getMedium(element_id);
+ checkRequiredProperties(medium, requiredMediumProperties);
+ checkRequiredProperties(medium.phase("Gas"), requiredGasProperties);
+ checkRequiredProperties(medium.phase("Solid"), requiredSolidProperties);
}
+ DBUG("Media properties verified.");
// Initial stress conditions
auto const initial_stress = ParameterLib::findOptionalTagParameter<double>(
--
GitLab