Begin rewrite of ScalarDiamonScheme

this is just running:
- did not check the validity of the results
- still requires **a lot** of clean-up and rework
- interpolation weights a still recomputed even if the operator is
discretized multiple times on the same mesh