Provide lifting type flux approximation

[nate-sime]

HDG flux functions incorporating an interior penalty flux function are less resilient to changes in the interior penalty parameter than their DG counterparts. Introducing a parameter-less lifting operator in its place will enable parameter free viscous flux approximation between elements. This requires computation of local problems on patches of the mesh neighbouring the interior facets. Cf. the methods of Bassi and Rebay.