Which play well with the latest change to not assemble on the overlap.

The old formula contained some 'dim's instead of '2's,
which obviously works for 2D....

Not all variants appear in YaspGrids. The number of variant
reduces from 4*d*d -> 2*d for skeleton integrals.