Bugfix: multiply residual by relaxation factor first

The residual needs to be multiplied by the relaxation factor before assigning it to the resulting field phi.