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

Only 6 orders of magnitude to go!

It caused quite a bit of trouble with reuse of large matrices.

That was a bad idea from the beginning! It lead to a situation,
where calls from within a generator inherited their context tags.
Should be fixed now.