wip: add option to replace displacements

Author: OmarDuran

src/coreComponents/integrationTests/mechanicTest/testFEMConsistency.cpp

@@ -131,9 +131,9 @@ class ConsistencyTest : public ::testing::TestWithParam< std::tuple< std::string
   </Tasks>  
   <Outputs>
   </Outputs>
-  <Events minTime="-1.0e11" maxTime="1.0">
-     <SoloEvent name="preFracture" target="/Solvers/SurfaceGen" targetTime="-1.0e11" beginTime="-1.0e11"/> 
-    <SoloEvent name="ELASTICITY.PRE.INIT.STEP" targetTime="-1.0e11" beginTime="-1.0e11" target="/Tasks/ELASTICITY.PRE.INIT.STEP"/>
+  <Events minTime="0.0" maxTime="1.0">
+     <SoloEvent name="preFracture" target="/Solvers/SurfaceGen" targetTime="0.0" beginTime="0.0"/> 
+    <SoloEvent name="ELASTICITY.PRE.INIT.STEP" targetTime="0.0" beginTime="0.0" target="/Tasks/ELASTICITY.PRE.INIT.STEP"/>
     <PeriodicEvent name="solverApplications" target="/Solvers/mechSolver" forceDt="1.0"/>
   </Events>
 </Problem>

src/coreComponents/physicsSolvers/solidMechanics/SolidMechanicsInitialization.cpp

@@ -53,6 +53,11 @@ SolidMechanicsInitialization< SOLID_SOLVER >::SolidMechanicsInitialization( cons
     setApplyDefaultValue( "" ).
     setDescription( "Name of the solid mechanics statistics task" );
 
+  registerWrapper( viewKeyStruct::resetDisplacementsString(), &m_resetDisplacements ).
+    setInputFlag( InputFlags::OPTIONAL ).
+    setApplyDefaultValue( 1 ).
+    setDescription( "Flag to reset displacements (1: reset, 0: do not reset)" );
+
   addLogLevel< logInfo::SolverInitialization >();
 }
 
@@ -106,7 +111,10 @@ bool SolidMechanicsInitialization< SOLID_SOLVER >::execute( real64 const time_n,
                                    getName(), time_n, m_solidSolverName ) );
 
   m_solidSolver->setStressInitialization( true );
-  m_solidMechanicsStateResetTask.execute( time_n, dt, cycleNumber, eventCounter, eventProgress, domain );
+  if( m_resetDisplacements )
+  {
+    m_solidMechanicsStateResetTask.execute( time_n, dt, cycleNumber, eventCounter, eventProgress, domain );
+  }
   m_solidSolver->execute( time_n, dt, cycleNumber, eventCounter, eventProgress, domain );
 
   GEOS_LOG_LEVEL_RANK_0( logInfo::SolverInitialization,
@@ -118,9 +126,11 @@ bool SolidMechanicsInitialization< SOLID_SOLVER >::execute( real64 const time_n,
   {
     m_solidMechanicsStatistics->execute( time_n, dt, cycleNumber, eventCounter, eventProgress, domain );
   }
-
-  m_solidMechanicsStateResetTask.execute( time_n, dt, cycleNumber, eventCounter, eventProgress, domain );
-
+  
+  if (m_resetDisplacements){
+    m_solidMechanicsStateResetTask.execute( time_n, dt, cycleNumber, eventCounter, eventProgress, domain );
+  }
+  
   return false;
 }
 

src/coreComponents/physicsSolvers/solidMechanics/SolidMechanicsInitialization.hpp

@@ -60,6 +60,7 @@ class SolidMechanicsInitialization : public TaskBase
   {
     constexpr static char const * solidSolverNameString() { return "solidSolverName"; }
     constexpr static char const * solidMechanicsStatisticsNameString() { return "solidMechanicsStatisticsName"; }
+    constexpr static char const * resetDisplacementsString() { return "reset_displacements"; }
   };
 
   void postInputInitialization() override;
@@ -71,6 +72,8 @@ class SolidMechanicsInitialization : public TaskBase
   SolidMechanicsStatistics * m_solidMechanicsStatistics;
 
   SolidMechanicsStateReset m_solidMechanicsStateResetTask;
+
+  integer m_resetDisplacements;
 };
 
 }