[PL/CT] pass process solutions of last times step.

ProcessLib/ComponentTransport/ComponentTransportProcess.cpp

@@ -212,7 +212,8 @@ void ComponentTransportProcess::
 }
 
 void ComponentTransportProcess::solveReactionEquation(
-    std::vector<GlobalVector*>& x, double const t, double const dt)
+    std::vector<GlobalVector*>& x, std::vector<GlobalVector*> const& x_prev,
+    double const t, double const dt)
 {
     if (!_chemical_solver_interface)
     {

ProcessLib/ComponentTransport/ComponentTransportProcess.h

@@ -124,8 +124,9 @@ public:
     void setCoupledTermForTheStaggeredSchemeToLocalAssemblers(
         int const process_id) override;
 
-    void solveReactionEquation(std::vector<GlobalVector*>& x, double const t,
-                               double const dt) override;
+    void solveReactionEquation(std::vector<GlobalVector*>& x,
+                               std::vector<GlobalVector*> const& x_prev,
+                               double const t, double const dt) override;
 
     void extrapolateIntegrationPointValuesToNodes(
         const double t,

ProcessLib/Process.h

@@ -173,8 +173,10 @@ public:
         return Eigen::Vector3d{};
     }
 
-    virtual void solveReactionEquation(std::vector<GlobalVector*>& /*x*/,
-                                       double const /*t*/, double const /*dt*/)
+    virtual void solveReactionEquation(
+        std::vector<GlobalVector*>& /*x*/,
+        std::vector<GlobalVector*> const& /*x_prev*/, double const /*t*/,
+        double const /*dt*/)
     {
     }
 

ProcessLib/TimeLoop.cpp

@@ -851,7 +851,8 @@ TimeLoop::solveCoupledEquationSystemsByStaggeredScheme(
 
     {
         auto& pcs = _per_process_data[0]->process;
-        pcs.solveReactionEquation(_process_solutions, t, dt);
+        pcs.solveReactionEquation(_process_solutions, _process_solutions_prev,
+                                  t, dt);
     }
 
     return nonlinear_solver_status;