function    value=T(D)
%  Operator T of D (psd)
n=size(D,1);
D=sparse(D);
J=sparse(speye(n)-ones(n)/n);
value=-sparse(J*D*J)/2;
values=sparse(value+value')/2;
