Change apply mesh movement to work with

c-p-schmidt · c-p-schmidt · commit 80b8cc4fd8ac · 2025-05-02T10:44:47.000+02:00
const reference instead of shared pointer
diff --git a/src/elch/4C_elch_moving_boundary_algorithm.cpp b/src/elch/4C_elch_moving_boundary_algorithm.cpp
@@ -81,7 +81,7 @@ void ElCh::MovingBoundaryAlgorithm::setup()
   }
 
   // transfer moving mesh data
-  scatra_field()->apply_mesh_movement(ale_field()->dispnp());
+  scatra_field()->apply_mesh_movement(*ale_field()->dispnp());
 
   // initialize the multivector for all possible cases
   fluxn_ = scatra_field()->calc_flux_at_boundary(false);
@@ -116,7 +116,7 @@ void ElCh::MovingBoundaryAlgorithm::time_loop()
   }
 
   // transfer moving mesh data
-  scatra_field()->apply_mesh_movement(ale_field()->dispnp());
+  scatra_field()->apply_mesh_movement(*ale_field()->dispnp());
 
   // time loop
   while (not_finished())
@@ -278,7 +278,7 @@ void ElCh::MovingBoundaryAlgorithm::solve_scatra()
   }
 
   // transfer moving mesh data
-  scatra_field()->apply_mesh_movement(ale_field()->dispnp());
+  scatra_field()->apply_mesh_movement(*ale_field()->dispnp());
 
   // solve coupled electrochemistry equations
   scatra_field()->solve();
diff --git a/src/fs3i/4C_fs3i_fps3i_partitioned.cpp b/src/fs3i/4C_fs3i_fps3i_partitioned.cpp
@@ -615,11 +615,11 @@ void FS3I::PartFPS3I::set_struct_scatra_solution()
 void FS3I::PartFPS3I::set_mesh_disp()
 {
   // fluid field
-  scatravec_[0]->scatra_field()->apply_mesh_movement(fpsi_->fluid_field()->dispnp());
+  scatravec_[0]->scatra_field()->apply_mesh_movement(*fpsi_->fluid_field()->dispnp());
 
   // Poro field
   scatravec_[1]->scatra_field()->apply_mesh_movement(
-      fpsi_->poro_field()->structure_field()->dispnp());
+      *fpsi_->poro_field()->structure_field()->dispnp());
 }
 
 
diff --git a/src/fs3i/4C_fs3i_partitioned.cpp b/src/fs3i/4C_fs3i_partitioned.cpp
@@ -683,12 +683,12 @@ void FS3I::PartFS3I::set_mesh_disp() const
   // fluid field
   std::shared_ptr<Adapter::ScaTraBaseAlgorithm> fluidscatra = scatravec_[0];
   fluidscatra->scatra_field()->apply_mesh_movement(
-      fluid_to_fluid_scalar(fsi_->fluid_field()->dispnp()));
+      *fluid_to_fluid_scalar(fsi_->fluid_field()->dispnp()));
 
   // structure field
   std::shared_ptr<Adapter::ScaTraBaseAlgorithm> structscatra = scatravec_[1];
   structscatra->scatra_field()->apply_mesh_movement(
-      structure_to_structure_scalar(fsi_->structure_field()->dispnp()));
+      *structure_to_structure_scalar(fsi_->structure_field()->dispnp()));
 }
 
 
diff --git a/src/poroelast_scatra/4C_poroelast_scatra_base.cpp b/src/poroelast_scatra/4C_poroelast_scatra_base.cpp
@@ -203,7 +203,7 @@ void PoroElastScaTra::PoroScatraBase::set_mesh_disp()
     dispnp = volcoupl_fluidscatra_->apply_vector_mapping21(*fluid_field()->dispnp());
   }
 
-  scatra_->scatra_field()->apply_mesh_movement(dispnp);
+  scatra_->scatra_field()->apply_mesh_movement(*dispnp);
 
   std::shared_ptr<const Core::LinAlg::Vector<double>> sdispnp = nullptr;
 
diff --git a/src/porofluid_pressure_based_elast_scatra/4C_porofluid_pressure_based_elast_scatra_base.cpp b/src/porofluid_pressure_based_elast_scatra/4C_porofluid_pressure_based_elast_scatra_base.cpp
@@ -292,7 +292,7 @@ void PoroPressureBased::PoroMultiPhaseScaTraBase::set_poro_solution()
   if (poroscatra == nullptr) FOUR_C_THROW("cast to ScaTraTimIntPoroMulti failed!");
 
   // set displacements
-  poroscatra->apply_mesh_movement(poromulti_->struct_dispnp());
+  poroscatra->apply_mesh_movement(*poromulti_->struct_dispnp());
 
   // set the fluid solution
   poroscatra->set_solution_field_of_multi_fluid(
diff --git a/src/scatra/4C_scatra_timint_implicit.cpp b/src/scatra/4C_scatra_timint_implicit.cpp
@@ -1644,8 +1644,7 @@ void ScaTra::ScaTraTimIntImpl::update()
 
 /*----------------------------------------------------------------------*
  *----------------------------------------------------------------------*/
-void ScaTra::ScaTraTimIntImpl::apply_mesh_movement(
-    std::shared_ptr<const Core::LinAlg::Vector<double>> dispnp)
+void ScaTra::ScaTraTimIntImpl::apply_mesh_movement(const Core::LinAlg::Vector<double>& dispnp) const
 {
   //---------------------------------------------------------------------------
   // only required in ALE case
@@ -1654,12 +1653,9 @@ void ScaTra::ScaTraTimIntImpl::apply_mesh_movement(
   {
     TEUCHOS_FUNC_TIME_MONITOR("SCATRA: apply mesh movement");
 
-    // check existence of displacement vector
-    if (dispnp == nullptr) FOUR_C_THROW("Got null pointer for displacements!");
-
     // provide scatra discretization with displacement field
-    discret_->set_state(nds_disp(), "dispnp", *dispnp);
-  }  // if (isale_)
+    discret_->set_state(nds_disp(), "dispnp", dispnp);
+  }
 }
 
 /*----------------------------------------------------------------------*
diff --git a/src/scatra/4C_scatra_timint_implicit.hpp b/src/scatra/4C_scatra_timint_implicit.hpp
@@ -383,7 +383,7 @@ namespace ScaTra
      *
      * @param[in] dispnp  displacement vector
      */
-    void apply_mesh_movement(std::shared_ptr<const Core::LinAlg::Vector<double>> dispnp);
+    void apply_mesh_movement(const Core::LinAlg::Vector<double>& dispnp) const;
 
     //! calculate fluxes inside domain and/or on boundary
     void calc_flux(const bool writetofile  //!< flag for writing flux info to file
diff --git a/src/ssi/4C_ssi_base.cpp b/src/ssi/4C_ssi_base.cpp
@@ -114,7 +114,7 @@ void SSI::SSIBase::setup()
   ssicoupling_->setup();
 
   // in case of an ssi  multi scale formulation we need to set the displacement here
-  auto dummy_vec = std::make_shared<Core::LinAlg::Vector<double>>(
+  auto dummy_vec = Core::LinAlg::Vector<double>(
       *Global::Problem::instance()->get_dis("structure")->dof_row_map(), true);
   ssicoupling_->set_mesh_disp(scatra_base_algorithm(), dummy_vec);
 
@@ -519,7 +519,7 @@ void SSI::SSIBase::test_results(MPI_Comm comm) const
 
 /*----------------------------------------------------------------------*/
 /*----------------------------------------------------------------------*/
-void SSI::SSIBase::set_struct_solution(std::shared_ptr<const Core::LinAlg::Vector<double>> disp,
+void SSI::SSIBase::set_struct_solution(const Core::LinAlg::Vector<double>& disp,
     std::shared_ptr<const Core::LinAlg::Vector<double>> vel, const bool set_mechanical_stress)
 {
   // safety checks
@@ -614,7 +614,7 @@ void SSI::SSIBase::set_mechanical_stress_state(
 
 /*----------------------------------------------------------------------*/
 /*----------------------------------------------------------------------*/
-void SSI::SSIBase::set_mesh_disp(std::shared_ptr<const Core::LinAlg::Vector<double>> disp)
+void SSI::SSIBase::set_mesh_disp(const Core::LinAlg::Vector<double>& disp)
 {
   // safety checks
   check_is_init();
diff --git a/src/ssi/4C_ssi_base.hpp b/src/ssi/4C_ssi_base.hpp
@@ -214,7 +214,7 @@ namespace SSI
     [[nodiscard]] std::shared_ptr<ScaTra::ScaTraTimIntImpl> scatra_manifold() const;
 
     /// set structure solution on other fields
-    void set_struct_solution(std::shared_ptr<const Core::LinAlg::Vector<double>> disp,
+    void set_struct_solution(const Core::LinAlg::Vector<double>& disp,
         std::shared_ptr<const Core::LinAlg::Vector<double>> vel, bool set_mechanical_stress);
 
     /// set scatra solution on other fields
@@ -377,7 +377,7 @@ namespace SSI
     [[nodiscard]] bool is_init() const { return isinit_; }
 
     /// set structure mesh displacement on scatra field
-    void set_mesh_disp(std::shared_ptr<const Core::LinAlg::Vector<double>> disp);
+    void set_mesh_disp(const Core::LinAlg::Vector<double>& disp);
 
     /// set structure velocity field on scatra field
     void set_velocity_fields(std::shared_ptr<const Core::LinAlg::Vector<double>> vel);
diff --git a/src/ssi/4C_ssi_coupling.cpp b/src/ssi/4C_ssi_coupling.cpp
@@ -129,8 +129,7 @@ void SSI::SSICouplingMatchingVolume::set_mechanical_stress_state(
 /*----------------------------------------------------------------------*/
 /*----------------------------------------------------------------------*/
 void SSI::SSICouplingMatchingVolume::set_mesh_disp(
-    std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,
-    std::shared_ptr<const Core::LinAlg::Vector<double>> disp)
+    std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra, const Core::LinAlg::Vector<double>& disp)
 {
   scatra->scatra_field()->apply_mesh_movement(disp);
 }
@@ -278,11 +277,10 @@ void SSI::SSICouplingNonMatchingBoundary::assign_material_pointers(
 /*----------------------------------------------------------------------*/
 /*----------------------------------------------------------------------*/
 void SSI::SSICouplingNonMatchingBoundary::set_mesh_disp(
-    std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,
-    std::shared_ptr<const Core::LinAlg::Vector<double>> disp)
+    std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra, const Core::LinAlg::Vector<double>& disp)
 {
   scatra->scatra_field()->apply_mesh_movement(
-      adaptermeshtying_->master_to_slave(*extractor_->extract_cond_vector(*disp)));
+      *adaptermeshtying_->master_to_slave(*extractor_->extract_cond_vector(disp)));
 }
 
 /*----------------------------------------------------------------------*/
@@ -402,11 +400,10 @@ void SSI::SSICouplingNonMatchingVolume::assign_material_pointers(
 /*----------------------------------------------------------------------*/
 /*----------------------------------------------------------------------*/
 void SSI::SSICouplingNonMatchingVolume::set_mesh_disp(
-    std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,
-    std::shared_ptr<const Core::LinAlg::Vector<double>> disp)
+    std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra, const Core::LinAlg::Vector<double>& disp)
 {
   scatra->scatra_field()->apply_mesh_movement(
-      volcoupl_structurescatra_->apply_vector_mapping21(*disp));
+      *volcoupl_structurescatra_->apply_vector_mapping21(disp));
 }
 
 /*----------------------------------------------------------------------*/
@@ -604,8 +601,7 @@ void SSI::SSICouplingMatchingVolumeAndBoundary::assign_material_pointers(
 /*----------------------------------------------------------------------*/
 /*----------------------------------------------------------------------*/
 void SSI::SSICouplingMatchingVolumeAndBoundary::set_mesh_disp(
-    std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,
-    std::shared_ptr<const Core::LinAlg::Vector<double>> disp)
+    std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra, const Core::LinAlg::Vector<double>& disp)
 {
   scatra->scatra_field()->apply_mesh_movement(disp);
 }
diff --git a/src/ssi/4C_ssi_coupling.hpp b/src/ssi/4C_ssi_coupling.hpp
@@ -69,7 +69,7 @@ namespace SSI
     //! \param scatra    underlying scatra problem of the SSI problem
     //! \param disp      displacement field to set
     virtual void set_mesh_disp(std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,
-        std::shared_ptr<const Core::LinAlg::Vector<double>> disp) = 0;
+        const Core::LinAlg::Vector<double>& disp) = 0;
 
     //! \brief set structure velocity fields on other field
     //!
@@ -119,7 +119,7 @@ namespace SSI
         unsigned nds) override;
 
     void set_mesh_disp(std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,
-        std::shared_ptr<const Core::LinAlg::Vector<double>> disp) override;
+        const Core::LinAlg::Vector<double>& disp) override;
 
     void set_velocity_fields(std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,
         std::shared_ptr<const Core::LinAlg::Vector<double>> convvel,
@@ -172,8 +172,8 @@ namespace SSI
   class SSICouplingNonMatchingBoundary : public SSICouplingBase
   {
    public:
-    SSICouplingNonMatchingBoundary()
-        : adaptermeshtying_(nullptr), extractor_(nullptr), issetup_(false), isinit_(false) {};
+    SSICouplingNonMatchingBoundary() = default;
+
     void init(const int ndim, std::shared_ptr<Core::FE::Discretization> structdis,
         std::shared_ptr<SSI::SSIBase> ssi_base) override;
 
@@ -189,7 +189,7 @@ namespace SSI
     }
 
     void set_mesh_disp(std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,
-        std::shared_ptr<const Core::LinAlg::Vector<double>> disp) override;
+        const Core::LinAlg::Vector<double>& disp) override;
 
     void set_velocity_fields(std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,
         std::shared_ptr<const Core::LinAlg::Vector<double>> convvel,
@@ -278,7 +278,7 @@ namespace SSI
     }
 
     void set_mesh_disp(std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,
-        std::shared_ptr<const Core::LinAlg::Vector<double>> disp) override;
+        const Core::LinAlg::Vector<double>& disp) override;
 
     void set_velocity_fields(std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,
         std::shared_ptr<const Core::LinAlg::Vector<double>> convvel,
@@ -355,7 +355,7 @@ namespace SSI
     }
 
     void set_mesh_disp(std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,
-        std::shared_ptr<const Core::LinAlg::Vector<double>> disp) override;
+        const Core::LinAlg::Vector<double>& disp) override;
 
     void set_velocity_fields(std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,
         std::shared_ptr<const Core::LinAlg::Vector<double>> convvel,
diff --git a/src/ssi/4C_ssi_monolithic.cpp b/src/ssi/4C_ssi_monolithic.cpp
@@ -568,7 +568,7 @@ void SSI::SsiMono::read_restart(int restart)
  *--------------------------------------------------------------------------*/
 void SSI::SsiMono::prepare_time_loop()
 {
-  set_struct_solution(structure_field()->dispnp(), structure_field()->velnp(),
+  set_struct_solution(*structure_field()->dispnp(), structure_field()->velnp(),
       is_s2i_kinetics_with_pseudo_contact());
 
   // calculate initial potential field if needed
@@ -589,7 +589,7 @@ void SSI::SsiMono::prepare_time_step()
   increment_time_and_step();
 
   // pass structural degrees of freedom to scalar transport discretization
-  set_struct_solution(structure_field()->dispnp(), structure_field()->velnp(),
+  set_struct_solution(*structure_field()->dispnp(), structure_field()->velnp(),
       is_s2i_kinetics_with_pseudo_contact());
 
   // prepare time step for scalar transport field
@@ -1155,7 +1155,7 @@ void SSI::SsiMono::distribute_solution_all_fields(const bool restore_velocity)
   }
 
   // distribute states to other fields
-  set_struct_solution(structure_field()->dispnp(), structure_field()->velnp(),
+  set_struct_solution(*structure_field()->dispnp(), structure_field()->velnp(),
       is_s2i_kinetics_with_pseudo_contact());
   set_scatra_solution(scatra_field()->phinp());
   if (is_scatra_manifold()) set_scatra_manifold_solution(*scatra_manifold()->phinp());
diff --git a/src/ssi/4C_ssi_partitioned_1wc.cpp b/src/ssi/4C_ssi_partitioned_1wc.cpp
@@ -167,7 +167,7 @@ void SSI::SSIPart1WCSolidToScatra::prepare_time_step(bool printheader)
   if (structure_field()->step() % diffsteps == 0)
   {
     if (is_s2i_kinetics_with_pseudo_contact()) structure_field()->determine_stress_strain();
-    set_struct_solution(structure_field()->dispn(), structure_field()->veln(),
+    set_struct_solution(*structure_field()->dispn(), structure_field()->veln(),
         is_s2i_kinetics_with_pseudo_contact());
     scatra_field()->prepare_time_step();
   }
@@ -238,7 +238,7 @@ void SSI::SSIPart1WCSolidToScatra::timeloop()
     if (structure_field()->step() % diffsteps == 0)
     {
       if (is_s2i_kinetics_with_pseudo_contact()) structure_field()->determine_stress_strain();
-      set_struct_solution(structure_field()->dispnp(), structure_field()->velnp(),
+      set_struct_solution(*structure_field()->dispnp(), structure_field()->velnp(),
           is_s2i_kinetics_with_pseudo_contact());
       do_scatra_step();  // It has its own time and timestep variables, and it increments them by
                          // itself.
@@ -289,7 +289,7 @@ void SSI::SSIPart1WCScatraToSolid::timeloop()
 
   // set zero velocity and displacement field for scatra
   auto zeros_structure = Core::LinAlg::create_vector(*structure_field()->dof_row_map(), true);
-  set_struct_solution(zeros_structure, zeros_structure, false);
+  set_struct_solution(*zeros_structure, zeros_structure, false);
 
   scatra_field()->prepare_time_loop();
 
diff --git a/src/ssi/4C_ssi_partitioned_2wc.cpp b/src/ssi/4C_ssi_partitioned_2wc.cpp
@@ -145,7 +145,7 @@ void SSI::SSIPart2WC::do_struct_step()
   if (is_s2i_kinetics_with_pseudo_contact()) structure_field()->determine_stress_strain();
 
   //  set mesh displacement and velocity fields
-  return set_struct_solution(structure_field()->dispnp(), structure_field()->velnp(),
+  return set_struct_solution(*structure_field()->dispnp(), structure_field()->velnp(),
       is_s2i_kinetics_with_pseudo_contact());
 }
 
@@ -196,7 +196,7 @@ void SSI::SSIPart2WC::prepare_time_loop()
   constexpr bool force_prepare = true;
   structure_field()->prepare_output(force_prepare);
   structure_field()->output();
-  set_struct_solution(structure_field()->dispnp(), structure_field()->velnp(), false);
+  set_struct_solution(*structure_field()->dispnp(), structure_field()->velnp(), false);
   scatra_field()->prepare_time_loop();
 }
 
@@ -207,7 +207,7 @@ void SSI::SSIPart2WC::prepare_time_step(bool printheader)
 {
   increment_time_and_step();
 
-  set_struct_solution(structure_field()->dispnp(), structure_field()->velnp(), false);
+  set_struct_solution(*structure_field()->dispnp(), structure_field()->velnp(), false);
   scatra_field()->prepare_time_step();
 
   // if adaptive time stepping and different time step size: calculate time step in scatra
@@ -503,7 +503,7 @@ void SSI::SSIPart2WCSolidToScatraRelax::outer_loop()
     // begin nonlinear solver / outer iteration ***************************
 
     // set relaxed mesh displacements and velocity field
-    set_struct_solution(dispnp, velnp, is_s2i_kinetics_with_pseudo_contact());
+    set_struct_solution(*dispnp, velnp, is_s2i_kinetics_with_pseudo_contact());
 
     // solve scalar transport equation
     do_scatra_step();
diff --git a/src/ssti/4C_ssti_algorithm.cpp b/src/ssti/4C_ssti_algorithm.cpp
@@ -281,10 +281,10 @@ void SSTI::SSTIAlgorithm::test_results(MPI_Comm comm) const
 
 /*----------------------------------------------------------------------*/
 /*----------------------------------------------------------------------*/
-void SSTI::SSTIAlgorithm::distribute_structure_solution()
+void SSTI::SSTIAlgorithm::distribute_structure_solution() const
 {
-  scatra_field()->apply_mesh_movement(structure_->dispnp());
-  thermo_field()->apply_mesh_movement(structure_->dispnp());
+  scatra_field()->apply_mesh_movement(*structure_->dispnp());
+  thermo_field()->apply_mesh_movement(*structure_->dispnp());
 
   // convective velocity is set to zero
   const auto convective_velocity = Core::LinAlg::create_vector(*structure_->dof_row_map());
@@ -295,7 +295,7 @@ void SSTI::SSTIAlgorithm::distribute_structure_solution()
 
 /*----------------------------------------------------------------------*/
 /*----------------------------------------------------------------------*/
-void SSTI::SSTIAlgorithm::distribute_scatra_solution()
+void SSTI::SSTIAlgorithm::distribute_scatra_solution() const
 {
   structure_field()->discretization()->set_state(
       1, "scalarfield", *scatra_->scatra_field()->phinp());
diff --git a/src/ssti/4C_ssti_algorithm.hpp b/src/ssti/4C_ssti_algorithm.hpp
@@ -133,8 +133,8 @@ namespace SSTI
     //! distribute states between subproblems
     //@{
     void distribute_solution_all_fields();
-    void distribute_scatra_solution();
-    void distribute_structure_solution();
+    void distribute_scatra_solution() const;
+    void distribute_structure_solution() const;
     void distribute_thermo_solution();
     //@}
 

Original file line number	Diff line number	Diff line change
`@@ -81,7 +81,7 @@ void ElCh::MovingBoundaryAlgorithm::setup()`
`81`	`81`	`}`
`82`	`82`
`83`	`83`	`// transfer moving mesh data`
`84`		`- scatra_field()->apply_mesh_movement(ale_field()->dispnp());`
	`84`	`+ scatra_field()->apply_mesh_movement(*ale_field()->dispnp());`
`85`	`85`
`86`	`86`	`// initialize the multivector for all possible cases`
`87`	`87`	`fluxn_ = scatra_field()->calc_flux_at_boundary(false);`
`@@ -116,7 +116,7 @@ void ElCh::MovingBoundaryAlgorithm::time_loop()`
`116`	`116`	`}`
`117`	`117`
`118`	`118`	`// transfer moving mesh data`
`119`		`- scatra_field()->apply_mesh_movement(ale_field()->dispnp());`
	`119`	`+ scatra_field()->apply_mesh_movement(*ale_field()->dispnp());`
`120`	`120`
`121`	`121`	`// time loop`
`122`	`122`	`while (not_finished())`
`@@ -278,7 +278,7 @@ void ElCh::MovingBoundaryAlgorithm::solve_scatra()`
`278`	`278`	`}`
`279`	`279`
`280`	`280`	`// transfer moving mesh data`
`281`		`- scatra_field()->apply_mesh_movement(ale_field()->dispnp());`
	`281`	`+ scatra_field()->apply_mesh_movement(*ale_field()->dispnp());`
`282`	`282`
`283`	`283`	`// solve coupled electrochemistry equations`
`284`	`284`	`scatra_field()->solve();`
Original file line number	Diff line number	Diff line change
`@@ -615,11 +615,11 @@ void FS3I::PartFPS3I::set_struct_scatra_solution()`
`615`	`615`	`void FS3I::PartFPS3I::set_mesh_disp()`
`616`	`616`	`{`
`617`	`617`	`// fluid field`
`618`		`- scatravec_[0]->scatra_field()->apply_mesh_movement(fpsi_->fluid_field()->dispnp());`
	`618`	`+ scatravec_[0]->scatra_field()->apply_mesh_movement(*fpsi_->fluid_field()->dispnp());`
`619`	`619`
`620`	`620`	`// Poro field`
`621`	`621`	`scatravec_[1]->scatra_field()->apply_mesh_movement(`
`622`		`- fpsi_->poro_field()->structure_field()->dispnp());`
	`622`	`+ *fpsi_->poro_field()->structure_field()->dispnp());`
`623`	`623`	`}`
`624`	`624`
`625`	`625`
Original file line number	Diff line number	Diff line change
`@@ -203,7 +203,7 @@ void PoroElastScaTra::PoroScatraBase::set_mesh_disp()`
`203`	`203`	`dispnp = volcoupl_fluidscatra_->apply_vector_mapping21(*fluid_field()->dispnp());`
`204`	`204`	`}`
`205`	`205`
`206`		`- scatra_->scatra_field()->apply_mesh_movement(dispnp);`
	`206`	`+ scatra_->scatra_field()->apply_mesh_movement(*dispnp);`
`207`	`207`
`208`	`208`	`std::shared_ptr<const Core::LinAlg::Vector<double>> sdispnp = nullptr;`
`209`	`209`
Original file line number	Diff line number	Diff line change
`@@ -383,7 +383,7 @@ namespace ScaTra`
`383`	`383`	`*`
`384`	`384`	`* @param[in] dispnp displacement vector`
`385`	`385`	`*/`
`386`		`- void apply_mesh_movement(std::shared_ptr<const Core::LinAlg::Vector<double>> dispnp);`
	`386`	`+ void apply_mesh_movement(const Core::LinAlg::Vector<double>& dispnp) const;`
`387`	`387`
`388`	`388`	`//! calculate fluxes inside domain and/or on boundary`
`389`	`389`	`void calc_flux(const bool writetofile //!< flag for writing flux info to file`
Original file line number	Diff line number	Diff line change
`@@ -129,8 +129,7 @@ void SSI::SSICouplingMatchingVolume::set_mechanical_stress_state(`
`129`	`129`	`/----------------------------------------------------------------------/`
`130`	`130`	`/----------------------------------------------------------------------/`
`131`	`131`	`void SSI::SSICouplingMatchingVolume::set_mesh_disp(`
`132`		`- std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,`
`133`		`- std::shared_ptr<const Core::LinAlg::Vector<double>> disp)`
	`132`	`+ std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra, const Core::LinAlg::Vector<double>& disp)`
`134`	`133`	`{`
`135`	`134`	`scatra->scatra_field()->apply_mesh_movement(disp);`
`136`	`135`	`}`
`@@ -278,11 +277,10 @@ void SSI::SSICouplingNonMatchingBoundary::assign_material_pointers(`
`278`	`277`	`/----------------------------------------------------------------------/`
`279`	`278`	`/----------------------------------------------------------------------/`
`280`	`279`	`void SSI::SSICouplingNonMatchingBoundary::set_mesh_disp(`
`281`		`- std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,`
`282`		`- std::shared_ptr<const Core::LinAlg::Vector<double>> disp)`
	`280`	`+ std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra, const Core::LinAlg::Vector<double>& disp)`
`283`	`281`	`{`
`284`	`282`	`scatra->scatra_field()->apply_mesh_movement(`
`285`		`- adaptermeshtying_->master_to_slave(extractor_->extract_cond_vector(disp)));`
	`283`	`+ adaptermeshtying_->master_to_slave(extractor_->extract_cond_vector(disp)));`
`286`	`284`	`}`
`287`	`285`
`288`	`286`	`/----------------------------------------------------------------------/`
`@@ -402,11 +400,10 @@ void SSI::SSICouplingNonMatchingVolume::assign_material_pointers(`
`402`	`400`	`/----------------------------------------------------------------------/`
`403`	`401`	`/----------------------------------------------------------------------/`
`404`	`402`	`void SSI::SSICouplingNonMatchingVolume::set_mesh_disp(`
`405`		`- std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,`
`406`		`- std::shared_ptr<const Core::LinAlg::Vector<double>> disp)`
	`403`	`+ std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra, const Core::LinAlg::Vector<double>& disp)`
`407`	`404`	`{`
`408`	`405`	`scatra->scatra_field()->apply_mesh_movement(`
`409`		`- volcoupl_structurescatra_->apply_vector_mapping21(*disp));`
	`406`	`+ *volcoupl_structurescatra_->apply_vector_mapping21(disp));`
`410`	`407`	`}`
`411`	`408`
`412`	`409`	`/----------------------------------------------------------------------/`
`@@ -604,8 +601,7 @@ void SSI::SSICouplingMatchingVolumeAndBoundary::assign_material_pointers(`
`604`	`601`	`/----------------------------------------------------------------------/`
`605`	`602`	`/----------------------------------------------------------------------/`
`606`	`603`	`void SSI::SSICouplingMatchingVolumeAndBoundary::set_mesh_disp(`
`607`		`- std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra,`
`608`		`- std::shared_ptr<const Core::LinAlg::Vector<double>> disp)`
	`604`	`+ std::shared_ptr<Adapter::ScaTraBaseAlgorithm> scatra, const Core::LinAlg::Vector<double>& disp)`
`609`	`605`	`{`
`610`	`606`	`scatra->scatra_field()->apply_mesh_movement(disp);`
`611`	`607`	`}`