OptAPP adding message to adjoint variable check

applications/OptimizationApplication/custom_python/add_custom_utilities_to_python.cpp

@@ -192,6 +192,7 @@ void  AddCustomUtilitiesToPython(pybind11::module& m)
         .def("GetVariableDimension", &OptimizationUtils::GetVariableDimension<array_1d<double, 3>>)
         .def("SetSolutionStepVariablesList", &OptimizationUtils::SetSolutionStepVariablesList, py::arg("destination_model_part"), py::arg("origin_model_part"))
         .def("IsSolutionStepVariablesListASubSet", &OptimizationUtils::IsSolutionStepVariablesListASubSet, py::arg("main_set_model_part"), py::arg("sub_set_model_part"))
+        .def("GetSolutionStepVariableNamesList", &OptimizationUtils::GetSolutionStepVariableNamesList, py::arg("model_part"))
         .def("GetComponentWiseModelParts", &OptimizationUtils::GetComponentWiseModelParts,
             py::arg("model"),
             py::arg("parameters"),

applications/OptimizationApplication/custom_utilities/optimization_utils.cpp

@@ -150,6 +150,19 @@ bool OptimizationUtils::IsSolutionStepVariablesListASubSet(
     return true;
 }
 
+std::vector<std::string> OptimizationUtils::GetSolutionStepVariableNamesList(const ModelPart& rModelPart)
+{
+    KRATOS_TRY
+
+    std::vector<std::string> result;
+    for (const auto& r_var : rModelPart.GetNodalSolutionStepVariablesList()) {
+        result.push_back(r_var.Name());
+    }
+    return result;
+
+    KRATOS_CATCH("");
+}
+
 std::vector<std::vector<ModelPart*>> OptimizationUtils::GetComponentWiseModelParts(
     Model& rModel,
     Parameters Settings)

applications/OptimizationApplication/custom_utilities/optimization_utils.h

@@ -13,6 +13,8 @@
 #pragma once
 
 // System includes
+#include <string>
+#include <vector>
 
 // Project includes
 #include "includes/define.h"
@@ -97,6 +99,8 @@ class KRATOS_API(OPTIMIZATION_APPLICATION) OptimizationUtils
         const ModelPart& rMainSetModelPart,
         const ModelPart& rSubSetModelPart);
 
+    static std::vector<std::string> GetSolutionStepVariableNamesList(const ModelPart& rModelPart);
+
     static std::vector<std::vector<ModelPart*>> GetComponentWiseModelParts(
         Model& rModel,
         Parameters Settings);

applications/OptimizationApplication/python_scripts/model_part_controllers/connectivity_preserving_model_part_controller.py

@@ -42,8 +42,13 @@ def ImportModelPart(self) -> None:
                 else:
                     # if it is existing, then check whether it has a subset of the variables list.
                     destination_root_model_part = self.model[destination_root_model_part_name]
+                    print(source_model_part)
+                    print(destination_root_model_part)
                     if not KratosOA.OptimizationUtils.IsSolutionStepVariablesListASubSet(source_model_part, destination_root_model_part):
-                        raise RuntimeError(f"The destination root model part solution step variables at {destination_root_model_part.FullName()} is not a sub set of {source_model_part.FullName()}")
+                        msg = f"The destination root model part solution step variables at {destination_root_model_part.FullName()} is not a sub set of {source_model_part.FullName()}"
+                        msg += f"\n\t{destination_root_model_part.FullName()} model part variables list:\n\t\t" + "\n\t\t".join(KratosOA.OptimizationUtils.GetSolutionStepVariableNamesList(destination_root_model_part))
+                        msg += f"\n\t{source_model_part.FullName()} model part variables list:\n\t\t" + "\n\t\t".join(KratosOA.OptimizationUtils.GetSolutionStepVariableNamesList(source_model_part))
+                        raise RuntimeError(msg)
                 KratosOA.OptimizationUtils.SetSolutionStepVariablesList(destination_root_model_part, source_model_part)
                 destination_model_part = self.model.CreateModelPart(destination_model_part_name)
             else: