scipy.sparse.dia_matrix.mean
scipy.sparse.dia_matrix.nonzero
Point-wise multiplication by another matrix