Cache the mass matrix and its inverse for the dual <-> primal maps.

Makes almost no difference speedwise