Write
Pf=proxf and
Ph=proxh. Since
wk=xk+zk−1,
yk=Ph(wk),zk=wk−Ph(wk),
and therefore
wk+1=T(wk),T=I−Ph+Pf(2Ph−I).
With reflected proximal
maps Rf=2Pf−I and
Rh=2Ph−I,
Firm nonexpansiveness of each proximal
map is equivalent to nonexpansiveness of its reflection. Thus
RfRh is nonexpansive, and its
average with the identity is firmly nonexpansive. This is the
Douglas–Rachford method.