Bugfix: Set phi correctly in the first smoother iteration.

If the flag initial_phi_is_zero, then in the first smoother iteration phi should just be set to 2*mu*D_diag_inv*r . Before fixing this bug, this was added to phi (which might be undefined or have wrong values), so gave incorrect results.