Panzer MiniEM: Fix bug in debug build

trilinos · Aug 17, 2023 · 844586d · 844586d
1 parent 4b126e5
commit 844586d
Showing 1 changed file with 12 additions and 0 deletions.
diff --git a/packages/panzer/mini-em/src/solvers/MiniEM_Interpolation.cpp b/packages/panzer/mini-em/src/solvers/MiniEM_Interpolation.cpp
@@ -155,6 +155,18 @@ Teko::LinearOp buildInterpolation(const Teuchos::RCP<const panzer::LinearObjFact
     ep_rowmap    = global_eloc->getMapForBlock(hoBlockIndex);
     ep_colmap    = ghosted_eloc->getMapForBlock(loBlockIndex);
 
+    {
+      // loop over element blocks
+      std::vector<std::string> elementBlockIds;
+      blockedDOFMngr->getElementBlockIds(elementBlockIds);
+      for(std::size_t blockIter = 0; blockIter < elementBlockIds.size(); ++blockIter) {
+
+        // loop over elements
+        std::vector<int> elementIds = ho_ugi->getElementBlock(elementBlockIds[blockIter]);
+        maxNumElementsPerBlock = std::max(maxNumElementsPerBlock, elementIds.size());
+      }
+    }
+
     // TODO: Fix this.
     size_t nnzPerRowEstimate = 25*loCardinality;