Correcting a ug in residual_criterion.hpp. Before this change, the in…

…itial residual was not changing between coupling iterations (See Issue #12394)
KratosMultiphysics · Sep 18, 2024 · d640871 · d640871
1 parent 3338e1c
commit d640871
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/...ions/SolidMechanicsApplication/custom_solvers/convergence_criteria/residual_criterion.hpp b/...ions/SolidMechanicsApplication/custom_solvers/convergence_criteria/residual_criterion.hpp
@@ -136,6 +136,10 @@ class ResidualCriterion : public  ConvergenceCriterion< TSparseSpace, TDenseSpac
                     const SystemVectorType& rb
                     ) override
   {
+    if (rModelPart.GetProcessInfo()[NL_ITERATION_NUMBER] <= 1) {
+      this->Set(LocalFlagType::INITIALIZED, false);
+    }
+
     if (TSparseSpace::Size(rb) != 0) //if we are solving for something
     {