[localoperator] Add wrapper that extends volume-only local operator block-diagonal fashion
This MR adds a new class BlockDiagonalLocalOperatorFullCoupling, which takes an existing, volume-only operator and applies it to all leaf spaces in pair of (compatible) ansatz and test function space trees.
As the name implies, it will create a fully coupled matrix (including all cell-local off-diagonal entries).