SciPy

scipy.sparse.dia_matrix.dot

dia_matrix.dot(other)[source]