function    value=KVs(D,V,H)
% adjoint operator of K_V
V=sparse(V);
D=sparse(D);
value=sparse(V'*(2*diag(sum(D))-2*D)*V);
 
