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