spdiags
