We currently defined the reduction/accum functions on arrays themselves, but this is kind of non-intuitive: what I really want is for the reduction to be attached to the operation, so I could say:
x = zeros()
scatter(src, target=x map_fn=..., accum_fn=np.add)