function    value=Ks(D)
% adjoint operator of K
D=sparse(D);
value=sparse(2*diag(sum(D))-2*D);
value=2*diag(sum(D))-2*D;
 
