From e8960ca6032292e673d59317a542b476a4f863e8 Mon Sep 17 00:00:00 2001
From: Wenqing Wang <wenqing.wang@ufz.de>
Date: Wed, 5 Dec 2018 15:51:51 +0100
Subject: [PATCH] [TimeLoop] added a functionality to check element
deactivation
---
ProcessLib/UncoupledProcessesTimeLoop.cpp | 10 +++++++++-
1 file changed, 9 insertions(+), 1 deletion(-)
diff --git a/ProcessLib/UncoupledProcessesTimeLoop.cpp b/ProcessLib/UncoupledProcessesTimeLoop.cpp
index 5d18bc31ad6..22ed4a62feb 100644
--- a/ProcessLib/UncoupledProcessesTimeLoop.cpp
+++ b/ProcessLib/UncoupledProcessesTimeLoop.cpp
@@ -542,6 +542,14 @@ bool UncoupledProcessesTimeLoop::loop()
INFO("=== Time stepping at step #%u and time %g with step size %g",
timesteps, t, dt);
+ // Check element deactivation:
+ int process_id = 0;
+ for (auto& process_data : _per_process_data)
+ {
+ process_data->process.checkElementDeactivation(t, process_id);
+ ++process_id;
+ }
+
if (is_staggered_coupling)
{
nonlinear_solver_succeeded =
@@ -803,7 +811,7 @@ bool UncoupledProcessesTimeLoop::solveCoupledEquationSystemsByStaggeredScheme(
timestep_id, t);
}
- unsigned process_id = 0;
+ int process_id = 0;
for (auto& process_data : _per_process_data)
{
if (process_data->skip_time_stepping)
--
GitLab