WIP: Add a Multilevel schwarz preconditioner

[pratikvn]

This PR uses the distributed coarse level generation from PGM to use as a coarse level for the additive Schwarz preconditioner. The only requirement is that the galerkin product generator ($RAP$) generate a multigrid::MultigridLevel which is a triplet of restrict, prolong and coarse operators. The user additionally also needs to set the solver for the coarse level.

TODO

• Add options to switch between multiplicative and additive
• Options for arbitrary number of levels.

Possible issues

1. The coarse level solve significantly reduces the number of iterations, but can be extremely expensive. Maybe makes sense to apply coarse solve, only a few times instead of every preconditioner apply (need to be careful of consistency of preconditioner apply) ?
2. Weighting between the local solver solution and the coarse solution is unclear.
3. Though the coarse solve and the local solve are independent, we have no way of performing them in parallel. This will probably require a rewrite of the apply interface to enable asynchronicity, in addition to support for multiple streams/queues.