The mpfa-l method currently only works on triangles when the volume variables and the fvGeometry caches are enabled. This could be resolved by using a specialized assembly map for the L-scheme.
TODO:
-
make the l-method work on triangular grids without caches