bug fix: numpy has changed sparse.identity, returns a dia_matrix and not a csr anymore