Use the operator (and not the matrix) to calculate defect. In parallel the former

procedure is not correct.

As the pre method of the smoother might (and in parallel actually does) change the
right hand and left hand side, we now copy that information to the outside world.

[[Imported from SVN: r515]]