[Coupling] Move a member of class Process to the derived class

ProcessLib/ComponentTransport/ComponentTransportProcess.cpp

@@ -63,8 +63,6 @@ void ComponentTransportProcess::assembleConcreteProcess(
     GlobalVector& b)
 {
     DBUG("Assemble ComponentTransportProcess.");
-    if (!_use_monolithic_scheme)
-        setCoupledSolutionsOfPreviousTimeStep();
 
     // Call global assembler for each local assembly item.
     GlobalExecutor::executeMemberDereferenced(
@@ -79,9 +77,6 @@ void ComponentTransportProcess::assembleWithJacobianConcreteProcess(
 {
     DBUG("AssembleWithJacobian ComponentTransportProcess.");
 
-    if (!_use_monolithic_scheme)
-        setCoupledSolutionsOfPreviousTimeStep();
-
     // Call global assembler for each local assembly item.
     GlobalExecutor::executeMemberDereferenced(
         _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,

ProcessLib/HT/HTProcess.cpp

@@ -189,5 +189,29 @@ NumLib::LocalToGlobalIndexMap* HTProcess::getDOFTableForExtrapolatorData(
         NumLib::ComponentOrder::BY_LOCATION);
 }
 
+void HTProcess::setCoupledSolutionsOfPreviousTimeStep()
+{
+    const auto number_of_coupled_solutions =
+        _coupled_solutions->coupled_xs.size();
+    _coupled_solutions->coupled_xs_t0.reserve(number_of_coupled_solutions);
+    const int process_id = _coupled_solutions->process_id;
+    for (std::size_t i = 0; i < number_of_coupled_solutions; i++)
+    {
+        const auto x_t0 = _xs_previous_timestep[process_id].get();
+        if (x_t0 == nullptr)
+        {
+            OGS_FATAL(
+                "Memory is not allocated for the global vector "
+                "of the solution of the previous time step for the ."
+                "staggered scheme.\n It can be done by overloading "
+                "Process::preTimestepConcreteProcess"
+                "(ref. HTProcess::preTimestepConcreteProcess) ");
+        }
+
+        MathLib::LinAlg::setLocalAccessibleVector(*x_t0);
+        _coupled_solutions->coupled_xs_t0.emplace_back(x_t0);
+    }
+}
+
 }  // namespace HT
 }  // namespace ProcessLib

ProcessLib/HT/HTProcess.h

@@ -67,13 +67,6 @@ public:
     bool isLinear() const override { return false; }
     //! @}
 
-    // Get the solution of the previous time step.
-    GlobalVector* getPreviousTimeStepSolution(
-        const int process_id) const override
-    {
-        return _xs_previous_timestep[process_id].get();
-    }
-
     void setCoupledTermForTheStaggeredSchemeToLocalAssemblers() override;
 
 private:
@@ -95,6 +88,10 @@ private:
                                     double const dt,
                                     const int process_id) override;
 
+    /// Set the solutions of the previous time step to the coupled term.
+    /// It only performs for the staggered scheme.
+    void setCoupledSolutionsOfPreviousTimeStep();
+
     NumLib::LocalToGlobalIndexMap* getDOFTableForExtrapolatorData(
         bool& manage_storage) const override;
 

ProcessLib/HydroMechanics/HydroMechanicsProcess-impl.h

@@ -289,7 +289,6 @@ void HydroMechanicsProcess<DisplacementDim>::
     }
 
     // For the staggered scheme
-    setCoupledSolutionsOfPreviousTimeStep();
     // For the equations of displacement
     if (_coupled_solutions->process_id == 1)
     {

ProcessLib/Process.cpp

@@ -375,27 +375,4 @@ NumLib::IterationResult Process::postIteration(const GlobalVector& x)
     return postIterationConcreteProcess(x);
 }
 
-void Process::setCoupledSolutionsOfPreviousTimeStep()
-{
-    const auto number_of_coupled_solutions =
-        _coupled_solutions->coupled_xs.size();
-    _coupled_solutions->coupled_xs_t0.reserve(number_of_coupled_solutions);
-    for (std::size_t i = 0; i < number_of_coupled_solutions; i++)
-    {
-        const auto x_t0 = getPreviousTimeStepSolution(i);
-        if (x_t0 == nullptr)
-        {
-            OGS_FATAL(
-                "Memory is not allocated for the global vector "
-                "of the solution of the previous time step for the ."
-                "staggered scheme.\n It can be done by overloading "
-                "Process::preTimestepConcreteProcess"
-                "(ref. HTProcess::preTimestepConcreteProcess) ");
-        }
-
-        MathLib::LinAlg::setLocalAccessibleVector(*x_t0);
-        _coupled_solutions->coupled_xs_t0.emplace_back(x_t0);
-    }
-}
-
 }  // namespace ProcessLib

ProcessLib/Process.h

@@ -118,14 +118,6 @@ public:
         return _secondary_variables;
     }
 
-    // Get the solution of the previous time step.
-
-    virtual GlobalVector* getPreviousTimeStepSolution(
-        const int /*process_id*/) const
-    {
-        return nullptr;
-    }
-
     // Used as a call back for CalculateSurfaceFlux process.
 
     virtual std::vector<double> getFlux(std::size_t /*element_id*/,
@@ -238,10 +230,6 @@ protected:
     /// references to the solutions of the coupled processes.
     CoupledSolutionsForStaggeredScheme* _coupled_solutions;
 
-    /// Set the solutions of the previous time step to the coupled term.
-    /// It only performs for the staggered scheme.
-    void setCoupledSolutionsOfPreviousTimeStep();
-
     /// Order of the integration method for element-wise integration.
     /// The Gauss-Legendre integration method and available orders is
     /// implemented in MathLib::GaussLegendre.