Just like the wilson one, this algorithm
• is currently only a 2-level method since I don't have correct implementations
for Mdir and Mdiag in CoarsenedMatrix yet (needed for further coarsening)
• needs levelization and refactoring into a proper algorithm