[PL/DS] Use unordered set for faster isDeactivated

ProcessLib/CreateDeactivatedSubdomain.cpp

@@ -110,8 +110,8 @@ static DeactivatedSubdomainMesh createDeactivatedSubdomainMesh(
     ranges::copy_if(elements, back_inserter(deactivated_elements), is_active,
                     [](auto const e) { return e->getID(); });
 
-    auto bulk_element_ids =
-        deactivated_elements | MeshLib::views::ids | ranges::to<std::vector>;
+    auto bulk_element_ids = deactivated_elements | MeshLib::views::ids |
+                            ranges::to<std::unordered_set>;
 
     static int mesh_number = 0;
     // Subdomain mesh consisting of deactivated elements.

ProcessLib/DeactivatedSubdomain.cpp

@@ -34,7 +34,7 @@ bool DeactivatedSubdomain::isDeactivated(MeshLib::Element const& element,
                                          double const time) const
 {
     auto const& bulk_element_ids = deactivated_subdomain_mesh.bulk_element_ids;
-    if (!ranges::contains(bulk_element_ids, element.getID()))
+    if (!bulk_element_ids.contains(element.getID()))
     {
         return false;
     }

ProcessLib/DeactivatedSubdomain.h

@@ -16,6 +16,7 @@
 #include <memory>
 #include <optional>
 #include <string>
+#include <unordered_set>
 #include <vector>
 
 #include "MathLib/InterpolationAlgorithms/PiecewiseLinearInterpolation.h"
@@ -40,7 +41,7 @@ struct DeactivatedSubdomainMesh
     /// A mesh created from material ids (independent of time) for the
     /// deactivated subdomain.
     MeshLib::Mesh mesh;
-    std::vector<std::size_t> bulk_element_ids;
+    std::unordered_set<std::size_t> bulk_element_ids;
 
     /// Inner nodes owned only by elements of the deactivated subdomain.
     /// \see ProcessLib::createDeactivatedSubdomainMesh()