Add debug check

Wells-Group · Jun 18, 2024 · 355e8a2 · 355e8a2
1 parent 0726c33
commit 355e8a2
Showing 1 changed file with 10 additions and 0 deletions.
diff --git a/cpp/Contact.cpp b/cpp/Contact.cpp
@@ -129,6 +129,16 @@ Contact::Contact(
 {
   assert(_mesh);
 
+  if (markers.size() != (std::size_t)surfaces.num_nodes())
+    throw std::runtime_error("maker and surfaces have different sizes.");
+
+  for (int i = 0; i < surfaces.num_nodes(); ++i)
+  {
+    if (surfaces.num_links(i) != 1)
+      throw std::runtime_error("Num surfces links !=0 : "
+                               + std::to_string(surfaces.num_links(i)));
+  }
+
   // TODO: remove the below code, compute num_local directly from
   // _cell_facet_pairs
   for (std::size_t s = 0; s < markers.size(); ++s)