From 844586df5a6ce7520669d1f7a1c76621ba99ed88 Mon Sep 17 00:00:00 2001
From: Christian Glusa <caglusa@sandia.gov>
Date: Thu, 17 Aug 2023 16:10:59 -0600
Subject: [PATCH] Panzer MiniEM: Fix bug in debug build

---
 .../mini-em/src/solvers/MiniEM_Interpolation.cpp     | 12 ++++++++++++
 1 file changed, 12 insertions(+)

diff --git a/packages/panzer/mini-em/src/solvers/MiniEM_Interpolation.cpp b/packages/panzer/mini-em/src/solvers/MiniEM_Interpolation.cpp
index f2af0035db57..dc95851fa840 100644
--- 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;