Move setting of wall shear stresses

to dedicated method and use const reference instead of shared pointer

Author: c-p-schmidt

src/fs3i/4C_fs3i_fps3i_partitioned.cpp

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

src/fs3i/4C_fs3i_partitioned.cpp

@@ -762,7 +762,7 @@ void FS3I::PartFS3I::set_wall_shear_stresses() const
   for (unsigned i = 0; i < scatravec_.size(); ++i)
   {
     std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra = scatravec_[i];
-    scatra->scatra_field()->set_wall_shear_stresses(vol_mortar_master_to_slavei(i, wss[i]));
+    scatra->scatra_field()->set_wall_shear_stresses(*vol_mortar_master_to_slavei(i, wss[i]));
   }
 }
 

src/scatra/4C_scatra_timint_implicit.cpp

@@ -1366,22 +1366,18 @@ void ScaTra::ScaTraTimIntImpl::set_external_force() const
 /*----------------------------------------------------------------------*
  *----------------------------------------------------------------------*/
 void ScaTra::ScaTraTimIntImpl::set_wall_shear_stresses(
-    std::shared_ptr<const Core::LinAlg::Vector<double>> wss)
+    const Core::LinAlg::Vector<double>& wall_shear_stress) const
 {
-  if (wss == nullptr) FOUR_C_THROW("WSS 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 wss->get_map().SameAs(*discret_->dof_row_map(nds_wall_shear_stress())))
-    FOUR_C_THROW("Maps are NOT identical. Emergency!");
-#endif
+  FOUR_C_ASSERT(wall_shear_stress.get_map().SameAs(*discret_->dof_row_map(nds_wall_shear_stress())),
+      "Maps are NOT identical. Emergency!");
 
-  discret_->set_state(nds_wall_shear_stress(), "WallShearStress", *wss);
+  discret_->set_state(nds_wall_shear_stress(), "WallShearStress", wall_shear_stress);
 }
 
 /*----------------------------------------------------------------------*

src/scatra/4C_scatra_timint_implicit.hpp

@@ -259,7 +259,8 @@ namespace ScaTra
     //! set the @velocity vector to the scalar transport discretization
     void set_velocity_field(const Core::LinAlg::Vector<double>& velocity);
 
-    void set_wall_shear_stresses(std::shared_ptr<const Core::LinAlg::Vector<double>> wss);
+    //! 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);