Make RHS BC application consistent

Zero the BC values in the residual rather than setting them to
bc-value. This makes the LHS and RHS consistent.
Actually zero the BC nodes.
Re-implement in terms of other methods
Re-implement in terms of other methods