From 04a20e3b1d1eae51f46e1ac37a85b694a384132e Mon Sep 17 00:00:00 2001
From: Christoph Lehmann <christoph.lehmann@ufz.de>
Date: Tue, 16 Aug 2016 16:50:28 +0200
Subject: [PATCH] pass time to setICs().
---
Applications/ApplicationsLib/UncoupledProcessesTimeLoop.cpp | 2 +-
ProcessLib/Process.cpp | 4 ++--
ProcessLib/Process.h | 2 +-
3 files changed, 4 insertions(+), 4 deletions(-)
diff --git a/Applications/ApplicationsLib/UncoupledProcessesTimeLoop.cpp b/Applications/ApplicationsLib/UncoupledProcessesTimeLoop.cpp
index 4eae2e6234e..cdfdecec8e2 100644
--- a/Applications/ApplicationsLib/UncoupledProcessesTimeLoop.cpp
+++ b/Applications/ApplicationsLib/UncoupledProcessesTimeLoop.cpp
@@ -89,7 +89,7 @@ void UncoupledProcessesTimeLoop::setInitialConditions(
ode_sys.getMatrixSpecifications()));
auto& x0 = *_process_solutions[pcs_idx];
- pcs.setInitialConditions(x0);
+ pcs.setInitialConditions(t0, x0);
MathLib::LinAlg::finalizeAssembly(x0);
time_disc.setInitialState(t0, x0); // push IC
diff --git a/ProcessLib/Process.cpp b/ProcessLib/Process.cpp
index 84d1c17eb88..c72224d7af1 100644
--- a/ProcessLib/Process.cpp
+++ b/ProcessLib/Process.cpp
@@ -69,7 +69,7 @@ void Process::initialize()
_process_variables, *_local_to_global_index_map, _integration_order);
}
-void Process::setInitialConditions(GlobalVector& x)
+void Process::setInitialConditions(double const t, GlobalVector& x)
{
DBUG("Set initial conditions.");
std::size_t const n_nodes = _mesh.getNumberOfNodes();
@@ -93,7 +93,7 @@ void Process::setInitialConditions(GlobalVector& x)
MeshLib::MeshItemType::Node, node_id);
pos.setNodeID(node_id);
- auto const& tup = ic->getTuple(0.0, pos); // 0.0 is t!
+ auto const& tup = ic->getTuple(t, pos);
for (int comp_id = 0; comp_id < num_comp; ++comp_id)
{
diff --git a/ProcessLib/Process.h b/ProcessLib/Process.h
index 7879580d639..b94198a31af 100644
--- a/ProcessLib/Process.h
+++ b/ProcessLib/Process.h
@@ -71,7 +71,7 @@ public:
void initialize();
- void setInitialConditions(GlobalVector& x);
+ void setInitialConditions(const double t, GlobalVector& x);
MathLib::MatrixSpecifications getMatrixSpecifications()
const override final;
--
GitLab