From d0b3e9ea08552d1ba14f2a861ad807010b2e7a1c Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Fri, 23 Feb 2024 17:30:53 +0100
Subject: [PATCH] [PL/TH2] Add global matrix output
---
.../ThermalTwoPhaseFlowWithPPProcess.cpp | 5 +++++
.../ThermalTwoPhaseFlowWithPPProcess.h | 2 ++
2 files changed, 7 insertions(+)
diff --git a/ProcessLib/ThermalTwoPhaseFlowWithPP/ThermalTwoPhaseFlowWithPPProcess.cpp b/ProcessLib/ThermalTwoPhaseFlowWithPP/ThermalTwoPhaseFlowWithPPProcess.cpp
index ca2a4a197b3..7b405f56903 100644
--- a/ProcessLib/ThermalTwoPhaseFlowWithPP/ThermalTwoPhaseFlowWithPPProcess.cpp
+++ b/ProcessLib/ThermalTwoPhaseFlowWithPP/ThermalTwoPhaseFlowWithPPProcess.cpp
@@ -96,6 +96,8 @@ void ThermalTwoPhaseFlowWithPPProcess::assembleConcreteProcess(
_global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
pv.getActiveElementIDs(), dof_table, t, dt, x, x_prev, process_id, M, K,
b);
+
+ _global_output(t, process_id, M, K, b);
}
void ThermalTwoPhaseFlowWithPPProcess::assembleWithJacobianConcreteProcess(
@@ -114,7 +116,10 @@ void ThermalTwoPhaseFlowWithPPProcess::assembleWithJacobianConcreteProcess(
_global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
_local_assemblers, pv.getActiveElementIDs(), dof_table, t, dt, x,
x_prev, process_id, M, K, b, Jac);
+
+ _global_output(t, process_id, M, K, b, &Jac);
}
+
void ThermalTwoPhaseFlowWithPPProcess::preTimestepConcreteProcess(
std::vector<GlobalVector*> const& x, double const t, double const delta_t,
const int process_id)
diff --git a/ProcessLib/ThermalTwoPhaseFlowWithPP/ThermalTwoPhaseFlowWithPPProcess.h b/ProcessLib/ThermalTwoPhaseFlowWithPP/ThermalTwoPhaseFlowWithPPProcess.h
index d887c8bb16a..d58e6fc1405 100644
--- a/ProcessLib/ThermalTwoPhaseFlowWithPP/ThermalTwoPhaseFlowWithPPProcess.h
+++ b/ProcessLib/ThermalTwoPhaseFlowWithPP/ThermalTwoPhaseFlowWithPPProcess.h
@@ -83,6 +83,8 @@ private:
std::vector<
std::unique_ptr<ThermalTwoPhaseFlowWithPPLocalAssemblerInterface>>
_local_assemblers;
+
+ Assembly::GlobalMatrixOutput _global_output;
};
} // namespace ThermalTwoPhaseFlowWithPP
--
GitLab