function A=spdiag(H)
% spdiag - Return elementwise multiplier as sparse diagonal matrix
%
% spdiag :: [[N]] -> [[N,N]].

A=diag(sparse(H));
