weak assembly of Dirichlet boundary conditions
In the context of finite volume methods incorporating Dirichlet boundary conditions strongly into the FE-space doesnt make much sense. Therefore, we usually want to assemble them weakly via an additional residual term. Do we want to support this in this module? If so, this needs additional functionality in the Subface
-class, primarily some kind of insideSubvolume()
function