Classic aggregation may not be handling Dirichlet BCs correctly. https://github.com/jhux2/Trilinos/blob/main/MueLu_AggregationPhase1Algorithm_def.hpp#L80-L100