MueLu CoalesceDropFactory: Avoid overflow

Signed-off-by: Christian Glusa <[email protected]>
trilinos · Jan 23, 2025 · fd2554a · fd2554a
1 parent 47d1a5e
commit fd2554a
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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
           {