WIP: Feature/new localresidual generalize direction
@kweis Same as for the assembly, I'm not sure if you always have all contributions on the boundary (should be visible in the boundary flux helper). If not the better solution might be to expect a scalar. On the other hand these special boudnary conditions also wouldn't be realizable with Crouzeix Raviart elements so maybe this problem doesn't occur.