This improvement is due to the use of `x[indices] += y` instead of `np.add.at(x,
indices, y)`, but comes at the cost of not handling duplicate `indices` values.