scipy.sparse.spmatrix.setdiag
scipy.sparse.spmatrix.toarray
Sum the matrix over the given axis. If the axis is None, sum over both rows and columns, returning a scalar.