From fd2554ae33e84df8c7b2004d32c6dc71f4b9a4d5 Mon Sep 17 00:00:00 2001
From: Christian Glusa <caglusa@sandia.gov>
Date: Thu, 23 Jan 2025 14:23:05 -0700
Subject: [PATCH] MueLu CoalesceDropFactory: Avoid overflow

Signed-off-by: Christian Glusa <caglusa@sandia.gov>
---
 .../MatrixTransformation/MueLu_CoalesceDropFactory_def.hpp      | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/packages/muelu/src/Graph/MatrixTransformation/MueLu_CoalesceDropFactory_def.hpp b/packages/muelu/src/Graph/MatrixTransformation/MueLu_CoalesceDropFactory_def.hpp
index d848a823743c..505d857a08ed 100644
--- a/packages/muelu/src/Graph/MatrixTransformation/MueLu_CoalesceDropFactory_def.hpp
+++ b/packages/muelu/src/Graph/MatrixTransformation/MueLu_CoalesceDropFactory_def.hpp
@@ -1192,7 +1192,7 @@ void CoalesceDropFactory<Scalar, LocalOrdinal, GlobalOrdinal, Node>::Build(Level
               // Deal with the situation where boundary conditions have only been enforced on rows, but not on columns.
               // We enforce dropping of these entries by assigning a very large number to the diagonal entries corresponding to BCs.
               if (!haveAddedToDiag)
-                localLaplDiagData[row] = STS::rmax();
+                localLaplDiagData[row] = STS::squareroot(STS::rmax());
             }
           }  // subtimer
           {