From b08d3bd8d84a1ce81d9bca4b98a7020fc1377d59 Mon Sep 17 00:00:00 2001
From: "laura.stengel" <laura.stengel@kit.edu>
Date: Thu, 26 Jan 2023 15:04:56 +0100
Subject: [PATCH] Documented SegregatedSolver
---
src/coupled/solvers/SegregatedSolver.cpp | 1 +
src/coupled/solvers/SegregatedSolver.hpp | 3 +++
2 files changed, 4 insertions(+)
diff --git a/src/coupled/solvers/SegregatedSolver.cpp b/src/coupled/solvers/SegregatedSolver.cpp
index a7a7d844b..f1cddf17b 100644
--- a/src/coupled/solvers/SegregatedSolver.cpp
+++ b/src/coupled/solvers/SegregatedSolver.cpp
@@ -61,6 +61,7 @@ SegregatedSolver::EvaluateIteration(const IElphyAssemble &elphyAssemble, const I
evaluationMetrics[metrics[19]].emplace_back(elphyAssemble.L2Norm(mechValues[1]));
evaluationMetrics[metrics[20]].emplace_back(elphyAssemble.H1Norm(mechValues[1]));
+ // Get some metrices from src/elasticity/problem/ElasticityProblem.hpp (for the mechanics part)
auto problemMetrics = mechAssemble.GetElasticityProblem().EvaluationMetrics();
auto problemEntries = mechAssemble.GetElasticityProblem().EvaluationResults(mechValues[0]);
for (int i = 0; i < problemMetrics.size(); ++i) {
diff --git a/src/coupled/solvers/SegregatedSolver.hpp b/src/coupled/solvers/SegregatedSolver.hpp
index 477fae62e..14a950553 100644
--- a/src/coupled/solvers/SegregatedSolver.hpp
+++ b/src/coupled/solvers/SegregatedSolver.hpp
@@ -23,6 +23,7 @@ class SegregatedSolver : public CardiacSolver {
std::unique_ptr<ElphySolver> elphySolver;
std::unique_ptr<ElastodynamicTimeIntegrator> mechSolver;
+ /// Creates some output for each step.
void printStep();
void printFullEvaluation();
void initEvaluation(const IElphyAssemble &elphyAssemble,
@@ -62,8 +63,10 @@ public:
void Initialize(IElphyAssemble &elphyAssemble, IElasticity &mechAssemble, Vector &potential, Vector &displacement);
+ /// Solves coupled equation from mechAssemble.firstTimeStep to mechAssemble.lastTimeStep
void Method(IElphyAssemble &elphyAssemble, IElasticity &mechAssemble, Vector &potential, Vector &displacement);
+ /// Calculates next timestep for the mechanics as well as for the electrophysiology.
void Step(IElphyAssemble &elphyAssemble, IElasticity &mechAssemble, Vectors &elphyValues, Vectors &mechValues);
void Finalize(Vectors &values) {};
--
GitLab