[PL] Infrastructure for setInitialConditions().

Some processes might require setup of the initial
conditions dependent variables in local assemblers.
The new interface allows such implementations.

ProcessLib/LocalAssemblerInterface.cpp

@@ -78,6 +78,16 @@ void LocalAssemblerInterface::computeSecondaryVariable(
     computeSecondaryVariableConcrete(t, local_x);
 }
 
+void LocalAssemblerInterface::setInitialConditions(
+    std::size_t const mesh_item_id,
+    NumLib::LocalToGlobalIndexMap const& dof_table, GlobalVector const& x)
+{
+    auto const indices = NumLib::getIndices(mesh_item_id, dof_table);
+    auto const local_x = x.get(indices);
+
+    setInitialConditionsConcrete(local_x);
+}
+
 void LocalAssemblerInterface::preTimestep(
     std::size_t const mesh_item_id,
     NumLib::LocalToGlobalIndexMap const& dof_table, GlobalVector const& x,

ProcessLib/LocalAssemblerInterface.h

@@ -34,6 +34,10 @@ class LocalAssemblerInterface
 public:
     virtual ~LocalAssemblerInterface() = default;
 
+    virtual void setInitialConditions(
+        std::size_t const mesh_item_id,
+        NumLib::LocalToGlobalIndexMap const& dof_table, GlobalVector const& x);
+
     virtual void preAssemble(double const /*t*/,
                              std::vector<double> const& /*local_x*/){};
 
@@ -96,6 +100,11 @@ public:
     }
 
 private:
+    virtual void setInitialConditionsConcrete(
+        std::vector<double> const& /*local_x*/)
+    {
+    }
+
     virtual void preTimestepConcrete(std::vector<double> const& /*local_x*/,
                                      double const /*t*/, double const /*dt*/)
     {

ProcessLib/Process.cpp

@@ -121,6 +121,7 @@ void Process::setInitialConditions(const int process_id, double const t,
          variable_id < per_process_variables.size();
          variable_id++)
     {
+        MathLib::LinAlg::setLocalAccessibleVector(x);
         SpatialPosition pos;
 
         auto const& pv = per_process_variables[variable_id];
@@ -163,6 +164,7 @@ void Process::setInitialConditions(const int process_id, double const t,
             }
         }
     }
+    setInitialConditionsConcreteProcess(x);
 }
 
 MathLib::MatrixSpecifications Process::getMatrixSpecifications(

ProcessLib/Process.h

@@ -176,6 +176,10 @@ private:
     /// processes. It is called by initialize().
     virtual void initializeBoundaryConditions();
 
+    virtual void setInitialConditionsConcreteProcess(GlobalVector const& /*x*/)
+    {
+    }
+
     virtual void preAssembleConcreteProcess(const double /*t*/,
                                             GlobalVector const& /*x*/)
     {