accumulation should be defined by operations, not the arrays

[rjpower]

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)