From 3c8fd1ca9f86b88552898b9cc85e8683a821d9d4 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Tue, 13 Aug 2019 18:39:19 +0200
Subject: [PATCH] [PL] TMPhF; split push back btw. init and postTS.
---
.../ThermoMechanicalPhaseFieldFEM.h | 16 +++++++++++++---
.../ThermoMechanicalPhaseFieldProcess.cpp | 8 ++++++++
2 files changed, 21 insertions(+), 3 deletions(-)
diff --git a/ProcessLib/ThermoMechanicalPhaseField/ThermoMechanicalPhaseFieldFEM.h b/ProcessLib/ThermoMechanicalPhaseField/ThermoMechanicalPhaseFieldFEM.h
index 7d6f68728df..a0d42131320 100644
--- a/ProcessLib/ThermoMechanicalPhaseField/ThermoMechanicalPhaseFieldFEM.h
+++ b/ProcessLib/ThermoMechanicalPhaseField/ThermoMechanicalPhaseFieldFEM.h
@@ -220,9 +220,19 @@ public:
std::vector<double>& local_b_data, std::vector<double>& local_Jac_data,
LocalCoupledSolutions const& local_coupled_solutions) override;
- void preTimestepConcrete(std::vector<double> const& /*local_x*/,
- double const /*t*/,
- double const /*delta_t*/) override
+ void initializeConcrete() override
+ {
+ unsigned const n_integration_points =
+ _integration_method.getNumberOfPoints();
+
+ for (unsigned ip = 0; ip < n_integration_points; ip++)
+ {
+ _ip_data[ip].pushBackState();
+ }
+ }
+
+ void postTimestepConcrete(std::vector<double> const& /*local_x*/
+ ) override
{
unsigned const n_integration_points =
_integration_method.getNumberOfPoints();
diff --git a/ProcessLib/ThermoMechanicalPhaseField/ThermoMechanicalPhaseFieldProcess.cpp b/ProcessLib/ThermoMechanicalPhaseField/ThermoMechanicalPhaseFieldProcess.cpp
index 69be0419107..bb006dc5c26 100644
--- a/ProcessLib/ThermoMechanicalPhaseField/ThermoMechanicalPhaseFieldProcess.cpp
+++ b/ProcessLib/ThermoMechanicalPhaseField/ThermoMechanicalPhaseFieldProcess.cpp
@@ -161,6 +161,14 @@ void ThermoMechanicalPhaseFieldProcess<DisplacementDim>::
_local_assemblers,
&ThermoMechanicalPhaseFieldLocalAssemblerInterface::
getIntPtHeatFlux));
+
+ // Initialize local assemblers after all variables have been set.
+ const int process_id = 0;
+ ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
+
+ GlobalExecutor::executeSelectedMemberOnDereferenced(
+ &LocalAssemblerInterface::initialize, _local_assemblers,
+ pv.getActiveElementIDs(), *_local_to_global_index_map);
}
template <int DisplacementDim>
--
GitLab