scipy.sparse.bsr_matrix.sum

bsr_matrix.sum(axis=None)

Sum the matrix over the given axis. If the axis is None, sum over both rows and columns, returning a scalar.