[staggered][localResidual] Allow resizing of faceResidual vector

* method evalFace() has now an optional bool parameter to specify whether the face
  residual should be resized
* this is important to avoid a segfault when only calling this method without
  having called eval() before, which is done e.g. in the multidimension coupling
  localJacobian