Move setting of pressures

to dedicated method and use const reference instead of shared pointer

Author: c-p-schmidt

src/fs3i/4C_fs3i_fps3i_partitioned.cpp

@@ -685,7 +685,7 @@ void FS3I::PartFPS3I::set_pressure_fields()
   for (unsigned i = 0; i < scatravec_.size(); ++i)
   {
     std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra = scatravec_[i];
-    scatra->scatra_field()->set_pressure_field(pressure[i]);
+    scatra->scatra_field()->set_pressure_field(*pressure[i]);
   }
 }
 

src/scatra/4C_scatra_timint_implicit.cpp

@@ -1391,22 +1391,18 @@ void ScaTra::ScaTraTimIntImpl::set_old_part_of_righthandside()
 /*----------------------------------------------------------------------*
  *----------------------------------------------------------------------*/
 void ScaTra::ScaTraTimIntImpl::set_pressure_field(
-    std::shared_ptr<const Core::LinAlg::Vector<double>> pressure)
+    const Core::LinAlg::Vector<double>& pressure) const
 {
-  if (pressure == nullptr) FOUR_C_THROW("Pressure state is nullptr");
-
-#ifdef FOUR_C_ENABLE_ASSERTIONS
   // We rely on the fact, that the nodal distribution of both fields is the same.
   // Although Scatra discretization was constructed as a clone of the fluid or
   // structure mesh, respectively, at the beginning, the nodal distribution may
   // have changed meanwhile (e.g., due to periodic boundary conditions applied only
   // to the fluid field)!
   // We have to be sure that everything is still matching.
-  if (not pressure->get_map().SameAs(*discret_->dof_row_map(nds_pressure())))
-    FOUR_C_THROW("Maps are NOT identical. Emergency!");
-#endif
+  FOUR_C_ASSERT(pressure.get_map().SameAs(*discret_->dof_row_map(nds_pressure())),
+      "Maps are NOT identical. Emergency!");
 
-  discret_->set_state(nds_pressure(), "Pressure", *pressure);
+  discret_->set_state(nds_pressure(), "Pressure", pressure);
 }
 
 /*----------------------------------------------------------------------*

src/scatra/4C_scatra_timint_implicit.hpp

@@ -262,7 +262,8 @@ namespace ScaTra
     //! set the @wall_shear_stress vector to the scalar transport discretization
     void set_wall_shear_stresses(const Core::LinAlg::Vector<double>& wall_shear_stress) const;
 
-    void set_pressure_field(std::shared_ptr<const Core::LinAlg::Vector<double>> pressure);
+    //! set the @pressure vector to the scalar transport discretization
+    void set_pressure_field(const Core::LinAlg::Vector<double>& pressure) const;
 
     void set_membrane_concentration(
         std::shared_ptr<const Core::LinAlg::Vector<double>> MembraneConc);