wolffd@0: function W = dualWlinear(alpha, Z, U, RHO, K)
wolffd@0: 
wolffd@0:     global PsiR;
wolffd@0:     d = size(K, 1);
wolffd@0:     m = length(alpha);
wolffd@0: 
wolffd@0:     W = Z - U - eye(d) /(d * RHO);
wolffd@0:     for i = 1:m
wolffd@0:         W = W + alpha(i) * PsiR{i} / RHO;
wolffd@0:     end
wolffd@0: 
wolffd@0: end