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::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 {