= Solution
Write $P_f=\operatorname{prox}_f$ and $P_h=\operatorname{prox}_h$. Since $w_k=x_k+z_{k-1}$,
$$
y_k=P_h(w_k),
\qquad z_k=w_k-P_h(w_k),
$$
and therefore
$$
\boxed{w_{k+1}=T(w_k),
\qquad T=I-P_h+P_f(2P_h-I).}
$$
With reflected proximal maps $R_f=2P_f-I$ and $R_h=2P_h-I$,
$$
T=\frac12(I+R_fR_h).
$$
Firm nonexpansiveness of each proximal map is equivalent to nonexpansiveness of its reflection. Thus $R_fR_h$ is nonexpansive, and its average with the identity is firmly nonexpansive. This is the <Douglas–Rachford method>.
Back to article page
