Hi, I just noted a small inconsistency in the constructors of the `P*VectorFunction`s compared to those of the scalar functions. All of - P1Function - P2Function - EdgeDoFFunction - FaceDoFFunction provide a constructor version that allows to set a `BoundaryCondition`. The constructor which does not allow this, in turn delegates to the other one using our default variant `BoundaryCondition::create0123BC()` to generate a B.C. object. - P1VectorFunction - P2VectorFunction both currently provide only a single constructor that - does not allow to set the boundary condition - and also does not set the default IMHO we should fix this by extending the two classes. Cheers Marcus