Add the ability to insert boundary condition values into a residual vector. This is achieved by adding the ability to evaluate a pointwise expression only on a subset.