Boundary conditions must somehow be copied when constructing temporary variables (functions) in solvers
For example this is not happening in the MinRes solver currently. Since it is heavily templated and boundary conditions cannot be set equally for example in scalar vs. StokesFunctions we need a clever solution here ;)
Anyhow, this should be considered during the badly needed refactoring of the solvers.