Fix for elements that might have unassigned partition via parallel De…

…launay meshgenerator
ecmwf · Jan 19, 2024 · b8f8a4d · b8f8a4d
1 parent dc63328
commit b8f8a4d
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/src/atlas/meshgenerator/detail/DelaunayMeshGenerator.cc b/src/atlas/meshgenerator/detail/DelaunayMeshGenerator.cc
@@ -197,6 +197,9 @@ void DelaunayMeshGenerator::generate(const Grid& grid, const grid::Distribution&
                 element_nodes_uncertainty.insert(g_node_connectivity(jelem,2));
                 element_uncertainty.insert(jelem);
             }
+            else if (elem_ownership != GHOST && elem_ownership != CERTAIN) {
+                ATLAS_NOTIMPLEMENTED;
+            }
         }
         // Log::info() << "element_uncertainty" << std::endl;
         // for( auto& jelem: element_uncertainty ) {
@@ -291,6 +294,9 @@ void DelaunayMeshGenerator::generate(const Grid& grid, const grid::Distribution&
                 else if (e2 != UNCERTAIN) {
                     elem_part = e2;
                 }
+                if (elem_part == UNCERTAIN) {
+                    elem_part = elem_node_partition(jelem,0);
+                }
                 if (elem_part == part_) {
                     collect_element(jelem);
                 }