Tweaks to sum across dimensions

- Implement sum(fn, A, region)
- Implement Base.sumabs across dimensions
- Use pairwise summation/BLAS for sums across first dimension