Consider some sort of assembly or diagonal extractor in LinSys_BDDC