“…Considering a patch of (sub)elements, supposing a piecewise affine, possibly nonconforming, potential approximation, imposing the normal flux continuity, and solving a local linear system on the patch is also the principle of the MPFA method [1,3,19], cf. also [28,41,14,56,42]. In the MPFA method, the fluxes on patches of (sub)elements are locally recovered from the potentials strongly by u h := −S∇p h .…”