Merge pull request #4143 from KratosMultiphysics/feature-core-mortar-mapping-reset-database

[Core] Mortar mapper ResetDatabase() implementation

Author: loumalouomega

kratos/processes/simple_mortar_mapper_process.cpp

@@ -121,6 +121,12 @@ void SimpleMortarMapperProcess<TDim, TNumNodes, TVarType, TNumNodesMaster>:: Exe
 {
     KRATOS_TRY;
 
+    // We reset the database if needed
+    const bool update_interface = mThisParameters["update_interface"].GetBool();
+    if (update_interface) {
+        this->UpdateInterface();
+    }
+
     if (mpThisLinearSolver == nullptr)
         ExecuteExplicitMapping();
     else
@@ -883,6 +889,28 @@ void SimpleMortarMapperProcess<TDim, TNumNodes, TVarType, TNumNodesMaster>::Exec
 /***********************************************************************************/
 /***********************************************************************************/
 
+template< SizeType TDim, SizeType TNumNodes, class TVarType, const SizeType TNumNodesMaster >
+void SimpleMortarMapperProcess<TDim, TNumNodes, TVarType, TNumNodesMaster>::UpdateInterface()
+{
+    // Iterate in the conditions
+    auto &r_destination_conditions_array = mDestinationModelPart.Conditions();
+    #pragma omp parallel for
+    for(int i = 0; i < static_cast<int>(r_destination_conditions_array.size()); ++i) {
+        auto it_cond = r_destination_conditions_array.begin() + i;
+        // Reset the index set
+        if (it_cond->Has(INDEX_SET)) {
+            (it_cond->GetValue(INDEX_SET))->clear();
+        }
+        // Reset the index set
+        if (it_cond->Has(INDEX_MAP)) {
+            (it_cond->GetValue(INDEX_SET))->clear();
+        }
+    }
+}
+
+/***********************************************************************************/
+/***********************************************************************************/
+
 template< SizeType TDim, SizeType TNumNodes, class TVarType, const SizeType TNumNodesMaster >
 Parameters SimpleMortarMapperProcess<TDim, TNumNodes, TVarType, TNumNodesMaster>::GetDefaultParameters()
 {
@@ -901,6 +929,7 @@ Parameters SimpleMortarMapperProcess<TDim, TNumNodes, TVarType, TNumNodesMaster>
         "destination_variable"             : "",
         "origin_variable_historical"       : true,
         "destination_variable_historical"  : true,
+        "update_interface"                 : false,
         "search_parameters"                : {
             "allocation_size"                  : 1000,
             "bucket_size"                      : 4,

kratos/processes/simple_mortar_mapper_process.h

@@ -760,6 +760,16 @@ class KRATOS_API(KRATOS_CORE) SimpleMortarMapperProcess
         }
     }
 
+    /**
+     * @brief Reset the interface database
+     * This method resets the mapping database saved in the destination database.
+     * Note that this needs to be done if such modelpart has changed its number
+     * of nodes or conditions. This needs to be done even though the mapping
+     * instance is deleted since such information is saved in the destination
+     * nodes and conditions.
+     */
+    void UpdateInterface();
+
     /**
      * @brief This method provides the defaults parameters to avoid conflicts between the different constructors
      */

kratos/processes/simple_mortar_mapper_wrapper_process.h

@@ -266,6 +266,7 @@ class SimpleMortarMapperProcessWrapper
             "destination_variable"             : "",
             "origin_variable_historical"       : true,
             "destination_variable_historical"  : true,
+            "update_interface"                 : false,
             "search_parameters"                : {
                 "allocation_size"                  : 1000,
                 "bucket_size"                      : 4,