Solution (source code)

= Solution

Projected gradient descent with unit step is
$$
\boxed{x_{k+1}=P_C[x_k-(x_k-P_Dx_k)]=P_C(P_Dx_k).}
$$
Thus this instance reduces to alternating Euclidean projections.