diff --git a/ProcessLib/ComponentTransport/ComponentTransportProcess.cpp b/ProcessLib/ComponentTransport/ComponentTransportProcess.cpp
index 95af9b5969b673efeb38e33cba3a262f16c9d9ea..fcec0ad9110000441047ee5e1269d4a6f435a865 100644
--- a/ProcessLib/ComponentTransport/ComponentTransportProcess.cpp
+++ b/ProcessLib/ComponentTransport/ComponentTransportProcess.cpp
@@ -199,36 +199,33 @@ void ComponentTransportProcess::assembleConcreteProcess(
MathLib::finalizeMatrixAssembly(M);
MathLib::finalizeMatrixAssembly(K);
MathLib::finalizeVectorAssembly(b);
- auto const residuum = computeResiduum(*x[0], *xdot[0], M, K, b);
- auto copyRhs = [&](int const variable_id, auto& output_vector)
- {
- if (_use_monolithic_scheme)
- {
- transformVariableFromGlobalVector(residuum, variable_id,
- dof_tables[0], output_vector,
- std::negate<double>());
- }
- else
- {
- transformVariableFromGlobalVector(
- residuum, 0, dof_tables[process_id], output_vector,
- std::negate<double>());
- }
- };
+
if (_use_monolithic_scheme || process_id == 0)
{
- copyRhs(0, *_hydraulic_flow);
- }
- if (_use_monolithic_scheme)
- {
- for (std::size_t c_idx = 0; c_idx < _molar_flow_rate.size(); ++c_idx)
+ auto const residuum = computeResiduum(*x[0], *xdot[0], M, K, b);
+ transformVariableFromGlobalVector(residuum, 0, dof_tables[0],
+ *_hydraulic_flow,
+ std::negate<double>());
+ if (_use_monolithic_scheme)
{
- copyRhs(c_idx + 1, *_molar_flow_rate[c_idx]);
+ for (std::size_t c_idx = 0; c_idx < _molar_flow_rate.size();
+ ++c_idx)
+ {
+ transformVariableFromGlobalVector(
+ residuum, c_idx + 1, dof_tables[0],
+ *_molar_flow_rate[c_idx], std::negate<double>());
+ }
}
+ return;
}
+
if (process_id > 0)
{
- copyRhs(process_id, *_molar_flow_rate[process_id - 1]);
+ auto const residuum =
+ computeResiduum(*x[process_id], *xdot[process_id], M, K, b);
+ transformVariableFromGlobalVector(residuum, 0, dof_tables[process_id],
+ *_molar_flow_rate[process_id - 1],
+ std::negate<double>());
}
}
diff --git a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay.prj b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay.prj
index 99dfd890fad5474304e4cf75468db957e7638088..044761213620a67f974a729b2eadc8143bf3f548 100644
--- a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay.prj
@@ -265,7 +265,7 @@
</vtkdiff>
<vtkdiff>
<regex>1D_AdvectionDiffusionSorptionDecay_ts_[0-9]*_t_[0-9]*.000000.vtu</regex>
- <field>MolarFlowRate</field>
+ <field>CsFlowRate</field>
<absolute_tolerance>1e-10</absolute_tolerance>
<relative_tolerance>1e-16</relative_tolerance>
</vtkdiff>
diff --git a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_0_t_0.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_0_t_0.000000.vtu
index f3ad438758cd75c486474624f19b13460a7d4c98..da1aa31304b65912c861fc7854c1933d2d2ba963 100644
--- a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_0_t_0.000000.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_0_t_0.000000.vtu
@@ -7,7 +7,7 @@
<Piece NumberOfPoints="2001" NumberOfCells="2000" >
<PointData>
<DataArray type="Float64" Name="Cs" format="appended" RangeMin="0" RangeMax="0" offset="84" />
- <DataArray type="Float64" Name="MolarFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="180" />
+ <DataArray type="Float64" Name="CsFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="180" />
<DataArray type="Float64" Name="pressure" format="appended" RangeMin="100000" RangeMax="100000" offset="276" />
</PointData>
<CellData>
diff --git a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_1000_t_31536000000000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_1000_t_31536000000000.000000.vtu
index 2a9d7cf400840917eaf5aaaa4a42f5f8bc437e5b..47e8cca6bdbda2e0cc09185bd8ef9067c85d073a 100644
--- a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_1000_t_31536000000000.000000.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_1000_t_31536000000000.000000.vtu
@@ -7,7 +7,7 @@
<Piece NumberOfPoints="2001" NumberOfCells="2000" >
<PointData>
<DataArray type="Float64" Name="Cs" format="appended" RangeMin="1.2774319966e-137" RangeMax="1" offset="84" />
- <DataArray type="Float64" Name="MolarFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="21356" />
+ <DataArray type="Float64" Name="CsFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="21356" />
<DataArray type="Float64" Name="pressure" format="appended" RangeMin="100000" RangeMax="499999.99999" offset="21452" />
</PointData>
<CellData>
diff --git a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_100_t_3153600000000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_100_t_3153600000000.000000.vtu
index 226895a230b2ba94366ecb34cab2aa2fb7e76f31..eff4fda3ae53f5dcdcca388cc225d63004c6e40c 100644
--- a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_100_t_3153600000000.000000.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_100_t_3153600000000.000000.vtu
@@ -7,7 +7,7 @@
<Piece NumberOfPoints="2001" NumberOfCells="2000" >
<PointData>
<DataArray type="Float64" Name="Cs" format="appended" RangeMin="-8.4963431577e-318" RangeMax="1" offset="84" />
- <DataArray type="Float64" Name="MolarFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="17948" />
+ <DataArray type="Float64" Name="CsFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="17948" />
<DataArray type="Float64" Name="pressure" format="appended" RangeMin="100000" RangeMax="499999.99999" offset="18044" />
</PointData>
<CellData>
diff --git a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_10_t_315360000000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_10_t_315360000000.000000.vtu
index 96b7e138157cae56be036b8d4eb10c361b3d5f94..cc375c66f5aa702496ef1471dcbc963b065df541 100644
--- a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_10_t_315360000000.000000.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_10_t_315360000000.000000.vtu
@@ -7,7 +7,7 @@
<Piece NumberOfPoints="2001" NumberOfCells="2000" >
<PointData>
<DataArray type="Float64" Name="Cs" format="appended" RangeMin="2.8191385752e-318" RangeMax="1" offset="84" />
- <DataArray type="Float64" Name="MolarFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="13680" />
+ <DataArray type="Float64" Name="CsFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="13680" />
<DataArray type="Float64" Name="pressure" format="appended" RangeMin="100000" RangeMax="499999.99999" offset="13776" />
</PointData>
<CellData>
diff --git a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_1_t_31536000000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_1_t_31536000000.000000.vtu
index 81421ad97fe237e120aa31d1c776a51a30e3d888..2df8dc3b16ea4b0329c1997da962e6052641db2c 100644
--- a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_1_t_31536000000.000000.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay_ts_1_t_31536000000.000000.vtu
@@ -7,7 +7,7 @@
<Piece NumberOfPoints="2001" NumberOfCells="2000" >
<PointData>
<DataArray type="Float64" Name="Cs" format="appended" RangeMin="8.6939545941e-318" RangeMax="1" offset="84" />
- <DataArray type="Float64" Name="MolarFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="12912" />
+ <DataArray type="Float64" Name="CsFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="12912" />
<DataArray type="Float64" Name="pressure" format="appended" RangeMin="100000" RangeMax="499999.99999" offset="13008" />
</PointData>
<CellData>
diff --git a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay.prj b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay.prj
index cc9e20d4633e2811a0e7e4cee3201b4013906bae..482e6778f7c2bb6148ee62a43e5fdc57c98f1dcd 100644
--- a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay.prj
@@ -264,7 +264,7 @@
</vtkdiff>
<vtkdiff>
<regex>1D_DiffusionSorptionDecay_ts_[0-9]*_t_[0-9]*.000000.vtu</regex>
- <field>MolarFlowRate</field>
+ <field>CsFlowRate</field>
<absolute_tolerance>1e-10</absolute_tolerance>
<relative_tolerance>1e-16</relative_tolerance>
</vtkdiff>
diff --git a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_0_t_0.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_0_t_0.000000.vtu
index e894b1f22712123aedd3b5ebef2a0720c705e684..cadf3dcad7c97a10d6706fc91f12c1680590cc3c 100644
--- a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_0_t_0.000000.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_0_t_0.000000.vtu
@@ -7,7 +7,7 @@
<Piece NumberOfPoints="161" NumberOfCells="160" >
<PointData>
<DataArray type="Float64" Name="Cs" format="appended" RangeMin="0" RangeMax="0" offset="84" />
- <DataArray type="Float64" Name="MolarFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="152" />
+ <DataArray type="Float64" Name="CsFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="152" />
<DataArray type="Float64" Name="pressure" format="appended" RangeMin="100000" RangeMax="100000" offset="220" />
</PointData>
<CellData>
diff --git a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_1000_t_31536000000000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_1000_t_31536000000000.000000.vtu
index e41aa48d4f5a6ea29816d6cf4849a371c65cb168..68781b1a84be4ed7c45243424d0eac9fa0dbbc48 100644
--- a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_1000_t_31536000000000.000000.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_1000_t_31536000000000.000000.vtu
@@ -7,7 +7,7 @@
<Piece NumberOfPoints="161" NumberOfCells="160" >
<PointData>
<DataArray type="Float64" Name="Cs" format="appended" RangeMin="-1.3032412715e-19" RangeMax="1" offset="84" />
- <DataArray type="Float64" Name="MolarFlowRate" format="appended" RangeMin="-2.1175823681e-22" RangeMax="1.4152178956e-08" offset="1860" />
+ <DataArray type="Float64" Name="CsFlowRate" format="appended" RangeMin="-2.1175823681e-22" RangeMax="1.4152178956e-08" offset="1860" />
<DataArray type="Float64" Name="pressure" format="appended" RangeMin="100000" RangeMax="100000" offset="2464" />
</PointData>
<CellData>
diff --git a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_100_t_3153600000000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_100_t_3153600000000.000000.vtu
index 44a5fa6da6db9a4dab70b4f6b3699ce24389cc60..118ea21493fec0a0bb4087abc8538bad7b8625b7 100644
--- a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_100_t_3153600000000.000000.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_100_t_3153600000000.000000.vtu
@@ -7,7 +7,7 @@
<Piece NumberOfPoints="161" NumberOfCells="160" >
<PointData>
<DataArray type="Float64" Name="Cs" format="appended" RangeMin="-8.2816973106e-22" RangeMax="1" offset="84" />
- <DataArray type="Float64" Name="MolarFlowRate" format="appended" RangeMin="-2.0514079191e-22" RangeMax="3.5930953041e-08" offset="1860" />
+ <DataArray type="Float64" Name="CsFlowRate" format="appended" RangeMin="-2.0514079191e-22" RangeMax="3.5930953041e-08" offset="1860" />
<DataArray type="Float64" Name="pressure" format="appended" RangeMin="100000" RangeMax="100000" offset="2416" />
</PointData>
<CellData>
diff --git a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_10_t_315360000000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_10_t_315360000000.000000.vtu
index 0003c65884801665d5b00fd0f936cd25f8477824..629ec7a281a908e54c697461e3a8c872f422e607 100644
--- a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_10_t_315360000000.000000.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_10_t_315360000000.000000.vtu
@@ -7,7 +7,7 @@
<Piece NumberOfPoints="161" NumberOfCells="160" >
<PointData>
<DataArray type="Float64" Name="Cs" format="appended" RangeMin="-2.3888639398e-34" RangeMax="1" offset="84" />
- <DataArray type="Float64" Name="MolarFlowRate" format="appended" RangeMin="-5.2939559203e-23" RangeMax="1.144954842e-07" offset="1860" />
+ <DataArray type="Float64" Name="CsFlowRate" format="appended" RangeMin="-5.2939559203e-23" RangeMax="1.144954842e-07" offset="1860" />
<DataArray type="Float64" Name="pressure" format="appended" RangeMin="100000" RangeMax="100000" offset="2408" />
</PointData>
<CellData>
diff --git a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_1_t_31536000000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_1_t_31536000000.000000.vtu
index 403a82973d56bfa1dcd0b24c40852de8ded3d818..1adbf9a07ca46a3d67753d7270bdfd8e5cbc869a 100644
--- a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_1_t_31536000000.000000.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay_ts_1_t_31536000000.000000.vtu
@@ -7,7 +7,7 @@
<Piece NumberOfPoints="161" NumberOfCells="160" >
<PointData>
<DataArray type="Float64" Name="Cs" format="appended" RangeMin="-1.3308255536e-45" RangeMax="1" offset="84" />
- <DataArray type="Float64" Name="MolarFlowRate" format="appended" RangeMin="-2.6469779602e-23" RangeMax="6.2627750966e-07" offset="1860" />
+ <DataArray type="Float64" Name="CsFlowRate" format="appended" RangeMin="-2.6469779602e-23" RangeMax="6.2627750966e-07" offset="1860" />
<DataArray type="Float64" Name="pressure" format="appended" RangeMin="100000" RangeMax="100000" offset="2452" />
</PointData>
<CellData>