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