[GeomechanicsApplication] Make quasi newton available for dsettlement

applications/GeoMechanicsApplication/custom_utilities/solving_strategy_factory.hpp

@@ -18,6 +18,7 @@
 #include "builder_and_solver_factory.hpp"
 #include "convergence_criteria_factory.hpp"
 #include "custom_strategies/strategies/geo_mechanics_newton_raphson_strategy.hpp"
+#include "custom_strategies/strategies/geo_mechanics_quasi_newton_strategy.hpp"
 #include "factories/standard_linear_solver_factory.h"
 #include "parameters_utilities.h"
 #include "scheme_factory.hpp"
@@ -76,6 +77,27 @@ class SolvingStrategyFactory
             return result;
         }
 
+        if (rSolverSettings[strategy_type].GetString() == "quasi_newton") {
+            const auto max_iterations    = rSolverSettings["max_iterations"].GetInt();
+            const auto compute_reactions = rSolverSettings["compute_reactions"].GetBool();
+            const auto reform_dof_set_at_each_step = rSolverSettings["reform_dofs_at_each_step"].GetBool();
+            const auto move_mesh_flag = rSolverSettings["move_mesh_flag"].GetBool();
+
+            auto strategy_parameters = Parameters{};
+
+            strategy_parameters.AddValue("quasi_newton_type", rSolverSettings["quasi_newton_type"]);
+            strategy_parameters.AddValue("quasi_newton_restart_interval",
+                                         rSolverSettings["quasi_newton_restart_interval"]);
+            strategy_parameters.AddValue("quasi_newton_max_rank",
+                                         rSolverSettings["quasi_newton_max_rank"]);
+            auto result =
+                std::make_unique<GeoMechanicsQuasiNewtonStrategy<TSparseSpace, TDenseSpace, TLinearSolver>>(
+                    rModelPart, scheme, criteria, builder_and_solver, strategy_parameters,
+                    max_iterations, compute_reactions, reform_dof_set_at_each_step, move_mesh_flag);
+            result->SetEchoLevel(echo_level);
+            return result;
+        }
+
         if (rSolverSettings[strategy_type].GetString() == "line_search") {
             const std::vector<std::string> strategy_entries = {"compute_reactions",
                                                                "max_line_search_iterations",

applications/GeoMechanicsApplication/tests/cpp_tests/custom_utilities/test_solving_strategy_factory.cpp

@@ -145,4 +145,22 @@ KRATOS_TEST_CASE_IN_SUITE(CreatedLineSearchStrategyUsesMaxIterationsInput, Krato
     KRATOS_EXPECT_EQ(static_cast<int>(p_line_search_strategy->GetMaxIterationNumber()), max_number_of_iterations);
 }
 
+KRATOS_TEST_CASE_IN_SUITE(Create_ReturnsSolvingStrategy_ForQuasiNewtonStrategy, KratosGeoMechanicsFastSuite)
+{
+    Model      model;
+    const int  buffer_size      = 2;
+    auto&      dummy_model_part = model.CreateModelPart("dummy", buffer_size);
+    Parameters parameters{testParameters};
+    parameters["strategy_type"].SetString("quasi_newton");
+    parameters.AddString("quasi_newton_type", "broyden");
+    parameters.AddInt("quasi_newton_restart_interval", 50);
+    parameters.AddInt("quasi_newton_max_rank", 10);
+
+    const auto created_strategy = SolvingStrategyFactoryType::Create(parameters, dummy_model_part);
+
+    KRATOS_EXPECT_TRUE(created_strategy->GetMoveMeshFlag()) // since it is set to true in the parameters
+    KRATOS_EXPECT_NE(created_strategy, nullptr);
+    KRATOS_EXPECT_EQ(created_strategy->Check(), 0);
+}
+
 } // namespace Kratos::Testing