[Core] Introducing proper insert method to PointerVectorSet

applications/MappingApplication/custom_utilities/closest_points.cpp

@@ -36,6 +36,11 @@ bool PointWithId::operator<(const PointWithId& rOther) const
     return mDistance < rOther.mDistance;
 }
 
+bool PointWithId::operator==(const PointWithId& rOther) const
+{
+    return Point::operator==(rOther);
+}
+
 void PointWithId::save(Serializer &rSerializer) const
 {
     KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, IndexedObject);

applications/MappingApplication/custom_utilities/closest_points.h

@@ -44,6 +44,8 @@ class KRATOS_API(MAPPING_APPLICATION) PointWithId : public IndexedObject, public
 
     bool operator<(const PointWithId& rOther) const;
 
+    bool operator==(const PointWithId& rOther) const;
+
     double GetDistance() const { return mDistance; }
 
 private:

applications/MeshingApplication/custom_processes/mmg/mmg_process.cpp

@@ -32,7 +32,7 @@
 /* The mappers includes */
 #include "spaces/ublas_space.h"
 #include "mappers/mapper_flags.h"
-#include "factories/mapper_factory.h" 
+#include "factories/mapper_factory.h"
 
 // NOTE: The following contains the license of the MMG library
 /* =============================================================================
@@ -366,6 +366,10 @@ void MmgProcess<TMMGLibrary>::InitializeMeshData()
     for (auto it_dof = mDofs.begin(); it_dof != mDofs.end(); ++it_dof)
         (**it_dof).FreeDof();
 
+    mrThisModelPart.Nodes().Unique();
+    mrThisModelPart.Conditions().Unique();
+    mrThisModelPart.Elements().Unique();
+
     // Generate the maps of reference
     mMmgUtilities.GenerateReferenceMaps(mrThisModelPart, aux_ref_cond, aux_ref_elem, mpRefCondition, mpRefElement);
 
@@ -590,7 +594,7 @@ void MmgProcess<TMMGLibrary>::ExecuteRemeshing()
             ModelPart& r_old_auxiliar_model_part = r_old_model_part.GetSubModelPart("AUXILIAR_COLLAPSED_PRISMS");
             ModelPart& r_auxiliary_model_part = mrThisModelPart.GetSubModelPart("AUXILIAR_COLLAPSED_PRISMS");
 
-            // Define mapper factory 
+            // Define mapper factory
             DEFINE_MAPPER_FACTORY_SERIAL
             if (MapperFactoryType::HasMapper("nearest_element") && mThisParameters["use_mapper_if_available"].GetBool()) {
                 KRATOS_INFO_IF("MmgProcess", mEchoLevel > 0) << "Using MappingApplication to interpolate values" << std::endl;
@@ -1501,7 +1505,7 @@ const Parameters MmgProcess<TMMGLibrary>::GetDefaultParameters() const
                 "max_num_search_iterations"     : 8,
                 "echo_level"                    : 0
             }
-        }, 
+        },
         "extrapolate_contour_values"           : true,
         "surface_elements"                     : false,
         "search_parameters"                    : {

applications/MeshingApplication/custom_utilities/mmg/mmg_utilities.cpp

@@ -2854,7 +2854,7 @@ void MmgUtilities<MMGLibrary::MMG2D>::SetLocalParameter(
     double HMin,
     double HMax,
     double HausdorffValue
-    ) 
+    )
 {
     if ( MMG2D_Set_localParameter(mMmgMesh, mMmgMet, MMG5_Edg, rColor, HMin, HMax, HausdorffValue) != 1)
         KRATOS_ERROR << "Unable to set local parameter" << std::endl;
@@ -2869,7 +2869,7 @@ void MmgUtilities<MMGLibrary::MMG3D>::SetLocalParameter(
     double HMin,
     double HMax,
     double HausdorffValue
-    ) 
+    )
 {
     if ( MMG3D_Set_localParameter(mMmgMesh, mMmgMet, MMG5_Triangle, rColor, HMin, HMax, HausdorffValue) != 1)
         KRATOS_ERROR << "Unable to set local parameter" << std::endl;
@@ -2884,7 +2884,7 @@ void MmgUtilities<MMGLibrary::MMGS>::SetLocalParameter(
     double HMin,
     double HMax,
     double HausdorffValue
-    ) 
+    )
 {
     if ( MMGS_Set_localParameter(mMmgMesh, mMmgMet, MMG5_Triangle, rColor, HMin, HMax, HausdorffValue) != 1)
         KRATOS_ERROR << "Unable to set local parameter" << std::endl;
@@ -3578,6 +3578,10 @@ void MmgUtilities<TMMGLibrary>::ReorderAllIds(ModelPart& rModelPart)
     for(IndexType i = 0; i < r_elements_array.size(); ++i)
         (it_elem_begin + i)->SetId(i + 1);
 
+    r_nodes_array.Unique();
+    r_conditions_array.Unique();
+    r_elements_array.Unique();
+
     KRATOS_CATCH("");
 }
 
@@ -3801,6 +3805,8 @@ void MmgUtilities<TMMGLibrary>::GenerateMeshDataFromModelPart(
         }
     }
 
+    r_nodes_array.Unique();
+
     /* Conditions */
     counter_to_remesh = block_for_each<SumReduction<IndexType>>(
     r_conditions_array,
@@ -3829,6 +3835,8 @@ void MmgUtilities<TMMGLibrary>::GenerateMeshDataFromModelPart(
         }
     }
 
+    r_conditions_array.Unique();
+
     /* Elements */
     counter_to_remesh = block_for_each<SumReduction<IndexType>>(
     r_elements_array,
@@ -3857,6 +3865,8 @@ void MmgUtilities<TMMGLibrary>::GenerateMeshDataFromModelPart(
         }
     }
 
+    r_elements_array.Unique();
+
     // Now we compute the colors
     rColors.clear();
     ColorsMapType nodes_colors, cond_colors, elem_colors;