feat: Do not allow linear solver fail feature

src/coreComponents/physicsSolvers/SolverBase.cpp

@@ -96,6 +96,12 @@
     setRestartFlags( RestartFlags::WRITE_AND_READ ).
     setDescription( "Write matrix, rhs, solution to screen ( = 1) or file ( = 2)." );
 
+  registerWrapper( viewKeyStruct::noLinearSolveFailString(), &m_noLinearSolveFail ).
+    setApplyDefaultValue( 0 ).
+    setInputFlag( InputFlags::OPTIONAL ).
+    setRestartFlags( RestartFlags::WRITE_AND_READ ).
+    setDescription( "Cut time step if linear solution fail without going until max nonlinear iterations." );
+
   registerGroup( groupKeyStruct::linearSolverParametersString(), &m_linearSolverParameters );
   registerGroup( groupKeyStruct::nonlinearSolverParametersString(), &m_nonlinearSolverParameters );
   registerGroup( groupKeyStruct::solverStatisticsString(), &m_solverStatistics );
@@ -1015,6 +1021,10 @@
 
       // Output the linear system solution for debugging purposes
       debugOutputSolution( time_n, cycleNumber, newtonIter, m_solution );
+
+      // Do not allow non converged linear solver - cut time step
+      if( m_noLinearSolveFail && m_linearSolverResult.status == LinearSolverResult::Status::NotConverged )
+        return false;
     }
 
     {

src/coreComponents/physicsSolvers/SolverBase.hpp

@@ -586,6 +586,7 @@ class SolverBase : public ExecutableGroup
     static constexpr char const * targetRegionsString() { return "targetRegions"; }
     static constexpr char const * meshTargetsString() { return "meshTargets"; }
     static constexpr char const * writeLinearSystemString() { return "writeLinearSystem"; }
+    static constexpr char const * noLinearSolveFailString() { return "noLinearSolveFail"; }
   };
 
   struct groupKeyStruct
@@ -786,6 +787,9 @@ class SolverBase : public ExecutableGroup
   real64 m_maxStableDt;
   real64 m_nextDt;
 
+  /// behavior in case of linear solver failure
+  integer m_noLinearSolveFail;
+
   /// name of the FV discretization object in the data repository
   string m_discretizationName;